Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,560,364 A→G 5.2% L189P (CTG→CCG)  ddpB ← D,D‑dipeptide ABC transporter permease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,560,3640AG 5.2% 241.9 / 3.0 77L189P (CTG→CCG) ddpBD,D‑dipeptide ABC transporter permease
Reads supporting (aligned to +/- strand):  ref base A (36/37);  new base G (2/2);  total (38/39)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.57e-01

GTACGAATGTAGTCTTCACTCAATTGTTCCAGCATCGCTGAGCGGATCTGGCGAGCGACAATTCCCAGGTGAACGAACGCCAGCGTTAATGCCGGTAAGATGAGATGTTGCAACGCATTGAAGAAGACTTCACCGTTGCCTTCAAGCAGCGCATCGAGCAGATAAAAGCCGGTAACGTGCGTTGGTGGATCCAGCCAGTCATCAAGCCTTCCGCCGCCGGGAAGAATTTGCAGATGACCATAAAACAGCACAATGACGCCCAGCCCAAGCCAGAACGCAGGTGTGGAGATTCC  >  NC_000913/1560217‑1560509
                                                                                                                                                   |                                                                                                                                                 
gTACGAATGTAGTCTTCACTCAATTGTTCCAGCATCGCTGAGCGGATCTGGCGAGCGACAATTCCCAGGTGAACGAACGCCAGCGTTAATGCCGGTAAGATGAGATGTTGCAACGCATTGAAGAAGACTGCACCGTTGCCTTCAagca                                                                                                                                                   <  4:440927/148‑1 (MQ=255)
 tACGAATGTAGTCTTCACTCAATTGTTCCAGCATCGCTGAGCGGATCTGGCGAGCGACATTTCCCAGGTGAACGAACGCCAGCGTTAATGCCGGTAAGATGAGATGTTGCAACGCATTGAAGAAGACTTCACCGTTGCCTTCAagcag                                                                                                                                                  >  1:148892/1‑148 (MQ=255)
     aaTGTAGTCTTCACTCAATTGTTCCAGCATCGCTGAGCGGATCTGGCGAGCGACAATTCCCAGGTGAACGAACGCCAGCGTTAATGCCGGTAAGATGAGATGTTGCAACGCATTGAAGAAGACTTCACCGTTGCCTTCAAGCAGCGCa                                                                                                                                              >  8:5885/1‑148 (MQ=255)
      aTGTAGTCTTCACTCAATTGTTCCAGCATCGCTGAGCGGATCTGGCGAGCGACAATTCCCAGGTGAACGAACGCCAGCGTTAATGCCGGTAAGATGAGATGTTGCAACGCATTGAAGAAGACTTCACCGTTGCCTTCAAGCAGCGCAt                                                                                                                                             >  4:222734/1‑148 (MQ=255)
                cACTCAATTGTTCCAGCATCGCTGAGCGGATCTGGCGAGCGACAATTCCCAGGTGAACGAACGCCAGCGTTAATGCCGGTAAGATGAGATGTTGCAACGCATTGAAGAAGACTTCACCGTTGCCTTCAAGCAGCGCATCGAGCAg                                                                                                                                      <  5:16084/145‑1 (MQ=255)
                cACTCAATTGTTCCAGCATCGCTGAGCGGATCTGGCGAGCGACAATTCCCAGGTGAACGAACGCCAGCGTTAATGCCGGTAAGATGAGATGTTGCAACGCATTGAAGAAGACTTCACCGTTGCCTTCAAGCAGCGCATCGAGCAGAt                                                                                                                                    <  6:148449/147‑1 (MQ=255)
                 aCTCAATTGTTCCAGCATCGCTGAGCGGATCTGGCGAGCGACAATTCCCAGGTGAACGAACGCCAGCGTTAATGCCGGTAAGATGAGATGTTGCAACGCATTGAAGAAGACTTCACCGTTGCCTTCAAGCAGCGCATCGAGCAGATa                                                                                                                                   >  5:249321/1‑147 (MQ=255)
                  cTCAATTGTTCCAGCATCGCTGAGCGGATCTGGCGAGCGACAATTCCCAGGTGAACGAACGCCAGCGTTAATGCCGGTAAGATGAGATGTTGCAACGCATTGAAGAAGACTTCACCGTTGCCTTCAAGCAGCGCATCGAGCAGATaaa                                                                                                                                 <  5:143508/148‑1 (MQ=255)
                   tCAATTGTTCCAGCATCGCTGAGCGGATCTGGCGAGCGACAATTCCCAGGTGAACGAACGCCAGCGTTAATGCCGGTAAGATGAGATGTTGCAACGCATTGAAGAAGACTTCACCGTTGCCTTCAAGCAGCGCATCGAGCAg                                                                                                                                      >  4:161099/1‑142 (MQ=255)
                   tCAATTGTTCCAGCATCGCTGAGCGGATCTGGCGAGCGACAATTCCCAGGTGAACGAACGCCAGCGTTAATGCCGGTAAGATGAGATGTTGCAACGCATTGAAGAAGACTTCACCGTTGCCTTCAAGCAGCGCATCGAGCAg                                                                                                                                      <  3:161099/142‑1 (MQ=255)
                    cAATTGTTCCAGCATCGCTGAGCGGATCTGGCGAGCGACAATTCCCAGGTGAACGAACGCCAGCGTTAATGCCGGTAAGATGAGATGTTGCAACGCATTGAAGAAGAATTCACCATTGCCGTAAAGCAGCTCATCTAGCAGATAAAAg                                                                                                                               >  6:347109/1‑148 (MQ=255)
                      aTTGTTCCAGCATCGCTGAGCGGATCTGGCGAGCGACAATTCCCAGGTGAACGAACGCCAGCGTTAATGCCGGTAAGATGAGATGTTGCAACGCATTGAAGAAGACTTCACCGTTGCCTTCAAGCAGCGCATCGAGCAGATAAAAGcc                                                                                                                             >  8:315110/1‑148 (MQ=255)
                      aTTGTTCCAGCATCGCTGAGCGGATCTGGCGAGCGACAATTCCCAGGTGAACGAACGCCAGCGTTAATGCCGGTAAGATGAGATGTTGCAACGCATTGAAGAAGACTTCACCGTTGCCTTCAAGCAGCGCATCGAGCAGATAAAAGcc                                                                                                                             <  2:322252/148‑1 (MQ=255)
                      aTTGGTCCAGCATCGCTGAGCGGATCTGGCGAGCGACAATTCCCAGGTGAACGAACGCCAGCGTTAATGCCGGTAAGATGAGATGTTGCAACGCATTGAAGAAGACTTCACCGTTGCGTTCAAGCAGCGCATCGAGCAGATAAAAGcc                                                                                                                             <  2:365352/148‑1 (MQ=255)
                          ttCCAGCATCGCTGAGCGGATCTGGCGAGCGACAATTCCCAGGTGAACGAACGCCAGCGTTAATGCCGGTAAGATGAGATGTTGCAACGCATTGAAGAAGACTTCACCGTTGCCTTCAAGCAGCGCATCGAGCAGATAAAAGCCGGTa                                                                                                                         <  7:198613/148‑1 (MQ=255)
                           tCCAGCATCGCTGAGCGGATCTGGCGAGCGACAATTCCCAGGTGAACGAACGCCAGCGTTAATGCCGGTAAGATGAGATGTTGCAACGCATTGAAGAAGACTTCACCGTTGCCTTCAAGCAGCGCATCGAGCAGATAAAAGCCGGTaa                                                                                                                        >  8:19672/1‑148 (MQ=255)
                              aGCATCGCTGAGCGGATCTGGCGAGCGACAATTCCCAGGTGAACGAACGCCAGCGTTAATGCCGGTAAGATGAGATGTTGCAACGCATTGAAGAAGACTTCACCGTTGCCTTCAAGCAGCGCATCGAGCAGATAAAAGCCGGTaa                                                                                                                        >  2:453703/1‑145 (MQ=255)
                              aGCATCGCTGAGCGGATCTGGCGAGCGACAATTCCCAGGTGAACGAACGCCAGCGTTAATGCCGGTAAGATGAGATGTTGCAACGCATTGAAGAAGACTTCACCGTTGCCTTCAAGCAGCGCATCGAGCAGATAAAAGCCGGTaa                                                                                                                        <  1:453703/145‑1 (MQ=255)
                                cATCGCTGAGCGGATCTGGCGAGCGACAATTCCCAGGTGAACGAACGCCAGCGTTAATGCCGGTAAGATGAGATGTTGCAACGCATTGAAGAAGACTTCACCGTTGCCTTCAAGCAGCGCATCGAGCAGATAAAAGCCGGTAACGTGc                                                                                                                   >  5:86310/1‑148 (MQ=255)
                                  tCGCTGAGCGGATCTGGCGAGCGACAATTCCCAGGTGAACGAACGCCAGCGTTAATGCCGGTAAGATGAGATGTTGCAACGCATTGAAGAAGACTTCACCGTTGCCTTCAAGCGGCGCATCGAGCAGAt                                                                                                                                    >  7:12838/1‑129 (MQ=255)
                                  tCGCTGAGCGGATCTGGCGAGCGACAATTCCCAGGTGAACGAACGCCAGCGTTAATGCCGGTAAGATGAGATGTTGCAACGCATTGAAGAAGACTTCACCGTTGCCTTCAAGCGGCGCATCGAGCAGAt                                                                                                                                    <  8:20005/129‑1 (MQ=255)
                                  tCGCTGAGCGGATCTGGCGAGCGACAATTCCCAGGTGAACGAACGCCAGCGTTAATGCCGGTAAGATGAGATGTTGCAACGCATTGAAGAAGACTTCACCGTTGCCTTCAAGCGGCGCATCGAGCAGAt                                                                                                                                    >  7:20005/1‑129 (MQ=255)
                                  tCGCTGAGCGGATCTGGCGAGCGACAATTCCCAGGTGAACGAACGCCAGCGTTAATGCCGGTAAGATGAGATGTTGCAACGCATTGAAGAAGACTTCACCGTTGCCTTCAAGCGGCGCATCGAGCAGAt                                                                                                                                    <  8:12838/129‑1 (MQ=255)
                                  tCGCTGAGCGGATCCGGCGAGCGACAATTCCCAGGTGAACGAACGCCAGCGTTAATGCCGGTAAGATGAGATGTTGCAACGCATTGAAGAAGACTTCACCGTTGCCTTCAAGCAGCGCATCGAGCAGATAAAAGCCGGTAACGTGCGt                                                                                                                 >  2:404359/1‑148 (MQ=255)
                                   cGCTGAGCGGATCTGGCGAGCGACAATTCCCAGGTGAACGAACGCCAGCGTTAATGCCGGTAAGATGAGATGTTGCAACGCATTGAAGAAGACTTCACCGTTGCCTTCAAGCAGCGCATCGAGCAGATAAAAGCCGGTAACGTGCGtt                                                                                                                >  2:277035/1‑148 (MQ=255)
                                    gCTGAGCGGATCTGGCGAGCGACAATTCCCAGCTGAACGAACGCCAGCGTTAATGCCGGTAAGATGAGATGTTGCAACGCATTGAAGAAGACTTCACCGTTGCCTTCAAGCAGCGCATCGAGCAGATAAAAGCCGGTAACGTGCGTtg                                                                                                               <  5:347109/148‑1 (MQ=255)
                                       gAGCGGATCTGGCGAGCGACAATTCCCAGGTGAACGAACGCCAGCGTTAATGCCGGTAAGATGAGATGTTGCAACGCATTGAAGAAGACTTCACCGTTGCCTTCAAGCAGCGCATCGAGCAGATAAAAGCCGGTAACGTGCGTtggtg                                                                                                            >  7:429842/1‑148 (MQ=255)
                                        aGCGGATCTGGCGAGCGACAATTCCCAGGTGAACGAACGCCAGCGTTAATGCCGGTAAGATGAGATGTTGCAACGCATTGAAGAAGACTTCACCGTTGCCTTCAAGCAGCGCATCGAGCAGATAAAAGCCGGTAACGTGCGTtggtg                                                                                                            <  4:303095/147‑1 (MQ=255)
                                         gCGGATCTGGCGAGCGACAATTCCCAGGTGAACGAACGCCAGCGTTAATGCCGGTAAGATGAGATGTTGCAACGCATTGAAGAAGACTTCACCGTTGCCTTCAAGCAGCGCATCGAGCAGATAAAAGCCGGTAACGTGCGTTGGTGGa                                                                                                          >  3:256633/1‑148 (MQ=255)
                                          cGGATCTGGCGAGCGACAATTCCCAGGTGAACGAACGCCAGCGTTAATGCCGGTAAGATGAGATGTTGCAACGCATTGAAGAAGACTTCACCGTTGCCTTCAAGCAGCGCATCGAGCAGATAAAAGCCGGTAACGTGCGTTGGTGGAt                                                                                                         <  4:400868/148‑1 (MQ=255)
                                          cGGATCTGGCGAGCGACAATTCCCAGGTGAACGAACGCCAGCGTTAATGCCGGTAAGATGAGATGTTGCAACGCATTGAAGAAGACTTCACCGTTGCCTTCAAGCAGCGCATCGAGCAGATAAAAGCCGGTAACGTGCGTTGGTGGAt                                                                                                         <  6:86310/148‑1 (MQ=255)
                                              tCTGGCGAGCGACAATTCCCAGGTGAACGAACGCCAGCGTTAATGCCGGTAAGATGAGATGTTGCAACGCATTGAAGAAGACTTCACCGTTGCCTTCAAGCAGCGCATCGAGCAGATAAAAGCCGGTAACGTGCGTTGGTGGATccag                                                                                                     <  2:419489/148‑1 (MQ=255)
                                                  gcgagcgaCAATTCCCAGGTGAACGAACGCCAGCGTTAATGCCGGTAAGATGAGATGTTGCAACGCATTGAAGAAGACTTCACCGTTGCCTTCAAGCAGCGCATCGAGCAGATAAAAGCCGGTAACGTGCGTTGGTGGATccagcc                                                                                                   <  1:132732/146‑1 (MQ=255)
                                                  gcgagcgaCAATTCCCAGGTGAACGAACGCCAGCGTTAATGCCGGTAAGATGAGATGTTGCAACGCATTGAAGAAGACTTCACCGTTGCCTTCAAGCAGCGCATCGAGCAGATAAAAGCCGGTAACGTGCGTTGGTGGATccagcc                                                                                                   <  2:334980/146‑1 (MQ=255)
                                                  gcgagcgaCAATTCCCAGGTGAACGAACGCCAGCGTTAATGCCGGTAAGATGAGATGTTGCAACGCATTGAAGAAGACTTCACCGTTGCCTTCAAGCAGCGCATCGAGCAGATAAAAGCCGGTAACGTGCGTTGGTGGATccag                                                                                                     <  7:315110/144‑1 (MQ=255)
                                                    gagcgaCAATTCCCAGGTGAACGAACGCCAGCGTTAATGCCGGTAAGATGAGATGTTGCAACGCATTGAAGAAGACTTCACCGTTGCCTTCAAGCAGCGCATCGAGCAGATAAAAGCCGGTAACGTGCGTTGGTGGATCCAGCCAGtc                                                                                               <  7:422097/148‑1 (MQ=255)
                                                       cgaCAATTCCCAGGTGAACGAACGCCAGCGTTAATGCCGGTAAGATGAGATGTTGCAACGCATTGAAGAAGACTTCACCGTTGCCTTCAAGCAGCGCATCGAGCAGATAAAAGCCGGTAACGTGCGTTGGTGGATCCAGCCAGtcat                                                                                             >  5:271476/1‑147 (MQ=255)
                                                         aCAATTCCCAGGTGAACGAACGCCAGCGTTAATGCCGGTAAGATGAGATGTTGCAACGCATTGAAGAAGACTTCACCGTTGCCTTCAAGCAGCGCATCGAGCAGATAAAAGCCGGTAACGTGCGTTGGTGGATCCAGCCAGtcatc                                                                                            >  4:208770/1‑146 (MQ=255)
                                                          cAATTCCCAGGTGAACGAACGCCAGCGTTAATGCCGGTAAGATGAGATGTTGCAACGCATTGAAGAAGACTTCACCGTTGCCTTCAAGCAGCGCATCGAGCAGATAAAAGCCGGTAACGTGCGTTGGTGGATCCAGCCAGTCATCAAg                                                                                         <  8:278563/148‑1 (MQ=255)
                                                             ttCCCAGGTGAACGAACGCCAGCGTTAATGCCGGTAAGATGAGATGTTGCAACGCATTGAAGAAGACTTCACCGTTGCCTTCAAGCAGCGCATCGAGCAGATAAAAGCCGGTAACGTGCGTTGGTGGATcca                                                                                                      <  6:434301/132‑1 (MQ=255)
                                                             ttCCCAGGTGAACGAACGCCAGCGTTAATGCCGGTAAGATGAGATGTTGCAACGCATTGAAGAAGACTTCACCGTTGCCTTCAAGCAGCGCATCGAGCAGATAAAAGCCGGTAACGTGCGTTGGTGGATcca                                                                                                      >  5:434301/1‑132 (MQ=255)
                                                             ttCCCAGGTGAACGAACGCCAGCGTTAATGCCGGTAAGATGAGATGTTGCAACGCATTGAAGAAGACTTCACCGTTGCCTTCAAGCAGCGCATCGAGCAGATAAAAGCCGGTAACGTGCGTTGGTGGATCCAGCCAGTCATCACGCCt                                                                                      >  2:239673/1‑148 (MQ=255)
                                                                ccAGGTGAACGAACGCCAGCGTTAATGCCGGTAAGATGAGATGTTGCAACGCATTGAAGAAGACTTCACCGTTGCCTTCAAGCAGCGCATCGAGCAGATAAAAGCCGGTAACGTGCGTTGGTGGATCCAGCCAGTCATCAAGCCTTc                                                                                    <  1:404359/147‑1 (MQ=255)
                                                                     tGAACGAACGCCAGCGTTAATGCCGGTAAGATGAGATGTTGCAACGCATTGAAGAAGACTTCACCGTTGCCTTCAAGCAGCGCATCGAGCAGATAAAAGCCGGTAACGTGCGTTGGTGGATCCAGCCAGTCATCAAGCCTTccgccgc                                                                              >  5:278753/1‑148 (MQ=255)
                                                                       aacgaacgCCAGCGTTAATGCCGGTAAGATGAGATGTTGCAACGCATTGAAGAAGACTTCACCGTTGCCTTCAAGCAGCGCATCGAGCAGATAAAAGCCGGTAACGTGCGTTGGTGGATCCAGCCAGTCATCAAGCCTTccgccgccg                                                                            <  6:271476/148‑1 (MQ=255)
                                                                         cgaacgCCAGCGTTAATGCCGGTAAGATGAGATGTTGCAACGCATTGAAGAAGACTTCACCGTTGCCTTCAAGCAGCGCATCGAGCAGATAAAAGCCGGTAACGTGCGTtgg                                                                                                              <  7:145694/112‑1 (MQ=255)
                                                                         cgaacgCCAGCGTTAATGCCGGTAAGATGAGATGTTGCAACGCATTGAAGAAGACTTCACCGTTGCCTTCAAGCAGCGCATCGAGCAGATAAAAGCCGGTAACGTGCGTtgg                                                                                                              >  8:145694/1‑112 (MQ=255)
                                                                         cgaacgCCAGCGTTAATGCCGGTAAGATGAGATGTTGCAACGCATTGAAGAAGACTTCACCGTTGCCTTCAAGCAGCGCATCGAGCAGATAAAAGCCGGTAACGTGCGTTGGTGGATCCAGCCAGTCATCAAGCCTTCCGCCGCCggg                                                                          <  8:429842/148‑1 (MQ=255)
                                                                          gaacgCCAGCGTTAATGCCGGTAAGATGAGATGTTGCAACGCATTGAAGAAGACTTCACCGTTGCCTTCAAGCAGCGCATCGAGCAGATAAAAGCCGGTAACGTGCGTTGGTGGATCCAGCCAGTCATCAAGCCTTCCGCCGCCggg                                                                          >  1:278140/1‑147 (MQ=255)
                                                                          gaacgCCAGCGTTAATGCCGGTAAGATGAGATGTTGCAACGCATTGAAGAAGACTTCACCGTTGCCTTCAAGCAGCGCATCGAGCAGATAAAAGCCGGTAACGTGCGTTGGTGGATCCAGCCAGTCATCAAGCCTTCCGCCGCCGGga                                                                         <  2:148892/148‑1 (MQ=255)
                                                                              gCCAGCGTTAATGCCGGTAAGATGAGATGTTGCAACGCATTGAAGAAGACTTCACCGTTGCCTTCAAGCAGCGCATCGAGCAGATAAAAGCCGGTAACGTGCGTTGGTGGATCCAGCCAGTCATCAAGCCTTCCGCCGCCGGgaagaa                                                                     <  7:5885/148‑1 (MQ=255)
                                                                                cAGCGTTAATGCCGGTAAGATGTGATGTTGCAACGCATTGAAGAAGACTTCACCGTTGCCTTCAAGCAGCGCATCGAGCAGATAAAAGCCGGTAACGTGCGTTGGTGGATCCAGCCAGTCATCAAGCCTTCCGCCGCCGGGAAGAAtt                                                                   <  5:2189/148‑1 (MQ=255)
                                                                                     ttAATGCCGGTAAGATGAGATGTTGCAACGCATTGAAGAAGACTTCACCGTTGCCTTCAAGCAGCGCATCGAGCAGATAAAAGCCGGTAACGTGCGTTGGTGGATCCAGCCAGTCATCAAGCCTTCCGCCGCCGGGAAGAATTTg                                                                 >  3:262964/1‑145 (MQ=255)
                                                                                     ttAATGCCGGTAAGATGAGATGTTGCAACCCATTGAAGAAGACTTCACCGTTGCCTTCAAGCAGCGCATCGAGCAGATAAAAGCCGGTAACGTGCGTTGGTGGATCCAGCCAGTCATCAAGCCTTCCGCCGCCGGGAAGAATTTGCAg                                                              <  6:278753/148‑1 (MQ=255)
                                                                                                                    aTTGAAGAAGACTTCACCGTTGCCTTCAAGCAGCGCATCGAGCAGATAAAAGCCGGTAACGTGCGTTGGTGGATCCAGCCAGTCATCAAGCCTTCCGCCGCCGGGAAGAATTTGCAGATGACCATAAAACAGCACAATGACgcccagc                               >  6:291123/1‑148 (MQ=255)
                                                                                                                      tGAAGAAGACTTCACCGTTGCCTTCAAGCAGCGCATCGAGCAGATAAAAGCCGGTAACGTGCGTTGGTGGATCCAGCCAGTCATCAAGCCTTCCGCCGCCGGGAAGAATTTGCAGATGACCATAAAACAGCACAATGACgcccagcc                              >  5:73075/1‑147 (MQ=255)
                                                                                                                        aagaagACTTCACCGTTGCCTTCAAGCAGCGCATCGAGCAGATAAAAGCCGGTAACGTGCGTTGGTGGATCCAGCCAGTCATCAAGCCTTCCGCCGCCGGGAAGAATTTGCAGATGACCATAAAACAGCACAATGACgcccagccca                            <  3:222734/147‑1 (MQ=255)
                                                                                                                          gaagaCTTCACCGTTGCCTTCAAGCAGCGCATCGAGCAGATAAAAGCCGGTAACGAGCGTTGCTGGATCCAGCCAGTCATCAAGCCTTCCGCCGCCGGCTACAATTTGCAGATGACCATAAAAAAGCACAATCACGCCCAGCCCAAGc                         >  6:451567/1‑148 (MQ=255)
                                                                                                                           aagaCTTCACCGTTGCCTTCAAGCAGCGCATCGAGCAGATAAAAGCCGGTAACGTGCGTTGGTGGATCCAGCCAGTCATCAAGCCTTCCGCCGCCGGGAAGAATTTGCAGATGACCATAAAACAGCACAATGACGCCCAGCCCAAGcc                        <  1:277035/148‑1 (MQ=255)
                                                                                                                              aCTTCACCGTTGCCTTCAAGCAGCGCATCGAGCAGATAAAAGCCGGTAACGTGCGTTGGTGGATCCAGCCAGTCATCAAGCCTTCCGCCGCCGGGAAGAATTTGCAGATGACCATAAAACAGcaca                                           <  5:100783/126‑1 (MQ=255)
                                                                                                                              aCTTCACCGTTGCCTTCAAGCAGCGCATCGAGCAGATAAAAGCCGGTAACGTGCGTTGGTGGATCCAGCCAGTCATCAAGCCTTCCGCCGCCGGGAAGAATTTGCAGATGACCATAAAACAGcaca                                           >  6:100783/1‑126 (MQ=255)
                                                                                                                                 tCACCGTTGCCTTCAAGCAGCGCATCGAGCAGATAAAAGCCGGTAACGTGCGTTGGTGGATCCAGCCAGTCATCAAGCCTTCCGCCGCCGGGAAGAATTTGc                                                                >  1:353004/1‑102 (MQ=255)
                                                                                                                                 tCACCGTTGCCTTCAAGCAGCGCATCGAGCAGATAAAAGCCGGTAACGTGCGTTGGTGGATCCAGCCAGTCATCAAGCCTTCCGCCGCCGGGAAGAATTTGc                                                                <  2:353004/102‑1 (MQ=255)
                                                                                                                                 tCACCGTTGCCTTCAAGCAGCGCATCGAGCAGATAAAAGCCGGTAACGTGCGTTGGTGGATCCAGCCAGTCATCAAGCCTTCCGCCGCCGGGAAGAATTTGCAGATGACCATAAAACAGCACAATGACGCCCAGCCCAAGCCAGAACg                  >  7:252822/1‑148 (MQ=255)
                                                                                                                                  cACCGTTGCCTTCAAGCAGCGCATCGAGCAGATAAAAGCCGGTAACGTGCGTTGGTGGATCCAGCCAGTCATCAAGCCTTCCGCCGCCGGGAAGAATTTGCAGATGACCATAAAACAGCACAATGACGCCCAGCCCAAGCCAGAACg                  >  5:9294/1‑147 (MQ=255)
                                                                                                                                   aCCGTTGCCTTCAAGCAGCGCATCGAGCAGATAAAAGCCGGTAACGTGCGTTGGTGGATCCAGCCAGTCATCAAGCCTTCCGCCGCCGGGAAGAATTTGCAGATGACCATAAAACAGCACAATGACGCCCAGCCCAAGCCAGAACg                  >  4:137680/1‑146 (MQ=255)
                                                                                                                                     cGTTGCCTTCAAGCAGCGCATCGAGCAGATAAAAGCCGGTAACGTGCGTTGGTGGATCCAGCCAGTCATCAAGCCTTCCGCCGCCGGGAAGAATTTGCAGATGACCATAAAACAGCACAATGACGCCCAGCCCAAGCCAGAACGCAgg              >  2:29307/1‑148 (MQ=255)
                                                                                                                                       ttGCCTTCAAGCAGCGCATCGAGCAGATAAAAGCCGGTAACGTGCGTTGGTGGATCCAGCCAGTCATCAAGCCTTCCGCCGCCGGGAAGAATTTGCAGATGACCATAAAACAGCACAATGACGCCCAGCCCAAGCCAGAACGCAGgtg            <  7:19672/148‑1 (MQ=255)
                                                                                                                                         gCCTTCAAGCAGCGCATCGAGCAGATAAAAGCCGGTAACGTGCGTTGGTGGATCCAGCCAGTCATCAAGCCTTCCGCCGCCGGGAAGAATTTGCAGATGACCATAAAACAGCACAATGACGCCCAGCCCAAGCCAGAACGCAGgtgt           <  1:29307/147‑1 (MQ=255)
                                                                                                                                          ccTTCAAGCAGCGCATCGAGCAGATAAAAGCCGGTAACGTGCGTTGGTGGATCCAGCCAGTCATCAAGCCTTCCGCCGCCGGGAAGAATTTGCAGATGACCATaaa                                                   <  2:449609/106‑1 (MQ=255)
                                                                                                                                          ccTTCAAGCAGCGCATCGAGCAGATAAAAGCCGGTAACGTGCGTTGGTGGATCCAGCCAGTCATCAAGCCTTCCGCCGCCGGGAAGAATTTGCAGATGACCATaaa                                                   >  1:449609/1‑106 (MQ=255)
                                                                                                                                          ccTTCAAGCAGCGCATCGAGCAGATAAAAGCCGGTAACGTGCGTTGGTGGATCCAGCCAGTCATCAAGCCTTCCGCCGCCGGGAAGAATTTGCAGATGACCATAAAACAGCACAATGACGCCCAGCCCAAGCCAGAACGCAGGtgtg          <  7:383486/147‑1 (MQ=255)
                                                                                                                                           cTTCAAGCAGCGCATCGAGCAGATAAAAGCCGGTAACGTGCGTTGGTGGATCCAGCCAGTCATCAAGCCTTCCGCCGCCGGGAAGAAt                                                                    >  7:318811/1‑88 (MQ=255)
                                                                                                                                           cTTCAAGCAGCGCATCGAGCAGATAAAAGCCGGTAACGTGCGTTGGTGGATCCAGCCAGTCATCAAGCCTTCCGCCGCCGGGAAGAAt                                                                    <  8:318811/88‑1 (MQ=255)
                                                                                                                                           cTTCAAGCAGCGCATCGAGCAGATAAAAGCCGGTAACGTGCGTTGGTGGATCCAGCCAGTCATCAAGCCTTCCGCCGCCGGGAAGAATTTGCAGATGACCATAAAACAGCACAATGACGCCCAGCCCAAGCCAGAACGCAGGtgtg          >  1:349247/1‑146 (MQ=255)
                                                                                                                                            ttCAAGCAGCGCATCGAGCAGATAAAAGCCGGTAACGTGAGTTGGTGGATCCAGCCAGTCATCAAGCCTTCCGCCGCCGGGAAGAATTTGCAGATGACCATAAAAGAGCACAATg                                        <  3:84023/115‑1 (MQ=255)
                                                                                                                                             tCAAGCAGCGCATCGAGCAGATAAAAGCCGGTAACGTGCGTTGGTGGATCCAGCCAGTCATCAAGCCTTCCGCCGCCGGGAAGAATTTGCAGATGACCATAAAACAGCACAATGACGCCCAGCCCAAGCCAGAACGCAGGtgtg          >  5:227846/1‑144 (MQ=255)
                                                                                                                                             tCAAGCAGCGCATCGAGCAGATAAAAGCCGGTAACGTGCGTTGGTGGATCCAGCCAGTCATCAAGCCTTCCGCCGCCGGGAAGAATTTGCAGATGACCATAAAACAGCACAATGACGCCCAGCCCAAGCCAGAACGCAGGtgtg          >  2:327879/1‑144 (MQ=255)
                                                                                                                                             tCAAGCAGCGCATAGAGCAGATAAAAGCCGGTAACGTGCGTTGGTGGATCCAGCCAGTCATCAAGCCATCCGCCGCCGGGAAGAATTAGCAGATGACCATAAAACATCACAATGa                                       >  4:84023/1‑115 (MQ=255)
                                                                                                                                                 gcagcGCATCGAGCAGATAAAAGCCGGTAACGTGCGTTGGTGGATCCAGCCAGTCATCAAGCCTTCCGCCGCCGGGAAGAATTTGCAGATGACCATAAAACAGCACAATGACGCCCAGCCCAAGCCAGAACGCAGGTGTGGAGATTcc  <  6:73075/148‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                 
GTACGAATGTAGTCTTCACTCAATTGTTCCAGCATCGCTGAGCGGATCTGGCGAGCGACAATTCCCAGGTGAACGAACGCCAGCGTTAATGCCGGTAAGATGAGATGTTGCAACGCATTGAAGAAGACTTCACCGTTGCCTTCAAGCAGCGCATCGAGCAGATAAAAGCCGGTAACGTGCGTTGGTGGATCCAGCCAGTCATCAAGCCTTCCGCCGCCGGGAAGAATTTGCAGATGACCATAAAACAGCACAATGACGCCCAGCCCAAGCCAGAACGCAGGTGTGGAGATTCC  >  NC_000913/1560217‑1560509

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: