Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,485,559 C→A 6.9% A145A (GCC→GCA PP_3972 → short‑chain dehydrogenase/reductase family oxidoreductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,485,5590CA 6.9% 1298.0 / 61.3 389A145A (GCC→GCAPP_3972short‑chain dehydrogenase/reductase family oxidoreductase
Reads supporting (aligned to +/- strand):  ref base C (174/188);  new base A (12/15);  total (186/203)
Fisher's exact test for biased strand distribution p-value = 8.42e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.86e-01

Only 100 of 401 total aligned reads displayed.

ATCATCGCACCCTTCGAACACCGTAGCAGCGGCACAGAGCAAGAAGCTCCCCGTAGGGGATATCCCGTCAGGGGGCCGCTGCGTGGCAAAGTTGTGGTAATCACAGGGGCATCCAGCGGCATAGGCCGCGCCGCCGCGCACGCCTTTGCCTGCAAAGGCGCACGCTTGGTGCTCGCAGCCCGCGATGAAGAAGCATTATTTGACGTTCTGGATGAATGCACCGATTGTGGTACGGATGCCATAGCGGTCACCACCGACGTCACCCACAGCGACCAGGTTCAAGCGCTGGCCGCCCAGGC  >  NC_002947/4485410‑4485708
                                                                                                                                                     |                                                                                                                                                     
 tcatcGCACCCTTCGAACACCGTAGCAGCGGCACAGAGCAAGAAGCTCCCCGTAGGGGATATCCCGTCAGGGGGCCGCTGCGTGGCAAAGTTGTGGTAATCACAGGGGCATCCAGCGGCATAGGCCGCGCCGCCGCGCACGCCTTTGCCt                                                                                                                                                      <  2:2522212/150‑1 (MQ=255)
      gCACCCTTCGAACACCGTAGCAGCGGCACAGAGCAAGAAGCTCCCCGTAGGGGATATCCCGTCAGGGGGCCGCTGCGTGGCAAAGTTGTGGTAATCACAGGGGCATCCAGCGGCATAGGCCGCGCCGCCGCGCACGCCTTTGCCTGCaaa                                                                                                                                                 <  2:5321683/150‑1 (MQ=255)
      gCACCCTTCGAACACCGTAGCAGCGGCACAGAGCAAGAAGCTCCCCGTAGGGGATATCCCGTCAGGGGGCCGCTGCGTGGCAAAGTTGTGGTAATCACAGGGGCATCCAGCGGCATAGGCCGCGCCGCCGCGCACGCCTTTGCCTGCaaa                                                                                                                                                 >  1:294181/1‑150 (MQ=255)
       cACCCGTCGAACACCGTAGCAGCGGCACAGAGCAAGAAGCTCCCCGTAGGGGATATCCCGTCAGGGGGCCGCTGCGTGGCAAAGTTGTGGTAATCACAGGGGCATCCAGCGGCATAGGCCGCGCCGCCGCGCACGCCTTTGCCTGCAAAg                                                                                                                                                <  2:7466442/150‑1 (MQ=255)
        aCCCTTCGAACACCGTAGCAGCGGCACAGAGCAAGAAGCTCCCCGTAGGGGATATCCCGTCAGGGGGCCGCTGCGTGGCAAAGTTGTGGTAATCACAGGGGCATCCAGCGGCATAGGCCGCGCCGCCGCGCACGCCTTTGCCTGCAAAgg                                                                                                                                               >  1:2708523/1‑150 (MQ=255)
        aCCCTTCGAACACCGTAGCAGCGGCACAGAGCAAGAAGCTCCCCGTAGGGGATATCCCGTCAGGGGGCCGCTGCGTGGCAAAGTTGTGGTAATCACAGGGGCATCCAGCGGCATAGGCCGCGCCGCCGCGCACGCCGTTGCCTGccatgg                                                                                                                                               >  1:4489708/1‑145 (MQ=255)
         cccTTCGAACACCGTAGCAGCGGCACAGAGCAAGAAGCTCCCCGTAGGGTATATCCCGTCAGGGGGCCGCTGCGTGGCAAAGTTGTGGTAATCACAGGGGCATCCAGCGGCATAGGCCGCGCCGCCGCGCACGCCTTTGCCTGCAAAGgc                                                                                                                                              <  1:3904789/150‑1 (MQ=255)
         cccTTCGAACACCGTAGCAGCGGCACAGAGCAAGAAGCTCCCCGTAGGGGATATCCCGTCAGGGGGCCGCTGCGTGGCAAAGTTGTGGTAATCACAGGGGCATCCAGCGGCATAGGCCGCGCCGCCGCGCACGCCTTTGCAGGCAAAGgc                                                                                                                                              <  2:8065589/150‑1 (MQ=255)
           cTTCGAACACCGTAGCAGCGGCACAGAGCAAGAAGCTCCCCGTAGGGGATATCCCGTCAGGGGGCCGCTGCGTGGCAAAGTTGTGGTAATCACAGGGGCATCCAGCGGCATAGGCCGCGCCGCCGCGCACGCCTTTGCCTGCAAAGgcgc                                                                                                                                            <  1:6515407/150‑1 (MQ=255)
           cTTCGAACACCGTAGCAGCGGCACAGAGCAAGAAGCTCCCCGTAGGGGATATCCCGTCAGGGGGCCGCTGCGTGGCAAAGTTGTGGTAATCACAGGGGCATCCAGCGGCATAGGCCGCGCCGCCGCGCACGCCTTTGCAGGCAAAGgcgt                                                                                                                                            <  1:741273/150‑2 (MQ=255)
           cTTCGAACACCGTAGCAGCGGCACAGAGCAAGAAGCTCCCCGTAGGGGATATCCCGTCAGGGGGCCGCTGCGTGGCAAAGTTGTGGTAATCACAGGGGCATCCAGCGGCATAGGCCGCGCCGCCGCGCACGCCTTTGCAGGCAAAGgcgt                                                                                                                                            <  1:6978768/150‑2 (MQ=255)
            ttCGAACACCGTAGCAGCGGCACAGAGCAAGAAGCTCCCCGTAGGGGATATCCCGTCAGGGGGCCGCTGCGTGGCAAAGTTGTGGTAATCACAGGGGCATCCAGCGGCATAGGCCGCGCCGCCGCGCACGCCTTTGCCTGCAAAGGCGCa                                                                                                                                           >  1:1771901/1‑150 (MQ=255)
            ttCGAACACCGTAGCAGCGGCACAGAGCAAGAAGCTCCCCGTAGGGGATATCCCGTCAGGGGGCCGCTGCGTGGCAAAGTTGTGGTAATCACAGGGGCATCCAGCGGCATAGGCCGCGCCGCCGCGCACGCCTTTGCCTGCAAAGGCGCa                                                                                                                                           >  1:5409034/1‑150 (MQ=255)
            ttCGAACACCGTAGCAGCGGCACAGAGCAAGAAGCTCCCCGTAGGGGATATCCCGTCAGGGGGCCGCTGCGTGGCAAAGTTGTGGTAATCACAGGGGCATCCAGCGGCATAGGCCGCGCCGCCGCGCACGCCTTTGCCTGCAAAGGCGCa                                                                                                                                           >  2:2837077/1‑150 (MQ=255)
            ttCGAACACCGTAGCAGCGGCACAGAGCAAGAAGCTCCCCGTAGGGGATATCCCGTCAGGGGGCCGCTGCGTGGCAAAGTTGTGGTAATCACAGGGGCATCCAGCGGCATAGGCCGCGCCGCCGCGCACGCCTTTGCCTGCAAAGGCGCa                                                                                                                                           >  2:8228206/1‑150 (MQ=255)
            ttCGAACACCGTAGCAGCGGCACAGAGCAAGAAGCTCCCCGTAGGGGATATCCCGTCAGGGGGCCGCTGCGTGGCAAAGTTGTGGTAATCACAGGGGCATCCAGCGGCATAGGCCGCGCCGCCGCGCACGCCTTTGCAGGCAAAGgcgtg                                                                                                                                           <  2:1386802/150‑3 (MQ=255)
               gAACACCGTAGCAGCGGCACAGAGCAAGAAGCTCCCCGTAGGGGATATCCCGTCAGGGGGCCGCTGCGTGGCAAAGTTGTGGTAATCACAGGGGCATCCAGCGGCATAGGCCGCGCCGCCGCGTGCGCCTTTGCAGGCAAAGgcgtgcgc                                                                                                                                        <  2:683785/150‑6 (MQ=255)
               gAACACCGTAGCAGCGGCACAGAGCAAGAAGCTCCCCGTAGGGGATATCCCGTCAGGGGGCCGCTGCGTGGCAAAGTTGTGGTAATCACAGGGGCATCCAGCGGCATAGGCCGCGCCGCCGCGCACGCCTTTGCCTGCAAAGGCGCACGc                                                                                                                                        >  2:1581861/1‑150 (MQ=255)
                aaCACCGTAGCGGCGGCACAGAGCAAGAAGCTCCCCGTAGGGGATATCCCGTCAGGGGGCCGCTGCGTGGCAAAGTTGTGGTAATCACAGGGGCATCCAGCGGCATAGGCCGCGCCGCCGCGCACGCCTTTGCCTGCAAAGGCGCACGCt                                                                                                                                       >  1:4932033/1‑150 (MQ=255)
                aaCACCGTAGCAGCGGCACAGAGCAAGAAGCTCCCCGTAGGGGATATCCCGTCAGGGGGCCGCTGCGTGGCAAAGTTGTGGTAATCACAGGGGCATCCAGCGGCATAGGCCGCGCCGCCGCGCGCGCCTTTGCAGGCAAAGgcgtg                                                                                                                                           <  1:8013167/146‑3 (MQ=255)
                 acacCGTAGCAGCGGCACAGAGCAAGAAGCTCCCCGTAGGGGATATCCCGTCAGGGGGCCGCTGCGTGGCAAAGTTGTGGTAATCACAGGGGCATCCAGCGGCATAGGCCGCGCCGCCGCGCACGCCTTTGCCTGCAAAGGCGCACGCtt                                                                                                                                      >  1:8770734/1‑150 (MQ=255)
                  cacCGTAGCAGCGGCACAGAGCAAGAAGCTCCCCGTAGGGGATATCCCGTCAGGGGGCCGCTGCGTGGCAAAGTTGTGGTAATCACAGGGGCATCCAGCGGCATAGGCCGCGCCGCCGCGCACGCCTTTGCCTGCAAAGGCGCACGCTTg                                                                                                                                     <  1:8608431/150‑1 (MQ=255)
                         gcagcGGCACAGAGCAAGAAGCTCCCCGTAGGGGATATCCCGTCAGGGGGCCGCTGCGTGGCAAAGTTGTGGTAATCACAGGGGCATCCAGCGGCATAGGCCGCGCCGCCGCGCACGCCTTTGCCTGCAAAGGCGCACGCTTGGTGCTCg                                                                                                                              <  1:6460894/150‑1 (MQ=255)
                          cagcGGCACAGAGCAAGAAGCTCCCCGTAGGGGATATCCCGTCAGGGGGCCGCTGCGTGGCAAAGTTGTGGTAATCACAGGGGCATCCAGCGGCATAGGCCGCGCCGCCGCGCACGCCTTTGCCTGCAAAGGCGCACGCTTGGTGCTCGc                                                                                                                             >  2:6171322/1‑150 (MQ=255)
                           agcGGCACAGAGCAAGAAGCTCCCCGTAGGGGATATCCCGTCAGGGGGCCGCTGCGTGGCAAAGTTGTGGTAATCACAGGGGCATCCAGCGGCATAGGCCGCGCCGCCGCGCACGCCTTTGCCTGCAAAGGCGCACGCTTGGTGCTCGCa                                                                                                                            >  1:1310824/1‑150 (MQ=255)
                              ggCACAGAGCAAGAAGCTCCCCGTAGGGGATATCCCGTCAGGGGGCCGCTGCGTGGCAAAGTTGTGGTAATCACAGGGGCATCCAGCGGCATAGGCCGCGCCGCCGCGCACGCCTTTGCCTGCAAAGGCGCACGCTTGGTGCTCGCAGcc                                                                                                                         >  1:7038767/1‑150 (MQ=255)
                                   agagCAAGAAGCTCCCCGTAGGGGATATCCCGTCAGGGGGCCGCTGCGTGGCAAAGTTGTGGTAATCACAGGGGCATCCAGCGGCATAGGCCGCGCCGCCGCGCACGCCTTTGCCTGCAAAGGCGCACGCTTGGTGCTCGCAGCCCGCGa                                                                                                                    <  1:5123762/150‑1 (MQ=255)
                                   agagCAAGAAGCTCCCCGTAGGGGATATCCCGTCAGGGGGCCGCTGCGTGGCAAAGTTGTGGTAATCACAGGGGCATCCAGCGGCATAGGCCGCGCCGCCGCGCACGCCTTTGCCTGCAAAGGCGCACGCTTGGTGCTCGCAGCCCGCGa                                                                                                                    <  1:3303241/150‑1 (MQ=255)
                                    gagCAAGAAGCTCCCCGTAGGGGATATCCCGTCAGGGGGCCGCTGCGTGGCAAAGTTGTGGTAATCACAGGGGCATCCAGCGGCATAGGCCGCGCCGCCGCGCACGCCTTTGCCTGCAAAGGCGCACGCTTGGTGCTCGCAGCCCGCGAt                                                                                                                   >  1:2014080/1‑150 (MQ=255)
                                     agCAAGAAGCTCCCCGTAGGGGATATCCCGTCAGGGGGCCGCTGCGTGGCAAAGTTGTGGTAATCACAGGGGCATCCAGCGGCATAGGCCGCGCCGCCGCGCACGCCTTTGCCTGCAAAGGCGCACGCTTGGTGCTCGCAGCCCGCGATg                                                                                                                  >  1:4110488/1‑150 (MQ=255)
                                       cAAGAAGCTCCCCGTAGGGGATATCCCGTCAGGGGGCCGCTGCGTGGCAAAGTTGTGGTAATCACAGGGGCATCCAGCGGCATAGGCCGCGCCGCCGCGCACGCCTTTGCCTGCAAAGGCGCACGCTTGGTGCTCGCAGCCCGCGATgaa                                                                                                                <  2:8760838/150‑1 (MQ=255)
                                        aagaagCTCCCCGTAGGGGATATCCCGTCAGGGGGCCGCTGCGTGGCAAAGTTGTGGTAATCACAGGGGCATCCAGCGGCATAGGCCGCGCCGCCGCGCACGCCTTTGCCTGCAAAGGCGCACGCTTGGTGCTCGCAGCCCGCGATgaag                                                                                                               >  1:7639820/1‑150 (MQ=255)
                                        aagaagCTCCCCGTAGGGGATATCCCGTCAGGGGGCCGCTGCGTGGCAAAGTTGTGGTAATCACAGGGGCATCCAGCGGCATAGGCCGCGCCGCCGCGCACGCCTTTGCCTGCAAAGGCGCACGCTTGGTGCTCGCAGCCCGCGATgaag                                                                                                               >  2:2180017/1‑150 (MQ=255)
                                         agaagCTCCCCGTAGGGGATATCCCGTCAGGGGGCCGCTGCGTGGCAATGTTGTGGTAATCACAGGGGCATCCAGCGGCATAGGCCGCGCCGCCGCGCACGCCTTTGCCTGCAAAGGCGCACGCTTGGTGCTCGCAGCCCGCGATgaaga                                                                                                              <  2:6114782/150‑1 (MQ=255)
                                            agCTCCCCGTAGGGGATATCCCGTCAGGGGGCCGCTGCGTGGCAAAGTTGTGGTAATCACAGGGGCATCCAGCGGCATAGGCCGCGCCGCCGCGCACGCCTTTGCCTGCAAAGGCGCACGCTTGGTGCTCGCAGCCCGCGATGAAGAAGc                                                                                                           <  1:3769780/150‑1 (MQ=255)
                                               tCCCCGTAGGGGATATCCCGTCAGGGGGCCGCTGCGTGGCAAAGTTGTGGTAATCACAGGGGCATCCAGCGGCATAGGCCGCGCCGCCGCGCACGCCTTTGCCTGCAAAGGCGCACGCTTGGTGCTCGCAGCCCGCGATGAAGAAGCatt                                                                                                        >  1:4290737/1‑150 (MQ=255)
                                               tCCCCGTAGGGGATATCCCGTCAGGGGGCCGCTGCGTGGCAAAGTTGTGGTAATCACAGGGGCATCCAGCGGCATAGGCCGCGCCGCCGCGCACGCCTTTGCCTGCAAAGGCGCACGCTTGGTGCTCGCAGCCCGCGATGAAGAAGCatt                                                                                                        >  1:153338/1‑150 (MQ=255)
                                                 cccGTAGGGGATATCCCGTCAGGGGGCCGCTGCGTGGCAACGTTGTGGTAATCACAGGGGCATCCAGCGGCATAGGCCGCGCCGCCGCGCACGCCTTTGCCTGCACAGGCGCACGCTTGGTGCTCGCAGCCCGCGATGAAGAAGCattct                                                                                                      >  2:7591834/1‑148 (MQ=255)
                                                 cccGTAGGGGATATCCCGTCAGGGGGCCGCTGCGTGGCAAAGTTGTGGTAATCACAGGGGCATCCAGCGGCATAGGCCGCGCCGCCGCGCACGCCTTTGCCTGCAAAGGCGCACGCTTGGTGCTCGCAGCCCGCGATGAAGAAGCattat                                                                                                      <  2:3682901/150‑1 (MQ=255)
                                                  ccGTAGGGGATATCCCGTCAGGGGGCCGCTGCGTGGCAAAGTTGTGGTAATCACAGGGGCATCCAGCGGCATAGGCCGCGCCGCCGCGCACGCCTTTGCCTGCAAAGGCGCACGCTTGGTGCTCGCAGCCCGCGATGAAGAAGCattatt                                                                                                     >  1:1403819/1‑150 (MQ=255)
                                                   cGTAGGGGATATCCCGTCAGGGGGCCGCTGCGTGGCAAAGTTGTGGTAATCACAGGGGCATCCAGCGGCATAGGCCGCGCCGCCGCGCACGCCTTTGCCTGCAAAGGCGCACGCTTGGTGCTCGCAGCCCGCGATGAAGAAGCATTAttt                                                                                                    <  2:6516497/150‑1 (MQ=255)
                                                          gATATCCCGTCAGGGGGCCGCTGCGTGGCAAAGTTGTGGTAATCACAGGGGCATCCAGCGGCATAGGCCGCGCCGCCGCGCACGCCTTTGCCTGCAAAGGCGCACGCTTGGTGCTCGCAGCCCGCGATGAAGAAGCATTATTTGACGTTc                                                                                             <  2:2061808/150‑1 (MQ=255)
                                                          gATATCCCGTCAGGGGGCCGCTGCGTGGCAAAGTTGTGGTAATCACAGGGGCATCCAGCGGCATAGGCCGCGCCCCCGCGCACGCCTTTGCCTGCAAGGGCGCACGCTTGGTGCTCGCAGCCCGCGATGAAGAAGCATTATTTGACGTTc                                                                                             <  2:4322686/150‑1 (MQ=255)
                                                           atatCCCGTCAGGGGGCCGCTGCGTGGCAAAGTTGTGGTAATCACAGGGGCATCCAGCGGCATAGGCCGCGCCGCCGCGCACGCCTTTGCCTGCAAAGGCGCACGCTTGGTGCTCGCAGCCCGCGATGAAGAAGCATTATTTGACGTTCt                                                                                            <  2:6312591/150‑1 (MQ=255)
                                                              tCCCGTCAGGGGGCCGCTGCGTGGCAAAGTTGTGGTAATCACAGGGGCATCCAGCGGCATAGGCCGCGCCGCCGCGCACGCCTTTGCCTGCAAAGGCGCACGCTTGGTGCTCGCAGCCCGCGATGAAGAAGCATTATTTGACGTTCTGGa                                                                                         <  1:3139109/150‑1 (MQ=255)
                                                                   tCAGGGGGCCGCTGCGTGGCAAAGTTGTGGTAATCACAGGGGCATCCAGCGGCATAGGCCGCGCCGCCGCGCACGCCTTTGCCTGCAAAGGCGCACGCTTGGTGCTCGCAGCCCGCGATGAAGAAGCATTATTTGACGTTCTGGATGAAt                                                                                    >  2:1780716/1‑150 (MQ=255)
                                                                       ggggCCGCTGCGTGGCAAAGTTGTGGTACTCACAGGGGCATCCAGCGGCATAGGCCGCGCCGCCGCGCACGCCTTTGCCTGCAAAGGCGCACGCTTGGTGCTCGCAGCCCGCGATGAAGAAGCATTATTTGACGTTCTGGATGAATGCAc                                                                                >  2:3232527/1‑150 (MQ=255)
                                                                       ggggCCGCTGCGTGGCAAAGTTGTGGTAATCACAGGGGCATCCAGCGGCATAGGCCGCGCCGCCGCGCACGCCTTTGCCTGCAAAGGCGCACGCTTGGTGCTCGCAGCCCGCGATGAAGAAGCATTATTTGACGTTCTGGATGAATGCAc                                                                                >  1:6862221/1‑150 (MQ=255)
                                                                       ggggCCGCTGCGTGGCAAAGTTGTGGTAATCACAGGGGCATCCAGCGGCATAGGCCGCGCCGCCGCGCACGCCTTTGCCTGCAAAGGCGCACGCTTGGTGCTCGCAGCCCGCGATGAAGAAGCATTATTTGACGTTCTGGATGAATGCAc                                                                                >  2:7967179/1‑150 (MQ=255)
                                                                       ggggCCGCTGCGTGGCAAAGTTGTGGTAATCACAGGGGCATCCAGCGGCATAGGCCGCGCCGCCGCGCACGCCTTTGCCTGCAAAGGCGCACGCTTGGTGCTCGCAGCCCGCGATGAAGAAGCATTATTTGACGTTCTGGATGAATGCAc                                                                                <  1:2870598/150‑1 (MQ=255)
                                                                            cGCTGCGTGGCAAAGTTGTGGTAATCACAGGGGCATCCAGCGGCATAGGCCGCGCCGCCGCGCACGCCTTTGCCTGCAAAGGCGCACGCTTGGTGCTCGCAGCCCGCGATGAAGAAGCATTATTTGACGTTCTGGATGAATGCACCGAtt                                                                           <  1:7163751/150‑1 (MQ=255)
                                                                              cTGCGTGGCAAAGTTGTGGTAATCACAGGGGCATCCAGCGGCATAGGCCGCGCCGCCGCGCACGCCTTTGCCTGCAAAGGCGCACGCTTGGTGCTCGCAGCCCGCGATGAAGAAGCATTATTTGACGTTCTGGATGAATGCACCGATtgt                                                                         <  2:5051677/150‑1 (MQ=255)
                                                                              cTGCGTGGCAAAGTTGTGGTAATCACAGGGGCATCCAGCGGCATAGGCCGCGCCGCCGCGCACGCCTTTGCCTGCAAAGGCGCACGCTTGGTGCTCGCAGCCCGCGATGAAGAAGCATTATTTGACGTTCTGGATGAATGCACCGATtgt                                                                         <  2:4482277/150‑1 (MQ=255)
                                                                               tGCGTGGCAAAGTTGTGGTAATCACAGGGGCATCCAGCGGCATAGGCCGCGCCGCCGCGCACGCCTTTGCCTGCAAAGGCGCACGCTTGGTGCTCGCAGCCCGCGATGAAGAAGCATTATTTGACGTTCTGGATGAATGCACCGATtgtg                                                                        <  2:8402550/150‑1 (MQ=255)
                                                                                 cGTGGCAAAGTTGTGGTAATCACAGGGGCATCCAGCGGCATAGGCCGCGCCGCCGCGCACGCCTTTGCCTGCAAAGGCGCACGCTTGGTGCTCGCAGCCCGCGATGAAGAAGCATTATTTGACGTTCTGGATGAATGCACCGATTGTGGt                                                                      <  2:4227654/150‑1 (MQ=255)
                                                                                  gTGGCAAAGTTGTGGTAATCACAGGGGCATCCAGCGGCATAGGCCGCGCCGCCGCGCACGCCTTTGCCTGCAAAGGCGCACGCTTGGTGCTCGCAGCCCGCGATGAAGAAGCATTATTTGACGTTCTGGATGAATGCACCGATTGTGGTa                                                                     <  1:7910033/150‑1 (MQ=255)
                                                                                  gTGGCAAAGTTGTGGTAATCACAGGGGCATCCAGCGGCATAGGCCGCGCCGCCGCGCACGCCTTTGCCTGCAAAGGCGCACGCTTGGTGCTCGCAGCCCGCGATGAAGAAGCATTATTTGACGTTCTGGATGAATGCACCGATTGTGGTa                                                                     <  1:5459385/150‑1 (MQ=255)
                                                                                             gtgGTAATCACAGGGGCATCCAGCGGCATAGGCCGCGCCGCCGCGCACGCCTTTGCCTGCAAAGGCGCACGCTTGGTGCTCGCAGCCCGCGATGAAGAAGCATTATTTGACGTTCTGGATGAATGCACCGATTGTGGTACGGATGCCATa                                                          >  1:3672853/1‑150 (MQ=255)
                                                                                                gTAATCACAGGGGCATCCAGCGGCATAGGCCGCGCCGCCGCGCACGCCTTTGCCTGCAAAGGCGCACGCTTGGTGCTCGCAGCCCGCGATGAAGAAGCATTATTTGACGTTCTGGATGAATGCACCGATTGTGGTACGGATGCCATAGCg                                                       >  2:3855456/1‑150 (MQ=255)
                                                                                                 tAATCACAGGGGCATCCAGCGGCATAGGCCGCGCCGCCGCGCACGCCTTTGCCTGCAAAGGCGCACGCTTGGTGCTCGCAGCCCGCGATGAAGAAGCATTATTTGACGTTCTGGATGAATGCACCGATTGTGGTACGGATGCCATAGCgg                                                      >  2:5197241/1‑150 (MQ=255)
                                                                                                    tCACAGGGGCATCCAGCGGCATAGGCCGCGCCGCCGCGCACGCCTTTGCCTGCAAAGGCGCACGCTTGGTGCTCGCAGCCCGCGATGAAGAAGCATTATTTGACGTTCTGGATGAATGCACCGATTGTGGTACGGATGCCATAGCGGTca                                                   >  1:6148250/1‑150 (MQ=255)
                                                                                                      acaGGGGCATCCAGCGGCATAGGCCGCGCCGCCGCGCACGCCTTTGCCTGCAAAGGCGCACGCTTGGTGCTCGCAGCCCGCGATGAAGAAGCATTATTTGACGTTCTGGATGAATGCACCGATTGTGGTACGGATGCCATAGCGGTcacc                                                 >  1:4808321/1‑150 (MQ=255)
                                                                                                        aGGGGCATCCAGCGGCATAGGCCGCGCCGCCGCGCACGCCTTTGCCTGCAAAGGCGCACGCTTGGTGCTCGCAGCCCGCGATGAAGAAGCATTATTTGACGTTCTGGATGAATGCACCGATTGTGGTACGGATGCCATAGCGGTcaccac                                               <  1:1170965/150‑1 (MQ=255)
                                                                                                         ggggCATCCAGCGGCATAGGCCGCGCCGCCGCGCACGCCTTTGCCTGCAAAGGCGCCCGCTTGGGGCTCGCAGCCCGCGATGAAGAAGCATTATTTGACGTTCTGGATGAATGCACCGATTGTGGTACGGATGCCATAGCGGTCaccacc                                              <  2:335149/150‑1 (MQ=255)
                                                                                                          gggCATCCAGCGGCATAGGCCGCGCCGCCGCGCACGCCTTTGCCTGCAAAGGCGCACGCTTGGTGCTCGCAGCCCGCGATGAAGAAGCATTATTTGACGTTCTGGATGAATGCACCGATTGTGGTACGGATGCCATAGCGGTCACCACCg                                             >  1:101929/1‑150 (MQ=255)
                                                                                                          gggCATCCAGCGGCATAGGCCGCGCCGCCGCGCACGCCTTTGCCTGCAAAGGCGCACGCTTGGTGCTCGCAGCCCGCGATGAAGAAGCATGATTTGACGTTCTGGATGAATGCACCGATTGTGGTACGGATGCCATAGCGGTCACCACCg                                             >  2:7744862/1‑150 (MQ=255)
                                                                                                            gCATCCAGCGGCATAGGCCGCGCCGCCGCGCACGCCTTTGCCTGCAAAGGCGCACGCTTGGTGCTCGCAGCCCGCGATGAAGAAGCATTATTTGACGTTCTGGATGAATGCACCGATTGTGGTACGGATGCCATAGCGGTCACCACCGAc                                           <  2:1973712/150‑1 (MQ=255)
                                                                                                            gCATCCAGCGGCATAGGCCGCGCCGCCGCGCACGCCTTTGCCTGCAAAGGCGCACGCTTGGTGCTCGCAGCCCGCGATGAAGAAGCATTATTTGACGTTCTGGATGAATGCACCGATTGTGGTACGGATGCCATAGCGGTCACCACCGAc                                           <  1:5947763/150‑1 (MQ=255)
                                                                                                              aTCCAGCGGCATAGGCCGCGCCGCCGCGCACGCCTTTGCCTGCAAAGGCGCACGCTTGGTGCTCGCAGCCCGCGATGAAGAAGCATTATTTGACGTTCTGGATGAATGCACCGATTGTGGTACGGATGCCATAGCGGTCACCACCGACGt                                         >  1:8692297/1‑150 (MQ=255)
                                                                                                               tCCAGCGGCATAGGCCGCGCCGCCGCGCACGCCTTTGCCTGCAAAGGCGCACGCTTGGTGCTCGCAGCCCGCGATGAAGAAGCATTATTTGACGTTCTGGATGAATGCACCGATTGTGGTACGGATGCCATAGCGGTCACCACCGACGTc                                        >  1:5830123/1‑150 (MQ=255)
                                                                                                                ccAGCGGCATAGGCCGCGCCGCCGCGCACGCCTTTGCCTGCAAAGGCGCACGCTTGGTGCTCGCAGCCCGCGATGAAGAAGCATTATTTGACGTTCTGGATGAATGCACCGATTGTGGTACGGATGCCATAGCGGTCACCACCGACGTCa                                       >  1:8523614/1‑150 (MQ=255)
                                                                                                                 cAGCGGCATAGGCCGCGCCGCCGCGCACGCCTTTGCCTGCAAAGGCGCACGCTTGGTGCTCGCAGCCCGCGATGAAGAAGCATTATTTGACGTTCTGGATGAATGCACCGATTGTGGTACGGATGCCATAGCGGTCACCACCGACGTCAc                                      <  1:6727345/150‑1 (MQ=255)
                                                                                                                  aGCGGCATAGGCCGCGCCGCCGCGCACGCCTTTGCCTGCAAAGGCGCACGCTTGGTGCTCGCAGCCCGCGATGAAGAAGCATTATTTGACGTTCTGGATGAATGCACCGATTGTGGTACGGATGCCATAGCGGTCACCACCGACGTCAcc                                     >  1:2507266/1‑150 (MQ=255)
                                                                                                                      gCATAGGCCGCGCCGCCGCGCACGCCTTTGCCTGCAAAGGCGCACGCTTGGTGCTCGCAGCCCGCGATGAAGAAGCATTATTTGACGTTCTGGATGAATGCACCGATTGTGGTACGGATGCCATAGCGGTCACCACCGACGTCACCcaca                                 <  1:1471312/150‑1 (MQ=255)
                                                                                                                       cATAGGCCGCGCCGCCGCGCACGCCTTTGCCTGCAAAGGCGCACGCTTGGTGCTCGCAGCCCGCGATGAAGAAGCATTATTTGACGTTCTGGATGAATGCACCGATTGTGGTACGGATGCCATAGCGGTCACCACCGACGTCACCCACAg                                <  2:2874149/150‑1 (MQ=255)
                                                                                                                        aTAGGCCGCGCCGCCGCGCACGCCTTTGCCTGCAAAGGCGCACGCTTGGTGCTCGCAGCCCGCGATGAAGAAGCATTATTTGACGTTCTGGATGAATGCACCGATTGTGGTACGGATGCCATAGCGGTCACCACCGACGTCACCCACAGc                               >  1:6553735/1‑150 (MQ=255)
                                                                                                                          aGGCCGCGCCGCCGCGCACGCCTTTGCCTGCAAAGGCGCACGCTTGGTTCTCGCAGCCCGCGATGAAGAAGCATTATTTGACGTTCTGGATGAATGCACCGATTGTGGTACGGATGCCATAGCGGTCACCACCGACGTCACCCACAGCGa                             <  2:660455/150‑1 (MQ=255)
                                                                                                                          aGGCCGCGCCGCCGCGCACGCCTTTGCCTGCAAAGGCGCACGCTTGGTGCTCGCAGCCCGCGATGAAGAAGCATTATTTGACGTTCTGGATGAATGCACCGATTGTGGTACGGATGCCATAGCGGTCACCACCGACGTCACCCACAGCGa                             <  2:584314/150‑1 (MQ=255)
                                                                                                                            gccgcgccgcCGCGCACGCCTTTGCCTGCAAAGGCGCACGCTTGGTGCTCGCAGCCCGCGATGAAGAAGCAGTATTTGACGTTCTGGATGAATGCACCGATTGTGGTCCGGATGCCATAGCGGTCACCACCGACGTCACCCACAGCGAcc                           >  2:1158142/1‑150 (MQ=255)
                                                                                                                              cgcgccgcCGCGCACGCCTTTGCCTGCAAAGGCGCACGCTTGGTGCTCGCAGCCCGCGATGAAGAAGCATTATTTGACGTTCTGGATGAATGCACCGATTGTGGTACGGATGCCATAGCGGTCACCACCGACGTCACCCACAGCGACCAg                         <  1:2240364/150‑1 (MQ=255)
                                                                                                                                   cgccgcGCACGCCTTTGCCTGCAAAGGCGCACGCTTGGTGCTCGCAGCCCGCGATGAAGAAGCATTATTTGACGTTCTGGATGAATGCACCGATTGTGGTACGGATGCCATAGCGGTCACCACCGACGTCACCCACAGCGACCAGGTTCa                    >  1:2081376/1‑150 (MQ=255)
                                                                                                                                    gccgcGCACGCCTTTGCCTGCAAAGGCGCACGCTTGGTGCTCGCAGCCCGCGATGAAGAAGCATTATTTGACGTTCTGGATGAATGCACCGATTGTGGTACGGATGCCATAGCGGTCACCACCGACGTCACCCACAGCGACCAGGTTCaa                   <  2:7206226/150‑1 (MQ=255)
                                                                                                                                     ccgcGCACGCCTTTGCCTGCAAAGGCGCACGCTTGGTGCTCGCAGCCCGCGATGAAGAAGCATTATTTGACGTTCTGGATGAATGCACCGATTGTGGTACGGATGCCATAGCGGTCACCACCGACGTCACCCACAGCGACCAGGTTCAAg                  <  2:6067744/150‑1 (MQ=255)
                                                                                                                                      cgcgCACGCCTTTGCCTGCAAAGGCGCACGCTTGGTGCTCGCAGCCCGCGATGAAGAAGCATTATTTGACGTTCTGGATGAATGCACCGATTGTGGTACGGATGCCATAGCGGTCACCACCGACGTCACCCACAGCGACCAgggagagct                 <  2:3509631/150‑8 (MQ=255)
                                                                                                                                      cgcgCACGCCTTTGCCTGCAAAGGCGCACGCTTGGTGCTCGCAGCCCGCGATGAAGAAGCATTATTTGACGTTCTGGATGAATGCACCGATTGTGGTACGGATGCCATAGCGGTCACCACCGACGTCACCCACAGCGACCAGGTTCAAgc                 >  1:7950520/1‑150 (MQ=255)
                                                                                                                                         gcACGCCTTTGCCTGCAAAGGCGCACGCTTGGTGCTCGCAGCCCGCGATGAAGAAGCATTATTTGACGTTCTGGATGAATGCACCGATTGTGGTACGGATGCCATAGCGGTCACCACCGACGTCACCCACAGCGACCAGGTTCAAGCGCt              <  1:8228206/150‑1 (MQ=255)
                                                                                                                                          tgcGCCTTTGCAGGCAAAGGCGTGCGCTTGGTGCTCGCAGCCCGCGATGAAGAAGCATTATTTGACGTTCTGGATGAATGCACCGATTGTGGTACGGATGCCATAGCGGTCACCACCGACGTCACCCACAGCGACCAGGTTCAAGCGCTg             >  2:1017882/3‑150 (MQ=255)
                                                                                                                                          agatccTTTGCAGGCAAAGGCGTGCGCTTGGTGCTCGAAGACATCGATTAAGAAGAATAAATTGACGTTCTGTATTAATGCACCGATTGTGGTACGGATTCCATAGCGGTCACCACCGACGTCACCCACAGCGAACAGGTACAAGCGCTg             >  2:8798681/5‑150 (MQ=255)
                                                                                                                                             gCCTTTGCCTGCAAAGGCGCACGCTTGGTGCTCGCAGCCCGCGATGAAGAAGCATTATTTGACGTTCTGGATGAATGCACCGATTGTGGTACGGATGCCATAGCGGTCACCACCGACGTCACCCACAGCGACCAGGTTCAAGCGCTGgcc          >  1:6254684/1‑150 (MQ=255)
                                                                                                                                               cTTTGCCTGCAAAGGCGCACGCTTGGTGCTCGCAGCCCGCGATGAAGAAGCATTATTTGACGTTCTGGATGAATGCACCGATTGTGGTACGGATGCCATAGCGGTCACCACCGACGTCACCCACAGCGACCAGGTTCAAGCGCTGgccgc        >  2:4925968/1‑150 (MQ=255)
                                                                                                                                               cTTTGCCTGCAAAGGCGCACGCTTGGTGCTCGCAGCCCGCGATGAAGAAGCATTATTTGACGTTCTGGATGAATGCACCGATTGTGGTACGGATGCCATAGCGGTCACCACCGACGTCACCCACAGCGACCAGGTTCAAGCGCTGgccgc        >  2:2447194/1‑150 (MQ=255)
                                                                                                                                               cTTTGCAGGCAAAGGCGCACGCTTGGTGCTCGCAGCCCGCGATGAAGAAGCATTATTTGACGTTCTGGATGAATGCACCGATTGTGGTACGGATGCCATAGCGGTCACCACCGACGTCACCCACAGCGACCAGGTTCAAGCGCTGgccgc        >  1:6341900/1‑150 (MQ=255)
                                                                                                                                               cTTTGCAGGCAAAGGCGCACGCTTGGTGCTCGCAGCCCGCGATGAAGAAGCATTATTTGACGTTCTGGATGAATGCACCGATTGTGGTACGGATGCCATAGCGGTCACCACCGACGTCACCCACAGCGACCAGGTTCAAGCGCTGgccgc        >  1:3744039/1‑150 (MQ=255)
                                                                                                                                                tttGCCTGCAAAGGCGCACGCTTGGTGCTCGCAGCCCGCGATGAAGAAGCATTATTTGACGTTCTGGATGAATGCACCGATTGTGGTACGGATGCCATAGCGGTCACCACCGACGTCACCCACAGCGACCAGGTTCAAGCGCTGgccgcc       >  2:8634747/1‑150 (MQ=255)
                                                                                                                                                 ttGCCTGCAAAGGCGCACGCTTGGTGCTCGCAGCCCGCGATGAAGAAGCATTATTTGACGTTCTGGATGAATGCACCGATTGTGGTACGGATGCCATAGCGGTCACCACCGACGTCACCCACAGCGACCAGGTTCAAGCGCTGGCCGccc      <  2:2708523/150‑1 (MQ=255)
                                                                                                                                                  tGCCTGCAAAGGCGCACGCTTGGTGCTCGCAGCCCGCGATGAAGAAGCATTATTTGACGTTCTGGATGAATGCACCGATTGTGGTACGGATGCCATAGCGGTCACCACCGACGTCACCCACAGCGACCAGGTTCAAGCGCTGGCCGCCCa     <  2:2977398/150‑1 (MQ=255)
                                                                                                                                                    ccTGCAAAGGCGCACGCTTGGTGCTCGCAGCCCGCGATGAAGAAGCATTATTTGACGTTCTGGATGAATGCACCGATTGTGGTACGGATGCCATAGCGGTCACCACCGACGTCACCCACAGCGACCAGGTTCAAGCGCTGGCCGCCCAgg   <  2:4600217/150‑1 (MQ=255)
                                                                                                                                                     cTGCAAAGGCGCACGCTTGGTGCTCGCAGCCCGCGATGAAGAAGCATTATTTGACGTTCTGGATGAATGCACCGATTGTGGTACGGATGCCATAGCGGTCACCACCGACGTCACCCACAGCGACCAGGTTCAAGCGCTGGCCGCCCAGGc  <  2:5654878/150‑1 (MQ=255)
                                                                                                                                                     cTGCAAAGGCGCACGCTTGGTGCTCGCAGCCCGCGATGAAGAAGCATTATTTGACGTTCTGGATGAATGCACCGATTGTGGTACGGATGCCATAGCGGTCACCACCGACGTCACCCACAGCGACCAGGTTCAAGCGCTGGCCGCCCAGGc  >  1:4008892/1‑150 (MQ=255)
                                                                                                                                                     cTGCAAAGGCGCACGCTTGGTGCTCGCAGCCCGCGATGAAGAAGCATTATTTGACGTTCTGGATGAATGCACCGATTGTGGTACGGATGCCATAGCGGTCACCACCGACGTCACCCACAGCGACCAGGTTCAAGCGCTGGCCGCCCAGGc  <  2:1261506/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
ATCATCGCACCCTTCGAACACCGTAGCAGCGGCACAGAGCAAGAAGCTCCCCGTAGGGGATATCCCGTCAGGGGGCCGCTGCGTGGCAAAGTTGTGGTAATCACAGGGGCATCCAGCGGCATAGGCCGCGCCGCCGCGCACGCCTTTGCCTGCAAAGGCGCACGCTTGGTGCTCGCAGCCCGCGATGAAGAAGCATTATTTGACGTTCTGGATGAATGCACCGATTGTGGTACGGATGCCATAGCGGTCACCACCGACGTCACCCACAGCGACCAGGTTCAAGCGCTGGCCGCCCAGGC  >  NC_002947/4485410‑4485708

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