Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,182,972 G→C 5.6% I270M (ATC→ATG)  PP_4563 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,182,9720GC 5.6% 285.2 / 7.0 89I270M (ATC→ATG) PP_4563hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (42/42);  new base C (2/3);  total (44/45)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.10e-01

TGGCAAGCGCCGCTCCACCATCAGCCTGGCCCTCAGCGCCGACCGGGACTCGCTGCTGCGCGACACCCGTCCGGGTGGCGAGCAGCTGGACTTCCGCAACTTCGTACAGGCGTAAGGCCGTCAGGCTTTCTCGATGCGATCGTCATGGATGACGATCAGGCCCTGCTTGAACAACCCGCCAATGGCCTTCTTGAAGTTGCCCTTGCTGACGTTGAACAGCTTGCTGATCAACGCAGGGTCACTCTTGTCGCACACCCCCAACACCCCGCCTTCGGCTTCCAGGCGCGCCATGAT  >  NC_002947/5182825‑5183118
                                                                                                                                                   |                                                                                                                                                  
tGGCAAGCGCCGCTCCACCATCAGCCTGGCCCTCAGCGCCGACCGGGACTCGCTGCTGCGCGACACCCGTCCGGGTGGCGAGCAGCTGGACTTCCGCAACTTCGTACAGGCGTAAGGCCGTCAGGCTTTCTCGATGCGATCGTCATGGAt                                                                                                                                                  <  2:1628394/150‑1 (MQ=255)
   cAAGCGCCGCTCCACCATCAGCCTGGCCCTCAGCGCCGACCGGGACTCGCTGCTGCGCGACACCCGTCCGGGTGGCGAGCAGCTGGACTTCCGCAACTTCGTACAGGCGTAAGGCCGTCAGGCTTTCTCGATGCGATCGTCATGGATGAc                                                                                                                                               <  1:358347/150‑1 (MQ=255)
      gcgcCGCTCCACCATCAGCCTGGCCCTCAGCGCCGACCGGGACTCGCTGCTGCGCGACACCCGTCCGGGTGGCGAGCAGCTGGACTTCCGCAACTTCGTACAGGCGTAAGGCCGTCAGGCTTTCTCGATGCGATCGTCATCCATGACGAt                                                                                                                                            <  2:1196563/150‑1 (MQ=255)
        gccgcTCCACCATCAGCCTGGCCCTCAGCGCCGACCGGGACTCGCTGCTGCGCGACACCCGTCCGGGTGGCGAGCAGCTGGACTTCCGCAACTTCTTAAAGGCGTAAGGCCGACCGGCTTTCCCGCTGAGATCGTCATGGATGACGATCa                                                                                                                                          <  1:54176/150‑1 (MQ=255)
          cgcTCCACCATCAGCCTGGCCCTCAGCGCCGACCGGGACTCGCTGCTGCGCGACACCCGTCCGGGTGGCGAGCAGCTGGACTTCCGCAACTTCGTACAGGCGTAAGGCCGTCAGGCTTTCTCGATGCGATCGTCATGGATGACGATCAgg                                                                                                                                        <  1:432345/150‑1 (MQ=255)
           tgttcaccaTCAGCCTGGCCCTCAGCGCCGACCGGGACTCGCTGCTGCGCGACACCCGTCCGGGTGGCGAGCAGCTGGACTTCCGCAACTTCGTACAGGCGTAAGGCCGTCAGGCTTTCTCGATGCGATCGTCATGGATGACGATCAGGc                                                                                                                                       >  2:1291994/5‑150 (MQ=255)
           gcTCCACCATCAGCCTGGCCCTCAGCGCCGACCGGGACTCGCTGCTGCGCGGCGCCCGTCCGGGTGGCGAGCAGCTGGACTTCCGCAACTTCGTACAGGCGTAAGGCCGTCAGGCTTTCTCGATGCGATCGTCATGGATGACGATCAGGc                                                                                                                                       >  1:2272951/1‑150 (MQ=255)
           gcTCCACCATCAGCCTGGCCCTCAGCGCCGACCGGGACTCGCTGCTGCGCGACACCCGTCCGGGTGGCGAGCAGCTGGCCTTCCGCAACTTCGTACAGGAGTAATGCCGTCAGGCTTTCTCGATGCGATCGTCATGGATGACGATCAGGc                                                                                                                                       <  1:1617185/150‑1 (MQ=255)
            cTCCACCATCAGCCTGGCCCTCAGCGCCGACCGGGACTCGCTGCTGCGCGACACCCGTCCGGGTGGCGAGCAGCTGGACTTCCGCAACTTCGTACAGGCGTAAGGCCGTCAGGCTTTCTCGATGCGATCGTCATGGATGACGATCAGGcc                                                                                                                                      <  1:1893590/150‑1 (MQ=255)
            cTCCACCATCAGCCTGGCCCTCAGCGCCGACCGGGACTCGCTGCTGCGCGACACCCGTCCGGGTGGCGAGCAGCTGGACTTCCGCAACTTCGTACAGGCGTAAGGCCGTCAGGCTTTCTCGATGCGATCGTCATCCATGACGATcgcatc                                                                                                                                      <  1:853112/150‑6 (MQ=255)
               caccaTCAGCCTGGCCCTCAGCGCCGACCGGGACTCGCTGCTGCGCGACACCCGTCCGGGTGGCGAGCAGCTGGACTTCCGCAACTTCGTACAGGCGTAAGGCCGTCAGGCTTTCTCGATGCGATCGTCATGGATGACGATCAGGCCCTg                                                                                                                                   >  2:2026703/1‑150 (MQ=255)
               caccaTCAGCCTGGCCCTCAGCGCCGACCGGGACTCGCTGCTGCGCGACACCCGTCCGGGTGGCGAGCAGCTGGACTTCCGCAACTTCGTACAGGCGTAAGGCCGTCAGGCTTTCTCGATGCGATCGTCATGGATGACGATCAGGCCCTg                                                                                                                                   >  2:2270612/1‑150 (MQ=255)
                 ccaTCAGCCTGGCCCTCAGCGCCGACCGGGACTCGCTGCTGCGCGACACCCGTCCGGGTGGCGAGCAGCTGGACTTCCGCAACTTCGTACAGGCGTAAGGCCGTCAGGCTTTCTCGATGCGATCGTCATGGATGa                                                                                                                                                <  2:1553041/135‑1 (MQ=255)
                 ccaTCAGCCTGGCCCTCAGCGCCGACCGGGACTCGCTGCTGCGCGACACCCGTCCGGGTGGCGAGCAGCTGGACTTCCGCAACTTCGTACAGGCGTAAGGCCGTCAGGCTTTCTCGATGCGATCGTCATGGATGa                                                                                                                                                >  1:1553041/1‑135 (MQ=255)
                 ccaTCAGCCTGGCCCTCAGCGCCGACCGGGACTCGCTGCTGCGCGACACCCGTCCGGGTGGCGAGCAGCTGGACTTCCGCAACTTCGTACAGGCGTAAGGCCGTCAGGCTTTCTCGATGCGATCGTCATGGATGACGATCAGGCCCTGc                                                                                                                                  <  1:249489/149‑1 (MQ=255)
                 ccaTCAGCCTGGCCCTCAGCGCCGACCGGGACTCGCTGCTGCGCGACACCCGTCCGGGTGGCGAGCAGCTGGACTTCCGCAACTTCGTACAGGCGTAAGGCCGTCAGGCTTTCTCGATGCGATCGTCATGGATGACGATCAGGCCCTGc                                                                                                                                  >  2:249489/1‑149 (MQ=255)
                    tCAGCCTGGCCCTCAGCGCCGACCGGGACTCGCTGCTGCGCGACACCCGTCCGGGTGGCGAGCAGCTGGACTTCCGCAACTTCGTACAGGCGTAAGGCCGTCAGGCTTTCTCGATGCGCTGGCGAGCGATGACGATCAGGCCCTGCTTGa                                                                                                                              >  1:1819292/1‑150 (MQ=255)
                     cAGCCTGGCCCTCAGCGCCGACCGGGACTCGCTGCTGCGCGACACCCGTCCGGGTGGCGAGCAGCTGGACTTCCGCAACTTCGTACAGGCGTAAGGCCGTCAGGCTTTCTCGATGCGATCGTCATGGATGACGATCAGGCCCTGCTTGaa                                                                                                                             <  1:1014463/150‑1 (MQ=255)
                     cAGCCTGGCCCTCAGCGCCGACCGGGACTCGCTGCTGCGCGACACCCGTCCGGGTGGCGAGCAGCTGGACTTCCGCAACTTCGTACAGGCGTAAGGCCGTCAGGCTTTCTCGATGCGATCGTCATGGATGACGATCAGGCCCTGCTTGaa                                                                                                                             <  2:956140/150‑1 (MQ=255)
                            gCCCTCAGCGCCGACCGGGACTCGCTGCTGCGCGACACCCGTCCGGGTGGCGAGCAGCTGGACTTCCGCAACTTCGTACAGGCGTAAGGCCGTCAGGCTTTCTCGATGCGATCGTCATGGATGACGATCAGGCCCTGCTTGAACAACCCg                                                                                                                      >  2:2193837/1‑150 (MQ=255)
                             cccTCAGCGCGGACCGGGACTCGCTGCTGCGCGACACCCGTCCGGGTGGCGAGCAGCTGGACTTCCGCAACTTCGTACAGGCGTAAGGCCGTCAGGCTTTCTCGATGCGATCGTCATGGATGACGATCAGGCCCTGCTTGAACAACCCGc                                                                                                                     <  1:1285256/150‑1 (MQ=255)
                              ccTCAGCGCCGACCGGGACTCGCTGCTGCGCGACACCCGTCCGGGTGGCGAGCAGCTGGACTTCCGCAACTTCGTACAGGCGTAAGGCCGTCAGGCTTTCTCGATGCGATCGTCATGGATGACGATCAGGCCCTGCTTGAACAACCCGcc                                                                                                                    >  2:1688232/1‑150 (MQ=255)
                                 cAGCGCCGACCGGGACTCGCTGCTGCGCGACACCCGTCCGGGTGGCGAGCAGCTGGACTTCCGCAACTTCGTACAGGCGTAAGGCCGTCAGGCTTTCTCGATGCGATCGTCATGGATGACGATCAGGCCCTGCTTGAACAACCCGCCAAt                                                                                                                 >  2:2087004/1‑150 (MQ=255)
                                   gcgcCGACCGGGAGTCGCTGCTGCGCGACAGCCGTGCGGGTGGCGCGCAGCTGGACTTCCGCAACTGCGTGCAGGCGTAAGGCCGTCAGGCTTTCTCGATGCGATCGTCATGGATGACGATCAGGCCCTGCTTGAACAACCCGCCAATgg                                                                                                               <  2:182181/150‑1 (MQ=255)
                                    cgcCGACCGGGACTCGCTGCTGCGCGACACCCGTCCGGGTGGCGAGCAGCTGGACTTCCGCAACTTCGTACAGGCGTAAGGCCGTCAGGCTTTCTCGATGCGATCGTCATGGATGACGATCAGGCCCTGCTTGAACAACCCGCCAATGGc                                                                                                              <  1:2026703/150‑1 (MQ=255)
                                     gcCGACCGGGACTCGCTGCTGCGCGACACCCGTCCGGGTGGCGAGCAGCTGGACTTCCGCAACTTCGTACAGGCGTAAGGCCGTCAGGCTTTCTCGATGCGATCGTCATGGATGACGATCAGGCCCTGCTTGAACAACCCGCCAATGGcc                                                                                                             <  1:2193837/150‑1 (MQ=255)
                                        gACCGGGACTCGCTGCTGCGCGACACCCGTCCGGGTGGCGAGCAGCTGGACTTCCGCAACTTCGTACAGGCGTAAGGCCGTCAGGCTTTCTCGATGTGATCGTCATCCATGACGATc                                                                                                                                           >  2:2066244/1‑117 (MQ=255)
                                        gACCGGGACTCGCTGCTGCGCGACACCCGTCCGGGTGGCGAGCAGCTGGACTTCCGCAACTTCGTACAGGCGTAAGGCCGTCAGGCTTTCTCGATGTGATCGTCATCCATGACGATc                                                                                                                                           <  1:2066244/117‑1 (MQ=255)
                                             ggACTCGCTGCTGCGCGACACCCGTCCGGGTGGCGAGCAGCTGGACTTCCGCAACTTCGTACAGGCGTACGCCAGTATCGCTTTCTCCAGGGGACTGTCATGGATGACGATCAGGCCCTGCTTGAACAACCCGCCAATGGCCTTCTTGaa                                                                                                     <  1:291011/150‑1 (MQ=255)
                                                cTCGCTGCTGCGCGACACCCGTCCGGGTGGCGAGCAGCTGGACTTCCGCAACTTCGTACAGGCGTAAGGCCGTCAGGCTTTCTCGATGCGATCGTCATGGATGACGATCAGGCCCTGCTTGaacaa                                                                                                                          <  2:1139983/126‑1 (MQ=255)
                                                cTCGCTGCTGCGCGACACCCGTCCGGGTGGCGAGCAGCTGGACTTCCGCAACTTCGTACAGGCGTAAGGCCGTCAGGCTTTCTCGATGCGATCGTCATGGATGACGATCAGGCCCTGCTTGaacaa                                                                                                                          >  1:1139983/1‑126 (MQ=255)
                                                      gctgcGCGACACCCGTCCGGGTGGCGAGCAGCTGGACTTCCGCAACTTCGTACAGGCGTAAGGCCGTCAGGCTTTCTCGATGCGATCGTCATGGATGACGATCAGGCCCTGCTTGAACAACCCGCCAATGGCCTTCTTGAAGTTGCCGtt                                                                                            >  1:503304/1‑150 (MQ=255)
                                                        tgcGCGACACCCGTCCGGGTGGCGAGCAGCTGGACTTCCGCAACTTCGTACAGGCGTAAGGCCGTCAGGCTTTCTCGATGCGATCGTCATGGATGACGATCAGGCCCTGCTTGAACAACCCGCCAATGGCCTTCTTGAAGTTGCCCTTGc                                                                                          <  1:678876/150‑1 (MQ=255)
                                                               cacCCGTCCGGGTGGCGAGCAGCTGGACTTCCGCAACTTCGTACAGGCGTAAGGCCGTCAGGCTTTCTCGATGCGATCGTCATGGATGACGATCAGGCCCTGCTTGAACAACCCGCCAATGGCCTTCTTGAAGTTGCCCTTGCTGACGtt                                                                                   <  2:55220/150‑1 (MQ=255)
                                                               cacCCGTCCGGGTGGCGAGCAGCTGGACTTCCGCAACTTCGTACAGGCGTAAGGCCGTCAGGCGTTCTCGATGCGATCGTCATGGATGACGATCAGGCCCTGCTTGAACAACCCGCCAATGGCCTTCTTGAAGTTGCGCTTGCTGACGtt                                                                                   >  1:813923/1‑150 (MQ=255)
                                                               cacCCGTCCGGGTGGCGAGCAGCTGGACTTCCGCAACTTCGTACAGGCGTAAGGCCGGCAGGCTTGCGCGATGCGATCGTCATGGATGACGATCAGGCCCTGCTTGAACAACCCGCCAATGGCATTCTTGAAGTTGCCCTTGCTGACGtt                                                                                   >  1:494830/1‑150 (MQ=255)
                                                                acCCGTCCGGGTGGTGAGCAGCTGGACTTCCGCAACTTCGTACAGGCGTAAGGCCGTCAGGCTTTCTCGATGCGATCGTCATGGATGACGATCAGGCCCTGCTTGAACAACCCGCCAATGGCCTTCTTGAAGtt                                                                                                  <  2:2065505/134‑1 (MQ=255)
                                                                acCCGTCCGGGTGGTGAGCAGCTGGACTTCCGCAACTTCGTACAGGCGTAAGGCCGTCAGGCTTTCTCGATGCGATCGTCATGGATGACGATCAGGCCCTGCTTGAACAACCCGCCAATGGCCTTCTTGAAGtt                                                                                                  >  1:2065505/1‑134 (MQ=255)
                                                                acCCGTCCGGGTGGCGAGCAGCTGGACTTCCGCAACTTCGTACAGGCGTAAGGCCGTCAGGCTTTCTCGATGCGATCGTCATGGATGACAGTCAGTCCCTGCTTGAACAACCCGCCAAGGGCCTTCTTGAAGTTGCCCTTGCTGACGTTg                                                                                  >  1:2231800/1‑150 (MQ=255)
                                                                 cccGTCCGGGTGGCGAGCAGCTGGACTTCCGCAACTTCGTACAGGCGTAAGGCCGTCAGGCTTTCTCGATGCGATCGTCATGGATGACGATCAGGCCCTGCTTGAACAACCCGCCAATGGCCTTCTTGAAGTTGCCCTTGCTGACGTTGa                                                                                 >  1:2356577/1‑150 (MQ=255)
                                                                  ccGTCCGGGTGGCGAGCAGCTGGACTTCCGCAACTTCGTACAGGCGTAAGGCCGTCAGGCTTTCTCGATGCGATCGTCATGGATGACGATCAGGCCCTGCTTGAACAACCCGCCAATGGCCTTCTTGaa                                                                                                     >  1:238606/1‑129 (MQ=255)
                                                                  ccGTCCGGGTGGCGAGCAGCTGGACTTCCGCAACTTCGTACAGGCGTAAGGCCGTCAGGCTTTCTCGATGCGATCGTCATGGATGACGATCAGGCCCTGCTTGAACAACCCGCCAATGGCCTTCTTGaa                                                                                                     <  2:238606/129‑1 (MQ=255)
                                                                  ccGTCCGGGTGGCGAGCAGCTGGACTTCCGCAACTTCGTACAGGCGTAAGGCCGTCAGGCTTTCTCGATGCGATCGTCATGGATGACGATCAGGCCCTGCTTGAACAACCCGCCAATGGCCTTCTTGAAGTTGCCCTTGCTGACGTTGaa                                                                                >  2:2120642/1‑150 (MQ=255)
                                                                  ccGTCCGGGTGGCGAGCAGCTGGACTTCCGCAACTTCGTACAGGCGTAAGGCCGTCAGGCTTTCTCGATGCGATCGTCATGGATGACGATCAGGCCCTGCTTAAACAACCCGCCAATGGCCTTCTTGAAGTTGCCCTTGCTGACGTTGaa                                                                                >  2:313958/1‑150 (MQ=255)
                                                                  ccGNCCGGGTGGCGAGCAGCTGGACTTCCGCAACTTCGTACAGGCGTAAGGCCGTCAGGCTTTCTCGATGCGATCGTCATGGATGACGATCAGGCCCTGCTTGAACAACCCGCCAATGGCCTTCTTGAAGTTGCCCTTGCTGACGTTGaa                                                                                >  2:1501490/1‑150 (MQ=255)
                                                                    gTCCGGGTGGCGAGCAGCTGGACTTCCGCAACTTCGTACAGGCGTAAGGCCGTCAGGCTTTCTCGATGCGATCGTCATGGATGACGATCAGGCCCTGCTTGAACAACCCGCCAATGGCCTTCTTGAAGTTGCCCTTGCTGACGTTGAACa                                                                              <  1:1291994/150‑1 (MQ=255)
                                                                      ccGGGTGGCGAGCAGCTGGACTTCCGCAACTTCGTACAGGCGTAAGGCCGTCAGGCTTTCTCGATGCGATCGTCATGGATGACGATCAGGCCCTGCTTAAACAACCCGCCAATGGCCTTCTTGAAGTTGCCCTTGCTGACGTTGAACAGc                                                                            <  1:299346/150‑1 (MQ=255)
                                                                        gggTGGCGAGCAGCTGGACTTCCGCAACTTCGTACAGGCGTAAGGCCGTCAGGCTTTCTCGATGCGATCGTCATGGATGACGATCAGGCCCTGCTTGAACAACCCGACAACGGCCTACTTTAATTTGGCGTTGTTGACTTTGAACAGCtt                                                                          >  2:623991/1‑150 (MQ=255)
                                                                               gAGCAGCTGGACTTCCGCAACTTCGTACAGGCGTAAGGCCGTCAGGCTTTCTCGATGCGATCGTCATGGATGACGATCAGGCCCTGCTTGAACAACCCGCCAATGGCCTTCTTGAAGTTGCCCTTGCTGACGTTGAACAGCTTGCTGATc                                                                   >  2:986615/1‑150 (MQ=255)
                                                                               gAGCAGCTGGACTTCCGCAACTTCGTACAGGCGTAAGGCCGTCAGGCTTTCTCGATGCGATCGTCATGGATGACGATCAGGCCCTGCTTGAACAACCCGCCAATGGCCTTCTTGAAGTTGCCCTTGCTGACGTTGAACAGCTTGCTGATc                                                                   >  1:1557222/1‑150 (MQ=255)
                                                                                     cTGGACTTCCGCAACTTCGTACAGGCGTAAGGCCGTCAGGCTTTCTCGATGCGATCGTCATGGATGACGATCAGGCCCTGCTTGaa                                                                                                                             >  1:621012/1‑86 (MQ=255)
                                                                                     cTGGACTTCCGCAACTTCGTACAGGCGTAAGGCCGTCAGGCTTTCTCGATGCGATCGTCATGGATGACGATCAGGCCCTGCTTGaa                                                                                                                             <  2:621012/86‑1 (MQ=255)
                                                                                      tGGACTTCCGCAACTTCGTACAGGCGTAAGGCCGTCAGGCTTTCTCGATGCGATCGTCATGGATGACGATCAGGCCCTGCTTGAACAACCCGCCAATGGCCTTCTTGAAGTTGCCCTTGCTGACGTTGAACAGCTTGCTGATCAACGCAg                                                            <  1:2270612/150‑1 (MQ=255)
                                                                                       ggACTTCCGCAACTTCGTACAGGCGTAAGGCCGTCAGGCTTTCTCGATGCGATCGTCATGGATGACGATCAGGCCCTGCTTGAACAACCCGCCAATGGCCTTCTTGAAGTTGCCCTTGCTGACGTTGAACAGCTTGCTGATCAACGCAgg                                                           <  1:2087004/150‑1 (MQ=255)
                                                                                        gACTTCCGCAACTTCGTACAGGCGTAAGGCCGTCAGGCTTTCTCGATGCGATCGTCATGGATGACGATCAGGTCCTGCTTGAACAACCCGCCAATg                                                                                                                <  2:1184843/96‑1 (MQ=255)
                                                                                        gACTTCCGCAACTTCGTACAGGCGTAAGGCCGTCAGGCTTTCTCGATGCGATCGTCATGGATGACGATCAGGCCCTGCTTGAACAACCCGCCAATg                                                                                                                >  1:1184843/1‑96 (MQ=255)
                                                                                          cTTCCGCAACTTCGTACAGGCGTAAGGCCGTCAGGCTTTCTCGATGCGATCGTCATGGATGACGATCAGGCCCTGCTTGAACAACCCGCCAATGGCCTTCTTGAAGTTGCCCTTGCTGACGTTGATCAGCAAGCTGTTCAACgtcagcaa                                                        <  2:474752/150‑8 (MQ=255)
                                                                                          cTTCCGCAACTTCGTACAGGCGTAAGGCCGTCAGGCTTTCTCGATGCGATCGTAATGGATGACGATCTGGCCCTGCTTGAACAACCCGCCAATGGCCTTCTTGAAGTTGCCCTTGCTGACGTTGATCAGCAAGCTGTTCAACgtcagcaa                                                        <  2:1305954/150‑8 (MQ=255)
                                                                                           ttCCGCAACTTCGTACAGGCGTAAGGCCGTCAGGCTTTCTCGATGCGATCGTCATGGATGACGATCAGGCCCTGCTTGAACAACCCGCCAATGGCCTTCTTGAAGTTGCCCTTGCTGACGTTGAACAGCTTGCTGATCAACGCAGGGTCa                                                       >  2:988795/1‑150 (MQ=255)
                                                                                            tCCGCAACTTCGTACAGGCGTAAGGCCGTCAGGCTTTCTCGATGCGATCGTCATGGATGACGATCAGGCCCTGCTTGAACAAcc                                                                                                                        >  1:1860403/1‑84 (MQ=255)
                                                                                            tCCGCAACTTCGTACAGGCGTAAGGCCGTCAGGCTTTCTCGATGCGATCGTCATGGATGACGATCAGGCCCTGCTTGAACAAcc                                                                                                                        <  2:1860403/84‑1 (MQ=255)
                                                                                              cGCAACTTCGTACAGGCGTAAGGCCGTCAGGCTTTCTCGATGCGATCGTCATGGATGACGATCAGGCCCTGCTTGAACAACCCGCCCATGGCCTTCTTGAAGTTGCCCTTGCTGACGTTGAACAGCTTGCTGATCAACGCAGGGTCActc                                                    >  1:400848/1‑150 (MQ=255)
                                                                                                cAACTTCGTACAGGCGTAAGGCCGTCAGGCTTTCTCGATGCGATCGTCATGGATGACGATCAGGCCCTGCTTGAACAACACGCCAATGGCCTTCTTGAAGTTGCCCTTGCTGACGTTGAACAGCTTGCTGATCAACGCAGGGTCACTCtt                                                  <  1:1792541/150‑1 (MQ=255)
                                                                                                  aCTTCGTACAGGCGTAAGGCCGTCAGGCTTTCTCGATGTGAGCGTCATGGATGACGATCAGGCCCTGCTTGAACAACCCGCCAATGGCCTTCTTGAAGTTGCCCTTGCTGACGTTGAACAGCTTGCTGATCAACGCAGGGTCACTCTTGt                                                <  2:564550/150‑1 (MQ=255)
                                                                                                  aCTTCGTACAGGCGTAAGGCCGTCAGGCTTTCTCGATGCGATCGTCATGGATGACGATCAGGCCCTGCTTGAACAACCCGCCAATGGCCTTCTTGAAGTTGCCCTTGCTGACGTTGAACAGCTTGCTGATCAACGCAGGGTCACTCTTGt                                                <  2:2082738/150‑1 (MQ=255)
                                                                                                   cTTCGTACAGGCGTAAGGCCGTCAGGCTTTCTCGATGCGATCGTCATGGATGACGATCAGGCCCTGCTTGAACAACCCGCCAATGGCCTTCTTGAAGTTGCCCTTGCTGTCGTTGAACAGCTTGCTGATCAACGCAGGGTCACTCTTGTc                                               >  2:144071/1‑150 (MQ=255)
                                                                                                   cTTCGTACAGGCGTAAGGCCGTCAGGCTTTCTCGATGCGATCGTCATGGATGACGATCAGGCCCTGCTTGAACAACCCGCCAATGGCCTTCTTGAAGTTGCCCTTGCTGACGTTGAACAGCTTGCTGATCAACGCAGGGTCATTCTTGTc                                               <  1:1528442/150‑1 (MQ=255)
                                                                                                   cTTCGTACAGGCGTAAGGCCGTCAGGCTTTCTCGATGCGATCGTCATGGATGACGATCAGGCCCTGCTTGAACAACCCGCCAATGGCCTTCTTGAAGTTGCCCTTGCTGACGTTGAACAGCTTGCTGATCAACGCAGGGTCACTCTTGTc                                               <  1:2075982/150‑1 (MQ=255)
                                                                                                        tACAGGCGTAAGGCCGTTAGGCTTTCTCGATGCGATCGTCATGGATGACGATCAGGCCCTGCTTGAACAACCCGCAAATGACCTTCTTGAAGTTGCCCTTGCTGACGTTGAACAGCTTGCTGATCAACGCAGGGTCACTCTTGTCGcaca                                          >  1:1980100/1‑150 (MQ=255)
                                                                                                          cAGGCGTAAGGCCGTCAGGCTTTCTCGATGCGATCGTCATGGATGACGATCAGGCCCTGCTTGAACAACCCGCCAATGGCCTTCTTGAAGTTGCCCTTGCTGACGTTGAACAGCTTGCTGATCAACGCAGGGTCCCTCTTGTCGCACAcc                                        >  2:152994/1‑150 (MQ=255)
                                                                                                          cAGGCGTAAGGCCGTCAGGCTTTCTCGATGCGATCGTCATGGATGACGATCAGGCCCTGCTTGAACAACCCGCCAATGGCCTTCTTGAAGTTGCCCTTGCTGACGTTGAACAGCTTGCTGATCAACGCAGGGTCACTCTTGTCGCACAcc                                        <  2:1557222/150‑1 (MQ=255)
                                                                                                           aGGCGTAAGGCCGTCAGGCTTTCTCGATGCGATCGTCATGGATGACGATCAGGCCCTGCTTGAACAACCCGCCAATGGCCTTCTTGAAGTTGCCCTTGCTGACGTTGAACAGCTTGCTGATCAACGCAGGGTCACTCTTGTCGCACAccc                                       >  1:899038/1‑150 (MQ=255)
                                                                                                            ggCGTAAGGCCGTCAGGCTTTCTCGATGCGATCGTCATGGATGACGATCAGGCCCTGCTTGAACAACCCGCCAATGGCCTTCTTGAAGTTGCCCTTGCTGACGTTGAACAGCTTGCTGATCAACGCAGGGTCACTCTTGTCGCACAcccc                                      <  2:880267/150‑1 (MQ=255)
                                                                                                             gCGTAAGGCCGTCAGGCTTTCTCGATGCGATCGTCATGGATGACGATCAGGCCCTGCTTGAACAACCCGCCAATGGCCTTCTTGAAGTTGCCCTTGCTGACGTTGAACAGCTTGCTGATCAACGCAGGGTCACTCTTGTCGCACAccccc                                     <  1:779982/150‑1 (MQ=255)
                                                                                                                tAAGGCCGTCAGGCTTTCTCGATGCGATCGTCATGGATGACGATCAGGCCCTGCTTGAACAACCCGCCAATGGCCTTCTTGAAGTTGCCCTTGCTGACGTTGAACAGCTTGCTGATCAACGCAGGGTCACTCTTGTCGCACACCCCCAac                                  >  1:2194607/1‑150 (MQ=255)
                                                                                                                     ccGTCAGGCTTTCTCGATGCGATCGTCATGGATGACGATCAGGCCCTGCTTGAACAACCCGCCAATGGCCTTCTTGAAGTTGCCCTTGCTGACGTTGAACAGCTTGCTGATCAACGCAGGGTCACTCTTGTCGCACACCCCCAACAcccc                             <  1:988795/150‑1 (MQ=255)
                                                                                                                       gTCAGGCTTTCTCGATGCGATCGTCATGGATGACGATCAGGCCCTGCTTGAACAACCCGCCAATGGCCTTCTTGAAGTTGCCCTTGCTGACGTTGAACAGCTTGCTGATCAACGCAGGGTCACTCTTGTCGCACACCCCCAACACCCCGc                           >  1:359269/1‑150 (MQ=255)
                                                                                                                          aGGCTTTCTCGATGCGATCGTCATGGATGACGATCAGGCCCTGCTTGAACAACCCGCCAATGGCCTTCTTGAAGTTGCCCTTGCTGACGTTGAACAGCTTGCTGATCAACGCAGGGTCACTCTTGTCGCACACCCCCAACACCCCGCCtt                        >  1:1411208/1‑150 (MQ=255)
                                                                                                                             ttttCTCGATGCGATCGTCATGGATGACGATCAGGCCCTGCTTGAACAACCCGCCAATGGCCTTCTTGAAGTTGCCCTTGCTGACGTTGAACAGCTTGCTGATCAACGCAGGGTCACTCTTGTCGCACACCCCCAACACCCCGCCTTCgg                     <  2:357036/149‑1 (MQ=255)
                                                                                                                              tttCTCGATGCGATCGTCATGGATGACGATCAGGCCCTGCTTGAACAACCCGCCAATGGCCTTCTTGAAGTTGCCCTTGCTGACGTTGAACAGCTTGCTGAt                                                                    >  1:1858832/1‑102 (MQ=255)
                                                                                                                              tttCTCGATGCGATCGTCATGGATGACGATCAGGCCCTGCTTGAACAACCCGCCAATGGCCTTCTTGAAGTTGCCCTTGCTGACGTTGAACAGCTTGCTGAt                                                                    <  2:1858832/102‑1 (MQ=255)
                                                                                                                                  tcGATGCGATCGTCATGGATGACGATCAGGCCCTGCTTGAACAACCCGCCAATGGCCTTCTTGAAGTTGCCCTTGCTGACGTTGAACAGCTTGCTGATCAACGCAGGGTCACTCTTGTCGCACACCCCCAACACCCCGCCTTCGGCTTcc                >  1:1744229/1‑150 (MQ=255)
                                                                                                                                    ggcctgATCGTCATCCATGACGATCAGGCCCTGCTTGAACAACCCGCCAATGGCCTTCTTGAAGTTGCCCTTGCTGACGTTGAACAGCTTGCTGATCAACGCAGGGTCACTCTTGTCGCACACCCCCAACACCCCGCCTTCGGCTTCCAg              >  2:193471/6‑150 (MQ=255)
                                                                                                                                        cGATCGTCATGGATGACGATCAGGCCCTGCTTGAACAACCCGCCAATGGCCTTCTTGAAGTTGCCCTTGCTGACGTTGAACAGCTTGCTGATCAACGCAGGGTCACTCTTGTCGCACACCCCCAACACCCCGCCTTCGGCTTCCAGgcgc          <  1:244706/150‑1 (MQ=255)
                                                                                                                                         gATCGTCATGGATGACGATCAGGCCCTGCTTGAACAACCCGCCAATGGCCTTCTTGAAGTTGc                                                                                                <  2:238058/63‑1 (MQ=255)
                                                                                                                                         gATCGTCATGGATGACGATCAGGCCCTGCTTGAACAACCCGCCAATGGCCTTCTTGAAGTTGc                                                                                                >  1:238058/1‑63 (MQ=255)
                                                                                                                                         gATCGTCATGGATGACGATCAGGCCCTGCTTGAACAACCCGCCAATGGCCTTCTTGAAGTTGCCCTTGCTGACGTTGAACAGCTTGCTGATCAACGCAGGGTCACTCTTGTCGCACACCCCCAACACCCCGCCTTCGGCTTCCAGgcgcg         >  1:211639/1‑150 (MQ=255)
                                                                                                                                              tCATGGATGACGATCAGGCCCTGCTTGAACAACCCGCCAATGGCCTTCTTGAAGTTGCCCTTGCTGACGTTGAACAGCTTGCTGATCAACGCAGGGTCACTCTTGTCGCACACCCCCAACACCCCGCCTTCGGCTTCCAGGCGCGCCATg    >  1:90871/1‑150 (MQ=255)
                                                                                                                                                aTGGATGACGATCAGGCCCTGCTTGAACAACCCGCCAATGGCCTTCTTGAAGTTGCCCTTGCTGACGTTGAACAGCTTGCTGATCAACGCAGGGTCACTCTTGTCGCACACCCCCAACACCCCGCCTTCGGCTTCCAGGCGCGCCATGAt  >  2:1891938/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
TGGCAAGCGCCGCTCCACCATCAGCCTGGCCCTCAGCGCCGACCGGGACTCGCTGCTGCGCGACACCCGTCCGGGTGGCGAGCAGCTGGACTTCCGCAACTTCGTACAGGCGTAAGGCCGTCAGGCTTTCTCGATGCGATCGTCATGGATGACGATCAGGCCCTGCTTGAACAACCCGCCAATGGCCTTCTTGAAGTTGCCCTTGCTGACGTTGAACAGCTTGCTGATCAACGCAGGGTCACTCTTGTCGCACACCCCCAACACCCCGCCTTCGGCTTCCAGGCGCGCCATGAT  >  NC_002947/5182825‑5183118

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