Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,836,906 C→T 5.7% S278S (AGC→AGT ynjB → putative ABC transporter periplasmic binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,836,9060CT 5.7% 398.0 / 8.2 121S278S (AGC→AGTynjBputative ABC transporter periplasmic binding protein
Reads supporting (aligned to +/- strand):  ref base C (51/63);  new base T (3/4);  total (54/67)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.96e-01

Only 100 of 126 total aligned reads displayed.

AATATCTTGATGTGCTGCATCCGTATTTGTGGCGCGAAGGAAAGGATTTCCCGCCTTCACCCGCGCGGATGGATGCTCTGCTGAAAGCCGGCACATTGCGCCTGTCGCTGACCTTTAACCCCGCGCATGCGCAGCAAAAAATCGCCAGCGGCGATTTGCCTGCAAGCAGTTACAGTTTTGGCTTTCGCGAGGGGATGATTGGCAACGTGCATTTCGTCACCATTCCTGCCAACGCGAATGCCAGTGCTGCGGCGAAGGTAGTTGCCAATTTCCTGCTCTCACCCGATGCGCAACTG  >  NC_000913/1836758‑1837053
                                                                                                                                                    |                                                                                                                                                   
aaTATCTTGATGTGCTGCATCCGTATTTGTGGCGCGAAGGAAAGGATTTCCCGCCTTCACCCGCGCGGATGGATGCTCTGCTGAAAGCCGGCACATTGCGCCTGTCGCTGACCTTTAACCCCGCGCATGCGCAGCAAAAAATCGCCAgc                                                                                                                                                     >  1:1032921/1‑149 (MQ=255)
aaTATCTTGATGTGCTGCATCCGTATTTGTGGCGCGAAGGAAAGGATTTCCCGCCTTCACCCGCGCGGATGGATGCTCTGCTGAAAGCCGGCACATTGCGCCTGTCGCTGACCTTTAACCCCGCGCATGCGCAGCAAAAAATCGCCAgc                                                                                                                                                     >  2:989143/1‑149 (MQ=255)
aaTATCTTGATGTGCTGCATCCGTATTTGTGGCGCGAAGGAAAGGATTTCCCGCCTTCACCCGCGCGGATGGATGCTCTGCTGAAAGCCGGCACATTGCGCCTGTCGCTGACCTTTAACCCCGCGCATGCGCAGCAAAAAATCGCCAgc                                                                                                                                                     >  4:492485/1‑149 (MQ=255)
 atatCTTGATGTGCTGCATCCGTATTTGTGGCGCGAAGGAAAGGATTTCCCGCCTTCACCCGCGCGGATGGATGCTCTGCTGAAAGCCGGCACATTGCGCCTGTCGCTGACCTTTAACCCCGCGAATGCGCAGCAAAAAATCGCCAgcg                                                                                                                                                    >  2:719878/1‑149 (MQ=255)
     cTTGATGTGCTGCATCCGTATTTGTGGCGCGAAGGAAAGGATTTCCCGCCTTCACCCGCGCGGATGGATGCTCTGCTGAAAGCCGGCACATTGCGCCTGTCGCTGACCTTTAACCCCGCGCATGCGCAGCAAAAAATCGCCAGCGGCGa                                                                                                                                                <  2:910046/149‑1 (MQ=255)
     cTTGATGTGCTGCATCCGTATTTGTGGCGCGAAGGAAAGGATTTCCCGCCTTCACCCGCGCGGATGGATGCTCTGCTGAAAGCCGGCACATTGCGCCTGTCGCTGACCTTTAACCCCGCGCATGCGCAGCAAAAAATCGCCAGCGGCGa                                                                                                                                                >  3:688311/1‑149 (MQ=255)
        gATGTGCTGCATCCGTATTTGTGGCGCGAAGGAAAGGATTTCCCGCCTTCACCCGCGCGGATGGATGCTCTGCTGAAAGCCGGCACATTGCGCCTGTCGCTGACCTTTAACCCCGCGCATGCGCAGCAAAAAATCGCCGCTGGCGAttt                                                                                                                                             <  1:214723/149‑1 (MQ=255)
        gATGTGCTGCATCCGTATTTGTGGCGCGAAGGAAAGGATTTCCCGCCTTCACCCGCGCGGATGGATGCTCTGCTGAAAGCCGGCACATTGCGCCTGTCGCTGACCTTTAACCCCGCGCATGCGCAGCAAAAAATCGCCGCTGGCGAttt                                                                                                                                             <  2:645652/149‑1 (MQ=255)
        gATGTGCTGCATCCGTATTTGTGGCGCGAAGGAAAGGATTTCCCGCCTTCACCCGCGCGGATGGATGCTCTGCTGAAAGCCGGCACATTGCGCCTGTCGCTGACCTTTAACCCCGCGCATGCGCAGCAAAAAATCGCCGCTGGCGAttt                                                                                                                                             <  3:1005099/149‑1 (MQ=255)
              ctgcATCCGTATTTGTGGCGCGAAGGAAAGGATTTCCCGCCTTCACCCGCGCGGATGGATGCTCTGCTGAAAGCCGGCACATTGCGCCTGTCGCTGACCTTTAACCCCGCGCATGCGCAGCAAAAAATCGCCAGCGGCGATTTGCCTGc                                                                                                                                       >  1:975214/1‑149 (MQ=255)
                gcATCCGTATTTGTGGCGCGAAGGAAAGGATTTCCCGCCTTCACCCGCGCGGATGTATGCTCTGCTGAAAGCCGGCACATTGCGCGTGTCGCTGACCTTTAACCCCGCGCATGCGCAGCAAAAAATCGCCAGCGGCGATTTGCCTGCaa                                                                                                                                     <  2:233606/149‑1 (MQ=255)
                gcATCCGTATTTGTGGCGCGAAGGAAAGGATTTCCCGCCTTCACCCGCGCGGATGGATGCTCTGCTGAAAGCCGGCACATTGCGCCTGTCGCTGACCTTTAACCCCGCGCATGCGCAGCAAAAAATCGCCAGCGGCGATTTGCCTGCaa                                                                                                                                     >  3:773085/1‑149 (MQ=255)
                  aTCCGTACCTGTGGCTCGAAGGAAAGGATTTCCCGCCTTCACCCGCGCGGCTGGATGCTCTGCTGAAAGCCGGCACATTGCGCCTGTCGCTGACCTTTAACCCCGCGCATGCGCAGCAAAAAATCGCCAGCGGCGATTTGCCTGCAAGc                                                                                                                                   >  2:1006355/1‑149 (MQ=255)
                   tCCGTATTTGTGGCGCGAAGGAAAGGATTTCCCGCCTTCACCCGCGCGGATGGATGCTCTGCTGAAATCCGGCACATTGCGCCTGTCGCTGACCTTTAACCCCGCGCATGCGCAGCAAAAAATCGCCAGCGGCGATTTGCCTGCAAGCa                                                                                                                                  >  1:390385/1‑149 (MQ=255)
                    ccGGATTTGGGGCGCGAAGGAAAGGATTTCCCGCCTTCCCCCGCGCGGATGGATGCTCTGCTGAAAGCCGGCACATTGCGCCTGTCGCTGACCTTTAACCCCGCGCATGCGCAGCAAAAAATCGCCAGCGGCGATTTGCCTGCAAGCAg                                                                                                                                 <  1:1027557/149‑1 (MQ=255)
                     cGTATTTGTGGCGCGAAGGAAAGGATTTCCCGCCTTCACCCGCGCGGATGGATGCTCTGCTGAAAGCCGGCACATTGCGCCTGTCGCTGACCTTTAACCCCGCGCATGCGCAGCAAAAAATCGCCAGCGGCGATTTGCCTGCAAGCAGt                                                                                                                                <  3:894464/149‑1 (MQ=255)
                     cGTATTTGTGGCGCGAAGGAAAGGATTTCCCGCCTTCACCCGCGCGGATGGATGCTCTGCTGAAAGCCGGCACATTGCGCCTGTCGCTGACCTTTAACCCCGCGCATGCGCAGCAAAAAATCGCCAGCGGCGATTTGCCTGCAAGCAGt                                                                                                                                >  1:606643/1‑149 (MQ=255)
                              ggCGCGAAGGAAAGGATTTCCCGGCTTCACCGGCGGGGATAGATGCTCTGCAGAAAGCGCGCACATTACGCCCGTCGTGCACCTTGCACCCATCCCATGGGCAGCAAAAAATCGCCAGCGGCGGTTGGCATGCAAGTAGTTACAGgttt                                                                                                                       <  2:31419/149‑1 (MQ=255)
                              ggCGCGAAGGAAAGGATTTCCCGCCTTCACCCGCGCGGATGGATGCTCTGCTGAAAGCCGGCACATTGCGCCTGTCGCTGACCTTTAACCCCGCGCATTCGCAGCAAAAAATCGCCAGCGGCGATTTGCCTGCAAGCAGTTACAGtttt                                                                                                                       <  4:726033/149‑1 (MQ=255)
                              ggCGCGAAGGAAAGGATTTCCCGCCTTCACCCGCGCGGATGGATGCTCTGCTGAAAGCCGGCACATTGCGCCTGTCGCTGACCTTTAACCCCGCGCATGCGCAGCAAAAAATCGCCAGCGGCGATTTGCCTGCAAGCAGTTACAGtttt                                                                                                                       <  1:727527/149‑1 (MQ=255)
                              ggCGCGAAGGAAAGGATTTCCCGCCTTCACCCGCGCGGATGGATGCTCTGCTGAAAGCCGGCACATTGCGCCTGTCGCTGACCTTTAACCCCGCGCATGCGCAGCAAAAAATCGCCAGCGGCGATTTGCCTGCAAGCAGTTACAGtttt                                                                                                                       >  4:163566/1‑149 (MQ=255)
                                 gcgAAGGAAAGGATTTCCCGCCTTCACCCGCGCGGATGGATGCTCTGCTGAAAGCCGGCACATTGCGCCTGTCGCTGACCTTTAACCCCGCGCATGCGCAGCAAAAAATCGCCAGCGGCGATTTGCCTGCAAGCAGTTACAGTTTTGGc                                                                                                                    <  2:142622/149‑1 (MQ=255)
                                     aggaaaggaTTTCCCGCCTTCACCCGCGCGGATGGATGCTCTGCTGAAAGCCGGCACATTGCGCCTGTCGCTGACCTTTAACCCCGCGCATGCGCAGCAAAAAATCGCCAGCGGCGATTTGCCTGCAAGCAGTTACAGTTTTGGCTTTc                                                                                                                <  2:341149/149‑1 (MQ=255)
                                       gaaaggaTTTCCCGCCTTCACCCGCGCGGATGGATGCTCTGCTGAAAGCCGGCACATTGCGCCTGTCGCTGACCTTTAACCCCGCGCATGCGCAGCAAAAAATCGCCAGCGGCGATTTGCCt                                                                                                                                         <  2:687926/122‑1 (MQ=255)
                                       gaaaggaTTTCCCGCCTTCACCCGCGCGGATGGATGCTCTGCTGAAAGCCGGCACATTGCGCCTGTCGCTGACCTTTAACCCCGCGCATGCGCAGCAAAAAATCGCCAGCGGCGATTTGCCTGCAAGCAGTTACAGTTTTGGCTTTcgc                                                                                                              >  1:903297/1‑149 (MQ=255)
                                        aaaggaTTTCCCGCCTTCACCCGCGCGGATGGATGCTCTGCTGAAAGCCGGCACATTGCGCCTGTCGCTGACCTTTAACCCCGCGCATGCGCAGCAAAAAATCGCCAGCGGCGATTTGCCTGCAAGCAGTTACAGTTTTGGCTTTcgcg                                                                                                             >  4:531508/1‑149 (MQ=255)
                                         aaggaTTTCCCGCCTTCACCCGCGCGGATGGATGCTCTGCTGAAAGCCGGCACATTGCGCCTGTCGCTGACCTTTAACCCCGCGCATGCGCAGCAAAAAATCGCCAGCGGCGATTTGCCTGCAAGCAGTTa                                                                                                                              >  1:177739/1‑131 (MQ=255)
                                         aaggaTTTCCCGCCTTCACCCGCGCGGATGGATGCTCTGCTGAAAGCCGGCACATTGCGCCTGTCGCTGACCTTTAACCCCGCGCATGCGCAGCAAAAAATCGCCAGCGGCGATTTGCCTGCAAGCAGTTa                                                                                                                              <  2:177739/131‑1 (MQ=255)
                                         aaggaTTTCCCGCCTTCACCCGCGCGGATGGATGCTCTGCTGAAAGCCGGCACATTGCGCCTGTCGCTGACCTTTAACCCCGCGCATGCGCAGCAAAAAATCGCCAGCGGCGATTTGCCTGCAAGCAGTTACAGTTTTGGCTTTCGCGa                                                                                                            <  2:140046/149‑1 (MQ=255)
                                         aaggaTTTCCCGCCTTCACCCGCGCGGATGGATGCTCTGCTGAAAGCCGGCACATTGCGCCTGTCGCTGACCTTTAACCCCGCGCATGCGCAGCAAAAAATCGCCAGCGGCGATTTGCCTGCAAGCAGTTACAGTTTTGGCTTTCGCGa                                                                                                            >  2:296118/1‑149 (MQ=255)
                                          aggaTTTCCCGCCTTCACCCGCGCGGATGGATGCTCTGCTGAAAGCCGGCACATTGCGCCTGTCGCTGACCTTTAACCCCGCGCATGCGCAGCAAAAAATCGCCAGCGGCGATTTGCCTGCAAGCAGTTACAGTTTTGGCTTTCGCGAg                                                                                                           <  1:989143/149‑1 (MQ=255)
                                          aggaTTTCCCGCCTTCACCCGCGCGGATGGATGCTCTGCTGAAAGCCGGCACATTGCGCCTGTCGCTGACCTTTAACCCCGCGCATGCGCAGCAAAAAATCGCCAGCGGCGATTTGCCTGCAAGCAGTTACAGTTTTGGCTTTCGCGAg                                                                                                           <  2:678099/149‑1 (MQ=255)
                                            gaTTTCCCGCCTTCACCCGCGCGGATGGATGCTCTGCTGAAAGCCGGCACATTGCGCCTGTCGCTGACCTTTAACCCCGCGCATGCGCAGCAAAAAATCGCCAGCGGCGATTTGCCTGCaa                                                                                                                                     >  3:60339/1‑121 (MQ=255)
                                            gaTTTCCCGCCTTCACCCGCGCGGATGGATGCTCTGCTGAAAGCCGGCACATTGCGCCTGTCGCTGACCTTTAACCCCGCGCATGCGCAGCAAAAAATCGCCAGCGGCGATTTGCCTGCaa                                                                                                                                     <  4:60339/121‑1 (MQ=255)
                                            gaTTTCCCGCCTTCACCCGCGCGGATGGATGCTCTGCTGAAAGCCGGCACATTGCGCCTGTCGCTGACCTTTAACCCCGCGCATGCGCAGCAAAAAATCGCCAGCGGCGATTTGCCTGCAAGCAGTTACAGTTTTGGCTTTCGCGAggg                                                                                                         <  3:243417/149‑1 (MQ=255)
                                            gaTTTCCCGCCTTCACCCGCGCGGATGGATGCTCTGCTGAAAGCCGGCACATTGCGCCTGTCGCTGACCTTTAACCCCGCGCATGCGCAGCAAAAAATCGCCAGCGGCGATTTGCCTGCAAGCAGTTACAGTTTTGGCTTTCGCGAggg                                                                                                         <  3:158719/149‑1 (MQ=255)
                                              tttCCCGCCTTCACCCGCGCGGATGGATGCTCTGCTGAAAGCCGGCACATTGCGCCTGTCGCTGACCTTTAACCCCGCGCATGCGCAGCAAAAAATCGCCAGCGGCGATTTGCCTGCAAGCAGTTACAGTTTTGGCTTTCGCGAGGGga                                                                                                       >  2:872258/1‑149 (MQ=255)
                                                 cccGCCTTCACCCGCGCGGATGGATGCTCTGCTGAAAGCCGGCACATTGCGCCTGTCGCTGACCTTTAACCCCGCGCATGCGCAGCAAAAAATCGCCAGCGGCGATTTGCCTGCAAGCAGTTACAGTTTTGGCTTTCGCGAGGGgatga                                                                                                    >  2:120824/1‑149 (MQ=255)
                                                 cccGCCTTCACCCGCGCGGATGGATGCTCTGCTGAAAGCCGGCACATTGCGCCTGTCGCTGACCTTTAACCCCGCGCATGCGCAGCAAAAAATCGCCAGCGGCGATTTGCCTGCAAGCAGTTACAGTTTTGGCTTTCGCGAGGGgatga                                                                                                    <  2:975214/149‑1 (MQ=255)
                                                  ccGCCTTCACCCGCGCGGATGGATGCTCTGCTGAAAGCCGGCACATTGCGCCTGTCGCTGACCTTTAACCCCGCGCATGCGCAGCAAAAAATCGCCAGCGGCGATTTGCCTGCAAGCAGTTACAGTTTTGGCTTTCGCGAGGGgatgat                                                                                                   >  3:667147/1‑149 (MQ=255)
                                                   cGCCTTCACCCGCGCGGATGGATGCTCTGCTGAAAGCCGGCACATTGCGCCTGTCGCTGACCTTTAACCCCGCGCATGCGCAGCAAAAAATCGCCAGCGGCGATTTGCCTGCAAGCAGTTACAGTTTTGGCTTTCGCGAGGGGATGAtt                                                                                                  <  1:872258/149‑1 (MQ=255)
                                                   cGCCTTCACCCGCGCGGATGGATGCTCTGCTGAAAGCCGGCACATTGCGCCTGTCGCTGACCTTTAACCCCGCGCATGCGCAGCAAAAAATCGCCAGCGGCGATTTGCCTGCAAGCAGTTACAGTTTTGGCTTTCGCGAGGGGATGAtt                                                                                                  <  3:73182/149‑1 (MQ=255)
                                                   cGCCTTCACCCGCGCGGATGGATGCTCTGCTGAAAGCCGGCACATTGCGCCTGTCGCTGACCTTTAACCCCGCGCATGCGCAGCAAAAAATCGCCAGCGGCGATTTGCCTGCAAGCAGTTACAGTTTTGGCTTTCGCGAGGGGATGAtt                                                                                                  <  3:785443/149‑1 (MQ=255)
                                                    gCCTTCACCCGCGCGGATGGATGCTCTGCTGAAAGCCGGCACATTGCGCCTGTCGCTGACCTTTAACCCCGCGCATGCGCAGCAAAAAATCGCCAGCGGCGATTTGCCTGCAAGCAGTTACAGTTTTGGCTTTCGCGAGGGGATGATTg                                                                                                 <  4:571923/149‑1 (MQ=255)
                                                    gCCTTCACCCGCGCGGATGGATGCTCTGCTGAAAGCCGGCACATTCCGCCTGTCGCTGACCTTTAGCCCCGCGCATGCGCAGCAAAAAATCGCCAGCGGCGATTTGCCTGCAAGCAGTTACAGTTTTGGCTTTCGCGAGGGGATGATTg                                                                                                 <  4:207397/149‑1 (MQ=255)
                                                     ccgtCACCCGCGCGGATGGATGCTCTGCTGAAAGCCGGCACATTGCGCCTGTCGCTGACCTTTAACCCCGCGCATGCGCAGCAAAAAATCGCCAGCGGCGATTTGCCTGCAAGCAGTTACAGTTTTGGCTTTCGCGAGGGGATGATTgg                                                                                                <  1:815056/146‑1 (MQ=255)
                                                     ccTTCACCCGCGCGGATGGATGCTCTGCTGAAAGCCGGCACATTGCGCCTGTCGCTGACCTTTAACCCCGCGCATGCGCAGCAAAAAATCGCCAGCGGCGATTTGCCTGCAAGCAGTTACAGTTTTGGCTTTCGCGAGGGGATGATTgg                                                                                                <  2:76852/149‑1 (MQ=255)
                                                     ccTTCACCCGCGCGGATGGATGCTCTGCTGAAAGCCGGCACATTGCGCCTGTCGCTGACCTTTAACCCCGCGCATGCGCAGCAAAAAATCGCCAGCGGCGATTTGCCTGCAAGCAGATACAGTTTTGGCTTTCGCGAGGGGATGATTgg                                                                                                <  4:924963/149‑1 (MQ=255)
                                                      cTTCACCCGCGCGGATGGATGCTCTGCTGAAAGCCGGCACATTGCGCCTGTCGCTGACCTTTAACCCCGCGCATGCGCAGCAAAAAATCGCCAGCGGCGATTTGCCTGCAAGCAGTTACAGTTTTGGCTTTCGCGAGGGGATGATTGGc                                                                                               >  2:1051173/1‑149 (MQ=255)
                                                       ttCACCCGCGCGGATGGATGCTCTGCTGAAAGCCGGCACATTGCGCCTGTCGCTGACCTTTAACCCCGCGCATGCGCAGCAAAAAATCGCCAGCGGCGATTTGCCTGCAAGCAGTTACAGTTTTGGCTTTCGCGAGGGGATGATTGGCa                                                                                              <  1:381956/149‑1 (MQ=255)
                                                                  ggatggatGCTCTGCTGAAAGCCGGCACATTGCGCCTGTCGCTGACCTTTAACCCCGCGCATGCGCAGCAAAAAATCGCCAGCGGCGATTTGCCTGCAAGCAGTTACAGTTTTGGCTTTCGCGAGGGGATGATTGGCAACGTGCATTTc                                                                                   <  2:606643/149‑1 (MQ=255)
                                                                   gatggatgCTCTGCTGAAAGCCGGCACATTGCGCCTGTCGCTGACCTTTAACCCCGCGCATGCGCAGCAAAAAATCGCCAGCGGCGATTTGCCTGCAAGCAGTTACAGTTTTGGCTTTCGCGAGGGGATGATTGGCAACGTGCATTTCg                                                                                  <  3:492485/149‑1 (MQ=255)
                                                                            tctGCTGAAAGCCGGCACATTGCGCCTGTCGCTGACCTTTAACCCCGCGCATGCGCAGCAAAAAATCGCCAGCGGCGATTTGCCTGCAAGCAGTTACAGTTTTGGCTTTCGCGAGGGGATGATTGGCAACGTGCATTTCGTCACCATTc                                                                         <  4:62103/149‑1 (MQ=255)
                                                                            tctGCTGAAAGCCGGCACATTGCGCCTGTCGCTGACCTTTAACCCCGCGCATGCGCAGCAAAAAATCGCCAGCGGCGATTTGCCTGCAAGCAGTTACAGTTTTGGCTTTCGCGAGGGGATGATTGGCAACGTGCATTTCGTCACCATTc                                                                         >  3:411871/1‑149 (MQ=255)
                                                                            tctGCTGAAAGCCGGCACATTGCGCCTGTCGCTGACCTTTAACCCCGCGCATGCGCAGCAAAAAATCGCCAGCGGCGATTTGCCTGCAAGCAGTTACAGTTTTGGCTTTCGCGAGGGGATGATTGGCAACGTGCATTTCGTCACCATTc                                                                         >  4:524148/1‑149 (MQ=255)
                                                                              tgctgAAAGCCGGCACATTGCGCCTGTCGCTGACCTTTAACCCCGCGCATGCGCAGCAAAAAATCGCCAGCGGCGATTTGCCTGCAAGCAGTTACAGTTTTGGCTTTCGCGAGGGGATGATTGGCAACGTGCATTTCGTCACCATTCCt                                                                       <  4:309274/149‑1 (MQ=255)
                                                                              tgctgAAAGCCGGCACATTGCGCCTGTCGCTGACCTTTAACCCCGCGCATGCGCAGCAAAAAATCGCCAGCGGCGATTTGCCTGCAAGCAGTTACAGTTTTGGCTTTCGCGAGGGGATGATTGGCAACGTGCATTTCGTCACCATTCCt                                                                       <  3:163566/149‑1 (MQ=255)
                                                                                ctgAAAGCCGGCACATTGCGCCTGTCGCTGACCTTTAACCCCGCGCATGCGCAGCAAAAAATCGCCAGCGGCGATTTGCCTGCAAGCAGTTACAGTTTTGGCTTTCGCGAGGGGATGATTGGCAACGTGCATTTCGTCACCATTCCTGc                                                                     >  3:629184/1‑149 (MQ=255)
                                                                                ctgAAAGCCGGCACATTGCGCCTGTCGCTGACCTTTAACCCCGCGCATGCGCAGCAAAAAATCGCCAGCGGCGATTTGCCTGCAAGCAGTTACAGTTTTGGCTTTCGCGAGGGGATGATTGGCAACGTGCATTTCGTCACCATTCCTGc                                                                     >  1:777723/1‑149 (MQ=255)
                                                                                 tgAAAGCCGGCACATTGCGCCTGTCGCTGACCTTTAACCCCGCGCATGCGCAGCAAAAAATCGCCAGCGGCGATTTGCCTGCAAGCAGTTACAGTTTTGGCTTTCGCGAGGGGATGATTGGCAACGTGCATTTCGTCACCATTCCTGcc                                                                    <  3:816407/149‑1 (MQ=255)
                                                                                  gAAAGCCGGCACATTGCGCCTGTCGCTGACCTTTAACCCCGCGCATGCGCAGCAAAAAATCGCCAGCGGCGATTTGCCTGCAAGCAGTTACAGTTTTGGCTTTCGCGAGGGGATGATTGGCAACGTGCATTTCGTCACCATTCCTGCCa                                                                   >  4:535356/1‑149 (MQ=255)
                                                                                      gCCGGCACATTGCGCCTGTCGCTGACCTTTAACCCCGCGCATGCGCAGCAAAAAATCGCCAGCGGCGATTTGCCTGCAAGCAGTTACAGTTTTGGCTTTCGCGAGGGGATGATTGGCAACGTGCATTTCGTCACCATTCCTGCCAAcgc                                                               <  1:296118/149‑1 (MQ=255)
                                                                                        cGGCACATTGCGCCTGTCGCTGACCTTTAACCCCGCGCATGCGCAGCAAAAAATCGCCAGCGGCGATTTGCCTGCAAGCAGTTACAGTTTTGGCTTTCGCGAGGGGATGATTGGCAACGTGCATTTCGTCACCATTCCTGCCAACGCGa                                                             >  3:681487/1‑149 (MQ=255)
                                                                                         ggCACATTGCGCCTGTCGCTGACCTTTAACCCCGCGCATGCGCAGCAAAAAATCGCCAGCGGCGATTTGCCTGCAAGCAGTTACAGTTTTGGCTTTCGCGAGGGGATGATTGGCAACGTGCATTTCGTCACCATTCCTGCCAACGCGaa                                                            <  3:40482/149‑1 (MQ=255)
                                                                                          gCACATTGCGCCTGTCGCTGACCTTTAACCCCGCGCATGCGCAGCAAAAAATCGCCAGCGGCGATTTGCCTGCAAGCAGTTACAGTTTTGGCTTTCGCGAGGGGATGATTGGCAACGTGCATTTCGTCACCATTCCTGCCAACGCGAAt                                                           <  1:63106/149‑1 (MQ=255)
                                                                                          gCACATTGCGCCTGTCGCTGACCTTTAACCCCGCGCATGCGCAGCAAAAAATCGCCAGCGGCGATTTGCCTGCAAGCAGTTACAGTTTTGGCTTTCGCGAGGGGATGATTGGCAACGTGCATTTCGTCACCATTCCTGCCAACGCGAAt                                                           <  1:120824/149‑1 (MQ=255)
                                                                                          gCACATTGCGCCTGTCGCTGACCTTTAACCCCGCGCATGCGCAGCAAAAAATCGCCAGCGGCGATTTGCCTGCAAGCAGTTACAGTTTTGGCTTTCGCGAGGGGATGATTGGCAACGTGCATTTCGTCACCATTCCTGCCAACGCGAAt                                                           <  4:773085/149‑1 (MQ=255)
                                                                                              aTTGCGTCTGTCGCTGACCTTTAACCCCGCGCATGCGCAGCAAAAAATCGCCAGCGGCGATTTGCCTGCAAGCAGTTACAGTTTTGGCTTTCGCGAGGGGATGATTGGCAACGTGCATTTCGTCACCATTCCTGCCAACGCGAATGCCa                                                       >  2:779421/1‑149 (MQ=255)
                                                                                              aTTGCGCCTGTCGCTGACCTTTAACCCCGCGCATGCGCAGCAAAAAATCGCCAGCGGCGATTTGCCTGCAAGCAGTTACAGTTTTGGCTTTCGCGAGGGGATGATTGGCAACGTGCATTTCGTCACCATTCCTGCCAACGCGAATGCCa                                                       <  4:667147/149‑1 (MQ=255)
                                                                                              aTTGCGCCTGTCGCTGACCTTTAACCCCGCGCATGCGCAGCAAAAAATCGCCAGCGGCGATTTGCCTGCAAGCAGTTACAGTTTTGGCTTTCGCGAGGGGATGATTGGCAACGTGCATTTCGTCACCATTCCTGCCAACGCGAATGCCa                                                       <  1:468270/149‑1 (MQ=255)
                                                                                               ttGCGCCTGTCGCTGACCTTTAACCCCGCGCATGCGCAGCAAAAAATCGCCAGCGGCGATTTGCCTGCAAGCAGTTACAGTTTTGGCTTTCGCGAGGGGATGATTGGCAACGTGCATTTCGTCACCATTCCTGCCAACGCGAATGCCAg                                                      >  4:79969/1‑149 (MQ=255)
                                                                                                   gcCTGTCGCTGACCTTTAACCCCGCGCATGCGCAGCAAAAAATCGCCAGCGGCGATTTGCCTGCAAGCAGTTACAGTTTTGGCTTTCGCGAGGGGATGATTGGCAACGTGCATTTCGTCACCATTCCTGCCAACGCGAATGCCAGtgct                                                  <  4:681487/149‑1 (MQ=255)
                                                                                                           cTGACCTTTAACCCCGCGCATGCGCAGCAAAAAATCGCCAGCGGCGATTTGCCTGCAAGCAGTTACa                                                                                                                            <  4:873001/67‑1 (MQ=255)
                                                                                                           cTGACCTTTAACCCCGCGCATGCGCAGCAAAAAATCGCCAGCGGCGATTTGCCTGCAAGCAGTTACa                                                                                                                            >  3:873001/1‑67 (MQ=255)
                                                                                                           cTGACCTTTAACCCCGCGCATGCGCAGCAAAAAATCGCCAGCGGCGATTTGCCTGCAAGCAGTTACAGTTTTGGCTTTCGCGAGGGGATGATTGGCAACGTGCATTTCGTCACCATTCCTGCCAACGCGAATGCCAGTGCTGCGGCGaa                                          <  1:126998/149‑1 (MQ=255)
                                                                                                           cTGACCTTTAACCCCGCGCATGCGCAGCAAAAAATCGCCAGCGGCGATTTGCCTGCAAGCAGTTACAGTTTTGGCTTTCGCGAGGGGATGATTGGCAACGTGCATTTCGTCACCATTCCTGCCAACGCGAATGCCAGTGCTGCGGCGaa                                          <  1:403068/149‑1 (MQ=255)
                                                                                                              aCCTTTAACCCCGCGCATGCGCAGCAAAAAATCGCCAGCGGCGATTTGCCTGCAAGCAGTTACAGTTTTGGCTTTCGCGAGGGGATGATTGGCAACGTGCATTTCGTCACCATTCCTGCCAACGCGAATGCCAGTGCTGCGGCGAAGGt                                       >  2:527115/1‑149 (MQ=255)
                                                                                                               ccTTTAACCCCGCGCATGCGCAGCAAAAAATCGCCAGCGGCGATTTGCCTGCAAGCAGTTACAGTTTTGGCTTTCGCGAGGGGATGATTGGCAACGTGCATTTCGTCACCATTCCTGCCAACGCGAATGCCAGTGCTGCGGCGAAGGTa                                      <  4:205757/149‑1 (MQ=255)
                                                                                                                cTTTAACCCCGCGCATGCGCAGCAAAAAATCGCCAGCGGCGATTTGCCTGCAAGCAGTTACAGTTTTGGCTTTCGCGAGGGGATGATTGGCAACGTGCATTTCGTCACCATTCCTGCCAACGCGAATGCCAGTGCTGCGGCGAAGGTAg                                     >  3:399787/1‑149 (MQ=255)
                                                                                                                 tttAACCCCGCGCATGCGCAGCAAAAAATCGCCAGCGGCGATTTGCCTGCAAGCAGTTACAGTTTTGGCTTTCGCGAGGGGATGATTGGCAACGTGCATTTCGTCACCATTCCTGCCAACGCGAATGCCAGTGCTGCGGCGAAGGTAGt                                    >  1:322442/1‑149 (MQ=255)
                                                                                                                  ttAACCCCGCGCATGCGCAGCAAAAAATCGCCAGCGGCGATTTGCCTGCAAGCAGTTACAGTTTTGGCTTTCGCGAGGGGATGATTGGCAACGTGCATTTCGTCACCATTCCTGCCAACGCGAATGCCAGTGCTGCGGCGAAGGTAGtt                                   <  3:637889/149‑1 (MQ=255)
                                                                                                                  ttAACCCCGCGCATGCGCAGCAAAAAATCGCCAGCGGCGATTTGCCTGCAAGCAGTTACAGTTTTGGCTTTCGCGAGGGGATGATTGGCAACGTGCATTTCGTCACCATTCCTGCCAACGCGAATGCCAGTGCTGCGGCGAAGGTAGtt                                   <  2:49160/149‑1 (MQ=255)
                                                                                                                       cccGCGCATGCGCAGCAAAAAATCGCCAGCGGCGATTTGCCTGCAAGCAGTTACAGTTTTGGCTTTCGCGAGGGGATGATTGGCAACGTGCATTTCGTCACCATTc                                                                         <  1:61857/106‑1 (MQ=255)
                                                                                                                       cccGCGCATGCGCAGCAAAAAATCGCCAGCGGCGATTTGCCTGCAAGCAGTTACAGTTTTGGCTTTCGCGAGGGGATGATTGGCAACGTGCATTTCGTCACCATTCCTGCCAACGCGAATGCCAGTGCTGCGGCGAAGGTAGTTGCCaa                              <  2:73659/149‑1 (MQ=255)
                                                                                                                         cgcgCATGCGCAGCAAAAAATCGCCAGCGGCGATTTGCCTGCAAGCAGTTACAGTTTTGGCTTTCGCGAGGGGATGATTGGCAACGTGCATTTCGTCACCATTCCTGCCAACTCGTATGCCATTGCTGCGGCGAAGGTAGTTGCCAAtt                            <  2:116775/149‑1 (MQ=255)
                                                                                                                            gcATGCGCAGCAAAAAATCGCCAGCGGCGATTTGCCTGCAAGCAGTTACAGTTTTGGCTTTCGCGAGGGGATGATTGGCAACGTGCATTTCGTCACCATTCCTGCCAACGCGAATGCCAGTGCTGCGGCGAAGGTAGTTGCCAATTTcc                         <  2:961881/149‑1 (MQ=255)
                                                                                                                             cATGCGCAGCAAAAAATCGCCAGCGGCGATTTGCCTGCAAGCAGTTACAGTTTTGGCTTTc                                                                                                                >  1:887950/1‑61 (MQ=255)
                                                                                                                             cATGCGCAGCAAAAAATCGCCAGCGGCGATTTGCCTGCAAGCAGTTACAGTTTTGGCTTTc                                                                                                                <  2:887950/61‑1 (MQ=255)
                                                                                                                               tGCGCAGCAAAAAATCGCCAGCGGCGATTTGCCTGCAAGCAGTTACAGTTTTGGCTTTCGCGAGGGGATGATTGGCAACGTGCATTTCGTCACCATTCCTGCCAACGCGAATGCCAGTGCTGCGGCGAAGGTAGTTGCCAATTTCCTGc                      <  4:189742/149‑1 (MQ=255)
                                                                                                                               tGCGCAGCAAAAAATCGCCAGCGGCGATTTGCCTGCAAGCAGTTACAGTTTTGGCTTTCGCGAGGGGATGATTGGCAACGTGCATTTCGTCACCATTCCTGCCAACGCGAATGCCAGTGCTGCGGCGAAGGTAGTTGCCAATTTCCTGc                      >  2:187219/1‑149 (MQ=255)
                                                                                                                                  gcagcaAAAAATCGCCAGCGGCGATTTGCCTGCAAGCAGTTACAGTTTTGGCTTTCGCGAGGGGATGATTGGCAACGTGCATTTCGTCACCATTCCTGCCAACGCGAATGCCAGTGCTGCGGCGAAGGTAGTTGCCAATTTCCTGctct                   >  1:842578/1‑149 (MQ=255)
                                                                                                                                        aaaaaTCGCCAGCGGCGATTTGCCTGCAAGCAGTTACAGTTTTGGCTTTCGCGAGGGGATGATTGGCAACGTGCATTTCGTCACCATTCCTGCCAACGCGAATGCCAGTGCTGCGGCGAAGGTAGTTGCCAATTTCCTGCTCTCACCCg             >  2:592288/1‑149 (MQ=255)
                                                                                                                                        aaaaaTCGCCAGCGGCGATTTGCCTGCAAGCAGTTACAGTTTTGGCTTTCGCGAGGGGATGATTGGCAACGTGCATTTCGTCACCATTCCTGCCAACGCGAATGCCAGTGCTGCGGCGAAGGTAGTTGCCAATTTCCTGCTCTCACCCg             >  3:128353/1‑149 (MQ=255)
                                                                                                                                          aaaTCGCCGCTGGCGATTTGCCTGCAAGCAGTTACAGTTTTGGCTTTCGCGAGGGGATGATTGGCAACGTGCATTTCGTCACCATTCCTGCCAACGCGAATGCCAGTGCTGCGGCGAAGGTAGTTGCCAATTTCCTGCTCTCACCCGAt           >  1:922485/1‑149 (MQ=255)
                                                                                                                                          aaaTCGCCAGCGGCGATTTGCCTGCAAGCAGTTACAGTTTTGGCTTTCGCGAGGGGATGATTGGCAACGTGCATTTCGTCACCATTCCTGCCAACGCGAATGCCAGTGCTGCGGCGAAGGTAGTTGCCAATTTCCTGCTCTCACCCGAt           >  2:39340/1‑149 (MQ=255)
                                                                                                                                          aaaTCGCCAGCGGCGATTTGCCTGCAAGCAGTTACAGTTTTGGCTTTCGCGAGGGGATGATTGGCAACGTGCATTTCGTCACCATTCCTGCCAACGCGAATGCCAGTGCTGCGGCGAAGGTAGTTGCCAATTTCCTGCTCTCACCCGAt           <  3:952354/149‑1 (MQ=255)
                                                                                                                                            aTCGCCAGCGGCGATTTGCCTGCAAGCAGTTACAGTTTTGGCTTTCGCGAGGGGATGATTGGCAACGTGCATTTCGTCACCATTCCTGCCAACGCGAATGCCAGTGCTGCGGCGAAGGTAGTTGCCAATTTCCTGCTCTCACCCGATgc         >  1:618776/1‑149 (MQ=255)
                                                                                                                                             tCGCCAGCGGCGATTTGCCTGCAAGCAGTTATAGTTTTGGCTTTCGCGAGGGGATGATTGGCAACGTGCATTTCGTCACCATTCCTGCCAACGCGAATGCCAGTGCTGCGGCGAAGGTAGTTGCCAACTTCCTGCTCTCACCCGATgcg        >  3:901265/1‑149 (MQ=255)
                                                                                                                                                  aGCGGCGATTTGCCTGCAAGCAGTTACAGTTTTGGCTTTCGCGAGGGGATGATTGGCAACGTGCATTTCGTCACCATTCCTGCCAACGCGAATGCCAGGGCTGCGGCGAAGGTAGTTGCCAATTTCCTGCTCTCACCCGATGCGCAACt   >  2:666216/1‑149 (MQ=255)
                                                                                                                                                   gcggcgATTTGCCTGCAAGCAGTTACAGTTTCGGCTTTCGCGAGGGGATGATTGGCAACGTGCATTTCGTCACCATTCCTGCCAACGCGAATGCCAGTGCTGCGGCGAAGGTAGTTGCCAATTTCCTGCTCTCACCCGATGCGCAACTg  >  3:598423/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
AATATCTTGATGTGCTGCATCCGTATTTGTGGCGCGAAGGAAAGGATTTCCCGCCTTCACCCGCGCGGATGGATGCTCTGCTGAAAGCCGGCACATTGCGCCTGTCGCTGACCTTTAACCCCGCGCATGCGCAGCAAAAAATCGCCAGCGGCGATTTGCCTGCAAGCAGTTACAGTTTTGGCTTTCGCGAGGGGATGATTGGCAACGTGCATTTCGTCACCATTCCTGCCAACGCGAATGCCAGTGCTGCGGCGAAGGTAGTTGCCAATTTCCTGCTCTCACCCGATGCGCAACTG  >  NC_000913/1836758‑1837053

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: