Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 733,648 A→T 10.9% intergenic (+19/‑28) PP_t09 → / → PP_t10 tRNA‑Lys/tRNA‑Pro

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947733,6480AT10.9% 210.4 / 16.4 74intergenic (+19/‑28)PP_t09/PP_t10tRNA‑Lys/tRNA‑Pro
Reads supporting (aligned to +/- strand):  ref base A (33/32);  new base T (4/4);  total (38/36)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.19e-01

TTTGGGAGCAGGATGTCAGGAGTTCGAATCCCCTCACCCCGACCATTTTCCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGGCTTTTAACCAATTGGTCGTAGGTTCGAATCCTACACGACCCACCATGTAGAAAAGGGCATCTCAAAGAGATGCCCTTTTTTTTGAAACACTCGGGGTGTAGCGCAGTCCGGTAGCGCGCCTGCTTTGGGAGCAGGATGTCAGGAGTTCGAATCCCCTCACCCCGACCATTTTTGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGGCTTTTAACC  >  NC_002947/733503‑733797
                                                                                                                                                 |                                                                                                                                                     
tttGGGAGCAGGATGTCAGGAGTTCGAATCCCCTCACCCCGACCATTTTCCGGGTCGTTAGTTAAGTTGGTAGAGCAGTGGCTTTTTAACCAATTGGTGTTAGGTTCGAATCTGACACGACCCACCAAGTAGAAAAGGGCTTCTCacggg                                                                                                                                                   >  2:656142/1‑146 (MQ=14)
tttGGGAGCAGGATGTCAGGAGTTCGAATCCCCTCACCCCGACCATTTTCCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGGCTTTTAACCAATTGGTCGTAGGTTCGAATCCTACACGACCCACCATGTAGAAAAGGGCATCTCAAaga                                                                                                                                                   <  2:696347/150‑1 (MQ=255)
  tGGGAGCAGGATGTCAGGAGTTCGAATCCCCTCACCCCGACCATTTTCCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGGCTTTTAACCAATTGGTCGTAGGTTCGAATCCTACACGACCCACCATGTAGAAAAGGGCATCTCAAAgaga                                                                                                                                                 <  1:1114646/150‑1 (MQ=255)
     gAGCAGGATGTCAGGAGTTCGAATCCCCTCACCCCGACCATTTTCCGTGTCGTTAGCTCAGTTGTGAGAGCAGTTGGCTGTTAACCAATTGGTCGTAGGTTCGAATCCTACACGACCTACCATGTAGAAAAGGGCATCTCAAGGAGGTGc                                                                                                                                              >  1:649428/1‑150 (MQ=17)
          ggATGTCAGGAGTTCGAATCCCCTCACCCCGACCATTTTCCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGGCTTTTAACCAATTGGTCGTAGGTTCGAATCCTACACGACCCACCATGTAGAAAAGGGCATCTCAAAGAGATGCCCttt                                                                                                                                         >  2:1432289/1‑150 (MQ=255)
          ggATGTCAGGAGTTCGAATCCCCTCACCCCGACCATTTTCCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGGCTTTTAACCAATTGGTCGTAGGTTCGAATCCTACACGACCCACCATGTAGAAAAGGGCATCTCAAAGAGATGCCCttt                                                                                                                                         >  2:1226568/1‑150 (MQ=255)
               tCAGGAGTTCGAATCCCCTCACCCCGACCATTTTCCGTGTCGTTCGTTCAGTTGGTAGAGCAGTTGGCTTTTAACCAATTGGTCGTAGGTTCGAATCCTACACGACCCACCATGTAGAAAAGGGCATCTCAAAGAGATGCCCttttgttt                                                                                                                                    <  1:1666220/150‑1 (MQ=255)
                      ttCGAATCCCCTCACCCCGACCATTTTCCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGGCTTTTAACCAATTGGTCGTAGGTTCGAATCCTACACGACCCACCATGTAGAAAAGGGCATCTCAAAGAGATGCCCTTTTTTTTGAAacac                                                                                                                             >  1:282750/1‑150 (MQ=255)
                         gAATCCCCTCACCCCGACCATTTTCCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGGCTTTTAACCAATTGGTCGTAGGTTCGAATCCTACACGACCCACCATGTAGAAAAGGGCATCTCAAAGAGATGCCCTTTTTTTTGAAACACTCg                                                                                                                          >  1:348719/1‑150 (MQ=255)
                          aaTCCCCTCACCCCGACCATTTTCCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGGCTTTTAACCAATTGGTCGTAGGTTCGAATCCTACACGACCCACCATGTAGAAAAGGGCATCTCAAaga                                                                                                                                                   <  2:555105/124‑1 (MQ=255)
                          aaTCCCCTCACCCCGACCATTTTCCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGGCTTTTAACCAATTGGTCGTAGGTTCGAATCCTACACGACCCACCATGTAGAAAAGGGCATCTCAAaga                                                                                                                                                   >  1:555105/1‑124 (MQ=255)
                          aaTCCCCTCACCCCGACCATTTTCCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGGCTTTTAACCAATTGGTCGTAGGTTCGAATCCTACACGACCCACCATGTAGAAAAGGGCATCTCAAAGAGATGCCCTTTTTTTTGAAACACTCgg                                                                                                                         <  2:1211353/150‑1 (MQ=255)
                                cTCACCCCGACCATTTTCCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGGCTTTTAACCAATTGGTCGTAGGTTCGAATCCTACACGACCCACCATGTAGAAAAGGGCATCTCAAAGAGATGCCCTTTTTTTTGAAACACTCGGGGTGTa                                                                                                                   <  2:1527444/150‑1 (MQ=255)
                                cTCACCCCGACCATTTTCCGGGTCGTTAGCTCAGTTGGCAGCATCGATGAACTTCAACCAATTGGTCGTAGGAGCGAATCATAAAAGACCCACCTTGTAGAGAAGGGCATCTCAAAGAGGTGCCCTTTTTTTTGAAACACTCGGGGTGTa                                                                                                                   >  1:1493153/1‑150 (MQ=21)
                                  cACCCCGACCATTTTCCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGGCTTTTAACCAATTGGTCGTAGGTTCGAATCCTACACGACCCACCATGTAGAAAAGGGCATCTCAAAGAGATGCCCTTTTTTTTGAAACACTCGGGGTGTAgc                                                                                                                 <  2:348719/150‑1 (MQ=255)
                                     cccGACCATTTTCCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGGCTTTTAACCAATTGGTCGTAGGTTCGAATCCTACACGACCCACCATGTAGAAAGGGGCACCTGCATGAGATGCCCTTTTTTTTTAACCTCTAAGGGTGTacccat                                                                                                              >  2:697685/1‑145 (MQ=14)
                                     cccGACCATTTTCCGGGTCGTTAGCGCAGTGGGTAGAGCAGTTGGCTTTTAACCAATTGGTCGTAGGTTCGAATCCTACACGACCCACCATGTAGAAAAGGGCATCTCAAAGAGACGCCCTTTTTTTTGAACCCCTCGGGGTGTaccccc                                                                                                              >  2:890762/1‑145 (MQ=21)
                                                   gggTCGTTAGCTCAGTTGGTAGAGCAGTTGGATTTTAACCAATTGGTCGTAGGGTCGAATCCTACACGACCCACCATGTAGAAAAGGGCATCTCAAAGAGATGCCCTTTTATTTGAAACACTTGGGGGGTAGCGCAGTCAGGTAgcgcgc                                                                                                >  1:697013/1‑150 (MQ=21)
                                                    ggTCGTTAGCTCAGTTTGTAGAGCAGTTGGGTTTTAACGTAGTGGTAGTAGGTTCGAATACTACACGACCCACCATGTAGAGAAGGGCATTTTAAAGAGATGCCCTTTTTTTGGAAACACTCGGGGTGTAGCGCAGTCCGGTAGCGCGcc                                                                                               <  2:649428/150‑1 (MQ=255)
                                                        gTTAGCTCAGTTGGTAGAGCAGTTGGCTTTTAACCAATTGGTCGTAGGTTCGAATCCTACACGACCCACCATGTAGAAAAGGGCATCTCAAAGAGATGCCCTTTTTTTTGAAACACTCGGGGTGTAGCGCAGTCCGGTAGCGCGCCTGCt                                                                                           <  1:546746/150‑1 (MQ=255)
                                                         ttAGCTCAGTTGGTAGAGCAGTTGGCTTTTAACCAATTGGTCGTAGGTTCGAATCCTACACGACCCACCATGTAGAAAAGGGCATCTCAAAGAGATGCCCTTTTTTTTGAAACACTCGGGGTGTAGCGCAGTCCGGTAGCGCGCCTGCtt                                                                                          <  1:1564380/150‑1 (MQ=255)
                                                         ttAGCTCAGTTGGTAGAGCAGTTGGCTTTTAACCAATTGGTCGTAGGTTCGAATCCTACACGACCCACCATGTAGAAAAGGGCATCTCAAAGAGATGCCCTTTTTTTTGAAACACTCGGGGTGTAGCGCAGTCCGGTAGCGCGCCTGCtt                                                                                          >  1:146233/1‑150 (MQ=255)
                                                             cTCAGTTGGTAGAGCAGTTGGCTTTTAACCAATTGGTCGTAGGTTCGAATCCTACACGACCCACCATGTAGAAAAGGGCATCTCAAAGAGATGCCCtttt                                                                                                                                        >  1:594280/1‑100 (MQ=255)
                                                             cTCAGTTGGTAGAGCAGTTGGCTTTTAACCAATTGGTCGTAGGTTCGAATCCTACACGACCCACCATGTAGAAAAGGGCATCTCAAAGAGATGCCCtttt                                                                                                                                        <  2:594280/100‑1 (MQ=255)
                                                              tccattGGTAGAGCGGTTGGCTTTTAACCATTTGGTCGTAGGTTCGAATCCTACACGACCCACCATGTAGAAAAGGGCATCTCAAAGAGATGCCCTTTTTTTTGAAACACTCGGGGTGTAGCGCAGTCCGGTAGCGCGCCTGCTTTGGGa                                                                                     <  1:1257446/146‑1 (MQ=255)
                                                              tCAGTTGGTAGAGCAGTTGGCTTTTAACCAATTGGTCGTAGGTTCGAATCCTACACGACCCACCATGTAGAAAAGGGCATCTCAAAGAGATGCCCTTTTTTTTGAAACACTCGGGGTGTAGCGCAGTc                                                                                                           <  2:144938/128‑1 (MQ=255)
                                                              tCAGTTGGTAGAGCAGTTGGCTTTTAACCAATTGGTCGTAGGTTCGAATCCTACACGACCCACCATGTAGAAAAGGGCATCTCAAAGAGATGCCCTTTTTTTTGAAACACTCGGGGTGTAGCGCAGTc                                                                                                           >  1:144938/1‑128 (MQ=255)
                                                              tCAGTTGGTAGAGCAGTTGGCTTTTAACCAATTGGTCGTAGGTTCGAATCCTACACGACCCACCATGTAGAAAAGGGCATCTCAAAGAGATGCCCTTTTTTTTGAAACACTCGGGGTGTAGCGCAGTCCGGTAGCGCGCCTGCTTTGGGa                                                                                     <  1:1432289/150‑1 (MQ=255)
                                                                 gTTGGTAGAGCAGTTGGCTTTTAACCAATTGGTCGTAGGTTCGAATCCTACACGACCCACCATGTAGAAAAGGGCATCTCAAAGAGATGCCCTTTTTTTTGAAACACTCGGGGTGTAGCGCAGTCCGGTAGCGCGCCTGCTTTGGGAGCa                                                                                  <  1:175202/150‑1 (MQ=255)
                                                                     gTAGAGCAGTTGGCTTTTAACCAATTGGTCGTAGGTTCGAATCCTACACGACCCACCATGTAGAAAAGGGCATCTCAAAGAGATGCCCTTTTTTTTGAAACACTCGGGGTGTAGCGCAGTCCGGTAGCGCGCCTGCTTTGGGAGCAGGAt                                                                              <  1:1226568/150‑1 (MQ=255)
                                                                      tAGAGCAGTTGGCTTTTAACGAATTGGTCGTAGGTTCGAATCCTACACGACCCACCATGTAGAAAAGGGCATCTCTTTGAGATGcc                                                                                                                                             >  1:365467/1‑86 (MQ=17)
                                                                      tAGAGCAGTTGGCTTTTAACCAATTGGTCGTAGGTTCGAATCCTACACGACCCACCATGTAGAAAAGGGCATCTCTTTGAGATGcc                                                                                                                                             <  2:365467/86‑1 (MQ=14)
                                                                      tAGAGCAGTTGGCTTTTAACCAATTGGTCGTAGGTTCGAATCCTACACGACCCACCATGTAGAAAAGGGCATCTCAAAGAGATGCCCTTTTTTTTTAAACACTCGGGGTGCAGCGCAGTCGGGTAGAGCGCCTGTTTTGGGAGCAGGCTg                                                                             >  2:1613600/1‑150 (MQ=25)
                                                                         agCAGTTGGCTTTTAACCAATTGGTCGTAGGTTCGAATCCTACACGACCCACCATGTAGAAAAGGGCATCTCAAAGAGATGCCCTTTTTTTTGAAACACTCGGGGTGTAGCGCAGTCCGGTAGCGCGCCTGCTTTGGGAGCAGGATGTCa                                                                          <  2:1493153/150‑1 (MQ=255)
                                                                          gCAGTTGGCTTCTAACCAATTGGTCGTAGGTTCGAATCCTACACGACCCACCATGTAGAAAAGGGCATCTCAAAGAGATGCCCTTTTTTTTGAAACACTCGGGGTGTAGCGCAGTCCGGTAGCGCGCCTGCTTTGGGAGCAGGATGTCAg                                                                         <  2:697013/150‑1 (MQ=255)
                                                                            aGTTGGCTTTTAACCAATTGGTCGTAGGTTCGAATCCTACACGACCCACCATGTAGAAAAGGGCATCTCAAAGAGATGCCCTTTTTTTTGAAACACTCGGGGTGTAGCGCAGTCCGGTAGCGCGCCTGCTTTGGGAGCAGGATGTCAGGa                                                                       >  2:44457/1‑150 (MQ=255)
                                                                             gTTGGCTTTTAACCAATTGGTCGTAGGTTCGAATCCTACACGACCCACCATGTAGAAAAGGGCATCTCAAAGAGATGCCCTTTTTTTTGAAACACTCGGGGTGTAGCGCAGTCCGGTAGCGCGCCTGCTTTGGGAGCAGGATGTCAGGAg                                                                      >  2:326245/1‑150 (MQ=255)
                                                                              ttGGCTTTTAACCAATTGGTCGTAGGTTCGAATCCTACACGACCCACCATGTAGAAAAGGGCATCTCTTTGAGATGCCCtttt                                                                                                                                        <  2:667967/83‑1 (MQ=17)
                                                                              ttGGCTTTTAACCAATTGGTCGTAGGTTCGAATCCTACACGACCCACCATGTAGAAAAGGGCATCTCTTTGAGATGCCCtttt                                                                                                                                        >  1:667967/1‑83 (MQ=17)
                                                                               tGGCTTTTAACCAATTGGTCGTAGGTTCGAATCCTACACGACCCACCATGTAGAAAAGGGCATCTCAAAGAGATGCCCTTTTTTTTGAAACACTCGGGGTGTAGCGCAGTCCGGTAGCGCGCCTGCTTTGGGAGCAGGATGTCAGGAGtt                                                                    <  1:890762/150‑1 (MQ=255)
                                                                                   ttttAACCAATTGGTCGTAGGTTCGAATCCTACACGACCCACCATGTAGAAAAGGGCATCTCAAAGAGATGCCCTTTTTTTTGAAACACTCGGGGTGTAGCGCAGTCCGGTAGCGCGCCTGCTTTGGGAGCAGGATGTCAGGAGTTCGaa                                                                <  1:326245/150‑1 (MQ=255)
                                                                                   ttttAACCAATTGGTCGTAGGTTCGAATCCTACACGACCCACCATGTAGAAAAGGGCATCTCAAAGAGATGCCCTTTTTTTTGAAACACTCGGGGTGTAGCGCAGTCCGGTAGCGCGCCTGCTTTGGGAGCAGGATGTCAGGAGTTCGaa                                                                >  1:1342627/1‑150 (MQ=255)
                                                                                   ttttAACCAATTGGTCGTAGGTTCGAATCCTACACGACCCACCATGTAGAAAAGGGCATCTCAAAGAGATGCCCTTTTTTTTGAAACACTCGGGGTGTAGCGCAGTCCGGTAGCGCGCCTGCTTTGGGAGCAGGATGTCAGGAGTTCGaa                                                                >  1:900340/1‑150 (MQ=255)
                                                                                   ttttAACCAATTGGTCGTAGGTTCGAATCCTACACGACCCACCATGTAGAAAAGGGCATCTCAAAGAGATGCCCTTTTTTTTGAAACACTCGGGGTGTAGCGCAGTCCGGTAGCGCGCCTGCTTTGGGAGCAGCATGTCAGGAGTTCGaa                                                                >  2:1651488/1‑150 (MQ=255)
                                                                                         ccAATTGGTCGTAGGTTCGAATCCTACACGACCCACCATGTAGAAAAGGGCATCTCAAAGAGATGCCCTTTTTTTTGAAACACTCGGGGTGTAGCGCAGTCCGGTAGCGCGCCTGCTTTGGGAg                                                                                    <  2:982165/124‑1 (MQ=255)
                                                                                         ccAATTGGTCGTAGGTTCGAATCCTACACGACCCACCATGTAGAAAAGGGCATCTCAAAGAGATGCCCTTTTTTTTGAAACACTCGGGGTGTAGCGCAGTCCGGTAGCGCGCCTGCTTTGGGAg                                                                                    >  1:982165/1‑124 (MQ=255)
                                                                                                gTCGTAGGTTCGAATCCTACACGACCCACCATGTAGAAAAGGGCATCTCAAAGAGATGCCCTTTTTTTTGAAACACTCGGGGTGTAGCGCAGTCCGGTAGCGCGCCTGCTTTGGGAGCAGGATGTCAGGAGTTCGAATCCCCTCACCCCg                                                   >  1:1483717/1‑150 (MQ=255)
                                                                                                   gTAGGTTCGAATCCTACACGACCCACCATGTAGAAAAGGGCATCTCTTTGAGATGCCCtttt                                                                                                                                        >  1:1201450/1‑62 (MQ=18)
                                                                                                   gTAGGTTCGAATCCTACACGACCCACCATGTAGAAAAGGGCATCTCTTTGAGATGCCCtttt                                                                                                                                        <  2:1201450/62‑1 (MQ=18)
                                                                                                        ttCGAATCCTACACGACCCACCATGTAGAAAAGGGCATCTCAAAGAGATGCCCTTTTTTTTGAAACACTCGGGGTGTAGCGCAGTCCGGTAGCGCGCCTGCTTTGGGAGCAGGATGTCAg                                                                         <  2:730816/120‑1 (MQ=255)
                                                                                                        ttCGAATCCTACACGACCCACCATGTAGAAAAGGGCATCTCAAAGAGATGCCCTTTTTTTTGAAACACTCGGGGTGTAGCGCAGTCCGGTAGCGCGCCTGCTTTGGGAGCAGGATGTCAg                                                                         >  1:730816/1‑120 (MQ=255)
                                                                                                        ttCGAATCCTACACGACCCACCATGTAGAAAAGGGCATCTCAAAGAGATGCCCTTTTTTTTGAAACACTCGGGGTGTAGCGCAGTCCGGTAGCGCGCCTGCTTTGGGAGCAGGATGTCAGGAGTTCGAATCCCCTCACCCCGACCAtttt                                           >  2:371731/1‑150 (MQ=255)
                                                                                                          cGAATCCTACACGACCCACCATGTAGAAAAGGGCATCTCAATGAGTTGTTCTTTTTTTTGCAACCCCCGGGATGTACCCCAGTACGGTCGCGCTCCCCCTTCAGGAGCATGTTGTGATGAGTTCGAATCCCCTCACCCCGACCATTTTTg                                         <  2:282750/150‑1 (MQ=19)
                                                                                                             aTCCTACACGACCCACCATGTAGAAAAGGGCATCTCAAAGAGATGCCCTTTTTTTTGAAACACTCGGGGTGTAGCGCAGTCCGGTAGCGCGCCTGCTTTGGGAGAAGGGTGTCAGGAGTGCGAATTCCCGCACCCCGACCAATTTTGGGt                                      >  2:452794/1‑150 (MQ=25)
                                                                                                              tCCTACACGACCCACCATGTAGAAAAGGGCATCTCAAAGAGATGCCCTTTTTTTTGAAac                                                                                                                               >  1:834613/1‑60 (MQ=255)
                                                                                                              tCCTACACGACCCACCATGTAGAAAAGGGCATCTCAAAGAGATGCCCTTTTTTTTGAAac                                                                                                                               <  2:834613/60‑1 (MQ=255)
                                                                                                              tCCTACACGACCCACCATGTAGAAAAGGGCATCTCAAAGAGATGCCCTTTTTTTTGAAACACTCGGGGTGTAGCGCAGTCCGGTAGCGCGCCTGCTTTGGGAGCAGGATGTCAGGAGTTCGAATCCCCTCACCCCGACCATTTTTGGGTc                                     >  2:826075/1‑150 (MQ=255)
                                                                                                                cTACACGACCCACCATGTAGAAAAGGGCATCTCAAAGAGATGCCCTTTTTTTTGAAACACTCGGGGTGTAGCGCAGTCCGGTAGCGCGCCTGCTTTGGGAGCAGGATGTCAGGAGTTCGAATCCCCTCACCCCGAcc                                                >  2:678836/1‑137 (MQ=255)
                                                                                                                cTACACGACCCACCATGTAGAAAAGGGCATCTCAAAGAGATGCCCTTTTTTTTGAAACACTCGGGGTGTAGCGCAGTCCGGTAGCGCGCCTGCTTTGGGAGCAGGATGTCAGGAGTTCGAATCCCCTCACCCCGAcc                                                <  1:678836/137‑1 (MQ=255)
                                                                                                                cTACACGACCCACCATGTAGAAAAGGGCATCTCAAAGAGATGCCCTTTTTTTTGAAACACTCGGGGTGTAGCGCAGTCCGGTAGCGCGCCTGCTTTGGGAGCAGGATGTCAGGAGTTCGAATCCCCTCACCCCGACCATTTTTGGGTCGt                                   <  1:44457/150‑1 (MQ=255)
                                                                                                                 tACACGACCCACCATGTAGAAAAGGGCATCTCAAAGAGATGCCCTTTTTTTTGAAACACTCGGGGTGTAACGCAGTCCGGTAGCGCGCCTGCTTTGGGAGCAGGATGTCAGGCATAAGAATCCCCTTACGCCCACCAAGTTTgtggtggt                                  >  2:1384383/1‑143 (MQ=25)
                                                                                                                 tACACGACCCACCATGTAGAAAAGGGCATCTCAAAGAGATGCCCTTTTTTTTGAAACACTCGGGGTGCAGCGCAGTCCGGTAGCGCGCCTGCTTTGGGAGCAGGATGTCAGGAGTTCGAATCCCCTCACCCCGACCATTTTTGGGTCGtt                                  >  2:191182/1‑150 (MQ=255)
                                                                                                                         ccaccaTGTAGAAAAGGGCATCTCAAAGAGATGCCCTTTTTTTTGAAACACTCGGGGTGTAGCGCAGTCCGGTAGCGCGCCTGCTTTGGGAGCAGGATGTCAGGAGTTCGAATCCCCTCACCCCGACCATTTTTGGGTc                                     <  2:390523/139‑1 (MQ=255)
                                                                                                                         ccaccaTGTAGAAAAGGGCATCTCAAAGAGATGCCCTTTTTTTTGAAACACTCGGGGTGTAGCGCAGTCCGGTAGCGCGCCTGCTTTGGGAGCAGGATGTCAGGAGTTCGAATCCCCTCACCCCGACCATTTTTGGGTc                                     >  1:390523/1‑139 (MQ=255)
                                                                                                                          caccaTGTAGAAAAGGGCATCTCAAAGAGATGCCCTTTTTTTTGAAACACTCGGGGTGTAGCGCAGTCCGGTAGCGCGCCTGCTTTGGGAGCAGGATGTCAGGAGTTCGAATCCCCTCACCCCGACCATTTTTg                                         <  1:1562111/134‑1 (MQ=255)
                                                                                                                          caccaTGTAGAAAAGGGCATCTCAAAGAGATGCCCTTTTTTTTGAAACACTCGGGGTGTAGCGCAGTCCGGTAGCGCGCCTGCTTTGGGAGCAGGATGTCAGGAGTTCGAATCCCCTCACCCCGACCATTTTTg                                         >  2:1562111/1‑134 (MQ=255)
                                                                                                                           accaTGTAGAAAAGGGCATCTCAAAGAGATGCCCTTTTTTTTGAAACACTCGGGGTGTAGCGCAGTCCGGTAGCGCGCCTGCTTTGGGAGCAGGATGTCAGGAGTTCGAATCCCCTCACCCCGACCATTTTTGGGTCGTTAGCTCAGTTg                        >  1:958476/1‑150 (MQ=255)
                                                                                                                           accaTGTAGAAAAGGGCATCTCAAAGAGATGCCCTTTTTTTTGAAACACTCGGGGTGTAGCGCAGTCCGGTAGCGCGCCTGCTTTGGGAGCAGGATGTCAGGAGTTCGAATCCCCTCACCCCGACCATTTTTGGGTCGTTAGCTCAGTTg                        <  1:826075/150‑1 (MQ=255)
                                                                                                                             caTGTAGAAAACGGCATCTCAAAGAAGTGCCCTTTTTTGTAAAACACTCGTGTTGTAGCGCAGTCCGGTAGCGCGCCTGCTTTGGGAGCAGGATGTCAGGAGTTCGAATCCCCTCACCCCGACCATTTTTGGGTCGTTAGCTCAGTTGGt                      <  1:240733/150‑1 (MQ=18)
                                                                                                                                gTAGAAAAGGGCATCTCAAAGAGATGCCCTTTTTTTTGAAACACTCGGGGTGTAGCGCAGTCCGGTAGCGCGCCTGCTTTGGGAGCAGGATGTCAGGAGTTCGAATCCCCTCACCCCGACCATTTTTGGGTCGTTAGCTCAGTTGGTaga                   <  2:1483717/150‑1 (MQ=255)
                                                                                                                                      aaGGGCATCTCTTTGAGATGCCCTTTTTTTTGAAACACTCGGGGTGTAGCGCAGTCCGGTAgcgcgc                                                                                                <  2:1326733/67‑1 (MQ=25)
                                                                                                                                      aaGGGCATCTCTTTGAGATGCCCTTTTTTTTGAAACACTCGGGGTGTAGCGCAGTCCGGTAgcgcgc                                                                                                >  1:1326733/1‑67 (MQ=25)
                                                                                                                                      aaGGGCATCTCAAAGAGATGCCCTTTTTTTTGAAACACTCGGGGTGTAACGCAGTCCGGTAGCGCGCCTGCTTTGGGAGCAGGATGTCAGGAGTTCGAATCCCCTCACCCCGACCATTTTTGGGTCGTTAGCTCAGTTGGTAGAGCAGtt             <  1:1384383/150‑1 (MQ=255)
                                                                                                                                        gggCATCTCAAAGAGATGCCCTTTTTTTTGAAACACTCGGGGTGTAGCGCAGTCCGGTAGCGCGCCTGCCTTGGGAGCAGGATGTCAGGAGTTCGAATCCCCTCACCCCGACCATTTTTGGGTCGTTAGCTCAGTTGGTAGAGCAGTTgg           >  1:378444/1‑150 (MQ=255)
                                                                                                                                           cATCTCAAAGAGATGCCCTTTTTTTTGAAACACTCGGGGTGTAGCGCAGTCCGGTAGCGCGCCTGCTTTGGGAGCAGGATGTCAGGAGTTCGAATCCCCTCACCCCGACCATTTTTGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGGCtt        >  1:343515/1‑150 (MQ=255)
                                                                                                                                                 aaaGAGATGCCCTTTTTTTTGAAACACTCGGGGTGTAGCGCAGTCCGGTAGCGCGCCTGCTTTGGGAGCAGGATGTCAGGAGTTCGAATCCCCTCACCCCGACCATTTTTGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGGCTTTTAAcc  >  2:1784512/1‑150 (MQ=255)
                                                                                                                                                 |                                                                                                                                                     
TTTGGGAGCAGGATGTCAGGAGTTCGAATCCCCTCACCCCGACCATTTTCCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGGCTTTTAACCAATTGGTCGTAGGTTCGAATCCTACACGACCCACCATGTAGAAAAGGGCATCTCAAAGAGATGCCCTTTTTTTTGAAACACTCGGGGTGTAGCGCAGTCCGGTAGCGCGCCTGCTTTGGGAGCAGGATGTCAGGAGTTCGAATCCCCTCACCCCGACCATTTTTGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGGCTTTTAACC  >  NC_002947/733503‑733797

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