Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 733,650 A→T 5.8% intergenic (+21/‑26) 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,6500AT 5.8% 273.7 / 6.5 86intergenic (+21/‑26)PP_t09/PP_t10tRNA‑Lys/tRNA‑Pro
Reads supporting (aligned to +/- strand):  ref base A (44/37);  new base T (2/3);  total (46/40)
Fisher's exact test for biased strand distribution p-value = 6.60e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.50e-01

TTTGGGAGCAGGATGTCAGGAGTTCGAATCCCCTCACCCCGACCATTTTCCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGGCTTTTAACCAATTGGTCGTAGGTTCGAATCCTACACGACCCACCATGTAGAAAAGGGCATCTCAAAGAGATGCCCTTTTTTTTGAAACACTCGGGGTGTAGCGCAGTCCGGTAGCGCGCCTGCTTTGGGAGCAGGATGTCAGGAGTTCGAATCCCCTCACCCCGACCATTTTTGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGGCTTTTAAC  >  NC_002947/733503‑733796
                                                                                                                                                   |                                                                                                                                                  
tttGGGAGCAGGATGTCAGGAGTTCGAATCCCCTCACCCCGACCATTTTCCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGGCTTTTAACCAATTGGTCGTAGGTTCGAATCCTACACGACCCACCATGTAGAAAAGGGCATCTCAAaga                                                                                                                                                  >  1:426637/1‑150 (MQ=255)
     gAGCAGGATGTCAGGAGTTCGAATCCCCTCACCCCGAGCATTTTCCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGGCTTTTAACCAATTGGTCGTAGGTTCGAATCCTACACGACCCACCATGTAGAAAAGGGCATCTCaaa                                                                                                                                                    >  1:1329613/1‑143 (MQ=255)
     gAGCAGGATGTCAGGAGTTCGAATCCCCTCACCCCGACCATTTTCCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGGCTTTTAACCAATTGGTCGTAGGTTCGAATCCTACACGACCCACCATGTAGAAAAGGGCATCTCaaa                                                                                                                                                    <  2:1329613/143‑1 (MQ=255)
      aGCAGGATGTCAGGAGTTCGAATCCACTCGCCCCGACCATTTTCCGGGTCGTTAGCTCCGTAGCTAGTCCGGCTTTCTTTTAACCAATTGGTCGTAGGTTCGAATCCTACACGACCCACCATGTAGAAAAGGGCATCTCTTTGAGATGcc                                                                                                                                            <  1:283383/150‑1 (MQ=17)
          ggATGTCAGGAGTTCGAATCCCCTCACCCCGACCATTTTCCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGGCTTTTAACCAATTGGTCGTAGGTTCGAATCCTACACGACCCACCATGTAGAAAAGGGCATCTCAAAGAGATg                                                                                                                                              >  1:308518/1‑144 (MQ=255)
          ggATGTCAGGAGTTCGAATCCCCTCACCCCGACCATTTTCCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGGCTTTTAACCAATTGGTCGTAGGTTCGAATCCTACACGACCCACCATGTAGAAAAGGGCATCTCAAAGAGATg                                                                                                                                              <  2:308518/144‑1 (MQ=255)
                cAGGAGTTCGAATCCCCTCACCCCGACCATTTTCCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGGCTTTTAACCAATTGGTCGTAGGTTCGAATCCTACACGACCCACCATGTAGAAAAGGGCATCTCAAAGAGATGCCCTTTTTTTTg                                                                                                                                  >  1:1446451/1‑150 (MQ=255)
                   gacttCGAATCCGCTCACCCCGACCATTTTCCGGCTCGTTAGCTCAGTTTGTAGAGCAGTTGGCTTTAAACCAATTGGTCGTAGGTTTGAAGCCTACGCGACCCACCATGTAGAAACGCGCATGTCAAAGAGATGCCCTTTTTTTTGaaa                                                                                                                               <  2:238188/147‑1 (MQ=18)
                           aTCCCCTCACCCCGACCATTTTCCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGGCTTTTAACCAATTGGTCGTAGGTTCGAATCCTACACGACCCACCATGTAGAAAAGGGCATCTCAAAGAGATGCCCTTTTTTTTGAAACACTCggg                                                                                                                       <  2:426637/150‑1 (MQ=255)
                             ccccTCACCCCGACCATTTTCCGGGTCGTTAGCTCAGTTGGTAGAGCCGTTGGCTTATAACCAATTGGTAGTCGGTTCGACTCCTACACGACCCACCATGTAGAAAAGGGCATCTCAAAGCGATGCCCTTTTTTTTGAAAAACTGggggc                                                                                                                     >  1:1485051/1‑149 (MQ=18)
                              cccTCACCCCGACCATTTTCCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGGCTTTTAACCAATTGGTCGTAGGTTCGAATCCTACACGACCCACCATGTAGAAAAGGGCATCTCTTTGAGATGccc                                                                                                                                           >  1:1094350/1‑127 (MQ=255)
                              cccTCACCCCGACCATTTTCCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGGCTTTTAACCAATTGGTCGTAGGTTCGAATCCTACACGACCCACCATGTAGAAAAGGGCATCTCTTTGAGATGccc                                                                                                                                           <  2:1094350/127‑1 (MQ=255)
                                cTCACCCCGACCATTTTCCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGGCTTTTAACCAATTGGTCGTAGGTTCGAATCCTACACGACCCACCATGTAGAAAAGGGCATCTCAAAGAGATGCCCTTTTTTTTGAAACACTCggg                                                                                                                       <  2:1576185/145‑1 (MQ=255)
                                cTCACCCCGACCATTTTCCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGGCTTTTAACCAATTGGTCGTAGGTTCGAATCCTACACGACCCACCATGTAGAAAAGGGCATCTCAAAGAGATGCCCTTTTTTTTGAAACACTCggg                                                                                                                       >  1:1576185/1‑145 (MQ=255)
                                       cGACCATTTTCCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGGCTTTTAACCAATTGGTCGTAGGTTCGAATCCTACACGACCCACCATGTAGAAAAGGGCATCTCAAAGAGATGCCCTTTTTTTTGAAACACTCGGGGTGTAGCGCAGt                                                                                                           >  1:48125/1‑150 (MQ=255)
                                         aCCATTTTCCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGGCTTTTAACCAATTGGTCGTAGGTTCGAATCCTACACGACCCACCATGTAGAAAAGGGCATCTCAAAGAGATGCCCTTTTTTTTGAAACACTCGGGGTGTAGCGCAGTcc                                                                                                         >  1:133564/1‑150 (MQ=255)
                                             ttttCCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGGCTTTTAACCAATTGGTCGTAGGTTCGAATCCTACACGACCCACCATGAAGAAAAGGGCATCTCAAAGAGATGCCCTTTTTTTTGAAACACTCGGGGTGTa                                                                                                                  >  1:1432214/1‑137 (MQ=255)
                                             ttttCCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGGCTTTTAACCAATTGGTCGTAGGTTCGAATCCTACACGACCCACCATGAAGAAAAGGGCATCTCAAAGAGATGCCCTTTTTTTTGAAACACTCGGGGTGTa                                                                                                                  <  2:1432214/137‑1 (MQ=255)
                                              tttCCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGGCTTTTAACCAATTGGTCGTAGGTTCGAATCCTACACGACCCACCATGTAGAAAAGGGCATCTCAAAGAGATGCCCTTTTTTTTGAAACACTCGGGGTGTAGCGCAGTCCGGTAg                                                                                                    >  1:565126/1‑150 (MQ=255)
                                                          tttttCAGTTGTTAGAGCGGTTGGCTCTTGATCAATTCGTGGTTTGTTCGTATTTCACACGACCCAGCTTGTACAAAAGGGCATCTGTAAGAGATGCCCTTTTCGTGACACTCTCGGGGTGTAGCGCAGTCCGGTAGCGCGCCTGCTTTg                                                                                       <  2:57039/146‑1 (MQ=16)
                                                           aGCTCAGTTGGTAGAGCATTTGGCTTTTAACCAATTGGTCGTAGGTTCGAATCCTACACGACCCACCATGTAGAAAAGGGCATCTCAAAGAGATGCCCTTTTTTTTGAAACACTCGGGGTGTAGCGCAGTCCGGTAGCGCGCCTGCTTTg                                                                                       <  2:1485051/150‑1 (MQ=255)
                                                             cTCAGTTGGTAGAGCAGTTGGCTTTTAACCAATTGGTCGTAGGTTCGAATCCTACAAGACCCACCACGTAGAAAAGGGCATCTCAAAGAGATGCCCTTTTTTTTTCAACACTCGGGGTATCGCGCAATCCGGTAGCGCGCCTGCTTTggg                                                                                     >  2:855670/1‑150 (MQ=21)
                                                              tCAGTTGGTAGAGCAGTTGGCTTTTAACCAATTGGTCGTAGGTTCGAATCCTACACTACCCACCATGTAGAAAAGGGCATCTCAAAGAGATGCCCTTTTTTTTGAAACACTCgggg                                                                                                                      >  2:127836/1‑116 (MQ=255)
                                                              tCAGTTGGTAGAGCAGTTGGCTTTTAACCAATTGGTCGTAGGTTCGAATCCTACACTACCCACCATGTAGAAAAGGGCATCTCAAAGAGATGCCCTTTTTTTTGAAACACTCgggg                                                                                                                      <  1:127836/116‑1 (MQ=255)
                                                              tCAGTTGGTAGAGCAGTTGGCTTTTAACCAATTGGTCGTAGGTTCGAATCCTACACGACCCACCATGTAGAAAAGGGCATCTCAAAGAGATGCCCtttttttt                                                                                                                                   >  1:530439/1‑103 (MQ=255)
                                                              tCAGTTGGTAGAGCAGTTGGCTTTTAACCAATTGGTCGTAGGTTCGAATCCTACACGACCCACCATGTAGAAAAGGGCATCTCAAAGAGATGCCCtttttttt                                                                                                                                   <  2:530439/103‑1 (MQ=255)
                                                                 gTTGGTAGAGCAGTTGGCTTTTAACCAATTGGTCGTAGGTTCGAATCCTACACGACCCACCATGTAGAAAAGGGCATCTCAAAGAGATGCCCTTTTTTTTGAAACACTCGGGGTGTAGCGCAGTCCGGTAGCGCGCCTGCTTTGGGAGCa                                                                                 <  2:1446451/150‑1 (MQ=255)
                                                                     gTAGAGCAGTTGGCTTTTAACCAATTGGTCGTAGGTTCGAATCCTACACGACCCACCATGTAGAAAAGGGCATCTCAAAGAGATGCCCTTTTTTTTGAAACACTCGGGGTGTAGCGCAGTCCGGTAGCGCGCCTGCTTTGGGAg                                                                                   <  1:1288376/144‑1 (MQ=255)
                                                                     gTAGAGCAGTTGGCTTTTAACCAATTGGTCGTAGGTTCGAATCCTACACGACCCACCATGTAGAAAAGGGCATCTCAAAGAGATGCCCTTTTTTTTGAAACACTCGGGGTGTAGCGCAGTCCGGTAGCGCGCCTGCTTTGGGAg                                                                                   >  2:1288376/1‑144 (MQ=255)
                                                                     gTAGAGCAGTTGGCTTTTAACCAATTGGTCGTAGGTTCGAATCCTACACGACCCACCATGTAGAAAAGGGCATCTCAAAGAGATGCCCTTTTTTTTGAAACACTCGGGGTGTAGCGCAGTCCGGTAGCGCGCCTGCTTTGGGAGCAGGAt                                                                             >  2:455299/1‑150 (MQ=255)
                                                                     gTAGAGCAGTTGGCTTTTAACCAATTGGTCGTAGGTTCGAATCCTACACGACCCACCATGTAAAAAAGGGCATCTCTTTGAGATGCCCttttc                                                                                                                                      >  1:661662/1‑92 (MQ=14)
                                                                     gTAGAGCAGTTGGCTTTTAACCAATTGGTCGTAGGTTCGAATCCTACACGACCCACCATGTAAAAAAGGGCATCTCTTTGAGATGCCCttttc                                                                                                                                      <  2:661662/93‑2 (MQ=14)
                                                                        gagCAGTTGGCTTTTAACCAATTGGTCGTAGGTTCGAATCCTACACGACCCACCATGTAGAAAAGGGCATCTCAAAGAGATGCCCTTTTTTTTGAAACACTCGGGGTGTAGCGCAGTCCGGTAGCGCGCCTGCTTTGGGAGCAGGATGTc                                                                          >  2:70545/1‑150 (MQ=255)
                                                                        gagCAGTTGGCTTTTAACCAATTGGTCGTAGGTTCGAATCCTACACGACCCACCATGTAGAAAAGGGCATCTCAAAGAGATGCCCTTTTTTTTGAAACACTCGGGGTGTAGCGCAGTCCGGTAGCGCGCCTGCTTTGGGAGCAGGATGTc                                                                          <  1:1240926/150‑1 (MQ=255)
                                                                              ttGGCTTTTAACCAATTGGTCGTAGGTTCGAATCCTACACGACCCACCATGTAGAAAAGGGCATCTCAAAGAGATGCCCTTTTTTTTGAAACACTCGGGGTGTAGCGCAGTCCGGTAGCGCGCCTGCTTTgg                                                                                      >  1:1167864/1‑132 (MQ=255)
                                                                              ttGGCTTTTAACCAATTGGTCGTAGGTTCGAATCCTACACGACCCACCATGTAGAAAAGGGCATCTCAAAGAGATGCCCTTTTTTTTGAAACACTCGGGGTGTAGCGCAGTCCGGTAGCGCGCCTGCTTTgg                                                                                      <  2:1167864/132‑1 (MQ=255)
                                                                              ttGGCTTTTAACCAATTGGTCGTAGGTTCGAATCCTACACGACCCACCATGTAGAAAAGGGCATCTCAAAGAGATGCCCTTTTTTTTGAAACACTCGGGGTGTAGCGCAGTCCGGTAGCGCGCCTGCTTTg                                                                                       <  2:574292/131‑1 (MQ=255)
                                                                              ttGGCTTTTAACCAATTGGTCGTAGGTTCGAATCCTACACGACCCACCATGTAGAAAAGGGCATCTCAAAGAGATGCCCTTTTTTTTGAAACACTCGGGGTGTAGCGCAGTCCGGTAGCGCGCCTGCTTTg                                                                                       >  1:574292/1‑131 (MQ=255)
                                                                                   ttttAACCAATTGGTCGTAGGTTCGAATCCTACACGACCCACCATGTAGAAAAGGGCATCTCAAAGAGATGCCCTTTTTTTTGAAACACTCGGGGTGTAg                                                                                                                 <  1:1203785/100‑1 (MQ=255)
                                                                                   ttttAACCAATTGGTCGTAGGTTCGAATCCTACACGACCCACCATGTAGAAAAGGGCATCTCAAAGAGATGCCCTTTTTTTTGAAACACTCGGGGTGTAg                                                                                                                 >  2:1203785/1‑100 (MQ=255)
                                                                                       aaCCAATTGGTCGTAGGTTCGAATCCTACACGACCCACCATGTAGAAAAGGGCATCTCAAAGAGATGCCCtttttt                                                                                                                                     <  1:27354/76‑1 (MQ=255)
                                                                                       aaCCAATTGGTCGTAGGTTCGAATCCTACACGACCCACCATGTAGAAAAGGGCATCTCAAAGAGATGCCCtttttt                                                                                                                                     >  2:27354/1‑76 (MQ=255)
                                                                                       aaCCAATTGGTCGTAGGTTCGAATCCTACACGACCCACCATGTAGAAAAGGGCATCTCAAAGAGATGCCCTTTTTTTTGAAACACTCGGGGTGTAGCGCAGTCCGGTAGCGCGCCTGCTTTGTGAGCa                                                                                 <  1:923138/128‑1 (MQ=255)
                                                                                       aaCCAATTGGTCGTAGGTTCGAATCCTACACGACCCACCATGTAGAAAAGGGCATCTCAAAGAGATGCCCTTTTTTTTGAAACACTCGGGGTGTAGCGCAGTCCGGTAGCGCGCCTGCTTTGTGAGCa                                                                                 >  2:923138/1‑128 (MQ=255)
                                                                                       aaCCAATTGGTCGTAGGTTCGAATCCTACACGACCCACCATGTAGAAAAGGGCATCTCAAAGAGATGCCCTTTTTTTTGAAACACTCGGGGTGTAGCGCAGTCCGGTAGCGCGCCTGCTTTGGGAGCa                                                                                 >  2:1053455/1‑128 (MQ=255)
                                                                                       aaCCAATTGGTCGTAGGTTCGAATCCTACACGACCCACCATGTAGAAAAGGGCATCTCAAAGAGATGCCCTTTTTTTTGAAACACTCGGGGTGTAGCGCAGTCCGGTAGCGCGCCTGCTTTGGGAGCa                                                                                 <  1:1053455/128‑1 (MQ=255)
                                                                                       aaCCAATTGGTCGTAGGTTCGAATCCTACACGACCCACCATGTAGAAAAGGGCATCTCAAAGAGATGCCCTTTTTTTTGAAACACTCGGGGTGTAGCGCAGTCCGGTAGCGCGCCTGCTTTGGGAGCAGGATGTCAGGAGTTCGAATccc                                                           <  2:48125/150‑1 (MQ=255)
                                                                                          cAATTGGTCGTAGGTTCGAATCCTACACGACCCACCATGTAGAAAAGGGCATCTCAAAGAGATGCCCTTTTTTTTGAAACACTCGGGGTGTAGCGCAGTCCGGTAGCGCGCCTGCTTTGGGAGCAGGATGTCAg                                                                        >  2:1402054/1‑134 (MQ=255)
                                                                                          cAATTGGTCGTAGGTTCGAATCCTACACGACCCACCATGTAGAAAAGGGCATCTCAAAGAGATGCCCTTTTTTTTGAAACACTCGGGGTGTAGCGCAGTCCGGTAGCGCGCCTGCTTTGGGAGCAGGATGTCAg                                                                        <  1:1402054/134‑1 (MQ=255)
                                                                                                gTCGTAGGTTCGAATCCTACACGACCCACGATGTAGAAAAGGGCATCTCAAAGAGATGCCCTTTTTTTTGAAACACTCGGGGTGTAGCGCAGTCCGGTAGCGCGCCTGCTTTGGGAGCAGGATGAAAGGAGTTCGAATCCCCTCACCACg                                                  >  2:1235177/1‑150 (MQ=255)
                                                                                                    tAGGTTCGAATCCTACACGACCCACCATGTAGAAAAGGGCATCTCAAAGAGATGCCCt                                                                                                                                          <  2:1279931/58‑1 (MQ=255)
                                                                                                    tAGGTTCGAATCCTACACGACCCACCATGTAGAAAAGGGCATCTCAAAGAGATGCCCt                                                                                                                                          >  1:1279931/1‑58 (MQ=255)
                                                                                                       gTTCGAATCCTACACGACCCACCATGTAGAAAAGGGCATATCAAAGAGATGCCCTTTTTTTTGAAACACTCGGGGTGTAGCGCAGTCCGGTAGCGCGCCTGCTTTGGGAGCAGGATGTCAGGAGTTCGAATCCCCTCACCCCGACCAttt                                           >  2:188477/1‑150 (MQ=255)
                                                                                                        ttCGAATCCTACACGACCCACCATGTAGAAAAGGGCATCTCAAAGAGATGCCCTTTTTTTTGAAACACTCGGGGTGTAGCGCAGTCCGGTAGCGCGCCTGCTTTg                                                                                       >  1:320036/1‑105 (MQ=255)
                                                                                                        ttCGAATCCTACACGACCCACCATGTAGAAAAGGGCATCTCAAAGAGATGCCCTTTTTTTTGAAACACTCGGGGTGTAGCGCAGTCCGGTAGCGCGCCTGCTTTg                                                                                       <  2:320036/105‑1 (MQ=255)
                                                                                                        ttCGAATCCTACACGACCCACCATGTAGAAAAGGGCATCTCAAAGAGATGCCCTTTTTTTTGAAACACTCGGGGTGTAGCGCAGTCCGGTAGCGCGCCTGCTTTGGGAGCAGGATGTCAGGAGTTCGAATCCCCTCACCCCGACCAtttt                                          >  1:524653/1‑150 (MQ=255)
                                                                                                           gAATCCTACACGACCCACCATGTAGAAAAGGGCATCTCAAAGAGATGCCCTTTTTTTTGAAACACTCGGGGTGTAGCGCAGTCCGGTAGCGCGCCTGCTTTGGGAGCAGGATGTCAGGAGTTCGAATCCCCTCACCCCGACCATTTTTgg                                       >  2:1439390/1‑150 (MQ=255)
                                                                                                           gAATCCTACACGACCCACCATGTAGAAAAGGGCATCTCAAAGAGATGCCCTTTTTTTTGAAACACTCGGGGTGTAGCGCAGTCCGGTAGCGCGCCTGCTTTGGGAGCAGGATGTCAGGAGTTCGAATCCCCTCACCCCGACCATTTTTgg                                       <  2:133564/150‑1 (MQ=255)
                                                                                                           gAATCCTACACGACCCACCATGTAGAAAAGGGCATCTCAAAGAGATGCCCTTTTTTTTGAAACACTCGGGGTGTAGCGCAGTCCGGTAGCGCGCCTGCTTTGGGAGCAGGATGTCAGGAGTTCGAATCCCCTCACCCCGACCATTTTTgg                                       <  1:455299/150‑1 (MQ=255)
                                                                                                                cTACACGACCCACCATGTAGAAAAGGGCATCTCAAAGAGATGCCCTTTTTTTTGAAACACTCGGGGTGTAGCGCAGTCCGGTAGCGCGCCTGCTTTGGGAGCAGGATGTCAGGAGTTCGAATCCCCTCACCCCGACCATTTTTGGGTCGt                                  >  1:1431603/1‑150 (MQ=255)
                                                                                                                cTACACGACCCACCATGTAGAAAAGGGCATCTCAAAGAGATGCCCTTTTTTTTGAAACACTCGGGGTGTAGCGCAGTCCGGTAGCGCGCCTGCTTTGGGAGCAGGATGTCAGGAGTTCGAATCCCCTCACCCCGACCATTTTTGGGTCGt                                  <  2:1431603/150‑1 (MQ=255)
                                                                                                                 tACACGACCCACCATGTAGAAAAGGGCATCTCAAAGAGATGCCCTTTTTTTTGAAACACTCgg                                                                                                                        <  2:908438/63‑1 (MQ=255)
                                                                                                                 tACACGACCCACCATGTAGAAAAGGGCATCTCAAAGAGATGCCCTTTTTTTTGAAACACTCgg                                                                                                                        >  1:908438/1‑63 (MQ=255)
                                                                                                                  acacGACCCACCATGTAGAAAAGGGCATCTCAAAGAGATGCCCTTTTTTTTGAAACACTCGGGGTGTAGCGCAGTCCGGTAGCGCGCCTGCTTTGGGAGCAGGATGTCAGGAGTTCGAATCCCCTCACCCCGACCATTTTTGGGTCGTTa                                >  1:286811/1‑150 (MQ=255)
                                                                                                                      taCCCACCATGTAGAAAAGGGCATCTCAAAGAGATGCCCTTTTTTTTGAAACACTCGGGgt                                                                                                                     <  2:652131/60‑1 (MQ=255)
                                                                                                                      gACCCACCATGTAGAAAAGGGCATCTCAAAGAGATGCCCTTTTTTTTGAAACACTCGGGgt                                                                                                                     >  1:652131/1‑61 (MQ=255)
                                                                                                                      gACCCACCATGTAGAAAAGGGCATCTCAAAGAGATGCCCTTTTTTTTGAAACACTCGGGGTGTAGCGCAGTCCGGTAGCGCGCCTGCTTTGGGAGCAGGATGTCAGGAGTTCGAATCCCCTCACCCCGACCATTTTTGGGTCGTTAGCTc                            >  2:1049931/1‑150 (MQ=255)
                                                                                                                      gACCCACCATGTAGAAAAGGGCATCTCAAAGAGATGCCCTTTTTTTTGAAACACTCGGGGTGTAGCGCAGTCCGGTAGCGCGCCTGCTTTGGGAGCAGGATGTCAGGAGTTCGAATCCCCTCACCCCGACCATTTTTGGGTCGTTAGCTc                            >  2:1485076/1‑150 (MQ=255)
                                                                                                                         ccaccaTGTAGAAAAGGGCATCTCAAAGAGATGCCCTTTTTTTTGAAACACTCGGGGTGTAGCGCAGTCCGGTAGCGCGCCTGCTTTGGGAGCAGGATGTCAGGAGTTCGAATCCCCTCACCCCGACCATTTTTGGGTCGTTAGCt                             >  1:210267/1‑146 (MQ=255)
                                                                                                                         ccaccaTGTAGAAAAGGGCATCTCAAAGAGATGCCCTTTTTTTTGAAACACTCGGGGTGTAGCGCAGTCCGGTAGCGCGCCTGCTTTGGGAGCAGGATGTCAGGAGTTCGAATCCCCTCACCCCGACCATTTTTGGGTCGTTAGCt                             <  2:210267/146‑1 (MQ=255)
                                                                                                                            ccaTGTAGAAAAGGGCATCTCAAAGAGATGCCCTTTTTTTTGAAACACTCGGGGTGTAGCGCAGTCCGGTAGCGCGCCTGCTTTGGGAGCAGGATGTCAGGAGTTCGAATCCCCTCACCCCGACCATTTTTGGGTCGTTAGCTCAGTTgg                      <  1:1049931/150‑1 (MQ=255)
                                                                                                                               tGTAGAAAAGGGCATCTCAAAGAGATGCCCTTTTTTTTGAAACACTCgg                                                                                                                        >  2:1305379/1‑49 (MQ=255)
                                                                                                                               tGTAGAAAAGGGCATCTCAAAGAGATGCCCTTTTTTTTGAAACACTCgg                                                                                                                        <  1:1305379/49‑1 (MQ=255)
                                                                                                                                   gAAAAGGGCATCTCAAAGAGATGCCCTTTTTTTTGAAACACTCGGGGTGTAGCGCAGTCCGGTAGCGCGCCTGc                                                                                           >  2:1257873/1‑74 (MQ=255)
                                                                                                                                   gAAAAGGGCATCTCAAAGAGATGCCCTTTNTTTTGAAACACTCGGGGTGTAGCGCAGTCCGGTAGCGCGCCTGc                                                                                           <  1:1257873/74‑1 (MQ=255)
                                                                                                                                    aaaaGGGCATCTCAAAGAGATGCCCTTTTTTTTGAAACACTCGGGGTGTAGCGCAGTCCGGTAGCGCGCCTGCTTTGGGAGCAGGATGTCAGGAGTTCGAATCCCCTCACCCCGACCATTTTTGGGTCGTTAGCTCAGTTGGTAGAGCAg              <  1:70545/150‑1 (MQ=255)
                                                                                                                                    aaaaGGGCATCTCAAAGAGATGCCCTTTTTTTTGAAACACTCGGGGTGTAGCGCAGTCCGGTAGCGCGCCTGCTTTGGGAGCAGGATGTCAGGAGTTCGAATCCCCTCACCCCGACCATTTTTGGGTCGTTAGCTCAGTTGGTAGAGCAg              >  2:165754/1‑150 (MQ=255)
                                                                                                                                      aaGGGCATCTCAAAGAGATGCCCTTTTTTTTGAAACACTCGGGGTGTAGCGCAGTCCGGTAGCGCGCCTGCTTTGGGAGCAGGATGTCAGGAGTTCGAATCCCCTCACCCCGACCATTTTTGGGTCGTTAGCTCAGTTGGTAGAGCAGtt            >  1:967954/1‑150 (MQ=255)
                                                                                                                                        gggCATCTCAAAGAGATGCCCTTTTTTTTGAAACACTCGGGGTGTAGCGCAGTCCGGTAGCGCGCCTGCTTTGGGAGCAGGATGTCAGGAGtt                                                                   >  2:239568/1‑93 (MQ=255)
                                                                                                                                        gggCATCTCAAAGAGATGCCCTTTTTTTTGAAACACTCGGGGTGTAGCGCAGTCCGGTAGCGCGCCTGCTTTGGGAGCAGGATGTCAGGAGtt                                                                   <  1:239568/93‑1 (MQ=255)
                                                                                                                                        gggCATCTCAAAGAGATGCCCTTTTTTTTGAAACACTCGGGGTGTAGCGCAGTCCGGTAGCGCGCCTGCTTTGGGAGCAGGATGTCAGGAGTTCGAATCCCCTCACCCCGACCAttttt                                         <  1:1182131/119‑1 (MQ=255)
                                                                                                                                        gggCATCTCAAAGAGATGCCCTTTTTTTTGAAACACTCGGGGTGTAGCGCAGTCCGGTAGCGCGCCTGCTTTGGGAGCAGGATGTCAGGAGTTCGAATCCCCTCACCCCGACCATTTTTGGGTCGTTAGCTCAGTTGGTAGAGCAGTTgg          >  1:1221496/1‑150 (MQ=255)
                                                                                                                                        gggCATCTCAAAGAGATGCCCTTTTTTTTGAAACACTCGGGGTGTAGCGCAGTCCGGTAGCGCGCATGCTTTGGGAGCAGGATGTCAGGAGTTCGAATCCCCTCAGCCCGACCAttttt                                         >  2:1182131/1‑119 (MQ=255)
                                                                                                                                            aTCTCAAAGAGATGCCCTTTTTTTTGAAACACTCGGGGTGTAGCGCAGTCCGGTAGCGCGCCTGCTTTGGGAGCAGGATGTCAGGAGTTCGAATCCCCTCACCCCGACCATTTTTGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGGCttt      >  1:202544/1‑150 (MQ=255)
                                                                                                                                            aTCTCAAAGAGATGCCCTTTTTTTTGAAACACTCGGGGTGTAGCGCAGTCCGGTAGCGCGCCTGCTTTGGGAGCAGGATGTCAGGAGTTCGAATCCCCTCACCCCGACCATTTTTGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGGCttt      <  1:1439390/150‑1 (MQ=255)
                                                                                                                                            aTCTCAAAGAGATGCCCTTTTTTTTGAAACACTCGGGGTGTAGCGCAGTCCGGTAGCGCGCCTGCTTTGGGAGCAGGATGTCAGGAGTTCGAATCCCCTCACCCCGACCATTTTTGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGGCgtt      >  1:1491569/1‑150 (MQ=255)
                                                                                                                                             tctcAAAGAGATGCCCTTTTTTTTGAAACACTCGGGGTGTAGCGAAGTCAGGTAGGTGGGTTGCTTTGGGAGCAGGATGTCAGGAGTTCGAATCCCCTCACCCCGACCATTTTTGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGGCtttt     >  1:1020943/1‑150 (MQ=18)
                                                                                                                                               tcAAAGAGATGCCCTTTTTTTTGAAACACTCGGGGTGTAGCGCAGTCCGGTAGCGCGCCTGCTTTGGGAGCAGGATGTCAGGAGTTCGAATCCCCTCACCCCGACCATTTTTGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGGCTTTTaa   <  1:156732/150‑1 (MQ=255)
                                                                                                                                                cAAAGAGATGCCCTTTTTTTTGAAACACTCGGGGTGTAGCGCAGTCCGGTAGCGCGCCTGCTTTGGGAGCAGGATGTCAGGAGTTCGAATCCCCTCACCCCGACCATTTTTGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGGCTTTTAAc  <  1:165754/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
TTTGGGAGCAGGATGTCAGGAGTTCGAATCCCCTCACCCCGACCATTTTCCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGGCTTTTAACCAATTGGTCGTAGGTTCGAATCCTACACGACCCACCATGTAGAAAAGGGCATCTCAAAGAGATGCCCTTTTTTTTGAAACACTCGGGGTGTAGCGCAGTCCGGTAGCGCGCCTGCTTTGGGAGCAGGATGTCAGGAGTTCGAATCCCCTCACCCCGACCATTTTTGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGGCTTTTAAC  >  NC_002947/733503‑733796

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