Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,641,005 G→A S48S (TCG→TCA)  dtpB → dipeptide and tripeptide permease B

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,641,0050GA89.0% 434.6 / 2.7 139S48S (TCG→TCA) dtpBdipeptide and tripeptide permease B
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (42/80);  minor base C (15/0);  total (58/81)
Fisher's exact test for biased strand distribution p-value = 5.68e-07
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

Only 100 of 139 total aligned reads displayed.

AAATATGAATACAACAACACCCATGGGGATGCTGCAGCAACCTCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGC‑GG‑TTTTCTTCG‑TTAAACAGCTTGGATTCTCGCAAGAGCAGGC‑TTTTGTCACTTTTGGTGCTTTTGCTGCGCTGGTCTATGGCCTCATTTCCATTGGCGGCTATGTCGGC‑GACCACCTGCTGGGGACCAAACGCACCATTGTTCTTGGAGCACTTGTGCTGGCGATTGGCTACTTCAT  >  NC_000913/3640858‑3641151
                                                                                                                                                      |                                                                                                                                                    
aaataTGAATACAACAACACCCATGGGGATGCTGCAGCAACCTCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGC‑GG‑TTTTCTTCG‑TTAAACAGCTTGGATTCTCACaa                                                                                                                                                   <  2:553420/151‑1 (MQ=255)
  atatGAATACAACAACACCCATGGGGATGCTGCAGCAACCTCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGC‑GG‑TTTTCTTCG‑TTAAACAGCTTGGATTCTCACAaga                                                                                                                                                 <  4:722126/151‑1 (MQ=255)
   tatGAATACAACAACACCCATGGGGATGCTGCAGCAACCTCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGC‑GG‑GTTTTTTTG‑TTAAAAAGCTTGGGTTTTCACaa                                                                                                                                                   >  3:437416/1‑148 (MQ=255)
    atGAATACAACAACACCCATGGGGATGCTGCAGCAACCTCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGC‑GG‑TTTTCTTCG‑TTAAACAGCTTGGATTCTCACAagag                                                                                                                                                <  2:1673402/150‑1 (MQ=255)
    atGAATACAACAACACCCATGGGGATGCTGCAGCAACCTCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGC‑GG‑TTTTCTTCG‑TTAAACAGCTTGGATTCTCACAagag                                                                                                                                                <  4:480501/150‑1 (MQ=255)
      gAATACAACAACACCCATGGGGATGCTGCAGCAACCTCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGC‑GG‑TTTTCTTCG‑TTAAACAGCTTGGATTCTCACAAGAGCa                                                                                                                                              <  4:1382050/150‑1 (MQ=255)
      gAATACAACAACACCCATGGGGATGCTGCAGCAACCTCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGC‑GG‑TTTTCTTCG‑TTAAACAGCTTGGATTCTCACAAGAGCAg                                                                                                                                             <  3:336541/151‑1 (MQ=255)
       aaTACAACAACACCCATGGGGATGCTGCAGCAACCTCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGCGGG‑TTTTTTTCT‑TTAAA‑AGCTTTGGTTTTCCCAAAAGCAgg                                                                                                                                            >  4:693968/1‑151 (MQ=255)
       aaTACAACAACACCCATGGGGATGCTGCAGCAACCTCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGC‑GG‑TTTTCTTCG‑TTAAACAGCTTGGATTCTCACAAGAGCAg                                                                                                                                             <  2:223613/150‑1 (MQ=255)
          acaacaacaCCCATGGGGATGCTGCAGCAACCTCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGC‑GG‑TTTTCTTCG‑TTAAACAGCTTGGATTCTCACAAGAGCAGGC‑tt                                                                                                                                        <  2:1285398/151‑1 (MQ=255)
          acaacaacaCCCATGGGGATGCTGCAGCAACCTCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGC‑GG‑TTTTCTTCG‑TTAAACAGCTTGGATTCTCACAAGAGCAGGC‑tt                                                                                                                                        <  2:592045/151‑1 (MQ=255)
           caacaacaCCCATGGGGATGCTGCAGCAACCTCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGC‑GG‑GTTTTTTTG‑TTAAAAAAGTTGGGTTTTCACAAAAGCAggggtt                                                                                                                                        >  2:84824/1‑147 (MQ=255)
            aacaacaCCCATGGGGATGCTGCAGCAACCTCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGC‑GG‑TTTTTTTCT‑TTTAAAAACTTGGATTTTCCCAAAAGCAGGC‑tttt                                                                                                                                      >  2:610272/1‑151 (MQ=255)
            aacaacaCCCATGGGGATGCTGCAGCAACCTCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGC‑GG‑TTTTCTTCG‑TTAAACAGCTTGGATTCTCACAAGAGCAGGC‑tttt                                                                                                                                      <  4:721089/151‑1 (MQ=255)
              caacaCCCATGGGGATGCTGCAGCAACCTCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGC‑GG‑TTTTTTTTG‑TTAAAAAGCTTGGATTCTCACAAGAGCAGGC‑tttt                                                                                                                                      >  3:126401/1‑149 (MQ=255)
               aacaCCCATGGGGATGCTGCAGCAACCTCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGC‑GG‑GTTTTTTCG‑TTTAACAGCTTGGGTTTTCACAAAAGCAGGC‑gttttt                                                                                                                                    >  1:14674/1‑148 (MQ=255)
                acacCCATGGGGATGCTGCAGCAACCTCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGC‑GG‑TTTTTTTTG‑TTTAAAAACTTGGGTTCTCACAAGAGCAGGC‑tttttt                                                                                                                                    >  3:196733/1‑147 (MQ=255)
                    ccATGGGGATGCTGCAGCAACCTCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGC‑GG‑TTTTCTTCG‑TTAAACAGCTTGGATTCTCACAAGAGCAGGC‑TTTTGTCACttt                                                                                                                              <  4:13558/151‑1 (MQ=255)
                    ccATGGGGATGCTGCAGCAACCTCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGC‑GG‑TTTTCTTCG‑TTAAACAGCTTGGATTCTCACAAGAGCAGGC‑TTTTGTCACttt                                                                                                                              <  1:19399/151‑1 (MQ=255)
                     cATGGGGATGCTGCAGCAACCTCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGC‑GG‑TTTTCTTCG‑TTTAACAACTTGGGTTCTCCCAAAAGCAGGG‑CTTTTTCTCtttt                                                                                                                             >  2:927652/1‑151 (MQ=255)
                      aTGGGGATGCTGCAGCAACCTCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGC‑GG‑GTTTCTTCG‑TTTAACAGCCTTGATTCTCCCAAGAGGAGGG‑TTTTTTTACctttt                                                                                                                            >  3:342012/1‑150 (MQ=255)
                       tGGGGATGCTGCAGCAACCTCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGC‑GG‑TTTTTTTCC‑TTTAAAACCTTGGGTTTTCCCAAAAGCAGGC‑TTTTTTTCCtttt                                                                                                                             >  4:278148/1‑149 (MQ=255)
                        ggggATGCTGCAGCAACCTCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGC‑GG‑TTTTCTTCG‑TTAAACAGCTTGGATTCTCACAAGAGCAGGC‑TTTTGTCACTTTTgg                                                                                                                           <  2:1212720/150‑1 (MQ=255)
                         gggATGCTGCAGCAACCTCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGC‑GG‑TTTTCTTCG‑TTAAACAGCTTGGATTCTCACAAGAGCAGGC‑TTTTGTCACTTTTGGTg                                                                                                                         <  2:725684/151‑1 (MQ=255)
                           gATGCTGCAGCAACCTCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGC‑GG‑TTTTTTTCGTTTAAAAAGTTTGG‑TTTTCCCAAGAGGAGGC‑TTTTTTTACtttttggggt                                                                                                                       >  4:365572/1‑145 (MQ=255)
                            aTGCTGCAGCAACCTCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGC‑GG‑TTTTCTTCG‑TTAAACAGCTTGGATTCTCACAAGAGCAGGC‑TTTTGTCACTTTTGGTGCtt                                                                                                                      <  2:497397/151‑1 (MQ=255)
                                tgcAGCAACCTCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGC‑GG‑TTTTTTTCG‑TTTAACAGCTTGGGTTCTCACAAGAGCAGGG‑GTTTTTCACTTTTTGTGCtttt                                                                                                                    >  4:1004526/1‑149 (MQ=255)
                                 gcagcaACCTCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGC‑GG‑TTTTTTTTT‑TTTAAAAACTTTGAGTTTCACAAGAGAAGGG‑TTTTTTTTCTTTTTTTGTttttttt                                                                                                                 >  1:94604/1‑148 (MQ=255)
                                 gcagcaACCTCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGC‑GG‑TTTTCTTCG‑TTAAACAGCTTGGATTCTCACAAGAGCAGGC‑TTTTGTCACTTTTGGTGCTTTtg                                                                                                                   <  4:65890/149‑1 (MQ=255)
                                  cagcaACCTCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGC‑GGTTTTTCTTC‑‑TTAAACAAATTGGGGTTTCTCAAGAGGAGGC‑GTTTGTTACTTTTTGGGCtttttgt                                                                                                                 >  4:812790/1‑147 (MQ=255)
                                  cagcaACCTCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGC‑GG‑TTTTCTTCG‑TTAAACAGCTTGGATTCTCACAAGAGCAGGC‑TTTTGTCACTTTTGGTGCTTTtgctg                                                                                                                <  4:731513/151‑1 (MQ=255)
                                   agcaACCTCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGC‑GG‑TTTTCTTCG‑TTAAACAGCTTGGATTCTCACAAGAGCAGGC‑TTTTGTCACTTTTGGTGCTTTtgctgc                                                                                                               <  4:150227/151‑1 (MQ=255)
                                   agcaACCTCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGC‑GG‑TTTTCTTCG‑TTAAACAGCTTGGATTCTCACAAGAGCAGGC‑TTTTGTCACTTTTGGTGCTTTtgctgc                                                                                                               <  2:1092838/151‑1 (MQ=255)
                                    gcaACCTCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGC‑GG‑TTTTCTTCG‑TTAAACAGCTTGGATTCTCACAAGAGCAGGC‑TTTTGTCACTTTTGGTGCTTTTGCTgcg                                                                                                              <  2:1282859/151‑1 (MQ=255)
                                     caACCTCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGC‑GG‑TTTTCTTCG‑TTAAACAGCTTGGATTCTCACAAGAGCAGGC‑TTTTGTCACTTTTGGTGCTTTTGCTgcgc                                                                                                             <  4:963668/151‑1 (MQ=255)
                                      aaCCTCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGC‑GG‑TCTTCTTCG‑TTAAACAGCTTGGATTCTCACAAGAGCAGGC‑TTTTGTCACTTTTGGTGCTTTTGCTGCGCt                                                                                                            <  4:492260/151‑1 (MQ=255)
                                       aCCTCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGC‑GG‑GTTTTTTCG‑GTGAAAAACTTTGGTTTTTACAAAAAGAGGC‑CTTTTTTACTTTTTGGGCTTTtgct                                                                                                                 >  3:440887/1‑145 (MQ=255)
                                          tCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGC‑GG‑TTTTTTTTG‑TTTAAAAAATTTGGTTTTCACAAAAGGAGGG‑TTTTGTTTCTTTTTGGGCGTTTTTTGCGCTgg                                                                                                          >  2:1365244/1‑149 (MQ=255)
                                           cGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGC‑GG‑TTTTCTTCG‑TTAAACAGCTTGGATTCTCACAAGAGCAGGC‑TTTTGTCACTTTTGGTGCTTTTGCTGCGCTGGt                                                                                                         <  4:710381/149‑1 (MQ=255)
                                           cGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGC‑GG‑TTTTCTTCG‑TTAAACAGCTTGGATTCTCACAAGAGCAGGC‑TTTTGTCACTTTTGGTGCTTTTGCTGCGCTGGt                                                                                                         <  2:340166/149‑1 (MQ=255)
                                           cGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGC‑GG‑TTTTCTTCG‑TTAAACAGCTTGGATTCTCACAAGAGCAGGC‑TTTTGTCACTTTTGGTGCTTTTGCTGCGCTGGTc                                                                                                        <  2:622617/150‑1 (MQ=255)
                                             cccATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGC‑GG‑TTTTCTTCG‑TTAAACAGCTTGGATTCTCACAAGAGCAGGC‑TTTTGTCACTTTTGGTGCTTTTGCTGCGCTGGTCTAt                                                                                                     <  2:101183/151‑1 (MQ=255)
                                              ccATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGC‑GG‑TTTTTTTCG‑TTTAAAAACATGGGTTCTCACAAAAAGAAGG‑TTTTTTTTCTTTTTGTGTTTTTTCTGCGCTgggttat                                                                                                     >  2:1366624/1‑145 (MQ=255)
                                               cATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGC‑GG‑TTTTTTTTT‑TTTAACAGCTTGGGTTTTCCCAAGAGCAGGC‑TTTTTTCACTTTTTGTGTTTTTTCTGCGCTGGGCTTTgg                                                                                                   >  1:345555/1‑151 (MQ=255)
                                                    ttCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGC‑GG‑TTTTCTTCG‑TTAAACAGCTTGGATTCTCACAAGAGCAGGC‑TTTTGTCACTTTTGGTGCTTTTGCTGCGCTGGTCTATGGCCTc                                                                                               <  4:610842/150‑1 (MQ=255)
                                                     tCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGC‑GG‑TTTTCTTCG‑TTAAACAGCTTGGATTCTCACAAGAGCAGGC‑TTTTGTCACTTTTGGTGCTTTTGCTGCGCTGGTCTATGGCCTCa                                                                                              <  4:1409138/150‑1 (MQ=255)
                                                      cATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGC‑GG‑TTTTCTTCG‑TTAAACAGCTTGGATTCTCACAAGAGCAGGC‑TTTTGTCACTTTTGGTGCTTTTGCTGCGCTGGTCTATGGCCTCAt                                                                                             <  4:1558361/150‑1 (MQ=255)
                                                      cATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGC‑GG‑GTTTTTTTG‑TTTAACAGCTTTGGTTCTCACAAAAACAGGG‑TTTTTTGTATTTTTGTGCTTTTTTTGCGCTgggttttgggctcctt                                                                                            >  2:1369610/1‑137 (MQ=255)
                                                              tttttGTCGAGTTATGGGAGCTATTCGGCTACTACGGCGTTCAGGGCGTACTGGC‑GG‑TTTTCTTCG‑TTAAACAGCTTGTATTCTCTCAAGAGCAGGC‑TTTTGTCACTTTTGGTGCTTTTGCTGCGCTGGTCTATGGCCTCATTTCCAtt                                                                                      <  4:1207139/149‑1 (MQ=255)
                                                              tttttGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGC‑GG‑TTTTCTTCG‑TTAAACAGCTTGGATTCTCACAAGAGCAGGC‑TTTTGTCACTTTTGGTGCTTTTGCTGCGCTGGTCTATGGCCTCATTTCCAtt                                                                                      <  4:779562/149‑1 (MQ=255)
                                                              tttttGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGC‑GG‑TTTTCTTCG‑TTAAACAGCTTGGATTCTCACAAGAGCAGGC‑TTTTGTCACTTTTGGTGCTTTTGCTGCGCTGGTCTATGGCCTCATTTCCAt                                                                                       <  3:381585/148‑1 (MQ=255)
                                                                  tGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGC‑GG‑TTTTCTTTG‑TTTAAAAGCTTTGGTTCTCCCAAAAGCAGGCGTTTTTTCCCTTTTGG‑GCTTTTTCTGCGCTGGTCTATTGGCTCCTTTCCATTtgcgg                                                                                 >  4:1017990/1‑150 (MQ=255)
                                                                   gTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGC‑GG‑TTTTCTTCG‑TTAAACAGCTTGGATTCTCACAAGAGCAGGC‑TTTTGTCACTTTTGGTGCTTTTGCTGCGCTGGTCTATGGCCTCATTTCCATTggcggc                                                                                <  2:1527098/150‑1 (MQ=255)
                                                                      gAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGC‑GG‑CTTTTTTCG‑TTTAACAACTTGGGTTTTCACAAGAACCGGG‑GTTTTTTACTTTTTGTGATTTTTGTGCCCTGGGCTATGGCCCCCttttcctttgggg                                                                                 >  4:1105490/1‑136 (MQ=255)
                                                                       aGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGC‑GG‑TTTTCTTCG‑TTAAACAGCTTGGATTCTCACAAGAGCAGGC‑TTTTGTCACTTTTGGTGCTTTTGCTGCGCTGGTCTATGGCCTCATTTCCATTGGCGGCTATGt                                                                           <  4:1075229/151‑1 (MQ=255)
                                                                        gTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGC‑GG‑TTTTCTTCG‑TTAAACAGCTTGGATTCTCACAAGAGCAGGC‑TTTTGTCACTTTTGGTGCTTTTGCTGCGCTGGTCTATGGCCTCATTTCCATTGGCGGCTATGTc                                                                          <  2:1072901/151‑1 (MQ=255)
                                                                        gTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGC‑GG‑TTTTCTTCG‑TTAAACAGCTTGGATTCTCACAAGAGCAGGC‑TTTTGTCACTTTTGGTGCTTTTGCTGCGCTGGTCTATGGCCTCATTTCCATTGGCGGCTATGTc                                                                          <  3:55755/151‑1 (MQ=255)
                                                                          tATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGC‑GG‑TTTTCTTCG‑TTAAACAGCTTGGATTCTCACAAGAGCAGGC‑TTTTGTCACTTTTGGTGCTTTTGCTGCGCTGGTCTATGGCCTCATTTCCATTGGCGGCTATGTCg                                                                         <  4:1139775/150‑1 (MQ=255)
                                                                           aTGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGC‑GG‑TTTTTTTTG‑TTTAAAAACTTGGGTTTTCCCAAAAGCAGGG‑TTTTTTTACTTTTTGTGCTTTTTTTTCGCTGGTCTATTGCCTCCTTTTCATTtgcgggctttttgg                                                                        >  4:776516/1‑141 (MQ=255)
                                                                            tGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGC‑GG‑TTTTTTTTG‑TTTAACAGCTTGGAATCTCACAAAAGCAGGG‑TTTTTTCACTTTTTGTGGTTTTTCTGCGCTGGGCTTTGGCCTCCTTTTCATTggggggctttttggcg                                                                      >  3:453612/1‑140 (MQ=255)
                                                                             gggAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGC‑GG‑TTTTCTTCG‑TTAAACAGCTTGGATTCTCACAAGAGCAGGC‑TTTTGTCACTTTTGGTGCTTTTGCTGCGCTGGTCTATGGCCTCATTTCCATTGGCGGCTATGTCGGC‑g                                                                     <  2:1337490/150‑1 (MQ=255)
                                                                                            actacGGCGTGCAGGGCGTACTGGC‑GG‑TTTTCTTCG‑TTAAACAGCTTGGATTCTCACAAGAGCAGGC‑TTTTGTCACTTTTGGTGCTTTTGCTGCGCTGGTCTATGGCCTCATTTCCATTGGCGGCTATGTCGGC‑GACCACCTGCTgggg                                                       <  2:1353211/149‑1 (MQ=255)
                                                                                            actacGGCGTGCAGGGCGTACTGGC‑GG‑TTTTCTTCG‑TTAAACAGCTTGGATTCTCACAAGAGCAGGC‑TTTTGTCACTTTTGGTGCTTTTGCTGCGCTGGTCTATGGCCTCATTTCCATTGGCGGCTATGTCGGC‑GACCACCTGCTGGGGAc                                                     <  3:314543/151‑1 (MQ=255)
                                                                                               acGGCGTGCAGGGCGTACTGGC‑GG‑TTTTCTTCG‑TTAAACAGCTTGGATTCTCACAAGAGCAGGC‑TTTTGTCACTTTTGGTGCTTTTGCTGCGCTGGTCTATGGCCTCATTTCCATTGGCGGCTATGTCGGC‑GACCACCTGCTGGGGACCa                                                   <  2:496521/150‑1 (MQ=255)
                                                                                                 ggCGTGCAGGGCGTACTGGC‑GG‑TTTTTTTCG‑TTAAACAGCCTTGATTCTCACAAGAGCAGGC‑GTTTGTCACTTTTTGTGCTTTTTCTGCGCTGGGCTATGGCCTCATTTCCATTGGGGGGTATGTCGGC‑GACCACCCCCTGGGGGCCcaa                                                 >  3:467741/1‑150 (MQ=255)
                                                                                                       cAGGGCGTACTGGC‑GG‑TTTTTTTCG‑TTAAAAAGCTTGGATTCTCCCAAAAGCCAGC‑TTTTGTCACTTTTTGTGGTTTTTCTTCGCTGGTCTATGGCCCCCTTTTCCTTGGCGGGCTTTTCGGCGGCCCACCTCCGGGGGCCCaacccccc                                            >  2:1015515/1‑144 (MQ=255)
                                                                                                       cAGGGCGTACTGGC‑GG‑TTTTCTTCG‑TTAAACAGCTTGGATTCTCACAAGAGCAGGC‑TTTTGTCACTTTTGGTGCTTTTGCTGCGCTGGTCTATGGCCTCATTTCCATTGGCGGCGATGTCGGC‑GACCACCTGCTGGGGACCAAACGCACCa                                          <  1:213286/151‑1 (MQ=255)
                                                                                                           gCGTACTGGC‑GG‑TTTTCTTTG‑TTTAAAAGCATTGGTTCTCCCAAGAGCAGGC‑GTTTGTTACTTTTTGGGGTTTTGTTGCGCCGGGCTATGGCCCCATTTTCATTTGCGGGGATTTTGGC‑GGCCACCCGCTGGGGGGCAcacccacc                                           >  2:629689/1‑146 (MQ=255)
                                                                                                            cGTACTGGC‑GG‑TTTTCTTCG‑TTAAACAGCTTGGATTCTCACAAGAGCAGGC‑TTTTGTCACTTTTGGTGCTTTTGCTGCGCTGGTCTATGGCCTCATTTCCATTGGCGGCTATGTCGGC‑GACCACCTGCTGGGGACCAAACGCACCATTGtt                                     <  4:1560922/151‑1 (MQ=255)
                                                                                                            cGTACTGGC‑GG‑TTTTCTTCG‑TTAAACAGCTTGGATTCTCACAAGAGCAGGC‑TTTTGTCACTTTTGGTGCTTTTGCTGCGCTGGTCTATGGCCTCATTTCCATTGGCGGCTATGTCGGC‑GACCACCTGCTGGGGACCAAACGCACCATTGtt                                     <  2:1606566/151‑1 (MQ=255)
                                                                                                               aCTGGC‑GG‑TTTTCTTCG‑TTAAACAGCTTGGATTCTCACAAGAGCAGGC‑TTTTGTCACTTTTGGTGCTTTTGCTGCGCTGGTCTATGGCCTCATTTCCATTGGCGGCTATGTCGGC‑GACCACCTGCTGGGGACCAAACGCACCATTGTTCtt                                  <  4:1000527/151‑1 (MQ=255)
                                                                                                                cTGGC‑GG‑TTTTCTTCG‑TTAAACAGCTTGGATTCTCACAAGAGCAGGC‑TTTTGTCACTTTTGGTGCTTTTGCTGCGCTGGTCTATGGCCTCATTTCCATTGGCGGCTATGTCGGC‑GACCACCTGCTGGGGACCAAACGCACCATTGTTCtt                                  >  4:571993/1‑150 (MQ=255)
                                                                                                                 tGGC‑GG‑TTTTCTTCG‑TTAAACCGCTTGGATTCTCACCAGAGCAGGC‑TTTTGTCCCTTTTGGTGCTTTTGCTGCGCTGGTCTATGGCCTCATTTCCATTGGCGGGTCTGTCGGC‑GGCCCCCTTCTGGGGGCCCAACACACCaattgtctttg                                >  2:795929/1‑141 (MQ=255)
                                                                                                                  ggC‑GG‑TTTTCTTCG‑TTAAACAGCTTGGATTCTCACAAGAGCAGGC‑TTTTGTCACTTTTGGTGCTTTTGCTGCGCTGGTCTATGGCCTCATTTCCATTGGCGGCTATGTCGGC‑GACCACCTGCTGGGGACCAAATGCACCATTGTTCTTGGa                               <  4:864072/151‑1 (MQ=255)
                                                                                                                  ggC‑GG‑TTTTCTTCG‑TTAAACAGCTTGGATTCTCACAAGAGCAGGC‑TTTTGTCACTTTTGGTGCTTTTGCTGCGCTGGTCTATGGCCTCATTTCCATTGGCGGCTATGTCGGC‑GACCACCTGCTGGGGACCAAACGCACCATTGTTCTTg                                 <  3:41941/149‑1 (MQ=255)
                                                                                                                  ggC‑GG‑TTTTCTTCG‑TTAAACAGCTTGGATTCTCACAAGAGCAGGC‑TTTTGTCACTTTTGGTGCTTTTGCTGCGCTGGTCTATGGCCTCATTTCCATTGGCGGCTATGTCGGC‑GACCACCTGCTGGGGACCAAACGCACCATTGTTCTTGGa                               <  4:1186106/151‑1 (MQ=255)
                                                                                                                   gC‑GG‑TTTTCTTCG‑TTAAACAGCTTGGATTCTCACAAGAGCAGGC‑TTTTGTCACTTTTGGTGCTTTTGCTGCGCTGGTCTATGGCCTCATTTCCATTGGCGGCTATGTCGGC‑GACCACCTGCTGGGGACCAAACGCACCATTGTTCTTGGAg                              <  3:423440/151‑1 (MQ=255)
                                                                                                                    c‑GG‑TTTTCTTCG‑TTAAACAGCTTGGATTCTCACAAGAGCAGGC‑TTTTGTCACTTTTGGTGCTTTTGCTGCGCTGGTCTATGGCCTCATTTCCATTGGCGGCTATGTCGGC‑GACCACCTGCTGGGGACCAAACGCACCATTGTTCTTGGAGc                             >  4:1528412/1‑151 (MQ=255)
                                                                                                                       g‑TTTTCTTCG‑TTAAACAGCTTGGATTCTCACAAGAGCAGGC‑TTTTGTCACTTTTGGTGCTTTTGCTGCGCTGGTCTATGGCCTCATTTCCATTGGCGGCTATGTCGGC‑GACCACCTGCTGGGGACCAAACGCACCATTGTTCTTGGAGCAc                           <  4:460424/151‑1 (MQ=255)
                                                                                                                           ttcttaG‑TTAAACAGCTTGGATTCTCACAAGAGCAGGC‑TTTTGTCACTTTTGGTGCTTTTGCTGCGCTGGTCTATGGCCTCATTTCCATTGGCGGCTATGTCGGC‑GACCACCTGCTGGGGACCAAACGCACCATTGTTCTTGGAGCACTtg                        >  3:19291/1‑151 (MQ=255)
                                                                                                                                cG‑TTAAACAGCTTGGATTCTCACAAGAGCAGGC‑TTTTGTCACTTTTGGTGCTTTTGCTGCGCTGGTCTATGGCCTCATTTCCATTGGCGGCTATGTCGGC‑GACCACCTGCTGGGGACCAAACGCACCATTGTTCTTGGAGCACTTGTGCTg                   <  4:330366/151‑1 (MQ=255)
                                                                                                                                cG‑TTAAACAGCTTGGATTCTCACAAGAGCAGGC‑TTTTGTCACTTTTGGTGCTTTTGCTGCGCTGGTCTATGGCCTCATTTCCATTGGCGGCTATGTCGGC‑GACCACCTGCTGGGGACCAAACGCACCATTGTTCTTGGAGCACTTGTGCTg                   <  2:212477/151‑1 (MQ=255)
                                                                                                                                   ttAAACAGCTTGGATTCTCACAAGAGCAGGC‑TTTTGTCACTTTTGGTGCTTTTGCTGCGCTGGTCTATGGCCTCATTTCCATTGGCGGCTATGTCGGC‑GACCACCTGCTGGGGACCAAACGCACCATTGTTCTTGGAGCACTTGTGCt                    >  2:1418145/1‑148 (MQ=255)
                                                                                                                                    tAAACAGCTTGGATTCTCACAAGAGCAGGC‑TTTTGTCACTTTTGGTGCTTTTGCTGCGCTGGTCTATGGCCTCATTTCCATTGGCGGCTATGTCGGC‑GACCACCTGCTGGGGACCAAACGCACCATTGTTCTTGGAGCACTTGTGCTGGc                 <  4:1676134/150‑1 (MQ=255)
                                                                                                                                    tAAACAGCTTGGATTCTCACAAGAGCAGGC‑TTTTGTCACTTTTGGTGCTTTTGCTGCGCTGGTCTATGGCCTCATTTCCATTGGCGGCTATGTCGGC‑GACCACCTGCTGGGGACCAAACGCACCATTGTTCTTGGAGCACTTGTGCTGGCg                >  2:944432/1‑151 (MQ=255)
                                                                                                                                     aaaCAGCTTGGATTCTCACAAGAGCAGGC‑TTTTGTCACTTTTGGTGCTTTTGCTGCGCTGGTCTATGGCCTCATTTCCATTGGCGGCTATGTCGGC‑GACCACCTGCTGGGGACCAAACGCACCATTGTTCTTGGAGCACTTGTGCTGGCg                <  4:921605/150‑1 (MQ=255)
                                                                                                                                     aaaCAGCTTGGATTCTCACAAGAGCAGGC‑TTTTGTCACTTTTGGTGCTTTTGCTGCGCTGGTCTATGGCCTCATTTCCATTGGCGGCTATGTCGGC‑GACCACCTGCTGGGGACCAAACGCACCATTGTTCTTGGAGCACTTGTGCTGGCGa               <  3:321155/151‑1 (MQ=255)
                                                                                                                                     aaaCAGCTTGGATTCTCACAAGAGCAGGC‑TTTTGTCACTTTTGGTGCTTTTGCTGCGCTGGTCTATGGCCTCATTTCCATTGGCGGCTATGTCGGC‑CACCACCTGCTGGGGACCAAAAGCCCCACTGTTCTTGGAGCACTTGTGCTGGCg                >  2:709880/1‑150 (MQ=255)
                                                                                                                                       aCAGCTTGGATTCTCACAAGAGCAGGC‑TTTTGTCACTTTTGGTGCTTTTGCTGCGCTGGTCTATGGCCTCATTTCCATTGGCGGCTATGTCGGC‑GACCACCTGCTGGGGACCAAACGCACCATTGTTCTTGGAGCACTTGTGCTGGCGAtt             <  1:126387/151‑1 (MQ=255)
                                                                                                                                       aCAGCTTGGATTCTCACAAGAGCAGGC‑TTTTGTCACTTTTGGTGCTTTTGCTGCGCTGGTCTATGGCCTCATTTCCATTGGCGGCTATGTCGGC‑GACCACCTGCTGGGGACCAAACGCACCATTGTTCTTGGAGCACTTGTGCTGGCGAtt             >  4:1607554/1‑151 (MQ=255)
                                                                                                                                        cAGCTTGGATTCTCACAAGAGCAGGC‑TTTTGTCACTTTTGGTGCTTTTGCTGCGCTGGTCTATGGCCTCATTTCCATTGGCGGCTATGTCGGC‑GACCACCTGCTGGGGACCAAACGCACCATTGTTCTTGGAGCACTTGTGCTGGCGATTg            >  4:1299443/1‑151 (MQ=255)
                                                                                                                                         aGCTTGGATTCTCACAAGAGCAGGC‑TTTTGTCACTTTTGGTGCTTTTGCTGCGCTGGTCTATGGCCTCATTTCCATTGGCGGCTATGTCGGC‑GACCACCTGCTGGGGACCAAACGCACCATTGTTCTTGGAGCACTTGTGCTGGCGATTg            <  2:811713/150‑1 (MQ=255)
                                                                                                                                          gCTTGGATTCTCACAAGAGCAGGC‑TTTTGTCACTTTTGGTGCTTTTGCTGCGCTGGTCTATGGCCTCATTTCCATTGGCGGCTATGTCGGC‑GACCACCTGCTGGGGACCAAACGCACCATTGTTCTTGGAGCACTTGTGCTGGCGATTGGc          <  2:694749/151‑1 (MQ=255)
                                                                                                                                           cTTCGATTCTCACAAGAGCAGGC‑TTTTGTCACTTTTGGTGCTTTTGCTGCGCTGGTCTATGGCCTCATTTCCATTGGCGGCTATGTCGGC‑GACCACCTGCTGGGGACCAAACGCACCATTGTTCTTGGAGCACTTGTGCTGGCGATTgg           >  1:445521/1‑149 (MQ=255)
                                                                                                                                            ttGGATTCTCACAAGAGCAGGC‑TTTTGTCACTTTTGGTGCTTTTGCTGCGCTGGTCTATGGCCTCATTTCCATTGGCGGCTATGTCGGC‑GACCACCTGCTGGGGACCAAACGCACCATTGTTCTTGGAGCACTTGTGCTGGCGATTGGCt         >  1:228560/1‑150 (MQ=255)
                                                                                                                                            ttGGATTCTCACAAGAGCAGGC‑TTTTGTCACTTTTGGTGCTTTTGCTGCGCTGGTCTATGGCCTCATTTCCATTGGCGGCTATGTCGGC‑GACCACCTGCTGGGGACCAAACGCACCATTGTTCTTGGAGCACTTGTGCTGGCGATTGGCt         >  2:15764/1‑150 (MQ=255)
                                                                                                                                             tGGATTCTCACAAGAGCAGGC‑TTTTGTCACTTTTGGTGCTTTTGCTGCGCTGGTCTATGGCCTCATTTCCATTGGCGGCTATGTCGGC‑GACCACCTGCTGGGGACCAAACGCACCATTGTTCTTGGAGCACTTGTGCTGGCGATTGGCTa        <  4:619359/150‑1 (MQ=255)
                                                                                                                                             tGGATTCTCACAAGAGCAGGC‑TTTTGTCACTTTTGGTGCTTTTGCTGCGCTGGTCTATGGCCTCATTTCCATTGGCGGCTATGTCGGC‑GACCACCTGCTGGGGACCAAACGCACCATTGTTCTTGGAGCACTTGTGCTGGCGATTGGCTa        <  4:84980/150‑1 (MQ=255)
                                                                                                                                                 ttCTCACAAGAGCAGGC‑TTTTGTCACTTTTGGTGCTTTTGCTGCGCTGGTCTATGGCCTCATTTCCATTGGCGGCTATGTCGGC‑GACCACCTGCTGGGGACCAAACGCACCATTGTTCTTg                                 <  1:296262/121‑1 (MQ=255)
                                                                                                                                                  tctcACAAGAGCAGGC‑TTTTGTCACTTTTGGTGCTTTTGCTGCGCTGGTCTATGGCCTCATTTCCATTGGCGGCTATGTCGGC‑GACCACCTGCTGGGGACCAAACGCACCATTGTTCTTGGAGCACTTGTGCTGGCGATTGGCTACTTCAt  >  4:591167/1‑151 (MQ=255)
                                                                                                                                                      |                                                                                                                                                    
AAATATGAATACAACAACACCCATGGGGATGCTGCAGCAACCTCGCCCATTTTTCATGATCTTTTTTGTCGAGTTATGGGAGCGATTCGGCTACTACGGCGTGCAGGGCGTACTGGC‑GG‑TTTTCTTCG‑TTAAACAGCTTGGATTCTCGCAAGAGCAGGC‑TTTTGTCACTTTTGGTGCTTTTGCTGCGCTGGTCTATGGCCTCATTTCCATTGGCGGCTATGTCGGC‑GACCACCTGCTGGGGACCAAACGCACCATTGTTCTTGGAGCACTTGTGCTGGCGATTGGCTACTTCAT  >  NC_000913/3640858‑3641151

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