Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,965,396 C→T 5.3% D79N (GAT→AAT)  rhlB ← ATP‑dependent RNA helicase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,965,3960CT 5.3% 419.7 / 6.7 131D79N (GAT→AAT) rhlBATP‑dependent RNA helicase
Reads supporting (aligned to +/- strand):  ref base C (64/59);  new base T (3/4);  total (68/63)
Fisher's exact test for biased strand distribution p-value = 7.12e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.18e-01

Only 100 of 134 total aligned reads displayed.

CAGCTGTTTGTCGTAGCCATCACCACCGTAAGCCAGACCCAGCTTCAGGCCAGTAGCTTCCGCCAGCGGTTCTGCGTCGGCATGGATCTGCACGGCAAGTTCACGCGTCGGTGCCATAATTAAGGCACGCGGCTGATTCACCTTGCGATCGGCAATCGCAGGATGAGAGAGAAGATAATGAAACGTTGACGTAAGAAACGCCATCGTTTTCCCGGTACCGGTTTGCGCCTGCCCGGCTACGTCACGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAA  >  NC_000913/3965247‑3965545
                                                                                                                                                     |                                                                                                                                                     
  gCTGTTTGTCGTAGCCATCACCACCGTAAGCCAGACCCAGCTTCAGGCCAGTAGCTTCCGCCAGCGGTTCTGCGTCGGCATGGATCTGCACGGCAAGTTCACGCGTCGGTGCCATAATTAAGGCACGCGGCTGATTCACCTTGCGATCgg                                                                                                                                                     <  2:1539379/150‑1 (MQ=255)
         gTCGTAGCCATCACCACCGTAAGCCAGACCCAGCTTCAGGCCAGTAGCTTCCGCCAGCGGTTCTGCGTCGGCATGGATCTGCACGGCAAGTTCACGCGTCGGTGCCATAATTAAGGCACGCGGCTGATTCACCTTGCGATCGGCAATCGc                                                                                                                                              >  1:878532/1‑150 (MQ=255)
           cGTAGCCATCACCACCGTAAGCCAGACCCAGCTTCAGGCCAGTAGCTTCCGCCAGCGGTTCTGCGTCGGCATGGATCTGCACGGCAAGTTCACGCGTCGGTGCCATAATTAAGGCACGCGGCTGATTCACCTTGCGATCGGCAATCGCAg                                                                                                                                            >  2:867078/1‑150 (MQ=255)
           cGTAGCCATCACCACCGTAAGCCAGACCCAGCTTCAGGCCAGTAGCTTCCGCCAGCGGTTCTGCGTCGGCATGGATCTGCACGGCAAGTTCACGCGTCGGTGCCATAATTAAGGCACGCGGCTGATTCACCTTGCGATCGGCAATCGCAg                                                                                                                                            >  1:663952/1‑150 (MQ=255)
           cGTAGCCATCACCACCGTAAGCCAGACCCAGCTTCAGGCCAGTAGCTTCCGCCAGCGGTTCTGCGTCGGCATGGATCTGCACGGCAAGTTCACGCGTCGGTGCCATAATTAAGGCACGCGGCTGATTCACCTTGCGATCGGCAATCGCAg                                                                                                                                            <  1:1266292/150‑1 (MQ=255)
             tAGCCATCACCACCGTAAGCCAGACCCAGCTTCAGGCCAGAAGCTTCCGCCAGCGGTTCAGCGGCGGCAAGGATCTGCACGGCAAGTTCACGCGTCGGTGCCATAATTAAGGCACGCGGCTGATTCACCCTGCGATTGCCGATCGCaagg                                                                                                                                          <  2:826717/150‑4 (MQ=255)
              aGCCATCACCACCGTAAGCCAGACCCAGCTTCAGGCCAGTAGCTTCCGCCAGCGGTTCTGCGTCGGCATGGATCTGCACGGCAAGTTCACGCGTCGGTGCCATAATTAAGGCACGCGGCTGATTCACCTTGCGATTGCCGATCGCaaggt                                                                                                                                         <  1:1028170/150‑5 (MQ=255)
               gCCATCACCGCCGTAAGCCAGACCCAGCTTCAGGCCAGTAGCTTCCGCCAGCGGTTCTGCGTCGGCATGGATCTGCACGGCAAGTTCACGCGTCGGTGCCATAATTAAGGCACGCGGCTGATTCACCTTGCGATTGGCAATCGCAGGATg                                                                                                                                        <  2:739805/150‑1 (MQ=255)
               gCCATCACCACCGTAAGCCAGACCCAGCTTCAGGCCAGTAGCTTCCGCCAGCGGTTCTGCGTCGGCATGGATCTGCACGGCAAGTTCACGCGTCGGTGCCATAATTAAGGCACGCGGCTGATTCACCTTGCGATCGGCAATCGCAGGATg                                                                                                                                        >  2:1569739/1‑150 (MQ=255)
               gCCATCACCACCGTAAGCCAGACCCAGCTTCAGGCCAGTAGCTTCCGCCAGCGGTTCTGCGTCGGCATGGATCTGCACGGCAAGTTCACGCGTCGGTGCCATAATTAAGGCACGCGGCTGATTCACCTTGCGATCGGCAATCGCAGGATg                                                                                                                                        >  1:952646/1‑150 (MQ=255)
               gCCATCACCACCGTAAGCCAGACCCAGCTTCAGGCCAGTAGCTTCCGCCAGCGGTTCTGCGTCGGCATGGATCTGCACGGCAAGTTCACGCGTCGGTGCCATAATTAAGGCACGCGGCTGATTCACCTTGCGATCGGCAATCGCAGGATg                                                                                                                                        >  1:958808/1‑150 (MQ=255)
               gCCATCACCACCGTAAGCCAGACCCAGCTTCAGGCCAGTAGCTTCCGCCAGCGGTTCTGCGTCCGCATGGATCTGCACGGCAAGTTAACGCGTCGGTGCCATAAATAAGGCACGCGGCTGATTCACCTTGCGATCGGCAATCGCAGGATg                                                                                                                                        >  2:317842/1‑150 (MQ=255)
                   tCACCACCGTAAGCCAGACCCAGCTTCAGGCCAGTAGCTTCCGCCAGCGGTTCTGCGTCGGCATGGATCTGCACGGCAAGTTCACGCGTCGGTGCCATAATTAAGGCACGCGGCTGATTCACCTTGCGATCGGCAATCGCAGGATgagag                                                                                                                                    <  2:546918/150‑1 (MQ=255)
                    caccacCGTAAGCCAGACCCAGCTTCAGGCCAGTAGCTTCCGCCAGCGGTTCTGCGTCGGCATGGATCTGCACGGCTAGTTCACGCGTCGGTGACATAATTAAGGCCCGCTGCTGATTCACCTTGCGAGCGTCTATCTCAGGATgagaga                                                                                                                                   >  2:1089307/1‑150 (MQ=255)
                    caccacCGTAAGCCAGACCCAGCTTCAGGCCAGTAGCTTCCGCCAGCGGTTCTGCGTCGGCATGGATCTGCACGGCAAGTTCACGCGTCGGTGCCATAATTAAGGCACGCGGCTGATTCACCTTGCGATCGGCAATCGCAGGATgagaga                                                                                                                                   >  2:1712194/1‑150 (MQ=255)
                      ccaccGTAAGCCAGACCCAGGTTCAGGCCAGTAGCTTCCGCCAGCGGTTCTGCGTCGGCATGGATCTGCACGGCAAGTTCACGCGTGGGTGCCTTAATTAAGGCACGCGGCTGATTCACCTTGCGATCGGCAATCGCAGGATgagagaga                                                                                                                                 <  2:1616586/150‑1 (MQ=255)
                         ccGTAAGCCAGACCCAGCTTCAGGCCAGTAGCTTCCGCCAGCGGTTCTGCGTCGGCATGGATCTGCACGGCAAGTTCACGCGTCGGTGCCATAATTAAGGCACGCGGCTGATTCACCTTGCGATCGGCAATCGCAGGATGAGAGagaaga                                                                                                                              >  1:965757/1‑150 (MQ=255)
                          cGTAAGCCAGACCCAGCTTCAGGCCAGTAGCTTCCGCCAGCGGTTCTGCGTCGGCATGGATCTGCACGGCAAGTTCACGCGTCGGTGCCATAATTAAGGCACGCGGCTGATTCACCTTGCGATCGGCAATCGCAGGATGAGAGAGAAGAt                                                                                                                             <  1:669907/150‑1 (MQ=255)
                             aaGCCAGACCCAGCTTCAGGCCAGTAGCTTCCGCCAGCGGTTCTGCGTCGGCATGGATCTGCACGGCAAGTTCACGCGTCGGTGCCATAATTAAGGCACGCGGCTGATTCACCTTGCGATCGGCAATCGCAGGATGAGAGAGAAGATAAt                                                                                                                          <  1:548414/150‑1 (MQ=255)
                              cgCCAGCCCCAGCTTCAGGCCAGTCGCTTCCGCCAGCGGTTCTGCGTCGGCATGGATCTGCACGGGTAGTTCACGCGTCGGTGCCATAATTAAGGCACGCGGCTGATTCACCTTGCGATCGGCAATCGCAGGATGAGAGAGAAGATAATg                                                                                                                         <  2:1274975/149‑1 (MQ=255)
                              aGCCAGACCCAGCTTCAGGCCAGTAGCTTCCGCCAGCGGTTCTGCGTCGGCATGGATCTGCACGGCAAGTTCACGCGTCGGTGCCATAATTAAGGCACGCGGCTGATTCACCTTGCGATCGGCAATCGCAGGATGAGAGAGAAGATAATg                                                                                                                         <  2:967472/150‑1 (MQ=255)
                              aGCCAGACCCAGCTTCAGGCCAGTAGCTTCCGCCAGCGGTTCTGCGTCGGCATGGATCTGCACGGCAAGTTCACGCGTCGGTGCCATAATTAAGGCACGCGGCTGATTCACCTTGCGATCGGCAATCGCAGGATGAGAGAGAAGATAATg                                                                                                                         >  2:595178/1‑150 (MQ=255)
                              aGCCAGACCCAGCTTCAGGCCAGTAGCTTCCGCCAGCGGTTCTGCGTCGGCATGGATCTGCACGGCAAGTTCACGCGTCGGTGCCATAATTAAGGCACGCGGCTGATTCACCTTGCGATCGGCAATCGCAGGATGAGAGAGAAGATAATg                                                                                                                         <  1:1716369/150‑1 (MQ=255)
                              aGCCAGACCCAGCTTCAGGCCAGTAGCTTCCGCCAGCGGTTCTGCGTCGGCATGGATCTGCACGGCAAGTTCACGCGTCGGTGCCATAATTAAGGCACGCGGCTGATTCACCTTGCGATCGGCAATCGCAGGATGAGAGAGAAGATAATg                                                                                                                         >  1:1080280/1‑150 (MQ=255)
                               gCCAGACCCAGCTTCAGGCCAGTAGCTTCCGCCAGCGGTTCTGCGTCGGCATGGATCTGCACGGCAAGTTCACGCGTCGGTGCCATAATTAAGGCACGCGGCTGATTCACCTTGCGATCGGCAATCGCAGGATGAGAGAGAAGATAATGa                                                                                                                        >  1:1603219/1‑150 (MQ=255)
                                ccAGACCCAGCTTCAGGCCAGTAGCTTCCGCCAGCGGTTCTGCGTCGGCATGGATCTGCACGGCAAGTTCACGCGTCGGTGCCATAATTAAGGCACGCGGCTGATTCACCTTGCGATCGGCAATCGCAGGATGAGAGAGAAGATAATGaa                                                                                                                       >  2:1247169/1‑150 (MQ=255)
                                 ctgACGCAGCTTCAGGCCAGTAGCTTCCGCCAGCGGTTCTGCGTCGGCATGGATCTGCACGGCAAGTTCACGCGTCGGTGCCATAATTAAGGCACGCGGCTGATTCACCTTGCGATCGGCAATCGCAGGATGAGAGAGAAGATAATGaaa                                                                                                                      <  2:932430/148‑1 (MQ=255)
                                 cAGACCCAGCTTCAGGCCAGTAGCTTCCGCCAGCGGTTCTGCGTCGGCATGGATCTGCACGGCAAGTTCACGCGTCGGTGCCATAATTAAGGCACGCGGCTGATTCACCTTGCGATCGGCAATCGCAGGATGAGAGAGAAGATAATGaaa                                                                                                                      <  1:591325/150‑1 (MQ=255)
                                   gACCCAGCTTCAGGCCAGTAGCTTCCGCCAGCGGTTCTGCGTCGGCATGGATCTGCACGGCAAGTTCACGCGTCGGTGCCATAATTAAGGCACGCGGCTGATTCACCTTGCGATCGGCAATCGCAGGATGAGAGAGAAGATAATGAAACg                                                                                                                    <  1:1493833/150‑1 (MQ=255)
                                      ccAGCTTCAGGCCAGTAGCTTCCGCCAGCGGTTCTGCGTCGGCATGGATCTGCACGGCAAGTTCACGCGTCGGTGCCATAATTAAGGCACGCGGCTGATTCACCTTGCGATCGGCAATCGCAGGATGAGAGAGAAGATAATGAAACGTTg                                                                                                                 >  1:823873/1‑150 (MQ=255)
                                        aGCTTCAGGCCAGTAGCTTCCGCCAGCGGTTCTGCGTCGGCATGGATCTGCACGGCAAGTTCACGCGTCGGTGCCATAATTAAGGCACGCGGCTGATTCACCTTGCGATCGGCAATCGCAGGATGAGAGAGAAGATAATGAAACGTTGAc                                                                                                               >  2:1025097/1‑150 (MQ=255)
                                        aGCTTCAGGCCAGTAGCTTCCGCCAGCGGTTCTGCGTCGGCATGGATCTGCAAGGCAAGTTCACGCGTCGTTGCCATAATTAAGGCACGCGGCTGATTCACCTTGCGATCGGCAATCGCAGGATTAGAGAGAAGATAATGAAACGTTGAc                                                                                                               >  2:253553/1‑150 (MQ=255)
                                         gCTTCAGGCCAGTAGCTTCCGCCAGCGGTTCTGCGTCGGCATGGATCTGCACGGCAAGTTCACGCGGCGGTGCCAGAATTAAGGCACGCGGCAGCTTCACCTTGCGATCGGCAATCGCAGGATGAGAGAGAAGATAATGAAACGTTGACg                                                                                                              >  2:786906/1‑150 (MQ=255)
                                             cAGGCCAGTAGCTTCCGCCAGCGGTTCTGCGTCGGCATGGATCTGCCCGGCAAGTTCACGCGTCGGTGCCATAATTAAGGCACGCGGCTGATTCACCTTGCGATCGGCAATCGCAGGATGAGAGAGAAGATAATGAAACGTTGACGTAAg                                                                                                          >  2:1128497/1‑150 (MQ=255)
                                             cAGGCCAGTAGCTTCCGCCAGCGGTTCTGCGTCGGCATGGATCTGCACGGCAAGTTCACGCGTCGGTGCCATAATTAAGGCACGCGGCTGATTCACCTTGCGATCGGCAATCGCAGGATGAGAGAGAAGATAATGAAACGTTGACGTAAg                                                                                                          >  2:422590/1‑150 (MQ=255)
                                             cAGGCCAGTAGCTTCCGCCAGCGGTTCTGCGTCGGCATGGATCTGCACGGCAAGTTCACGCGTCGGTGCCATAATTAAGGCACGCGGCTGATTCACCTTGCGATCGGCAATCGCAGGATGAGAGAGAAGATAATGAAACGTTGACGTAAg                                                                                                          <  2:57837/150‑1 (MQ=255)
                                             cAGGCCAGTAGCTTCCGCCAGCGGTTCTGCGTCGGCATGGATCTGCACGGCAAGTTCACGCGTCGGTGCCATAATTAAGGCACGCGGCTGATTCACCTTGCGATCGGCAATCGCAGGATGAGAGAGAAGATAATGAAACGTTGACGTAAg                                                                                                          <  2:599830/150‑1 (MQ=255)
                                                 ccAGTAGCTTCCGCCAGCGGGTCTGCGTCGGCATGGATCTGCACGGCAAGTTCACGCGGCGGTGCCATAATTAAGGCACGCGGCTGATTCACCTTGCGATCGGCAATCGCAGGATGAGAGAGAAGATAATGAAACGTTGACGTAAGAAAc                                                                                                      <  2:1055046/150‑1 (MQ=255)
                                                  cAGTAGCTTCCGCCAGCGGTTCTGCGTCGGCATGGATCTGCACGGCAAGTTCACGCGTCGGTGCCATAATTAAGGCACGCGGCTGATTCACCTTGCGATCGGCAATCGCAGGATGAGAGAGAAGATAATGAAACGTTGACGTAAGAAACg                                                                                                     >  2:114074/1‑150 (MQ=255)
                                                  cAGTAGCTTCCGCCAGCGGTTCTGCGTCGGCATGGATCTGCACGGCAAGTTCACGCGTCGGTGCCATAATTAAGGCACGCGGCTGATTCACCTTGCGATCGGCAATCGCAGGATGAGAGAGAAGATAATGAAACGTTGACGTAAGAAACg                                                                                                     <  2:878532/150‑1 (MQ=255)
                                                   aGTAGCTTACGCAAGCGGTACTCCGTCGGCAAGGATCTGCACGGCCAGTTCACGCGTCGGTGAAATAATTAAGGCACGCGGCTGATTCACCTTGCGATAGGCATTCGGAGGATGAGAGAGAAGATAATGAAACGTTGACGAAAGTAACGc                                                                                                    >  2:1076216/1‑150 (MQ=255)
                                                            cGCCAGCGGTTCTGCGTCGGCATGGATCTGCACGGCAAGTTCACGCGTCGGTGCCATAATTAAGGCACGCGGCTGATTCACCTTGCGATCGGCAATCGCAGGATGAGAGAGAAGATAATGAAACGTTGACGTAAGAAACGCCATCGtttt                                                                                           >  1:503652/1‑150 (MQ=255)
                                                                aGCGGTTCTGCGTCGGCATGGATCTGCACGGCAAGTTCACGCGTCGGTGCCATAATTAAGGCACGCGGCTGATTCACCTTGCGATCGGCAATCGCAGGATGAGAGAGAAGATAATGAAACGTTGACGTAAGAAACGCCATCGTTTTCCCg                                                                                       >  1:1393358/1‑150 (MQ=255)
                                                                        tGCGTCGGCATGGATCTGCACGGCAAGTTCACGCGTCGGTGCCATAATTAAGGCACGCGGCTGATTCACCTTGCGATCGGCAATCGCAGGATGAGAGAGAAGATAATGAAACGTTGACGTAAGAAACGCCATCGTTTTCCCGGTACCGGt                                                                               >  1:1542507/1‑150 (MQ=255)
                                                                           gTCGGCATGGATCTGCACGGCAAGTTCACGCGTCGGTGCCATAATTAAGGCACGCGGCTGATTCACCTTGCGATCGGCAATCGCAGGATGAGAGAGAAGATAATGAAACGTTGACGTAAGAAACGCCATCGTTTTCCCGGTACCGGTTTg                                                                            <  1:226737/150‑1 (MQ=255)
                                                                               gCATGGATCTGCACGGCAAGTTCACGCGTCGGTGCCATAATTAAGGCACGCGGCTGATTCTCGTTGCGATCGGCAATCGCAGGATGAGAGAGAAGATAATGAAAGGTTGACGTTAGAAACGCCATCGTTTTCCCGGTACCGGTTTGCGcc                                                                        <  1:1179903/150‑1 (MQ=255)
                                                                               gCATGGATCTGCACGGCAAGTTCACGCGTCGGTGCCATAATTAAGGCACGCGGCTGATTCACCTTGCGATCGGCAATCGCAGGATGAGAGAGAAGATAATGAAACGTTGACGTAAGAAACGCCGTCGTTTTCCCGGTACCGGTTTGCGcc                                                                        <  1:44916/150‑1 (MQ=255)
                                                                               gCATGGATCTGCACGGCAAGTTCACGCGTCGGTGCCATAATTAAGGCACGCGGCTGATTCACCTTGCGATCGGCAATCGCAGGATGAGAGAGAAGATAATGAAACGTTGACGTAAGAAACGCCATCGTTTTCCCGGTACCGGTTTGCGcc                                                                        <  1:888929/150‑1 (MQ=255)
                                                                               gCATGGATCTGCACGGCAAGTTCACGCGTCGGTGCCATAATTAAGGCACGCGGCTGATTCACCTTGCGATCGGCAATCGCAGGATGAGAGAGAAGATAATGAAACGTTGACGTAAGAAACGCCATCGTTTTCCCGGTACCGGTTTGCGcc                                                                        <  1:132170/150‑1 (MQ=255)
                                                                                  tGGATCTGCACGGCAAGTTCACGCGTCGGTGCCATAATTAAGGCACGCGGCTGATTCACCTTGCGATCGGCAATCGCAGGATGAGAGAGAAGATAATGAAACGTTGACGTAAGAAACGCCATCGTTTTCCCGGTACCGGTTTGCGCCTGc                                                                     >  1:78348/1‑150 (MQ=255)
                                                                                  tGGATCTGCACGGCAAGTTCACGCGTCGGTGCCATAATTAAGGCACGCGGCTGATTCACCTTGCGATCGGCAATCGCAGGATGAGAGAGAAGATAATGAAACGTTGACGTAAGAAACGCCATCGTTTTCCCGGTACCGGTTTGCGCCTGc                                                                     <  1:789621/150‑1 (MQ=255)
                                                                                  tGGATCTGCACGGCAAGTTCACGCGTCGGTGCCATAATTAAGGCACGCGGCTGATTCACCTTGCGATCGGCAATCGCAGGATGAGAGAGAAGATAATGAAACGTTGACGTAAGAAACGCCATCGTTTTCCCGGTACCGGTTTGCGCCTGc                                                                     >  2:569474/1‑150 (MQ=255)
                                                                                  tGGATCTGCACGGCAAGTTCACGCGTCGGTGCCATAATTAAGGCACGCGGCTGATTCACCTTGCGATCGGCAATCGCAGGATGAGAGAGAAGATAATGAAACGTTGACGTAAGAAACGCCATCGTTTTCCCGGTACCGGTTTGCGCCTGc                                                                     >  1:1653189/1‑150 (MQ=255)
                                                                                  tGGATCTGCACGGCAAGTTCACGCGTCGGTGCCATAATTAAGGCACGCGGCTGATTCACCTTGCGATCGGCAATCGCAGGATGAGAGAGAAGATAATGAAACGTTGACGTAAGAAACGCCATCGTTTTCCCGGTACCGGTTTGCGCCTGc                                                                     <  1:428660/150‑1 (MQ=255)
                                                                                  tGGATCTGCACGGCAAGTTCACGCGTCGGTGCCATAATTAAGGCACGCGGCTGATTCACCTTGCGATCGGCAATCGCAGGATGAGAGAGAAGATAATGAAACGTTGACGTAAGAAACGCCATCGTTTTCCAGGTACCGGTTTGCGCCTGc                                                                     >  2:463830/1‑150 (MQ=255)
                                                                                    gATCTGCACGGCAAGTTCACGCGTCGGTGCCATAATTAAGGCACGCGGCTGATTCACCTTGCGATCGGCAATCGCAGGATGAGAGAGAAGATAATGAAACGTTGACGTAAGAAACGCCATCGTTTTCCCGGTACCGGTTTGCGCCTGccc                                                                   <  2:1028292/150‑1 (MQ=255)
                                                                                       cTGCACGGCAAGTTCACGCGGCGGTGCCATAATTAAGGCACGCGGCTGATTCACCTTGCGATCGGCAATCGCAGGATGAGAGAGAAGATAATGAAACGTTGACGTAAGAAACGCCATCGTTTTCCCGGTACCGGTTTGCGCCTGCCCGGc                                                                <  2:1201476/150‑1 (MQ=255)
                                                                                        tGCACGGCAAGTTCACGCGTCGGTGCCATAATTAAGGCACGCGGCTGATTCACCTTGCGATCGGCAATCGCAGGATGAGAGAGAAGATAATGAAACGTTGACGTAAGAAACGCCATCGTTTTCCCGGTACCGGTTTGCGCCTGCCCGGCt                                                               >  1:1254418/1‑150 (MQ=255)
                                                                                        tGCACGGCAAGTTCACGCGTCGGTGCCATAATTAAGGCACGCGGCTGATTCACCTTGCGATCGGCAATCGCAGGATGAGAGAGAAGATAATGAAACGTTGACGTAAGAAACGCCATCGTTTTCCCGGTACAGGTTTGCGCCTGCCCGGCt                                                               >  2:797469/1‑150 (MQ=255)
                                                                                        tGCACGGCAAGTTCACGCGTCGGTGCCATAATTAAGGCACGCGGCTGAATCACCTGGCGATCGGCAATCGCAGGATGAGAGAGAAGATAATGAAACGTTGACGAAAGAAACGCCATCGTTTTCCCGGTACCGGTTTGCGCCTGCCCGGCt                                                               <  2:831692/150‑1 (MQ=255)
                                                                                         gCACGGCAAGTTCACGCGTCGGTGCCATAATTAAGGCACGCGGCTGATTCACCTTGCGATCGGCAATCGCAGGATGAGAGAGAAGATAATGAAACGTTGACGTAAGAAACGCCATCGTTTTCCCGGTACCGGTTTGCGCCTGCCCGGCTa                                                              <  1:1712194/150‑1 (MQ=255)
                                                                                           aCGGCAAGTTCACGCGTCGGTGCCATAATTAAGGCACGCGGCTGATTCACCTTGCGATCGGCAATCGCAGGATGAGAGAGAAGATAATGAAACGTTGACGTAAGAAACGCCATCGTTTTCCCGGTACCGGTTTGCGCCTGCCCGGCTACg                                                            >  1:1343060/1‑150 (MQ=255)
                                                                                              gCAAGTTCACGCGTCGGTGCCATAATTAAGGCACGCGGCTGATTCACCTTGCGATCGGCAATCGCAGGATGAGAGAGAAGATAATGAAACGTTGACGTAAGAAACGCCATCGTTTTCCCGGTACCGGTTTGCGCCTGCCCGGCTACGTCa                                                         <  1:659820/150‑1 (MQ=255)
                                                                                                  gTTCACGCGTCGGTGCCATAATTAAGGCACGCGGCTGATTCACCTTGCGATCGGCAATCGCAGGATGAGAGAGAAGATAATGAAACGTTGACGTAAGAAACGCCATCGTTTTCCCGGTACCGGTTTGCGCCTGCCCGGCTACGTCACGAc                                                     >  2:610366/1‑150 (MQ=255)
                                                                                                     cACGCGTCGGTGCCATAATTAAGGCACGCGGCTGATTCACCTTGCGATCGGCAATCGCAGGATGAGAGAGAAGATAATGAAACGTTGACGTAAGAAACGCCATCGTTTTCCCGGTACCGGTTTGCGCCTGCCCGGCTACGTCACGACCCg                                                  >  2:1597540/1‑150 (MQ=255)
                                                                                                     cACGCGTCGGTGCCATAATTAAGGCACGCGGCTGATTCACCTTGCGATCGGCAATCGCAGGATGAGAGAGAAGATAATGAAACGTTGACGTAAGAAACGCCATCGTTTTCCCGGTACCGGTTTGCGCCTGCCCGGCTACGTCACGACCCg                                                  >  2:76856/1‑150 (MQ=255)
                                                                                                     cACGCGTCGGTGCCATAATTAAGGCACGCGGCTGATTCACCTTGCGATCGGCAATCGCAGGATGAGAGAGAAGATAATGAAACGTTGACGTAAGAAACGCCATCGTTTTCCCGGTACCGGTTTGCGCCTGCCCGGCTACGTCACGACCCg                                                  >  1:396668/1‑150 (MQ=255)
                                                                                                      aCGCGTCGGTGCCATAATTAAGGCACGCGGCTGATTCACCTTGCGATCGGCAATCGCAGGATGAGAGAGAAGATAATGAAACGTTGACGTAAGAAACGCCATCGTTTTCCCGGTACCGGTTTGCGCCTGCCCGGCTACGTCACGACCCGc                                                 >  1:1042064/1‑150 (MQ=255)
                                                                                                       cgcgTCGGTGCCATAATTAAGGCACGCGGCTGATTCACCTTGCGATCGGCAATCGCAGGATGAGAGAGAAGATAATGAAACGTTGACGTAAGAAACGCCAGCGTTTTCCCGGTACCGGTTTGCGCCTGCCCGGCTACGTCACGACCCGcc                                                <  1:867078/150‑1 (MQ=255)
                                                                                                       cgcgACGGCGCCATAATTAAGGCACGCGGCTGATTCACCTTGCGATCGGCAATCGCAGGATGAGAGAGAAGATAATGAAACGTTGACGTAAGAAACGCCATCGTTTTCCCGGTACCGGTTTGCGCCTGCCCGGCTACGTCACGACCCGcc                                                <  1:180678/150‑1 (MQ=255)
                                                                                                         cgTCGGTGCCATAATTAAGGCACGCGGCTGATTCACCTTGCGATCGGCAATCGCAGGATGAGAGAGAAGATAATGAAACGTTGACGTAAGAAACGCCATCGTTTTCCCGGTACCGGTTTGCGCCTGCCCGGCTACGTCACGACCCGCCAg                                              >  2:1531727/1‑150 (MQ=255)
                                                                                                          gTCGGTGCCATAATTAAGGCACGCGGCTGATTCACCTTGCGATCGGCAATCGCAGGATGAGAGAGAAGATAATGAAACGTTGACGTAAGAAACGCCATCGTTTTCCCGGTACCGGTTTGCGCCTGCCCGGCTACGTCACGACCCGCCAGc                                             >  2:1676313/1‑150 (MQ=255)
                                                                                                             ggTGCCATAATTAAGGCACGCGGCTGATTCACCTTGCGATCGGCAATCGCAGGATGAGAGAGAAGATAATGAAACGTTGACGTAAGAAACGCCATCGTTTTCCCGGTACCGGTTTGCGCCTGCCCGGCTACGTCACGACCCGCCAGCGTc                                          <  1:499302/150‑1 (MQ=255)
                                                                                                                gCCATAATTAAGGCACGCGGCTGATTCACCTTGCGATCGGCAATCGCAGGATGAGAGAGAAGATAATGAAACGTTGACGTAAGAAACGCCATCGTTTTCCCGGTACCGGTTTGCGCCTGCCCGGCTACGTCACGACCCGCCAGCGTCAGc                                       <  2:1393358/150‑1 (MQ=255)
                                                                                                                  cATAATTAAGGCACGCGGCTGATTCACCTTGCGATCGGCAATCGCAGGATGAGAGAGAAGATAATGAAACGTTGACGTAAGAAACGCCATCGTTTTCCCGGTACCGGTTTGCGCCTGCCCGGCTACGTCACGACCCGCCAGCGTCAGCgg                                     >  2:891229/1‑150 (MQ=255)
                                                                                                                     aaTTAAGGCACGCGGCTGATTCACCTTGCGATCGGCAATCGCAGGATGAGAGAGAAGATAATGAAACGTTGACGTAAGAAACGCCATCGTTTTCCCGGTACCGGTTTGCGCCTGCCCGGCTACGTCACGACCCGCCAGCGTCAGCGGAAg                                  >  1:1206063/1‑150 (MQ=255)
                                                                                                                      aTTAAGGCACGCGGCTGATTCACCTTGCGATCGGCAATCGCAGGATGAGAGAGAAGATAATGAAACGTTGACGTAAGAAACGCCATCGTTTTCCCGGTACCGGTTTGCGCCTGCCCGGCTACGTCACGACCCGCCAGCGTCAGCGGAAgg                                 >  2:378479/1‑150 (MQ=255)
                                                                                                                         aaGGCACGCGGCTGATTCACCTTGCGATCGGCAATCGCAGGATGAGAGAGAAGATAATGAAACGTTGACGTAAGAAACGCCATCGTTTTCCCGGTACCGGTTTGCGCCTGCCCGGCTACGTCACGACCCGCCAGCGTCAGCGGAAGGGcc                              <  1:317842/150‑1 (MQ=255)
                                                                                                                              aCGCGGCTGATTCCCCTTGCGATCGGCAATCGCAGGATGAGAGAGCAGATCATGAAACGGTGACGTAAGAAACGCCATCGGTTTCCCGGTACCGGTTTGCGCCTGCCCGGCTACGTCCCGACCCGCCAGCGTCAGCGGAAGGGCCAGTGc                         <  1:1129429/150‑1 (MQ=255)
                                                                                                                              aCGCGGCTGATTCACCTTGCGATCGGCAATCGCAGGATGAGAGAGAAGATAATGAAACGTTGACGTAAGAAACGCCATCGTTTTCCCGGTACCGGTTTGCGCCTGCCCGGCTACGTCACGACCCGCCAGCGTCAGCGGAAGGGCCAGTGc                         <  2:1603219/150‑1 (MQ=255)
                                                                                                                                   gCTGATTCACCTTGCGATCGGCAATCGCAGGATGAGAGAGAAGATAATGAAACGTTGACGTAAGAAACGCCATCGTTTTCCCGGTACCGGTTTGCGCCTGCCCGGCTACGTCACGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGaa                    >  1:1673926/1‑150 (MQ=255)
                                                                                                                                      gATTCACCTTGCGATCGGCAATCGCAGGATGAGAGAGAAGATAATGAAACGTTGACGTAAGAAACGCCATCGTTTTCCCGGTACCGGTTTGCGCCTGCCCGGCTACGTCACGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATgg                 >  1:573921/1‑150 (MQ=255)
                                                                                                                                        ttCACCTTGCGATCGGCTATCGCAGGATGAGAGAGAAGATAATGAACCTCTGACGTAAGAAACGCCATCGTTTTCCCGGTTCCGGTGAGCGCCGGCCCGGCTACGTCACGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGc               <  1:1142542/150‑1 (MQ=255)
                                                                                                                                        ttCACCTTGCGATCGGCAATCGCAGGATGAGAGAGAAGATAATGAAACGTTGACGTAAGAAACGCCATCGTTTTCCCGGTACCGGTTTGCGCCTGCCCGGCTACGTCACGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGc               <  1:1521627/150‑1 (MQ=255)
                                                                                                                                          catcctGCGATTGCCGATCGCAGGATGAGAGAGAAGATAATGAAACGTTGACGTAAGAAACGCCATCGTTTTCCCGGTACCGGTTTGCGCCTCCCCGGCTACGACACGAACCGCCAGCGTCCGCGGAAGGGCCAGAGCCCGAATCGGcca             >  1:259329/6‑148 (MQ=255)
                                                                                                                                           atcctGCGATTGCCGATCGCAGGATGAGAGAGAAGATAATGAAACGTTGACGTAAGAAACGCCATCGTTTTCCCGGTACCGGTTTGCGCCTGCCCGGCTACGTCACGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTa            >  1:1697770/5‑150 (MQ=255)
                                                                                                                                           atcctGCGATTGCCGATCGCAGGATGAGAGAGAAGATAATGAAACGTTGACGTAAGAAACGCCATCGTTTTCCCGGTACCGGTTTGCGCCTGCCCGGCTACGTCACGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTa            >  1:1604554/5‑150 (MQ=255)
                                                                                                                                            ccTTGCGATCGGCAATCGCAGGATGAGAGAGAAGATAATGAAACGTTGACGTAAGAAACGCCATCGTTTTCCCGGTACCGGTTTGCGCCTGCCCGGCTACGTCACGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTAc           >  2:1651828/1‑150 (MQ=255)
                                                                                                                                            ccTTGCGATCGGCAATCGCAGGATGAGAGAGAAGATAATGAAACGTTGACGTAAGAAACGCCATAGTTTTCCCGGTACCGGTTTGCGCCTGCCCGGCTACGTCACGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTAc           <  1:1128497/150‑1 (MQ=255)
                                                                                                                                              ttGCGATCGGCAATCGCAGGATGAGAGAGAAGATAATGAAACGTTGACGTAAGAAACGCCATCGTTTTCCCGGTACCGGTTTGCGCCTGCCCGGCTACGTCACGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAg         <  2:1299948/150‑1 (MQ=255)
                                                                                                                                              ttGCGATCGGCAATCGCAGGATGAGAGAGAAGATAATGAAACGTTGACGTAAGAAACGCCATCGTTTTCCCGGTACCGGTTTGCGCCTGCCCGGCTACGTCACGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAg         <  2:1132212/150‑1 (MQ=255)
                                                                                                                                              ttGCGATCGGCAATCGCAGGATGAGAGAGAAGATAATGAAACGTTGACGTAAGAAACGCCATCGTTTTCCCGGTACCGGTTTGCGCCTGCCCGGCTACGTCACGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAg         >  2:205074/1‑150 (MQ=255)
                                                                                                                                               tGCGATCGGCAATCGCAGGATGAGAGAGAAGATAATGAAACGTTGACGTAAGAAACGCCATCGTTTTCCCGGTACCGGTTTGCGCCTGCCCGGCTACGTCACGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGt        <  1:1386359/150‑1 (MQ=255)
                                                                                                                                               tGCGATCGGCAATCGCAGGATGAGAGAGAAGAAAATGAAACGTTGACGTAAGAAACGCCATCGTTTTCCCGGTACCGGTTTGCGCCTGCCCGGCTACGTCACGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGt        <  1:1185852/150‑1 (MQ=255)
                                                                                                                                                 cGATCGGCAATCGCAGGATGAGAGAGAAGATAATGAAACGTTGACGTAAGAAACGCCATCGTTTTCCCGGTACCGGTTTGCGCCTGCCCGGCTACGTCACGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTa      <  1:1531727/150‑1 (MQ=255)
                                                                                                                                                 cGATCGGCAATCGCAGGATGAGAGAGAAGATAATGAAACGTTGACGTAAGAAACGCCATCGTTTTCCCGGTACCGGTTTGCGCCTGCCCGGCTACGTCACGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTa      <  2:319383/150‑1 (MQ=255)
                                                                                                                                                  gATCGGCAATCGCAGGATGAGAGAGAAGATAATGAAACGTTGACGTAAGAACCGCCATCGGTTTCCCGGTACCGGTTTGCGCATGCCCGGCTACGTCACGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCATGAATGGGAGTACAGTTAt     >  2:263144/1‑150 (MQ=255)
                                                                                                                                                  gATCGGCAATCGCAGGATGAGAGAGAAGATAATGAAACGTTGACGTAAGAAACGCCATCGTTTTCCCGGTACCGGTTTGCGCCTGCCCGGCTACGTCACGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTAt     >  2:332497/1‑150 (MQ=255)
                                                                                                                                                   aTCGGCAATCGCAGGATGAGAGAGAAGATAATGAAACGTTGACGTAAGAACCGCCATCGTTTTCCCGGTACCGGTTTGCGCCTGCCCGGCTACGTCACGACCCGCCAGCTTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATg    <  1:32434/150‑1 (MQ=255)
                                                                                                                                                     cGGCAATCGCAGGATGAGAGAGAAGATAATGAAACTTAGACGTAAGAAACGCCATAGTTTTCCCGGTACCGGTTTGCGCCTGCCCGGCTACGTCACGACCCGCCAGCGTCCGCGGAAGGGCCAGTGCATGAATGGGCGTACAGTTATGaa  >  2:1162801/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
CAGCTGTTTGTCGTAGCCATCACCACCGTAAGCCAGACCCAGCTTCAGGCCAGTAGCTTCCGCCAGCGGTTCTGCGTCGGCATGGATCTGCACGGCAAGTTCACGCGTCGGTGCCATAATTAAGGCACGCGGCTGATTCACCTTGCGATCGGCAATCGCAGGATGAGAGAGAAGATAATGAAACGTTGACGTAAGAAACGCCATCGTTTTCCCGGTACCGGTTTGCGCCTGCCCGGCTACGTCACGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAA  >  NC_000913/3965247‑3965545

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