Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,560,207 G→T 5.1% A241A (GCC→GCA 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,2070GT 5.1% 393.3 / 7.5 126A241A (GCC→GCAddpBD,D‑dipeptide ABC transporter permease
Reads supporting (aligned to +/- strand):  ref base G (59/60);  new base T (2/5);  total (61/65)
Fisher's exact test for biased strand distribution p-value = 4.41e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.10e-01

Only 100 of 128 total aligned reads displayed.

CACCCATTCCGGGCCAGGCAAAAACGGTTTCGGTGAGCACTGCGCCATACAACAAATCGCCCAGCGCCAAACCCAATACGGTAATCGATGGGATCAACGCATTGGGTAGCGCATAACATAAAACGATATACCAGCCGGGCAAGCCGCTGGCCCGGGCGGTACGAATGTAGTCTTCACTCAATTGTTCCAGCATCGCTGAGCGGATCTGGCGAGCGACAATTCCCAGGTGAACGAACGCCAGCGTTAATGCCGGTAAGATGAGATGTTGCAACGCATTGAAGAAGACTTCACCG  >  NC_000913/1560059‑1560351
                                                                                                                                                    |                                                                                                                                                
cACCCATTCCGGGCTAGGCAAAAACGGTTTCGGTGAGCACTGCGCCATACAAAAAATCGCCCAGCGCCAAACCCAATACGGTAATCGATGGGATCAACGCATTGGGTAGCGCATAACATAAAACGATATACCAGCCGGGCAAGCCGCTgg                                                                                                                                                 <  2:1352282/150‑1 (MQ=255)
  cccATTCCGGGCCAGGCAAAAACGGTTTCGGTGAGCACTGCGCCATACAACAAATCGCCCAGCGCCAAACCCAATACGGTAATCGATGGGATCAACGCATTGGGTAGCGCATAACATAAAACGATATACCAGCCGGGCAAGCCGCTGGcc                                                                                                                                               >  2:1732675/1‑150 (MQ=255)
   ccATTCCGGGCCAGGCAAAAACGGATTCGGTGAGCACTGCGCCATACAACAAATCGCCCAGCGCCAAACCCACTACGGTAATCGATGGGATCAACGCATTGGGTAGCGCATAACATAAAACGATATACCAGCCGGGCAAGCCGCTGGccc                                                                                                                                              <  1:317448/150‑1 (MQ=255)
     aTTCCGGGCCAGGCAAAAACGGTTTCGGTGAGCACTGCGCCATACAACAAATCGCCCAGCGCCAAACCCAATACGGTTATCGATGGGATCAACGCATTGGGTAGCGCATAACATAAAACGATATACCAGCCGGGCAAGCCGCTGGCCCgg                                                                                                                                            <  1:471419/150‑1 (MQ=255)
      ttCCGGGCCAGGCAAAAACGGTTTCGGTGAGCACTGCGCCATACAACAAATCGCCCAGCGCCAAACCCAATACGGTAATCGATGGGATCAACGCATTGGGTAGCGCATAACATAAAACGATATACCAGCCGGGCAAGCCGCTGGCCCggg                                                                                                                                           <  2:1484873/150‑1 (MQ=255)
       tCCGGGCCAGGCAAAAACGGTTTCGGTGAGCACTGCGCCATACAACAAATCGCCCAGCGCCAAACCCAATACGGTAATCGATGGGATCAACGCATTGGGTAGCGCATAACATAAAACGATATACCAGCCGGGCCAGCGGCTTGCCCggct                                                                                                                                          <  1:900650/150‑3 (MQ=255)
         cGGGCCAGGCAAAAACGGTTTCGGTGAGCACTGCGCCATACAACAAATCGCCCAGCGCCAAACCCAATACGGTAATCGATGGGATCAACGCATTGGGTAGCGCATAACATAAAACGATATACCAGCCGGGCCAGCGGCTTGCCCggctgg                                                                                                                                        <  2:687279/150‑5 (MQ=255)
         cGGGCCAGGCAAAAACGGTTTCGGTGAGCACTGCGCCATACAACAAATCGCCCAGCGCCAAACCCAATACGGTAATCGATGGGATCAACGCATTGGGTAGCGCATAACATAAAACGATATACCAGCCGGGCCAGCGGCTTGCCCggctgg                                                                                                                                        <  2:509803/150‑5 (MQ=255)
          gggCCAGGCAAAAACGGTTTCGGTGAGCACTGCGCCATACAACAAATCGCCCAGCGCCAAACCCAATACGGTAATCGATGGGATCAACGCATTGGGTAGCGCATAACATAAAACGATATACCAGCCGGGCAAGCCGCTGGCCCGGGCGGt                                                                                                                                       >  1:2121887/1‑150 (MQ=255)
          gggCCAGGCAAAAACGGTTTCGGTGAGCACTGCGCCATACAACAAATCGCCCAGCGCCAAACCCAATACGGTAATCGATGGGATCAACGCATTGGGTAGCGCATAACATAAAACGATATACCAGCCGGGCAAGCCGCTGGCCCGGGCGGt                                                                                                                                       >  2:2007228/1‑150 (MQ=255)
          gggCCAGGCAAAAACGGTTTCGGTGAGCACTGCGCCATACAACAAATCGCCCAGCGCCAAACCCAATACGGTAATCGATGGGATCAACGCATTGGGTAGCGCATAACATAAAACGATATACCAGCCGGGCAAGCCGCTGGCCCGGGCGGt                                                                                                                                       <  1:741207/150‑1 (MQ=255)
           ggCCAGGCAAAAACGGTTTCGGTGAGCACTGCGCCATACAACAAATCGCCCAGCGCCAAACCCAATACGGTAATCGATGGGATCAACGCATTGGGTAGCGCATAACATAAAACGATATACCAGCCGGGCAAGCCGCTGGCCCGGGCGGTa                                                                                                                                      <  1:889805/150‑1 (MQ=255)
             ccAGGCAAAAACGGTTTCGGTGAGCACTGCGCCATACAACAAATCGCCCAGCGCCAAACCCAATACGGTAATCGATGGGATCAACGCATTGGGTAGCGCATAACATAAAACGATATACCAGCCGGGCAAGCCGCTGGCCCGGGCGGTACg                                                                                                                                    <  1:1271688/150‑1 (MQ=255)
              cAGGCAAAAACGGTTTCGGTGAGCACTGCGCCATACAACAAATCGCCCAGCGCCAAACCCAATACGGTAATCGATGGGATCAACGCATTGGGTAGCGCATAACATAAAACGATATACCAGCCGGGCAAGCCGCTGGCCCGGGCGGTACGa                                                                                                                                   >  2:2213794/1‑150 (MQ=255)
                ggCAAAAACGGTTTCGGTGAGCACTGCGCCATACAACAAATCGCCCAGCGCCAAACCCAATACGGTAATCGATGGGATCAACGCATTGGGTAGCGCATAACATAAAACGATATACCAGCCGGGCAAGCCGCTGGCCCGGGCGGTACGAAt                                                                                                                                 <  2:511020/150‑1 (MQ=255)
                 gCAAAAACGGTTTCGGTGAGCACTGCGCCATACAACAAATCGCCCAGCGCCAACCCCAATACGGTAATCGATGGGATCAACGCAGTGGGGAGCGCATAACATAAAACGATATACCAGCCGGGCAAGCCGCTGGCCCGGGCGGTACGAATg                                                                                                                                >  1:145357/1‑150 (MQ=255)
                 gCAAAAACGGTTTCGGTGAGCACTGCGCCATACAACAAATCGCCCAGCGCCAAACCCAATACGGTAATCGATGGGATCAACGCATTGGGTAGCGCATAACATAAAACGATATACCAGCCGGGCAAGCCGCTGGCCCGGGCGGTACGAATg                                                                                                                                >  2:654138/1‑150 (MQ=255)
                  cAAAAACGGTTTCGGTGAGCACTGCGCCATACAACAAATCGCCCAGCGCCAAACCCAATACGGTAATCGATGGGATCAACGCATTGGGTAGCGCATAACATAAAACGATATACCAGCCGGGCAAGCCGCTGGCCCGGGCGGTACGAATGt                                                                                                                               <  2:1293893/150‑1 (MQ=255)
                    aaaaCGGTTTCGGTGAGCACTGCGCCATACAACAAATCGCCCAGCGCCAAACCCAATACGGTAATCGATGGGATCAACGCATTGGGTAGCGCATAACATAAAACGATATACCAGCCGGGCAAGCCGCTGGCCCGGGCGGTACGAATGTAg                                                                                                                             >  2:1677587/1‑150 (MQ=255)
                    aaaaCGGTTTCGGTGAGCACTGCGCCATACAACAAATCGCCCAGCGCCAAACCCAATACGGTAATCGATGGGATCAACGCATTGGGTAGCGCATAACATAAAACGATATACCAGCCGGGCAAGCCGCTGGCCCGGGCGGTACGAATGTAg                                                                                                                             <  1:1206334/150‑1 (MQ=255)
                          gTTTCGGTGAGCACTGCGCCATACAACAAATCGCCCAGCGCCAAACCCAATACGGTAATCGATGGGATAAACGCATTGGGTAGCGCATAACATAAAACGATATACCAGCCGGGCAAGCCGCTGGCCCGGGCGGTACGAATGTAGTCTTCa                                                                                                                       >  1:1104849/1‑150 (MQ=255)
                                gTGAGCACTGCGCCATACAACAAATCGCCCAGCGCCAAACCCAATACGGTAATCGATGGGATCAACGCATTGGGTAGCGCATAACATAAAACGATATACCAGCCGGGCAAGCCGCTGGCCCGGGCGGTACGAATGTAGTCTTCACTCAAt                                                                                                                 <  2:541503/150‑1 (MQ=255)
                                gTGAGCACTGCGCCATACAACAAATCGCCCAGCGCCAAACCCAATACGGTAATCGATGGGATCAACGCATTGGGTAGCGCATAACATAAAACGATATACCAGCCGGGCAAGCCGCTGGCCCGGGCGGTACGAATGTAGTCTTCACTCAAt                                                                                                                 <  2:257014/150‑1 (MQ=255)
                                gTGAGCACTGCGCCATACAACAAATCGCCCAGCGCAAAACCCAATACGGTAATCGCTGGGAGCAACGCATTGGGTAGCGCTTAACATAAAACGATATACCAGCCGGGCAAGCCGCTGGCCCGGGCGGTACGAATGTAGTCTTCACTCAAt                                                                                                                 <  1:293830/150‑1 (MQ=255)
                                 tGAGCACTGCGCCATACAACAAATCGCCCAGCGCCAAACCCAATACGGTAATCGATGGGATCAACGCATTGGGTAGCGCATAACATAAAACGATATACCAGCCGGGCAAGCCGCTGGCCCGGGCGGTACGAATGTAGTCTTCACTCAAtt                                                                                                                >  2:847835/1‑150 (MQ=255)
                                 ggAGCACTGCGCCATACAACAACTCGCCCAGCGCCAAACCCACTACGGTAATCGATGGGATCAACGCATTGGGTAGCGCATAACATAAAACGATATACCAGCCGGGCAAGCCGCTGGCCCGGGCGGTACGAATGTAGTCTTCACTCAAtt                                                                                                                <  1:111443/149‑1 (MQ=255)
                                    gCACTGCGCCATACAACAAATCGCCCAGCGCCAAACCCAATACGGTAATCGATGGGATCAACGCATTGGGTAGCGCATAACATAAAACGATATACCAGCCGGGCAAGCCGCTGGCCCGGGCGGTACGAATGTAGTCTTCACTCAATTGtt                                                                                                             >  2:1725606/1‑150 (MQ=255)
                                      aCTGCGCCATACAACAAATCGCCCAGCGCCAAACCCAATACGGTAATCGATGGGATCAACGCATTGGGTAGCGCATAACATAAAACGATATACCAGCCGGGCAAGCCGCTGGCCCGGGCGGTACGAATGTAGTCTTCACTCAATTGTTcc                                                                                                           <  2:2121088/150‑1 (MQ=255)
                                        tGCGCCATACAAGAAATCGCCCAGCGACAAACCCAATACGGTAATCGATGGGATCAACGCATTGGGTATCGCATAACATAAAACGATATACCAGCCGGGCAAGCCGCTTGCCCGGGCGGTACGAATGTAGTCTTCACTCAATTGTTCCAg                                                                                                         <  1:1089703/150‑1 (MQ=255)
                                         gcgcCATACAACAAATCGCCCAGCGCCAAACCCAATACGGTAATCGATGGGATCAACGCATTGGGTAGCGCATAACATAAAACGATATACCAGCCGGGCAAGCCGCTGGCCCGGGCGGTACGAATGTAGTCTTCACTCAATTGTTCCAGc                                                                                                        >  1:1404906/1‑150 (MQ=255)
                                                acaacaAATCGAGCAGCGCCAAAACCAATACGGTAATCGATGTGATCAACGCATTGGGTAGCGCATAACATAAAACGATATACCAGCCTTGCAAGCCGCTGGCCCGGGCGGTACGAATGTAGTCTTCACTCAATTGTTCCAGCATCGCTg                                                                                                 <  1:413316/150‑1 (MQ=255)
                                                aaaacaaATCGCCCAGCGACAAAACAAATACGTTAATCGATTGGATCCACGCAGTGGGTAGCGCATAACATAAAACGAAATAAAAGCCTGGCAATACGCTGGCCCGGGCGTTACTAAAGTAGGCTTAACACATTGGTTCATGCAGCGCtt                                                                                                 >  2:2015667/3‑149 (MQ=255)
                                                 caacaaATCGCCCAGCGCCAAACCCAATACGGTAATCGATGGGATCAACGCATTGGGTAGCGCATAACATAAAACGATATACCAGCCGGGCAAGCCGCTGGCCCGGGCGGTACGAATGTAGTCTTCACTCAATTGTTCCAGCATCGCTGa                                                                                                >  1:1855155/1‑150 (MQ=255)
                                                 caacaaATCGCCCAGCGCCAAACCCAATACGGTAATCGATGGGATCAACGCATTGGGTAGCGCATAACATAAAACGATATACCAGCCGGGCAAGCCGCTGGCCCGGGCGGTACGAAAGTAGTCTTCACTCAATTGTTCCAGCATCGCTgc                                                                                                >  1:567790/1‑149 (MQ=255)
                                                  aacaaATCGCCCAGCGCCAAACCCAATACGGTAATCGATGGGATCAACGCATTGGGTAGCGCATAACATAAAACGATATACCAGCCGGGCAAGCCGCTGGCCCGGGCGGGACGCATGTAGTCTTCACTCAATTGTTCCAGCATCGCTGAg                                                                                               >  1:1978299/1‑150 (MQ=255)
                                                   acaaATCGCCCAGCGCCAAACCCAATACGGTAATCGATGGGATCAACGCATTGGGTAGCGCATAACATAAAACGATATACCAGCCGGGCAAGCCGCTGGCCCGGGCGGTACGAATGTAGTCTTCACTCAATTGTTCCAGCATCGCTGAGc                                                                                              <  2:223145/150‑1 (MQ=255)
                                                   acaaATCGCCCAGCGCCAAACCCAATACGGTAATCGATGGGATCAACGCATTGGGTAGCGCATAACATAAAACGATATACCAGCCGGGCAAGCCGCTGGCCCGGGCGGTACGAATGTAGTCTTCACTCAATTGTTCCAGCATCGCTGAGc                                                                                              <  1:1870388/150‑1 (MQ=255)
                                                    caaATCGCCCAGCGCCAAACCCAATACGGTAATCGATGGGATCAACGCATTGGGTAGCGCATAACATAAAACGATATACCAGCCGGGCAAGCCGCTGGCCCGGGCGGTACGAATGTAGTCTTCACTCAATTGTTCCAGCATCGCTGAGCg                                                                                             >  1:2102994/1‑150 (MQ=255)
                                                    caaATCGCCCAGCGCCAAACCCAATACGGTAATCGATGGGATCAACGCATTGGGTAGCGCATAACATAAAACGATATACCAGCCGGGCAAGCCGCTGGCCCGGGCGGTACGAATGTAGTCTTCACTCAATTGTTCCAGCATCGCTGAGCg                                                                                             <  2:2063291/150‑1 (MQ=255)
                                                         cGCCCAGCGCCAAACCCAATACGGTAATCGATGGGATCAACGCATTGGGTAGCGCATAACATAAAACGATATACCAGCCGGGCAAGCCGCTGGCCCGGGCGGTACGAATGTAGTCTTCACTCAATTGTTCCAGCATCGCTGAGCGGATCt                                                                                        >  2:2010959/1‑150 (MQ=255)
                                                           cccAGCGCCAAACCCAATACGGTAAGCGATGGGATCAACGCATTGGGTAGCGCATAACATAAAACGATATACCAGCCGGGCAAGCCGCTGGCCCGGGCGGTACGAATGTAGTCTTCACTCAATTGTTCCAGCATCGCTGAGCGGATCTgg                                                                                      <  2:2139757/150‑1 (MQ=255)
                                                                cgcCAAACCCAATACGGTAATCGATGGGATCAACGCATTGGGTAGCGCATAACATAAAACGATATACCAGCCGGGCAAGCCGCTGGCCCGGGCGGTACGAATGTAGTCTTCACTCAATTGTTCCAGCATCGCTGAGCGGATCTGgcgagc                                                                                 <  1:3517/150‑1 (MQ=255)
                                                                 gcCAAACCCAATACGGTAATCGATGGGATCAACGCATTGGGTAGCGCATAACATAAAACGATATACCAGCCGGGCAAGCCGCTGGCCCGGGCGGTACGAATGTAGTCTTCACTCAATTGTTCCAGCATCGCTGAGCGGATCTGgcgagcg                                                                                >  2:1798402/1‑150 (MQ=255)
                                                                 gcCAAACCCAATACGGTAATCGATGGGATCAACGCATTGGGTAGCGCATAACATAAAACGATATACCAGCCGGGCAAGCCGCTGGCCCGGGCGGTACGAATGTAGTCTTCACTCAATTGTTCCAGCATCGCTGAGCGGATCTGgcgagcg                                                                                >  2:1110344/1‑150 (MQ=255)
                                                                   cAAACCCAATACGGTAATCGATGGGATCAACGCATTGGGTAGCGCATAACATAAAACGATATACCAGCCGGGCAAGCCGCTGGCCCGGGCGGTACGAATGTAGTCTTCACTCAATTGTTCCAGCATCGCTGAGCGGATCTGGCGAGCGAc                                                                              <  2:898631/150‑1 (MQ=255)
                                                                   cAAACCCAATACGGTAATCGATGGGATCAACGCATTGGGTAGCGCATAACATAAAACGATATACCAGCCGGGCAAGCCGCTGGCCCGGGCGGTACGAATGTAGTCTTCACTCAATTGTTCCAGCATCGCTGAGCGGATCTGGCGAGCGAc                                                                              <  2:758417/150‑1 (MQ=255)
                                                                      aCCCAATACGGTAATCGATGGGATCAACGCATTGGGTAGCGCATAACATAAAACGATATACCAGCCGGGCAAGCCGCTGGCCCGGGCGGTACGAATGTAGTCTTCACTCAAATGTTCCAGCATCGCTGAGCGGATCTGGCGCGCGCCAAt                                                                           >  2:196893/1‑150 (MQ=255)
                                                                        ccAATACGGTAATCGATGGGATCATCGCATTGGGTAGCGCATAACATATAACGATATACCAGCCGGTCTAGCCGCTGGCCCGGGCGGTACGAATGTAGTCTTCACTCAATTGTTCCAGCATCGCTGAGCGGATCTGGCGAGCGACAATTc                                                                         <  2:164617/150‑1 (MQ=255)
                                                                        ccAATACGGTAATCGATGGGATCAACGCATTGGGTAGCGCATAACATAAAACGATATACCAGCCGGGCAAGCCGCTGGCCCGGGCGGTACGAATGTAGTCTTCACTCAATTGTTCCAGCATCGCTGAGCGGATCTGGCGAGCGACAATTc                                                                         <  2:1106296/150‑1 (MQ=255)
                                                                         cAATACGGTAATCGATGGGATCAACGCATTGGGTAGCGCATAACATAAAACGATATACCAGCCGGGCAAGCCGCTGGCCCGGGCGGTACGAATGTAGTCTTCACTCAAATGTTCCAGCATCGCTGAGCGGATCTGGCGAGCGACAATTcc                                                                        >  1:1927915/1‑150 (MQ=255)
                                                                          aaTACTGTAATCTATGGGATCAACGCATTGGGTAGCGCATAACATAAAACGATATACCAGCCGGGCAAGCCGCTGGCCCGGGCGGTACGAATGTAGTCTTCACTCAATTGTTCCAGCATCGCTGAGCGGATCTGGCGAGCGACAATTccc                                                                       <  2:332333/150‑1 (MQ=255)
                                                                             aCGGTAATCGATGGGATCAACGCATTGGGTAGCGCATAACATAAAACGATATACCAGCCGGGCAAGCCGCTGGCCCGGGCGGTACGAATGTAGTCTTCACTCAATTGTTCCAGCATCGCTGAGCGGATCTGGCGAGCGACAATTCCCAgg                                                                    <  2:206583/150‑1 (MQ=255)
                                                                                  aaTCGATGGGATCAACGCATTGGGTAGCGCATAACATAAAACGATATACCAGCCGGGCAAGCCGCGGGCCCGGGCGGTACGAATGTAGTCTTCACTCAATTGTTCCAGCATCGCTGAGCGGATCTGGCGAGCGACAATTCCCAGGTgaac                                                               <  1:38250/150‑1 (MQ=255)
                                                                                   aTCGATGGGATCAACGCATTGGGTAGCGCATAACATAAAACGATATACCAGCCGGGCAAGCCGCTGGCCCGGGCGGTACGAATGTAGTCTTCACTCAATTGTTCCAGCATCGCTGAGCGGATCTGGCGAGCGACAATTCCCAGGTgaacg                                                              >  1:667910/1‑150 (MQ=255)
                                                                                      gATGGGATCAACGCATTGGGTAGCGCATAACATAACACGATATACCAGCCGGGCAAGCCGCTGGCCCGGGCGGTACGAATGTAGTCTTCACTCAATTGTTCCAGCATCGCTGAGCGGATCTGGCGAGCGACAATTCCCAGGTgaacgaac                                                           <  1:1747188/150‑1 (MQ=255)
                                                                                      gATGGGATCAACGCATTGGGTAGCGCATAACATAAAACGATATACCAGCCGGGCAAGCCGCTGGCCCGGGCGGTACGAATGTAGTCTTCACTCAATTGTTCCAGCATCGCTGAGCGGATCTGGCGAGCGACAATTCCCAGGTgaacgaac                                                           >  1:681687/1‑150 (MQ=255)
                                                                                             tCAACGCATTGGGTAGCGCATAACATAAAACGATATACCAGCCGGGCAAGCCGCTGGCCCGGGCGGTACGAATGTAGTCTTCACTCAATTGTTCCAGCATCGCTGAGCGGATATGGCGAGCGACAATTCCCAGGTGAACGAACGCAAGCg                                                    >  2:1110664/1‑150 (MQ=255)
                                                                                               aaCGCATTGGGTAGCGCATAACATAAAACGATATACCAGCCGGGCAAGCCGCTGGCCCGGGCGGTACGAATGTAGTCTTCACTCAATTGTTCCAGCATCGCTGTGCGGATCTGGCGAGCGACAATTCCCACGTGAACGAACGCCAGCGtt                                                  >  2:609501/1‑150 (MQ=255)
                                                                                                  gCATTGGGTAGCGCATAACATAAAACGATATACCAGCCGGGCAAGCCGCTGGCCCGGGCGGTACGAATGTAGTCTTCACTCAATTGTTCCAGCATCGCTGAGCGGATCTGGCGAGCGACAATTCCCAGGTGAACGCACGCCAGCGTTAAt                                               >  1:1508387/1‑150 (MQ=255)
                                                                                                   cATTGGGTAGCGCATAACATAAAACGATATACCAGCCGGGCAAGCCGCTGGCCCGGGCGGTACGAATGTAGTCTTCACTCAATTGTTCCAGCATCGCTGAGCGGATCTGGCGAGCGACAATTCCCAGGTGAACGAACGCCAGCGTTAATg                                              <  1:1655478/150‑1 (MQ=255)
                                                                                                      tGGGTAGCGCATAACATAAAACGATATACCAGCCGGGCAAGCCGCTGGCCCGGGCGGTACGAATGTAGTCTTCACTCAATTGTTCCAGCATCGCTGAGCGGATCTGGCGAGCGACAATTCCCAGGTGAACGAACGCCAGCGTTAATGCCg                                           <  1:968343/150‑1 (MQ=255)
                                                                                                       gggTAGCGCATAACATAAAACGATATACCAGCCGGGCAAGCCGCTGGCCCGGGCGGTACGAATGTAGTCTTCACTCAATTGTTCCAGCATCGCTGAGCGGATCTGGCGAGCGACAATTCCCAGGTGAACGAACGCCAGCGTTAATGCCgg                                          >  2:1204287/1‑150 (MQ=255)
                                                                                                        ggTAGCGCATAACATAAAACGATATACCAGCCGGGCAAGCCGCTGGCCCGGGCGGTACGAATGTAGTCTTCACTCTATTGTTCCAGCATCGCTGAGCGGATCTGGCGAGCGACAATTCCCAGGTGAACGAACGCCAGCGTTAATGCCGGt                                         >  2:1250160/1‑150 (MQ=255)
                                                                                                              gcATAACATAAAACGATATACCAGCCGGGCAAGCCGCTGGCCCGGGCGTTACGAATGTAGTCTTCACTCAAATGTTCCATCATCTCTGAGAGGATCTGGCGAGCGACAATTCCCAGGTGAACGAACGCCAGAGTTAATGCCGGGAagatg                                   >  2:1450272/1‑150 (MQ=255)
                                                                                                                ataacataaAACGATATACCAGCCGGGCAAGCCGCTGGCCCGGGCGGTACGAATGTAGTCTTCACTCAATTGTTCCAGCATCGCTGAGCGGATCTGGCGAGCGACAATTCCCAGGTGAACGAACGCCAGCGTTAATGCCGGTAagatgag                                 <  1:89658/150‑1 (MQ=255)
                                                                                                                ataacataaAACGATATACCAGCCGGGCAAGCCGCTGGCCCGGGCGGTACGAATGTAGTCTTCACTCAATTGTTCCAGCATCGCTGAGCGGATCTGGCGAGCGACAATTCCCAGGTGAACGAACGCCAGCGTTAATGCCGGTAagatgag                                 <  1:952080/150‑1 (MQ=255)
                                                                                                                 taacataaAACGATATACCAGCCGGGCAAGCCGCTGGCCCGGGCGGTACGAATGTAGTCTTCACTCAATTGTTCCAGCATCGCTGAGCGGATCTGGCGAGCGACAATTCCCAGGTGACCGAACGCCAGCGTTAATGCCGGTAtgatgaga                                >  2:1931671/1‑150 (MQ=255)
                                                                                                                    cataaAACGATATACCAGCCGGGCAAGCCGCTGGCCCGGGCGGTACGAATGTAGTCTTCACTCAATTGTTCCAGCATCGCTGAGCGGATCTGGCGCGCGACAATTCCCAGGTGAACGAACGCCAGCGTTAATGCCGGTAAGATGAGATGt                             >  1:1649795/1‑150 (MQ=255)
                                                                                                                    cataaAACGATATACCAGCCGGGCAAGCCGCTGGCCCGGGCGGTACGAATGTAGTCTTCACTCAATTGTTCCAGCATCGCTGAGCGGATCTGGCGAGCGACAATTCCCAGGTGAACGAACGCCAGCGTTAATGCCGGTAAGATGAGATGt                             <  2:1650457/150‑1 (MQ=255)
                                                                                                                      taaAACGATATACCAGCCGGGCAAGCCGCTGGCCCGGGCGGTACGAATGTAGTCTTCACTCAATTGTTCCAGCATCGCTGAGCGGATCTGGCGAGCGACAATTCCCAGGTGAACGCACGCCAGCGTTAATGCCGGTAAGATGAGATGTTg                           >  2:289923/1‑150 (MQ=255)
                                                                                                                      taaAACGATATACCAGCCGGGCAAGCCGCTGGCCCGGGCGGTACGAATGTAGTCTTCACTCAATTGTTCCAGCATCGCTGAGCGGATCTGGCGAGCGACAATTCCCAGGTGAACGAACGCCAGCGTTAATGCCGGTAAGATGAGATGTTg                           <  2:1404906/150‑1 (MQ=255)
                                                                                                                        aaaCGATATACCAGCCGGGCAAGCCGCTGGCCCGGGCGGTACGAATGTAGTCTTCACTCAATTGTTCCAGCATCGCTGAGCGGATCTGGCGAGCGACAATTCACAGGTGAACGAACGCCAGCGTTAATGCCGGTAAGATGAGATGTTGCa                         >  2:1214620/1‑150 (MQ=255)
                                                                                                                         aaCGATATACCAGCCGGGCAAGCCGCTGGCCCGGGCGGTACGTATGTAGGCTTCACTCAATTGTTCCACCATCGCTGAGCGGATCTGGCGAGCGACAATTCCCAGGTGAACGAACGCCAGCGTTAATGCCGGTAAGATGAGATGTTGCaa                        <  2:856580/150‑1 (MQ=255)
                                                                                                                          aCGATATACCAGCCGGGCAAGCCGCTGGCCCGGGCGGTACGAATGTAGTCTTCACTCAATTGTTCCAGCATCGCTGAGCGGATCTGGCGAGCGACAATTCCCAGGTGAACGAACGCCAGCGTTAATGCCGGTAAGATGAGATGTTGCAAc                       >  1:1049487/1‑150 (MQ=255)
                                                                                                                           cGATATACCAGCCGTGCAAGCCGCGGGCCCGGGCGGGACGAATGTAGTCTTCACTCAATTTTTCCATCATAGCTGAGCGGATCTGGCGAGCGACAATTCCCAGGTGAACGAACGCCAGCGTTAATGCCGGTAAGATGAGATGTTGCAACg                      <  2:1427953/150‑1 (MQ=255)
                                                                                                                           cGATATACCAGCCGGGCAAGCCGCTGGCCCGGGCGGTACGAATGTAGTCTTCACTCAATTGTTCCAGCATCGCTGAGCGGATCTGGCGAGCGACAATTCCCAGGTGAACGAACGCCAGCGTTAATGCCGGTAAGATGAGATGTTGCAACg                      <  1:847835/150‑1 (MQ=255)
                                                                                                                           cGATATACCAGCCGGGCAAGCCGCTGGCCCGGGCGGTACGAATGTAGTCTTCACTCAATTGTTCCAGCATCGCTGAGCGGATCTGGCGAGCGACAATTCCCAGGTGAACGAACGCCAGCGTTAATGCCGGTAAGATGAGATGTTGCAACg                      <  2:721532/150‑1 (MQ=255)
                                                                                                                            gATATACCAGCCGGGCAAGCCGCTGGCCCGGGCGGTACGAATGTAGTCTTCACTCAATTGTTCCAGCATCGCTGAGCGGATCTGGCGAGCGACAATTCCCAGGTGAACGAACGCCAGCGTTAATGCCGGTAAGATGAGATGTTGCAACGc                     <  2:928454/150‑1 (MQ=255)
                                                                                                                            gATATACCAGCCGGGCAAGCCGCTGGCCCGGGCGGTACGAATGTAGTCTGCACTCACTTGTTCCAGCATCGCTGAGCGGATCTGGCGAGCGACAATTCCCAGGTGAACGCACGCCAGCGTTAAGGCCTGTATGATGAGATGTTGCAACGc                     >  1:412753/1‑150 (MQ=255)
                                                                                                                                 aCCAGCCGGGCAAGCCGCTGGCCCGGGCGGTACGAATGTAGTCTTCACTCAATTGTTCCAGCATCGCTGAGCGGATCTGGCGAGCGACAATTCCCAGGTGAACGAACGCCAGCGTTAATGCCGGTAAGATGAGATGTTGCAACGCATTga                <  1:2007228/150‑1 (MQ=255)
                                                                                                                                  ccAGCCGGGCAAGCCGCTGGCCCGGGCGGTACGAATGTAGTCTTCACTCAATTGTTCCAGCATCGCTGAGCGGATCTGGCGAGCGACAATTCCCAGGTGAACGAACGCCAGCGTTAATGCCGGTAAGATGAGATGTTGCAACGCATTgaa               >  2:1291237/1‑150 (MQ=255)
                                                                                                                                   caggagtgCAAGCCGCTGGCCAGGGCGGTACGAATGTAGTATTCACTCAATTGTTCCAGCATCGCTGAGCGGATCTGGCGAGCGACAATTCCCAGGTGAACGAACGCCAGCGTTAATGCCGGTAAGATGAGATGTTGCAACGCATTgaag              <  1:1223539/143‑1 (MQ=255)
                                                                                                                                   cAGCCGGGCAAGCCGCTGGCCCGGGCGGTACGAATGTAGTCTTCACTCAATTGTTCCAGCATCGCTGAGCGGATCTGGCGAGCGACAATTCCCAGGTGAACGAACGCCAGCGTTAATGCCGGTAAGATGAGATGTTGCAACGCATTgaag              >  2:1933674/1‑150 (MQ=255)
                                                                                                                                    gcccGGGCCAGCGGCTTGCCCGGGCGGTACGAATGTAGTCTTCACTCAATTGTTCCAGCATTGCTGAGCGGCTCTGGCGAGCGACAATTCCCAGGTGAACGAACGCCAGCGTTAATGCCGGTAAGATGAGATGTTGCAACGCATTgaaga             >  1:1733910/3‑150 (MQ=255)
                                                                                                                                    gcccGGGCCAGCGGCTTGCCCGGGCGGTACGAATGTAGTCTTCACTCAATTGTTCCAGCATCGCTGAGCGGATCTGGCGAGCGACAATTCCCAGGTGAACGAACGCCAGCGTTAATGCCGGTAAGATGAGATGTTGCAACGCATTgaaga             >  1:492888/3‑150 (MQ=255)
                                                                                                                                    aGCCGGGCAAGCCGCTGGCCCGGGCGGTACGAATGTAGTCTTCACTCAATTGTTCCAGCATCGCTGAGCGGATCTGGCGAGCGACAATTCCCAGGTGAACGAACGCCAGCGTTAATGACGGTAAGATGAGATGTTGCAACGCATTgaaga             >  1:2164703/1‑150 (MQ=255)
                                                                                                                                     gCCGGTCAAGCCGCTGGCCCGGGCGGTACGAATGTAGTCTTCACTCAATTGGTCCAGCATCGCTGAGCGGATCTGGCGAGCGACAATTCCCAGGTGAACGAACGCCAGCGTTAATGCCGGTAAGATGAGATGTTGCAACGCATTgaagaa            <  1:418582/150‑1 (MQ=255)
                                                                                                                                      ccGGGCAAGCCGCTGGCCCGGGCGGTACGAATGTAGTCTTCACTCAATTGTTCCAGCATCGCTGAGCGGATCTGGCGAGCGACAATTCCCAGGTGAACGAACGCCAGCGTTAATGCCGGTAAGATGAGATGTTGCAACGCATTGaagaag           <  2:2069233/150‑1 (MQ=255)
                                                                                                                                       cGGGCAAGCCGCTGGCCCGGGCGGTACGAATGTAGTCTTCACTCAATTGTTCCAGCATCGCTGAGCGGATCTGGCGCGCGACAATTCCCCGGTGAACGAACGCCAGCGTTAATGCCGGTAAGATGAGATGTTGCACCGCATTGAagaaga          >  1:942960/1‑150 (MQ=255)
                                                                                                                                        gggCAAGCCGCTGGCCCGGGCGGTACGAATGTAGTCTTCACTCAATTGTTCCAGCATCGCTGAGCGGATCTGGCGAGCGACAATTCCCAGGTGAACGAACGCCAGCGTTAATGCCGGTAAGATGAGATGTTGCAACGCATTGAAGAAGAc         <  2:659553/150‑1 (MQ=255)
                                                                                                                                        gggCAAGCCGCTGGCCCGGGCGGTACGAATGTAGTCTTCACTCAATTGTTCCAGCATCGCTGAGCGGATCTGGCGAGCGACAATTCCCAGGTGAACGAACGCCAGCGTTAATGCCGGTAAGATGAGATGTTGCAACGCATTGAAGAAGAc         >  2:12915/1‑150 (MQ=255)
                                                                                                                                           cAAGCCGCTGGCCCGGGCGGTACGAATGTAGTCTTCACTCAATTGTTCCAGCATCGCTGAGCGGATCTGGCGAGCGACAATTCCCAGGTGAACGAACGCCAGCGTTAATGCCGGTAAGATGAGATGTTGCAACGCATTGAAGAAGACTTc      >  2:2012476/1‑150 (MQ=255)
                                                                                                                                           cAAGCCGCTGGCCCGGGCGGTACGAATGTAGTCTTCACTCAATTGTTCCAGCATCGCTGAGCGGATCTGGCGAGCGACAATTCCCAGGTGAACGAACGCCAGCGTTAATGCCGGTAAGATGAGATGTTGCAACGCATTGAAGAAGACTTc      >  1:2025364/1‑150 (MQ=255)
                                                                                                                                           cAAGCCGCTGGCCCGGGCGGTACGAATGTAGTCTTCACTCAATTGTTCCAGCATCGCTGAGCGGATCTGGCGAGCGACAATTCCCAGGTGAACGAACGCCAGCGTTAATGCCGGTAAGATGAGATGTTGCAACGCATTGAAGAAGACTTc      >  2:759736/1‑150 (MQ=255)
                                                                                                                                           cAAGCCGCTGGCCCGGGCGGTACGAATGTAGTCTTCACTCAATTGTTCCAGCATCGCTGAGCGGATCTGGCGAGCGACAATTCCCAGGTGAACGAACGCCAGCGTTAAAGCCGGTAAGATGAGATGAAGCACCGCATTGATGAAGCCTTc      >  2:117288/1‑150 (MQ=255)
                                                                                                                                             aGCCGCTGGCCCGGGCGGTACGAATGTAGTCTTCACTCAATTGTTCCAGCATCGCTGAGCGGATCTGGCGAGCGACAATTCCCAGGTGAACGAACGCCAGCGTTAATGCCGGTAAGATGAGATGTTGCAACGCATTGAAGAAGACTTccc    >  2:1002110/1‑148 (MQ=255)
                                                                                                                                             aGCCGCTGGCCCGGGCGGTACGAATGTAGTCTTCACTCAATTGTTCCAGCATCGCTGAGCGGATCTGGCGAGCGACAATTCCCAGGTGAACGAACGCCAGCGTGAGTGCCGGTATGATGAGATGTTGGGGCGCAGTGAAGATGAGttgag    >  2:2045599/1‑147 (MQ=255)
                                                                                                                                              gCCGCTGGCCCGGGCGGTACGAATGTAGTATTCACTCAATTGGTCCAGGATCGCTGAGAGGATCTGGCGAGCGACATATCCCAGGTCAACGAACTCCAGCTTTAATGACGGTAAGATGAGATGTTGCAACGCATTGAAGAAGACTTCAcc   <  1:336855/150‑1 (MQ=255)
                                                                                                                                              cccGCTGGCCCGGGCGGTACGAATGTAGTCTTCACTCAATTGTTCCAGCATCGCTGAGCGGATCTGGCGAGCGACAATTCCCAGGTGAACGAACGCCAGCGTTAATGCCGGTAAGATGAGATGTTGCAACGCATTGAAGAAGACTTCAcc   <  1:2045599/149‑1 (MQ=255)
                                                                                                                                               ccGCTGGCCCGGGCGGTACGAATGTAGTCTTCACTCAATTGTTCCAGCATCGCTGAGCGGATCTGGCGAGCGACAATTCCCAGGTGAACGAACGCCAGCGTTAATGCCGGTAAGATGAGATGTTGCAACGCATTGAAGAAGACTTCACCg  <  1:1798402/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                
CACCCATTCCGGGCCAGGCAAAAACGGTTTCGGTGAGCACTGCGCCATACAACAAATCGCCCAGCGCCAAACCCAATACGGTAATCGATGGGATCAACGCATTGGGTAGCGCATAACATAAAACGATATACCAGCCGGGCAAGCCGCTGGCCCGGGCGGTACGAATGTAGTCTTCACTCAATTGTTCCAGCATCGCTGAGCGGATCTGGCGAGCGACAATTCCCAGGTGAACGAACGCCAGCGTTAATGCCGGTAAGATGAGATGTTGCAACGCATTGAAGAAGACTTCACCG  >  NC_000913/1560059‑1560351

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