Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,836,905 G→C 5.1% S278T (AGC→ACC)  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,9050GC 5.1% 286.2 / 4.5 94S278T (AGC→ACC) ynjBputative ABC transporter periplasmic binding protein
Reads supporting (aligned to +/- strand):  ref base G (52/35);  new base C (3/2);  total (55/39)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.20e-01

GCAATATCTTGATGTGCTGCATCCGTATTTGTGGCGCGAAGGAAAGGATTTCCCGCCTTCACCCGCGCGGATGGATGCTCTGCTGAAAGCCGGCACATTGCGCCTGTCGCTGACCTTTAACCCCGCGCATGCGCAGCAAAAAATCGCCAGCGGCGATTTGCCTGCAAGCAGTTACAGTTTTGGCTTTCGCGAGGGGATGATTGGCAACGTGCATTTCGTCACCATTCCTGCCAACGCGAATGCCAGTGCTGCGGCGAAGGTAGTTGCCAATTTCCTGCTCTCACCCGATGCGCAA  >  NC_000913/1836756‑1837050
                                                                                                                                                     |                                                                                                                                                 
gCAATATCTTGATGTGCTGCATCCGTATTTGTGGCGCGAAGGAAAGGATTTCCCGCCTTCACCCGCGCGGATGGATGCTCTGCTGAAAGCCGGCACATTGCGCCTGTCGCTGACCTTTAACCCCGCGCATGCGCAGCAAAAAATCGCCAg                                                                                                                                                   >  2:585886/1‑150 (MQ=255)
  aaTATCTTGATGTGCTGCATCCGAATTTGTGGCGCGAAGGAAAGGATTTCCCGCCTTCAAACGCGCGGATGGATGCTCTGCTGAAAGCCGGCACATTGAGCCTGTCGCTGACCTTTAACCCCGCGCATGCGCAGCAAAAAATCGCCAgcg                                                                                                                                                 >  1:1253765/1‑150 (MQ=255)
   atatCTTGATGTGCTGCATCCGTATTTGTGGCGCGAAGGAAAGGATTTCCCGCCTTCACCCGCGCGGATGGATGCTCTGCTGAAAGCCGGCACATTGCGCCTGTCGCTGACCTTTAACCCCGCGCATGCGCAGCAAAAAATCGCCAgcgg                                                                                                                                                <  1:1291510/150‑1 (MQ=255)
      tCTTGATGTGCTGCATCCGTATTTGTGGCGCGAAGGAAAGGATTTCCCGCCTTCACCCGCGCGGATGGATGCTCTGCTGAAAGCCGGCACATTGCGCCTGTCGCTGACCTTTAACCCCGCGCATGCGCAGCAAAAAATCGCCAGCGGCGa                                                                                                                                             <  2:197266/150‑1 (MQ=255)
       cTTGATGTGCTGCATCCGTATTTGTGGCGCGAAGGAAAGGATTTCCCGCCTTCACCCGCGCGGATGGATGCTCTGCTGAACGCCGGCACATTGCGCCTGTCGCTGACCTTTAACCCCGCGCATGCGCAGCAAAAAATCGCCAGCGGCGAt                                                                                                                                            >  2:867410/1‑150 (MQ=255)
       cTTGATGTGCTGCATCCGTATTTGTGGCGCGAAGGAAAGGATTTCCCGCCTTCACCCGCGCGGATGGATGCTCTGCTGAAAGCCGGCACATTGCGCCTGTCGCTGACCTTTAACCCCGCGCATGCGCAGCAAAAAATCGCCAGCGGCGAt                                                                                                                                            >  1:1706678/1‑150 (MQ=255)
        ttGATGTGCTGCATCCGTATTTGTGGCGCGAAGGAAAGGATTTCCCGCCTTCACCCGCGCGGATGGATGCTCTGCTGAAAGCCGGCACATTGCGCCTGTCGCTGACCTTTAACCCCGCGCATGCGCAGCAAAAAATCGCCAGCGGCGAtt                                                                                                                                           <  1:1620897/150‑1 (MQ=255)
        ttGAGGTGCTGCATCCGTATTTGTGGCGCGAAGGAAAGGATTTCCCGCCTTCAACCGCGCGGATTGATGCTCTGCTGAAAGCCGGCACATTGCGCCTGTCGCTGAACTTTAACCCCGCGCATGCGCAGCAAAAAATCGCCATCGGCGAtt                                                                                                                                           <  2:1069777/150‑1 (MQ=255)
           aTGTGCTGCATCCGTATTTGTGGCGCGAAGGAAAGGATTTCCCGCCTTCACCCGCGCGGATGGATGCTCTGCTGAAAGCCGGCACATTGCGCCTGTCGCTGACCTTTAACCCCGCGCATGCGCAGCAAAAAATCGCCGCTGGCGAttttt                                                                                                                                        <  1:711524/150‑3 (MQ=255)
            tgtgCTGCATCCGTATTTGTGGCGCGAAGGAAAGGATTTCCCGCCTTCCCCCGCGCGGATGGATGCTCTGCTGAAAGCCGGCACATTGCGCCTGTCGCTGACCTTTAACCCCGCGCATGCGCAGCAAAAAATCGCCAGCGGCGATTTGcc                                                                                                                                       >  2:1274421/1‑150 (MQ=255)
             gtgCTGCATCCGTATTTGTGGCGCGAAGGATAGGATTTCCCGCCTTCACCCGCGCGGATGGATGCGCTGCTGAAAGCCGGCACATTGCGCCTGTCGCTGACCTTTAACCCCGCGCATGCGCAGCAAAAAATCGCCAGCGGCGATTTGCCt                                                                                                                                      <  1:28436/150‑1 (MQ=255)
              tgctgcATCCGTATTTGTGGCGCGAAGGAAAGGATTTCCCGCCTTCACCCGCGCGGATGGATGCTCTGCTGAAAGCCGGCACATTGCGCCTGTCGCTGACCTTTAACCCCGCGCATGCGCAGCAAAAAATCGCCAGCGGCGATTTGCCTg                                                                                                                                     <  1:1709854/150‑1 (MQ=255)
              ggctgcATCCGTATTTGTGGCGCGAAGGAAAGGATTTCCCGCCTTCACCCGCGCGGATGGATGCTCTGCTGAAAGCCGGCACATTGCGCCTGTCGCTGACCTTTAACCCCGCGCATGCGCAGCAAAAAATCGCCAGCGGCGATTTGCCTg                                                                                                                                     <  1:358861/149‑1 (MQ=255)
                  gcATCCGTATTTGTGGCGCGAAGGAAAGGATTTCCCGCCTTCACCCGCGCGGATGGATGCTCTGCTGAAAGCCGGCACATTGCGCCTGTCGCTGCCGTTTAACCCCGCGCATGCGCAGCAAAAAATCGCCAGCGGCGATTTGCCTGCAAg                                                                                                                                 >  2:824705/1‑150 (MQ=255)
                       cGTATTTGTGGCGCGAAGGAAAGGATTTCCCGCCTTCACCCGCGCGGATGGATGCTCTGCTGAAAGCCGGCACATTTCGCCTGTCGCTGACCTTTAACCACGAGCATGCGCATCAAAAAATAGCCAGCGGCGATTTGCCTCAAAGCAGtt                                                                                                                            >  2:126320/1‑150 (MQ=255)
                                ggCGCGAAGGAAAGGATTTCCCGCCTTCACCCGCGCGGATGGATGCTCTGCTGAAAGCCGGCACATTGCGCCTGTCGCTGACCTTTAACCCCGCGCATGCGCAGCAAAAAATCGCCAGCGGCGATTTGCCTGCAAGCAGTTACAGTTTTg                                                                                                                   <  1:1160200/150‑1 (MQ=255)
                                 gcgcGAAGGAAAGGATTTCCCGCCTTCACCCGCGCGGATGGATGCTCTGCTGAAAGCCGGCACATTGCGCCTGTCGCTGACCTTTAACCCCGCGCATGCGCAGCAAAAAATCGCCAGCGGCGATTTGCCTGCAAGCAGTTACAGTTTTgg                                                                                                                  <  2:657926/150‑1 (MQ=255)
                                  cgcgAAGGAAAGGATTTCCCGCCTTCACCCGCGCGGATGGATGCTCTGCTGAAAGCCGGCACATTGCGCCTGTCGCTGCCCTTTAACCCCGCGCATGCGCAGCAAAAAATCGCCAGCGGCGATTTGCCTGCAAg                                                                                                                                 >  1:135475/1‑134 (MQ=255)
                                  cgcgAAGGAAAGGATTTCCCGCCTTCACCCGCGCGGATGGATGCTCTGCTGAAAGCCGGCACATTGCGCCTGTCGCTGACCTTTAACCCCGCGCATGCGCAGCAAAAAATCGCCAGCGGCGATTTGCCTGCAAg                                                                                                                                 <  2:135475/134‑1 (MQ=255)
                                   gcgAAGGAAAGGATTTCCCGCCTTCACCCGCGCGGATGGATGCTCTGCTGAAAGCCGGCACATTGCGCCTGTCGCTGACCTTTAACCCCGCGCATGCGCAGCAAAAAATCGCCAGCGGCGATTTGCCTGCAAGCAGTTACAGTTTTGGCt                                                                                                                <  2:1475055/150‑1 (MQ=255)
                                    cgAAGGAAAGGATTTCCCGCCTTCACCCGCGCGGATGGATGCTCTGCTGAAAGCCGGCACATTGCGCCTGTCGCTGACCTTTAACCCCGCGCATGCGCAGCAAAAAATCGCCAGCGGCGATTTGCCTGCAAGCAGTTACAGTTTTGGCtt                                                                                                               >  1:908665/1‑150 (MQ=255)
                                    cgAAGGAAAGGATTTCCCGCCTTCACCCGCGCGGATGGATGCTCTGATGAAAGCCGGCACATTGCGCCTGTCGCTGACCTTTAACCCCGCGCATGCGCAGCAAAAAATCGCCAGCGGCGATTTGCCTGCAAGCAGTTACAGTTTTGGCtt                                                                                                               >  1:922435/1‑150 (MQ=255)
                                     gAAGGAAAGGATTTCCCGCCTTCACCCGCGCGGATGGATGCTCTGCTGAAAGCCGGCACATTGCGCCTGTCGCTGACCTTTAACCCCGCGCATGCGCAGCAAAAAATCGCCAGCGGCGATTTGCCTGCAAGCAGTTACAGTTTTGGCttt                                                                                                              >  2:1472431/1‑150 (MQ=255)
                                     gAAGGAAAGGATTTCCCGCCTTCACCCGCGCGGATGGATGCTCTGCTGAAAGCCGGCACATTGCGCCTGTCGCTGACCTTTAACCCCGCGCATGCGCAGCAAAAAATCGCCAGCGGCGATTTGCCTGCAAGCAGTTACAGTTTTGGCttt                                                                                                              >  2:1003369/1‑150 (MQ=255)
                                     gAAGGAAAGGATTTCCCGCCTTCACCCGCGCGGATGGATGCTCTGCTGAAAGCCGGCACATTGCGCCTGTCGCTGACCTTTAACCCCGCGCATGCGCAGCAAAAAATCGCCAGCGGCGATTTGCCTGCAAGCAGTTACAGTTTTGGCttt                                                                                                              >  2:1168268/1‑150 (MQ=255)
                                      aaggaaaggaTTTCCCGCCTTCACCCGCGCGGATGGATGCTCTGCTGAAAGCCGGCACATTGCGCCTGTCGCTGACCTTTAACCCCGCGCATGCGCAGCAAAAAATCGCCAGCGGCGATTTGCCTGCAAGCAGTTACAGTTTTGGCTTTc                                                                                                             <  2:1606464/150‑1 (MQ=255)
                                        ggaaaggaTTTCCCGCCTTCACCCGCGCGGATGGATGCTCTGCTGAAAGCCGGCACATTGCGCCTGTCGCTGACCTTTAACCCCGCGCATGCGCAGCAAAAAATCGCCAGCGGCGATTTGCCTGCAAGCAGTTACAGTTTTGGCTTTcgc                                                                                                           <  2:1202491/150‑1 (MQ=255)
                                         caaaggaTTACCCGACTACACCCGATCCGATCTATGCACTGCTGAAAGCCGGCACATTGCGAATGAAGCTCACCTTAAACCCCTCGCATGCGCAGCAAAAAATCGCAAGCGGCGATTTCCCTGCAAGCAGTAACAGTATAGGATTAcgcg                                                                                                          >  2:758858/2‑150 (MQ=255)
                                          aaaggaTTTCCCGCCTTCACCCGCGCGGATGGATGCTCTGCTGAAAGCCGGCACATTGCGCCTGTCGCTGACCTTTAACCCCGCGCATGCGCAGCAAAAAATCGCCAGCGGCGATTTGCCTGCAAGCAGTTACAGTTTTGGCTTTCGCGa                                                                                                         >  2:985304/1‑150 (MQ=255)
                                           aaggaTTTCCCGCCTTCACCCGCGCGGATGGATGCTCTGCTGAAAGCCGGCAGATTGCGCCTGTCGCTGACCTTTAACCCCGCGCATGCGCAGCAAAAAATCGCCAGCGGCGATTTGCCTGCAAGCAGTTACAGTTTTGGCTTTCGCGAg                                                                                                        <  1:167250/150‑1 (MQ=255)
                                             ggaTTTCCCGCCTTCACCCGCGCGGATGGATGCTCTGCTGAAAGCCGGCACATTGCGCCTGTCGCTGACCTTTAACCCCGCGCATGCGCAGCAAAAAATCGCCAGCGGCGATTTGCCTGCAAGCAGTTACAGTTTTGGCTTTCGCGAggg                                                                                                      <  2:695024/150‑1 (MQ=255)
                                                tttCCCGCCTTCACCCGCGCGGATGGATGCTCTGCTGAAAGCCGGCACATTGCGCCTGTCGCTGACCTTTAACCCCGCGCATGCGCAGCAAAAAATCGCCAGCGGCGATTTGCCTGCAAGCAGTTACAGTTTTGGCATTcgcgttggcag                                                                                                   >  2:1722766/1‑143 (MQ=255)
                                                   cccGCCTTCACCCGCGCGGATGGATGCTCTGCTGAAAGCCGGCACATTGCGCCTGTCGCTGACCTTTAACCCCGCGCATGCGCAGCAAAAAATCGCCCGCGGCGATTTGCCTGCAAGCAGTAACAGTTTTGGCTTTCGCGAGGGgatgat                                                                                                >  2:143181/1‑150 (MQ=255)
                                                    ccGCCTTCACCCGCGCGGATGGATGCTCTGCTGAAAGCCGGCACATTGCGCCTGTCGCTGACCTTTAACCCCGCGCATGCGCAGCAAAAAATCGCCAGCGGCGATTTGCCTGCAAGCAGTTACAGTTTTGGCTTTCGCGAGGGGATGAtt                                                                                               >  1:1519297/1‑150 (MQ=255)
                                                          tCACCCGCGCGGATGGATGCTCTGCTGAAAGCCGGCACATTGCGCCTGTCGCTGACCTTTAACCCCGCGCATGCGCAGCAAAAAATCGCCAGCGGCGATTTGCCTGCAAGCAGTTACAGTTTTGGCTTTCGCGAGGGGATGATTGCCAAc                                                                                         >  2:610912/1‑150 (MQ=255)
                                                          tCACCCGCGCGGATGGATGCTCTGCTGAAAGCCGGCACATTGCGCCTGTCGCTGACCATTAACCCCGCGCATGCGCAGCAAAAAATCGCCAGCGGCGATTTGCCTGCAAGCAGTTACAGTTTTGGCTTTCGCGAGGGGATGATTGGCAAc                                                                                         <  1:236521/150‑1 (MQ=255)
                                                                gcgcgGATGGATGCTCTGCTGAAAGCCGGCACATTGCGCCTGTCGCTGACCTTTAACCCCGCGCATGCGCAGCAAAAAATCGCCAGCGGCGATTTGCCTGCAAGCAGTTACAGTTTTGGCTTTCGCGAGGGGATGATTGGCAACGTGCAt                                                                                   >  2:320215/1‑150 (MQ=255)
                                                                   cgGATGGATGCTCTGCTGAAAGCCGGCACATTGCGCCTGTCGCTGACCTTTAACCCCGCGCATGCGCAGCAAAAAATCGCCCGCGGCGATTTGCCTGCAAGCAGTTACAGTTTTGGCTTTCGCGAGGGGATGATTGGCAACGTGCATTTc                                                                                >  2:1041974/1‑150 (MQ=255)
                                                                   cgGATGGATGCTCTGCTGAAAGCCGGCACATTGCGCCTGTCGCTGACCTTTAACCCCGCGCATGCGCAGCAAAAAATCGCCAGCGGCGATTTGCCTGCAAGCAGTTACAGTTTTGGCTTTCGCGAGGGGATGATTGGCAACGTGCATTTc                                                                                >  2:1047389/1‑150 (MQ=255)
                                                                   cgGATGGATGCTCTGCTGAAAGCCGGCACATTGCGCCTGTCGCTGACCTTTAACCCCGCGCATGCGCAGCAAAAAATCGCCAGCGGAGATTTGCCTGCAAGCAGTTACAGTTTTGGCTTTCGCGAGGGGATGATTGGCAACGTGCATTTc                                                                                >  2:1170611/1‑150 (MQ=255)
                                                                      ctggatgCGCTGCTGAAAGCCGGCACCTTGCGCCTGTCGCTGACCTTTAACCCCGCGCATGCGCAGCAAAAAATCGCCAGCGGCGATTTGCCTGCAAGCAGTTACAGTTGTGGCTTTCGCGAGGGGATGATTGGCAACGTGCATTTCGTc                                                                             <  1:1038089/149‑1 (MQ=255)
                                                                       tggatgCTCTGCTGAAAGCCGGCACATTGCGCCTGTCGCTGACCTTTAACCCCGCGCATGCGCAGCAAAAAATCGCCAGCGGCGATTTGCCTGCAAGCAGTTACAGTTTTGGCTTTCGCGAGGGGATGATTGGCAACGTGCATTTCGTCa                                                                            >  2:449800/1‑150 (MQ=255)
                                                                          atgCTCTGCTGAAAGCCGGCACATTGCGCCTGTCGCTGACCTTTAACCCCGCGCATGCGCAGCAAAAAATCGCCAGCGGCGATTTGCCTGCAAGCAGTTACAGTTTTGGCTTTCGCGAGGGGATGATTGGCAACGTGCATTTCGTCACCa                                                                         <  1:520025/150‑1 (MQ=255)
                                                                             ctctGCTGAAAGCCGGCACATTGCGCCTGTCGCTGACCTTTAACCCCGCGCCTGCGCAGCAAAAAATCGCCAGCGGCGGTTTGCCTGCAAGCAGTTTCAGTTTTGGCTTTCGCGTGGGGATGATTGGCAACGTGGATTTCGGCCCCATTc                                                                      >  2:618455/1‑150 (MQ=255)
                                                                               ctgctgAAAGCCGGCACATTGCGCCTGTCGCTGACCTTTAACCCCGCGCATGCGCAGCAAAAAATCGCCAGCGGCGATTTGCCTGCAAGCAGTTACAGTTTTGGCTTTCGCGAGGGGATGATTGGCAACGTGCATTTCGTCACCATTCCt                                                                    <  2:1317885/150‑1 (MQ=255)
                                                                               ctgctgAAAGCCGGCACATTGCGCCTGTCGCTGACCTTTAACCCCGCGCATGCGCAGCAAAAAATCGCCAGCGGCGATTTGCCTGCAAGCAGTTACAGTTTTGGCTTTCGCGAGGGGATGATTGGCAACGTGCATTTCGTCACCATTCCt                                                                    <  2:1310813/150‑1 (MQ=255)
                                                                                tgctgAAAGCCGGCACATTGCGCCTGTCGCTGACCTTTAACCCCGCGCATGCGCAGCAAAAAATCGCCAGCGGCGATTTGCCTGCAAGCAGTTACAGTTTTGGCTTTCGCGAGGGGATGATTGGCAACGTGCATTTCGTCACCATTCCTg                                                                   <  1:585886/150‑1 (MQ=255)
                                                                                tgctgAAAGCCGGCACATTGCGCCTGTCGCTGACCTTTAACCCCGCGCATGCGCAGCAAAAAATCGCCAGCGGCGATTTGCCTGCAAGCAGTTACAGTTTTGGCTTTCGCGAGGGGATGATTGGCAACGTGCATTTCGTCACCATTCCTg                                                                   >  1:1653303/1‑150 (MQ=255)
                                                                                    gAAAGCCGGCACATTGCGCCTGTCGCTGACCTTTAACCCCGCGCATGCGCAGCAAAAAATCGCCAGCGGCGATTTGCCTGCAAGCAGTAACAGATTTGGCTTTCGCGAGGGGATGATTGGCAACGTGCATTTCGTCCCCATTCCTGCCaa                                                               >  2:844668/1‑150 (MQ=255)
                                                                                     aaaGCCGGCACATTGCGCCTGTCGCTGACCTTTAACCCCGCGCATGCGCAGCTAAAAATCGCCAGCGGCGATTTGCCTGCAAGCAGTTACAGTTTTGGCTTTCGCGAGGGGATGATTGGCAACGTGCATTTCGTCACCATTCCTGCCAAc                                                              <  1:106564/150‑1 (MQ=255)
                                                                                     aaaGCCGGCACATTGCGCCTGTCGCTGACCTTTAACCCCGCGCATGCGCAGCAAAAAATCGCCAGCGGCGATTTGCCTGCAAGCAGTTACAGTTTTGGCTTTCGCGAGGGGATGATTGGCAACGTGCATTTCGTCACCATTCCTGCCAAc                                                              <  2:1646485/150‑1 (MQ=255)
                                                                                      aaGCCGGCACATTGCGCCTGTCGCTGACCTTTAACCCCGCGCATGCGCAGCAAAAAATCGCCAGCGGCGATTTGCCTGCAAGCAGTTACAGTTTTGGCTTTCGCGAGGGGATGATTGGCAACGTGCATTTCGTCACCATTCCTGCCAAcg                                                             >  1:976369/1‑150 (MQ=255)
                                                                                          cGGCACATTGCGCCTGTCGCTGACCTTTAACCCCGCGCATGCGCAGCAAAAAATCGCCAGCGGCGATTTGCCTGCAAGCAGTTACAGTTTTGGCTTTCGCGAGGGGATGATTGGCAACGTGCATTTCGTCACCATTCCTGCCAACGCGaa                                                         >  1:1689887/1‑150 (MQ=255)
                                                                                          cGGCACATTGCGCCTGTCGCTGACCTTTAACCCCGCGCATGCGCAGCAAAAAATCGCCAGCGGCGATTTGCCTGCAAGCAGTTACAGTTTTGGCTTTCGCGAGGGGATGATTGGCAACGTGCATTTCGTCACCATTCCTGCCAACGCGaa                                                         >  1:1594022/1‑150 (MQ=255)
                                                                                          cGGCACATTGCGCCTGTCGCTGACCTTTAACCCCGCGCATGCGCAGCAAAAAATCGCCAGCGGCGATTTGCCTGCAAGCAGTTACAGTTTTGGCTTTCGCGAGGGGATGATTGGCAACGTGCATTTCGTCACCATTCCTGCCAACGCGaa                                                         >  2:479870/1‑150 (MQ=255)
                                                                                                aTTGCGCCTGTCGCTGACCTTTAACCCCGCGCATGCGCAGCAAAAAATCGCCAGCGGCGATTTGCCTGCAAGCAGTTACAGTTTTGGCTTTCGCGAGGGGATGATTGGCAACGTGCATTTCGTCACCATTCCTGCCAACGCGAATGCCAg                                                   <  1:424209/150‑1 (MQ=255)
                                                                                                aTTGCGCCTGTCGCTGACCTTTAACCCCGCGCATGCGCAGCAAAAAATCGCCAGCGGCGATTTGCCTGCAAGCAGTTACAGTTTTGGCTTTCGCGAGGGGATGATTGGCAACGTGCATTTCGTCACCATTCCTGCCAACGCGAATGCCAg                                                   >  2:535898/1‑150 (MQ=255)
                                                                                                  tGCGCCTGTCGCTGACCTTTAACCCCGCGCTTGAGCAGCAAAAAATCGCCAGCGGCGATTTGCCTGCAAGCAGTTACAGTTTTGGCTTTCGCGAGGGGATGATTGGCAACGTGCATTTCGTCACCATTCCTGCCATCGCGAATGCCAGtg                                                 <  2:1615094/150‑1 (MQ=255)
                                                                                                  tGCGCCTGTCGCTGACCTTTAACCCCGCGCATGCGCAGCAAAAAATCGCCAGCGGCGATTTGCCTGCAAGCAGTTACAGTTTTGGCTTTCGCGAGGGGATGATTGGCAACGTGCATTTCGTCACCATTCCTGCCAACGCGAATGCCAGtg                                                 >  2:518868/1‑150 (MQ=255)
                                                                                                   gcgcCTGTCGCTGACCTTTAACCCCGCGCATGCGCAGCAAAAAATCGCCAGCGGCGATTTGCCTGCAAGCAGTTACAGTTTTGGCTTTCGCGAGGGGATGATTGGCAACGTGCATTTCGTCACCATTCCTGCCAACGCGAATGCCAGtgc                                                >  1:218163/1‑150 (MQ=255)
                                                                                                    cgcCTGTCGCTGACCTTTAACCCCGCGCATGCGCAGCAAAAAATCGCCAGCGGCGATTTGCCTGCAAGCAGTTACAGTTTTGGCTTTCGCGAGGGGATGATTGGCAACGTGCATTTCGTCACCATTCCTGCCAACGCGAATGCCAGtgct                                               >  2:1149250/1‑150 (MQ=255)
                                                                                                      ccTGTCGCTGACCTTTAACCCCGCGCATGCGCAGCAAAAAATCGCCAGCGGCGATTTGCCTGCAAGCAGTTACAGTTTTGGCTTTCGCGAGGGGATGATTGGCAACGTGCATTTCGTCACCATTCCTGCCAACGCGAATGCCAGtgctgc                                             >  1:793983/1‑150 (MQ=255)
                                                                                                        tggcGCTGACCTTTAACCCCGCGCATGCGCAGCAAAAAATCGCCAGCGGCGATTTGCCTGCAAGCAGTTACTGTTTTGGCTTTCGCGAGGGGATGATTGGCAACGTGCATTTCGTCACCATTCCTGCCAACGCGAATGCCAGTGCTgcgg                                           <  1:1588646/147‑1 (MQ=255)
                                                                                                         gTCGCTGACCTTTAACCCCGCGCATGCGCAGCAAAAAATCGCCAGCGGCGATTTGCCTGCAAGCAGTTACAGTTTTGGCTTTCGCGAGGGGATGATTGGCAACGTGCATTTCGTCACCATTCCTGCCAACGCGAATGCCAGTGCTgcggc                                          <  1:1563993/150‑1 (MQ=255)
                                                                                                               gACCTTTAACCCCGCGCATGCGCAGCAAAAAATCGCCAGCGGCGATTAGCCTGCAAGCAGTTACAGTTTTGGCTTTCGCGAGGGGATGATTGGCAACGTGCATTTCGTCACCATTCCTGCCAACGCGAATGCCAGTGCTGCGGCGAAGGt                                    >  2:966065/1‑150 (MQ=255)
                                                                                                                 ccTTTAACCCCGCGCATGCGCAGCAAAAAATCGCCAGCGGCGATTTGCCTGCAAGCAGTTACAGTTTTGGCTTTCGCGAGGGGATGATTGGCAACGTGCATTTCGTCACCATTCCTGCCAACGCGAATGCCAGTGCTGCGGCGAAGGTAg                                  >  2:1484006/1‑150 (MQ=255)
                                                                                                                 ccTTTAACCCCGCGCATGCGCAGCAAAAAATCGCCAGCGGCGATTTGCCTGCAAGCAGTTACAGTTTTGGCTTTCGCGAGGGGATGATTGGCAACGTGCATTTCGTCACCATTCCTGCCAACGCGAATGCCAGTGCTGCGGCGAAGGTAg                                  >  2:855811/1‑150 (MQ=255)
                                                                                                                 ccTTTAACCCCGCGCATGCGCAGCAAAAAATCGCCAGCGGCGATTTGCCTGCAAGCAGTTACAGTTTTGGCTTTCGCGAGGGGATGATTGGCAACGTGCATTTCGTCACCATTCCTGCCAACGCGAATGCCAGTGCTGCGGCGAAGGTAg                                  >  2:432951/1‑150 (MQ=255)
                                                                                                                 ccTTTAACCCCGCGCATGCGCAGCAAAAAATCGCCAGCGGCGATTTGCCTGCAAGCAGTTACAGTTTTGGCTGTCGCGAGGGGATGATTGGCAACGTGCATTTCGTCACCATTCCTGCCAACGCGAATGCCAGTGCTGCGGCGAAGGTAg                                  >  2:1055471/1‑150 (MQ=255)
                                                                                                                  cTTTAACCCCGCGCATGCGCAGCAAAAAATCGCCAGCGGCGATTTGCCTGCAAGGAGTTACAGTTTTGGCTTTCGCGAGGGGATGATTGGCAACGTGCATTTCGTCACCATTCCTGCCAACGCGAATGCCAGTGCTGCGGCGAAGGTAGt                                 <  1:16942/150‑1 (MQ=255)
                                                                                                                  cTTTAACCCCGCGCATGCGCAGCAAAAAATCGCCAGCGGCGATTTGCCTGCAAGCAGTTACAGTTTTGGCTTTCGCGAGGGGATGATTGGCAACGTGCATTTCGTCACCATTCCTGCCAACGCGAATGCCAGTGCTGCGGCGAAGGTAGt                                 <  1:762303/150‑1 (MQ=255)
                                                                                                                    ttAACCCCGCGCATGCGCAGCAAAAAATCGCCAGCGGCGATTTGCCTGCAAGCAGTTACAGTTTTGGCTTTCGCGAGGGGATGATTGGCAACGTGCATTTCGTCACCATTCCTGCCAACGCGAATGCCAGTGCTGCGGCGAAGGTAGTTg                               <  2:632141/150‑1 (MQ=255)
                                                                                                                      aaCCCCGCGCATGCGCAGCAAAAAATCGCCAGCGGCGATTTGCCTGCAAGCAGTTACAGTTTTGGCTTTCGCGAGGGGATGATTGGCAACGTGCATTTCGTCACCATTCCTGCCAACGCGAATGCCAGTGCTGCGGCGAAGGTAGTTGcc                             <  1:1207589/150‑1 (MQ=255)
                                                                                                                      aaCCCCGCGCATGCGCAGCAAAAAATCGCCAGCGGCGATTTGCCTGCAAGCAGTTACAGTTTAGGCTTTCTCGAGGGGATGATTGGCAACGTGCATTTCGTCACCATTACTGCAAACGCGAATGCCAGTGCTTCGGCGAAGGTAGTTGcc                             >  2:1047812/1‑150 (MQ=255)
                                                                                                                       aCCCCGCGCATGCGCAGCAAAAAATCGCCAGCGGCGATTTGCCTGCAAGCAGTTACAGTTTTGGCTTTCGCGAGGGGATGATTGGCAACGTGCATTTCGTCACCATTCCTGCCAACGCGAATGCCAGTGCTGCGGCGAAGGTAGTTGCCa                            >  1:1547374/1‑150 (MQ=255)
                                                                                                                            gcgcATGCGCAGCAAAAAATCGCCATCGGCGATTTGCCTGCAAGCAGTTACAGTTTTGGCTTTCGCGAGGGGATGATTGGCAACGTGCATTTCGTCACCATTCCTGCCAACGCGAATGCCAGTGCTGCGGCGAAGGTAGTTGCCAATTTc                       <  2:1306749/150‑1 (MQ=255)
                                                                                                                             cgcATGCGCAGCAAAAAATCGCCAGCGGCGATTTGCCTGCAAGCAGTTACAGTTTTGGCTTTCGCGAGGGGATGATTGGCAACGTGCATTTCGTCACCATTCCTGCCAACGCGAATGCCAGTGCTGCGGCGAAGGTAGTTGCCAATTTcc                      <  1:1331385/150‑1 (MQ=255)
                                                                                                                               cATGCGCAGCAAAAAATCGCCAGCGGCGATTTGCCTGCAAGCAGTTACAGTTTTGGCTTTCGCGAGGGGATGATTGGCAACGTGCATTTCGTCACCATTCCTGCCAACGCGAATGCCAGTGCTGCGGCGAAGGTAGTTGCCAATTTCCTg                    >  1:686266/1‑150 (MQ=255)
                                                                                                                                 tgagcagcaAAAAATCGCCACCGGCGATTTGCCTGCAAGCAGTTACAGTTTTGGCTTTCGCGAGGGGATGATTGGCAACGTGCATTTCGTCACCATTCCTGCCAACGCGAATGCCAGTGCTGCGGCGAAGGTAGTTGCCAATTTCCTGct                  <  2:330369/147‑1 (MQ=255)
                                                                                                                                 tGCGCAGCAAAAAATCGCCAGCGGCGATTTGCCTGCAAGCAGTTACAGTTTTGGCTTTCGCGAGGGGATGATTGGCAACGTGCATTTCGTCACCATTCCTGCCAACGCGAATGCCAGTGCTGCGGCGAAGGTAGTTGCCAATTTCCTGct                  <  1:1472431/150‑1 (MQ=255)
                                                                                                                                    gcagcaAAAAATCGCCAGCGGCGATTTGCCTGCAAGCAGTTACAGTTTTGGCTTTCGCGAGGGGATGATTGGCAACGTGCATTTCGTCACCATTCCTGCCAACGCGAATGCCAGTGCTGCGGCGAAGGTAGTTGCCAATTTCCTGCtctc               >  1:1200411/1‑150 (MQ=255)
                                                                                                                                    gcagcaAAAAATCGCCAGCGGCGATTTGCCTGCAAGCAGTTACAGTTTTGGCTTTCGCGAGGGGATGATTGGCAACGTGCATTTCGTCACCATTCCTGCCAACGCGAATGCCAGTGCTGCGGCGAAGGTAGTTGCCAATTTCCTGCtctc               >  1:474953/1‑150 (MQ=255)
                                                                                                                                     cagcaAAAAATCGCCAGCGGCGATTTGCCTGCAAGCAGTTACAGTTTTGGCTTTCGCGAGGGGATGATTGGCAACGTGCATTTCGTCACCATTCCTGCCAACGCGAATGCCAGTGCTGCGGCGAAGGTAGTTGCCATTTTCCTGCTCTCa              >  1:393411/1‑150 (MQ=255)
                                                                                                                                       gcaAAAAATCGCCAGCGGCGATTTGCCTGCAAGCAGTTACAGTTTTGGCTTTCGCGAGGGGATGATTGGCAACGTGCATTTCGTCACCATTCCTGCCAACGCGAATGCCAGTGCTGCGGCGAAGGTAGTTGCCAATTTCCTGCTCTCAcc            >  1:485382/1‑150 (MQ=255)
                                                                                                                                       gcaAAAAATCGCCAGCGGCGATTTGCCTGCAAGCAGTTACAGTTTTGGCTTTCGCGAGGGGATGATTGGCAACGTGCATTTCGTCACCATTCCTGCCAACGCGAATGCCAGTGCTGCGGCGAAGGTAGTTGCCAATTTCCTGCTCTCAcc            >  1:1225998/1‑150 (MQ=255)
                                                                                                                                       caggcaaaTCGCCGCTGGCGATTTTCCTGCAAGCAGTTACAGTTTTGGCTTTCGCGAGGGGATGATTGGCAACGTGCATTTCGTCACCATTCCTGCCAACGCGAATGCCAGTGCTGCGGCGAAGGTAGTTGCCAATTTCCTGCTCTCAcc            >  2:549874/6‑150 (MQ=255)
                                                                                                                                        caAAAAATCGCCAGCGGCGATTTGCCTGCAAGCAGTTACAGTTTTGGCTTTCGCGAGGGGATGATTGGCAACGTGCATTTCGTCACCATTCCTGCCAACGCGAATGCCAGTGCTGCGGCGAAGGTAGTTGCCAATTTCCTGCTCTCAccc           >  1:447403/1‑150 (MQ=255)
                                                                                                                                        aaaaaaaTCGCCAGAGGCGATTTGCCTGCAAGCAGTTACAGTTTTGGCTTTCGCGAGGGGATGATTGGCAACGTGCATTTCGTCACCATTCCTGCCAACGCGAATGCCAGTGCTGCGGCGAAGGTAGTTGCCAATTTCCTTCTCTCAccc           <  2:1253765/149‑1 (MQ=255)
                                                                                                                                         ggcaaaTCGCCGCTGGCGATTTGCCTGCAAGCAGTTACAGTTTTGGCTTTCGCGAGGGGATGATTGGCAACGTGCATTTCGTCACCATTCCTGCCAACGCGAATGCCAGTGCTGCGGCGAAGGTAGTTGCCAATTTCCTGCTCTCACCCg          >  1:1045685/4‑150 (MQ=255)
                                                                                                                                         aaaaaaTCGCCAGCGGCGATTTGCCTGCAAGCAGTTACAGATTTGGCTTTCGCGAGGGGATGATTCGCAACGTGCATTACGTCACCATACCTGCCAACGCGAATGCCAGTGCTGCGGCGAAGGAAGTTGCCAAGTGCGTGCTCGCAcgcg          >  1:315691/1‑150 (MQ=255)
                                                                                                                                          gcaaaTCGCCGCTGGCGATTTGCCTGCAAGCAGTTACAGTTTTGGCTTTCGCGAGGGGATGATTGGCAACGTGCATTTCGTCACCATTCCTGCCAACGCGAATGCCAGTGCTGCGGCGAAGGTAGTTGCCAATTTCCTGCTCTCACCCGa         >  2:1592315/3‑150 (MQ=255)
                                                                                                                                           aaaaTCGCCAGCGGCGATTTGCCTGCAAGCAGTTACAGTTTTGGCTTTCGCGAGGGGATGATTGGCAACGTGCATTTCGTCACCATTCCTGCCAACGCGAATGCCAGTGCTGCGGCGAAGGTAGTTGCCATTTTCCTGCTCTCACCCGAt        >  1:354877/1‑150 (MQ=255)
                                                                                                                                            aaaTCGCCAGCGGCGATTTGCCTGCAAGCAGTTACAGTTTTGGCTTTCGCGAGGGGATGATTGGCAACGTGCATTTCGTCACCATTCCTGCCAACGCGAATGCCAGTGCTGCGGCGAAGGTAGTTGCCAATTTCCTGCTCTCACCCGATg       <  1:985304/150‑1 (MQ=255)
                                                                                                                                              aTCGCCAGCGGCGATTTGCCTGCAAGCAGTTACAGTTTTGGCTTTCGCGAGGGGATGATTGGCAACGTGCATTTCGTCACCATTCCTGCCAACGCGAATGCCAGTGCTGCGGCGAAGGTAGTTGCCAATTTCTTGCTCTCCCCCGATgcg     >  2:1151725/1‑150 (MQ=255)
                                                                                                                                               tCGCCAGCGGCGATTTGCCTGCAAGCTGTTACAGTTTTGGCTTTCGCGAGGGGATGATTGGCAACGTGCATTTCGTCACCATTCCTGCCAACGCGAATGCCAGTGCTGCGGCGAAGGTAGTTGCCAATTTCCTGCTCTCACCCGATgcgc    <  2:864659/150‑1 (MQ=255)
                                                                                                                                               tCGCCAGCGGCGATTTGCCTGCAAGCAGTTACAGTTTTGGCTTTCGCGAGGGGATGATTGGCAACGTGCATTTCGTCACCATTCCTGCCAACGCGAATGCCAGTGCTGCGGCGAAGGTAGTTGCCAATTTCCTGCTCTCACCCGATgcgc    <  1:867410/150‑1 (MQ=255)
                                                                                                                                                 gCCAGCGGCGATTTGCCTGCAAGCAGTTACAGTTTTGGCTTTCGCGAGGGGATGATTGGCAACGTGCATTTCGTCACCATTCCTGCCAACGCGAATGCCAGTGCTGCGGCGAAGGTAGTTGCCAATTTCCTGCTCTCACCCGATGCGCaa  >  1:1222963/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                 
GCAATATCTTGATGTGCTGCATCCGTATTTGTGGCGCGAAGGAAAGGATTTCCCGCCTTCACCCGCGCGGATGGATGCTCTGCTGAAAGCCGGCACATTGCGCCTGTCGCTGACCTTTAACCCCGCGCATGCGCAGCAAAAAATCGCCAGCGGCGATTTGCCTGCAAGCAGTTACAGTTTTGGCTTTCGCGAGGGGATGATTGGCAACGTGCATTTCGTCACCATTCCTGCCAACGCGAATGCCAGTGCTGCGGCGAAGGTAGTTGCCAATTTCCTGCTCTCACCCGATGCGCAA  >  NC_000913/1836756‑1837050

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: