Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,965,400 A→G 5.2% I77I (ATT→ATC)  rhlB ← ATP‑dependent RNA helicase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,965,4000AG 5.2% 409.7 / 7.5 135I77I (ATT→ATC) rhlBATP‑dependent RNA helicase
Reads supporting (aligned to +/- strand):  ref base A (69/57);  new base G (4/3);  total (74/61)
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.04e-01

Only 100 of 136 total aligned reads displayed.

TGTCGTAGCCATCACCACCGTAAGCCAGACCCAGCTTCAGGCCAGTAGCTTCCGCCAGCGGTTCTGCGTCGGCATGGATCTGCACGGCAAGTTCACGCGTCGGTGCCATAATTAAGGCACGCGGCTGATTCACCTTGCGATCGGCAATCGCAGGATGAGAGAGAAGATAATGAAACGTTGACGTAAGAAACGCCATCGTTTTCCCGGTACCGGTTTGCGCCTGCCCGGCTACGTCACGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCC  >  NC_000913/3965255‑3965549
                                                                                                                                                 |                                                                                                                                                     
tGTCGTAGCCCTCACCACCGTAAGCCAGACCCAGCTTCAGGCCAGTAGCTTCCGCCAGCGGTTCTGCGTCGGCATGGAGCTGCACGGCCAGATCACGCGTCGCTGCCATAATTAAGGCACGCGGCGGATTCCCCTTGCGATCGGAAATCg                                                                                                                                                   >  2:1522691/1‑150 (MQ=255)
 gTCGTAGCCATCACCACCGTAAGCCAGACCCAGCTTCAGGCCAGTAGCTTCCGCCAGCGGTTCTGCGTCGGCATGGATCTGCACGGCAAGTTCACGCGTCGGTGCCATAATTAAGGCACGCGGCTGATTCACCTTGCGATCGGCAATCGc                                                                                                                                                  >  1:878532/1‑150 (MQ=255)
   cGTAGCCATCACCACCGTAAGCCAGACCCAGCTTCAGGCCAGTAGCTTCCGCCAGCGGTTCTGCGTCGGCATGGATCTGCACGGCAAGTTCACGCGTCGGTGCCATAATTAAGGCACGCGGCTGATTCACCTTGCGATCGGCAATCGCAg                                                                                                                                                >  2:867078/1‑150 (MQ=255)
   cGTAGCCATCACCACCGTAAGCCAGACCCAGCTTCAGGCCAGTAGCTTCCGCCAGCGGTTCTGCGTCGGCATGGATCTGCACGGCAAGTTCACGCGTCGGTGCCATAATTAAGGCACGCGGCTGATTCACCTTGCGATCGGCAATCGCAg                                                                                                                                                <  1:1266292/150‑1 (MQ=255)
     tAGCCATCACCACCGTAAGCCAGACCCAGCTTCAGGCCAGTAGCATCCGCCAGCGGTTCTGCGTCGGCATGGATCTGCACGGCAAGTTCACGCGTCGGTGCCATAATTAAGGCACGCGGCTGATTCACCTTGCGATTGCCGATCGCaagg                                                                                                                                              <  2:351785/150‑4 (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)
         cATCACCACCGTAAGCCAGACCCAGCTTCAGGCCAGTAGCTTCCGCCAGCGGTTCTGCGTCGGCATGGATCTGCACGGCAAGTTCACGCGTCGGTGCCATAATTAAGGCACGCGGCTGATTCACCTTGCGATCGGCAATCGCAGGATgag                                                                                                                                          >  1:1621769/1‑150 (MQ=255)
          aTCACCACCGTAAGCCAGACCCAGCTTCAGGCCAGTAGCTTCCGCCAGCGGTTCTGCGTCGGCATGGATCTGCACGGCAAGTTCACGCGTCGGTGCCATAATGAAGGCACGCGGCTGATTCACCTTGCGATCGGCAATCGCAGGATgaga                                                                                                                                         >  1:52834/1‑150 (MQ=255)
            caccacCGTAAGCCAGACCCAGCTTCAGGCCAGTAGCTTCCGCCAGCGGTTCTGCGTCGGCATGGATCTGCACGGCAAGTTCACGCGTCGGTGCCATAATTAAGGCACGCTGCTGATTCACCTAGCGATCGGCAATCGCAGGATgagaga                                                                                                                                       >  2:1597450/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)
                    tAAGCCAGACCCAGCTTCAGGCCAGTAGCTTCCGCCCGCGGTTCTGCGGCGGCATGGATCTGCACGGCAAGTTCCCGCGGCGGTGCCATAATTAAGGCACGCGGCTGATTCACCTTGCGATCGGCAATCGCAGGATGAGAGAGAAGATaa                                                                                                                               <  2:301476/150‑1 (MQ=255)
                     aaGCCAGACCCAGCTTCAGGCCAGTAGCTTCCGCCAGCGGTTCTGCGTCGGCATGGATCTGCACGGCAAGTTCACGCGTCGGTGCCATAATTAAGGCACGCGGCTGATTCACCTTGCGATCGGCAATCGCAGGATGAGAGAGAAGATAAt                                                                                                                              <  1:548414/150‑1 (MQ=255)
                      cgCCAGCCCCAGCTTCAGGCCAGTCGCTTCCGCCAGCGGTTCTGCGTCGGCATGGATCTGCACGGGTAGTTCACGCGTCGGTGCCATAATTAAGGCACGCGGCTGATTCACCTTGCGATCGGCAATCGCAGGATGAGAGAGAAGATAATg                                                                                                                             <  2:1274975/149‑1 (MQ=255)
                      aGCCAGACCCAGCTTCAGGCCAGTAGCTTCCGCCAGCGGTTCTGCGTCGGCATGGATCTGCACGGCAAGTTCACGCGTCGGTGCCATAATTAAGGCACGCGGCTGATTCACCTTGCGATCGGCAATCGCAGGATGAGAGAGAAGATAATg                                                                                                                             >  1:319383/1‑150 (MQ=255)
                      aGCCAGACCCAGCTTCAGGCCAGTAGCTTCCGCCAGCGGTTCTGCGTCGGCATGGATCTGCACGGCAAGTTCACGCGTCGGTGCCATAATTAAGGCACGCGGCTGATTCACCTTGCGATCGGCAATCGCAGGATGAGAGAGAAGATAATg                                                                                                                             >  1:1080280/1‑150 (MQ=255)
                      aGCCAGACCCAGCTTCAGGCCAGTAGCTTCCGCCAGCGGTTCTGCGTCGGCATGGATCTGCACGGCAAGTTCACGCGTCGGTGCCATAATTAAGGCACGCGGCTGATTCACCTTGCGATCGGCAATCGCAGGATGAGAGAGAAGATAATg                                                                                                                             <  2:967472/150‑1 (MQ=255)
                      aGCCAGACCCAGCTTCAGGCCAGTAGCTTCCGCCAGCGGTTCTGCGTCGGCATGGATCTGCACGGCAAGTTCACGCGTCGGTGCCATAATTAAGGCACGCGGCTGATTCACCTTGCGATCGGCAATCGCAGGATGAGAGAGAAGATAATg                                                                                                                             <  1:1716369/150‑1 (MQ=255)
                       gCCAGACCCAGCTTCAGGCCAGTAGCTTCCGCCAGCGGTTCTGCGTCGGCATGGATCTGCACGGCAAGTTCACGCGTCGGTGCCATAATTAAGGCACGCGGCTGATTCACCTTGCGATCGGCAATCGCAGGATGAGAGAGAAGATAATGa                                                                                                                            >  1:1603219/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)
                                 gCTTCATGCCAGTAGCTTCCGCCAGCGGTTCTGCGTCGGAATGGATCTGCACGGAAAGTTCACGCGTCGGTGCCATAATTAAGGCACGCGGCTGATTCACCTTGCGATCGGCAATCGCAGGATGAGAGAGAAGATAATGAAACGTTGACg                                                                                                                  <  1:1542287/150‑1 (MQ=255)
                                     cAGGCCAGTAGCTTCCGCCAGCGGTTCTGCGTCGGCATGGATCTGCCCGGCAAGTTCACGCGTCGGTGCCATAATTAAGGCACGCGGCTGATTCACCTTGCGATCGGCAATCGCAGGATGAGAGAGAAGATAATGAAACGTTGACGTAAg                                                                                                              >  2:1128497/1‑150 (MQ=255)
                                     cAGGCCAGTAGCTTCCGCCAGCGGTTCTGCGTCGGCATGGATCTGCACGGCAAGTTCACGCGTCGGTGCCATAATTAAGGCACGCGGCTGATTCACCTTGCGATCGGCGATCGCAGGATGAGAGAGAAGATAATGAAACGTTGACGTAAg                                                                                                              >  2:589769/1‑150 (MQ=255)
                                     cAGGCCAGTAGCTTCCGCCAGCGGTTCTGCGTCGGCATGGATCTGCACGGCAAGTTCACGCGTCGGTGCCATAATTAAGGCACGCGGCTGATTCACCTTGCGATCGGCAATCGCAGGATGAGAGAGAAGATAATGAAACGTTGACGTAAg                                                                                                              <  2:599830/150‑1 (MQ=255)
                                     cAGGCCAGTAGCTTCCGCCAGCGGTTCTGCGTCGGCATGGATCTGCACGGCAAGTTCACGCGTCGGTGCCATAATTAAGGCACGCGGCTGATTCACCTTGCGATCGGCAATCGCAGGATGAGAGAGAAGATAATGAAACGTTGACGTAAg                                                                                                              >  2:422590/1‑150 (MQ=255)
                                     cAGGCCAGTAGCTTCCGCCAGCGGTTCTGCGTCGGCATGGATCTGCACGGCAAGTTCACGCGTCGGTGCCATAATTAAGGCACGCGGCTGATTCACCTTGCGATCGGCAATCGCAGGATGAGAGAGAAGATAATGAAACGTTGACGTAAg                                                                                                              <  2:57837/150‑1 (MQ=255)
                                     cAGGCCAGTAGCTTCCGCCAGCGGTTCTGCGTCGGCATGGATCTGCACGGAAAGTTCACGCGGCGGTGCCATAATTAAGGCACGCGGCTGATGCACCTTGCGATCGGCAATCGCAGGATGAGAGAGAAGATAATGAAACGTTGACGTAAg                                                                                                              <  1:1438571/150‑1 (MQ=255)
                                         ccAGTAGCTTCCGCCAGCGGGTCTGCGTCGGCATGGATCTGCACGGCAAGTTCACGCGGCGGTGCCATAATTAAGGCACGCGGCTGATTCACCTTGCGATCGGCAATCGCAGGATGAGAGAGAAGATAATGAAACGTTGACGTAAGAAAc                                                                                                          <  2:1055046/150‑1 (MQ=255)
                                          cAGTAGCTTCCGCCAGCGGTTCTGCGTCGGCATGGATCTGCACGGCAAGTTCACGCGTCGGTGCCATAATTAAGGCACGCGGCTGATTCACCTTGCGATCGGCAATCGCAGGATGAGAGAGAAGATAATGAAACGTTGACGTAAGAAACg                                                                                                         <  2:1270619/150‑1 (MQ=255)
                                          cAGTAGCTTCCGCCAGCGGTTCTGCGTCGGCATGGATCTGCACGGCAAGTTCACGCGTCGGTGCCATAATTAAGGCACGCGGCTGATTCACCTTGCGATCGGCAATCGCAGGATGAGAGAGAAGATAATGAAACGTTGACGTAAGAAACg                                                                                                         <  2:878532/150‑1 (MQ=255)
                                           aGTAGCTTACGCAAGCGGTACTCCGTCGGCAAGGATCTGCACGGCCAGTTCACGCGTCGGTGAAATAATTAAGGCACGCGGCTGATTCACCTTGCGATAGGCATTCGGAGGATGAGAGAGAAGATAATGAAACGTTGACGAAAGTAACGc                                                                                                        >  2:1076216/1‑150 (MQ=255)
                                             tAGCTTCCGCCAGCGGTTCTGCGTCGGCATGGATCTGCACGGCAAGTTCACGCGTCGGTGCCATAATTAAGGCACGCGGCTGATTCACCTTGCGATCGGCAATCGCAGGATGAGAGAGAAGATAATGAAACGTTGACGTAAGAAACGCCa                                                                                                      <  2:172222/150‑1 (MQ=255)
                                                     gCCAGCGGTTCTGCGTCGGCATGGATCTGCACGGCAAGTTCACGCGTCGGTGCCATAATTAAGGCACGCGGCTGATTCACCTTGCGATCGGCAATCGCAGGATGAGAGAGAAGATAATGAAACGTTGACGTAAGAACCGCCATCGTTTTc                                                                                              >  1:612642/1‑150 (MQ=255)
                                                        aGCGGTTCTGCGTCGGCATGGATCTGCACGGCAAGTTCACGCGTCGGTGCCATAATTAAGGCACGCGGCTGATTCACCTTGCGATCGGCAATCGCAGGATGAGAGAGAAGATAATGAAACGTTGACGTAAGAAACGCCATCGTTTTCCCg                                                                                           >  1:1393358/1‑150 (MQ=255)
                                                               cTGCGTCGGCATGGATCTGCACGGCAAGTTGACGCGTCGGTGCCATAATTAAGGCACGCGGCTGATTCACCTTGCGATCGTCAATCGCAGGATGAGAGAGAAGATAATGAAACGTTGACGTAAGAAACGCCATCGTTTTCCCGGTACCgg                                                                                    <  2:1094094/150‑1 (MQ=255)
                                                                   gTCGGCATGGATCTGCACGGCAAGTTCACGCGTCGGTGCCATAATTAAGGCACGCGGCTGATTCACCTTGCGATCGGCAATCGCAGGATGAGAGAGAAGATAATGAAACGTTGACGTAAGAAACGCCATCGTTTTCCCGGTACCGGTTTg                                                                                <  1:226737/150‑1 (MQ=255)
                                                                       gCATGGATCTGCACGGCAAGTTCACGCGTCGGTGCCATAATTAAGGCACGCGGCTGATTCTCGTTGCGATCGGCAATCGCAGGATGAGAGAGAAGATAATGAAAGGTTGACGTTAGAAACGCCATCGTTTTCCCGGTACCGGTTTGCGcc                                                                            <  1:1179903/150‑1 (MQ=255)
                                                                       gCATGGATCTGCACGGCAAGTTCACGCGTCGGTGCCATAATTAAGGCACGCGGCTGATTCACCTTGCGATCGGCAATCGCAGGATGAGAGAGAAGATAATGAAACGTTGACGTAAGAAACGCCGTCGTTTTCCCGGTACCGGTTTGCGcc                                                                            <  1:44916/150‑1 (MQ=255)
                                                                          tGGATCTGCACGGCAAGTTCACGCGTCGGTGCCATAATTAAGGCACGCGGCTGATTCACCTTGCGATCGGCAATCGCAGGATGAGAGAGAAGATAATGAAACGTTGACGTAAGAAACGCCATCGTTTTCCCGGTACCGGTTTGCGCCTGc                                                                         >  2:569474/1‑150 (MQ=255)
                                                                          tGGATCTGCACGGCAAGTTCACGCGTCGGTGCCATAATTAAGGCACGCGGCTGATTCACCTTGCGATCGGCAATCGCAGGATGAGAGAGAAGATAATGAAACGTTGACGTAAGAAACGCCATCGTTTTCCCGGTACCGGTTTGCGCCTGc                                                                         >  1:78348/1‑150 (MQ=255)
                                                                          tGGATCTGCACGGCAAGTTCACGCGTCGGTGCCATAATTAAGGCACGCGGCTGATTCACCTTGCGATCGGCAATCGCAGGATGAGAGAGAAGATAATGAAACGTTGACGTAAGAAACGCCATCGTTTTCCAGGTACCGGTTTGCGCCTGc                                                                         >  2:463830/1‑150 (MQ=255)
                                                                               cTGCACGGCAAGTTCACGCGGCGGTGCCATAATTAAGGCACGCGGCTGATTCACCTTGCGATCGGCAATCGCAGGATGAGAGAGAAGATAATGAAACGTTGACGTAAGAAACGCCATCGTTTTCCCGGTACCGGTTTGCGCCTGCCCGGc                                                                    <  2:1201476/150‑1 (MQ=255)
                                                                                tGCACGGCAAGTTCACGCGTCGGTGCCATAATTAAGGCACGCGGCTGATTCACCTTGCGATCGGCAATCGCAGGATGAGAGAGAAGATAATGAAACGTTGACGTAAGAAACGCCATCGTTTTCCCGGTACCGGTTTGCGCCTGCCCGGCt                                                                   >  1:440277/1‑150 (MQ=255)
                                                                                tGCACGGCAAGTTCACGCGTCGGTGCCATAATTAAGGCACGCGGCTGATTCACCTTGCGATCGGCAATCGCAGGATGAGAGAGAAGATAATGAAACGTTGACGTAAGAAACGCCATCGTTTTCCCGGTACCGGTTTGCGCCTGCCCGGCt                                                                   <  2:692365/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:890962/1‑150 (MQ=255)
                                                                                      gCAAGTTCACGCGTCGGTGCCATAATTAAGGCACGCGGCTGATTCACCTTGCGATCGGCAATCGCAGGATGAGAGAGAAGATAATGAAACGTTGACGTAAGAAACGCCATCGTTTTCCCGGTACCGGTTTGCGCCTGCCCGGCTACGTCa                                                             <  1:659820/150‑1 (MQ=255)
                                                                                          gTTCACGCGTCGGTGCCATAATTAAGGCACGCGGCTGATTCACCTTGCGATCGGCAATCGCAGGATGAGAGAGAAGATAATGAAACGTTGACGTAAGAAACGCCATCGTTTTCCCGGTACCGGTTTGCGCCTGCCCGGCTACGTCACGAc                                                         >  2:610366/1‑150 (MQ=255)
                                                                                             cACGCGTCGGTGCCATAATTAAGGCACGCGGCTGATTCACCTTGCGATCGGCAATCGCAGGATGAGAGAGAAGATAATGAAACGTTGACGTAAGAAACGCCATCGTTTTCCCGGTACCGGTTTGCGCCTGCCCGGCTACGTCACGACCCg                                                      >  1:396668/1‑150 (MQ=255)
                                                                                             cACGCGTCGGTGCCATAATTAAGGCACGCGGCTGATTCACCTTGCGATCGGCAATCGCAGGATGAGAGAGAAGATAATGAAACGTTGACGTAAGAAACGCCATCGTTTTCCCGGTACCGGTTTGCGCCTGCCCGGCTACGTCACGACCCg                                                      >  2:537715/1‑150 (MQ=255)
                                                                                             cACGCGTCGGTGCCATAATTAAGGCACGCGGCTGATTCACCTTGCGATCGGCAATCGCAGGATGAGAGAGAAGATAATGAAACGTTGACGTAAGAAACGCCATCGTTTTCCCGGTACCGGTTTGCGCCTGCCCGGCTACGTCACGACCCg                                                      >  2:1597540/1‑150 (MQ=255)
                                                                                             cACGCGTCGGTGCCATAATTAAGGCACGCGGCTGATTCACCTTGCGATCGGCAATCGCAGGATGAGAGAGAAGATAATGAAACGTTGACGTAAGAAACGCCATCGTTTTCCCGGTACCGGTTTGCGCCTGCCCGGCTACGTCACGACCCg                                                      >  2:76856/1‑150 (MQ=255)
                                                                                              aCGCGTCGGTGCCATAATTAAGGCACGCGGCTGATTCACCTTGCGATCGGCAATCGCAGGATGAGAGAGAAGATAATGAAACGTTGACGTAAGAAACGCCATCGTTTTCCCGGTACCGGTTTGCGCCTGCCCGGCTACGTCACGACCCGc                                                     >  1:1042064/1‑150 (MQ=255)
                                                                                               cgcgTCGGTGCCATAATTAAGGCACGCGGCTGATTCACCTTGCGATCGGCAATCGCAGGATGAGAGAGAAGATAATGAAACGTTGACGTAAGAAACGCCATCGTTTTCCCGGTACCGGTTTGCGCCTGCCCGGCTACGTCACGACCCGcc                                                    >  1:1335179/1‑150 (MQ=255)
                                                                                               cgcgTCGGTGCCATAATTAAGGCACGCGGCTGATTCACCTTGCGATCGGCAATCGCAGGATGAGAGAGAAGATAATGAAACGTTGACGTAAGAAACGCCAGCGTTTTCCCGGTACCGGTTTGCGCCTGCCCGGCTACGTCACGACCCGcc                                                    <  1:867078/150‑1 (MQ=255)
                                                                                               cgcgACGGCGCCATAATTAAGGCACGCGGCTGATTCACCTTGCGATCGGCAATCGCAGGATGAGAGAGAAGATAATGAAACGTTGACGTAAGAAACGCCATCGTTTTCCCGGTACCGGTTTGCGCCTGCCCGGCTACGTCACGACCCGcc                                                    <  1:180678/150‑1 (MQ=255)
                                                                                                agccgcgTGCCTTAATTATGGCAC‑CGGCTGATTCACCTTGCGATCGGCAATCGCAGGATGAGAGAGAAGATAATGAAACGTTGACGTAAGAAACGCCATCGTTTTCCCGGTACCGGTTTGCGCCTGCCCGGCTACGTCACGACCCGCCAg                                                  >  2:1148118/7‑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)
                                                                                                             aaTTAAGGCACGCGGCTGATTCACCTTGCGATCGGCAATCGCAGGATGAGAGAGAAGATAATGAAACGTTGACGTAAGAAACGCCATCGTTTTCCCGGTACCGGTTTGCGCCTGCCCGGCTACGTCACGACCCGCCAGCGTCAGCGGAAg                                      <  2:1621769/150‑1 (MQ=255)
                                                                                                                      aCGCGGCTGATTCCCCTTGCGATCGGCAATCGCAGGATGAGAGAGCAGATCATGAAACGGTGACGTAAGAAACGCCATCGGTTTCCCGGTACCGGTTTGCGCCTGCCCGGCTACGTCCCGACCCGCCAGCGTCAGCGGAAGGGCCAGTGc                             <  1:1129429/150‑1 (MQ=255)
                                                                                                                           gCTGATTCACCTTGCGATCGGCAATCGCAGGATGAGAGAGAAGATAATGAAACGTTGACGTAAGAAACGCCATCGTTTTCCCGGTACCGGTTTGCGCCTGCCCGGCTACGTCACGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGaa                        >  1:1673926/1‑150 (MQ=255)
                                                                                                                             tGATTCACCTTGCGATCGGCAATCGCAGGATGAGAGAGAAGATAATGAAACGTTGACGTAAGAAACGCCATCGTTTTCCCGGTACCGGTTTGCGCCTGCCCGGCTACGTCACGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATg                      >  1:963384/1‑150 (MQ=255)
                                                                                                                                ttCACCTTGCGATCGGCTATCGCAGGATGAGAGAGAAGATAATGAACCTCTGACGTAAGAAACGCCATCGTTTTCCCGGTTCCGGTGAGCGCCGGCCCGGCTACGTCACGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGc                   <  1:1142542/150‑1 (MQ=255)
                                                                                                                                  catcctGCGATTGCCGATCGCAGGATGAGAGAGAAGATAATGAAACGTTGACGTAAGAAACGCCATCGTTTTCCCGGTACCGGTTTGCGCCTCCCCGGCTACGACACGAACCGCCAGCGTCCGCGGAAGGGCCAGAGCCCGAATCGGcca                 >  1:259329/6‑148 (MQ=255)
                                                                                                                                   atcctGCGATTGCCGATCGCAGGATGAGAGAGAAGATAATGAAACGTTGACGTAAGAAACGCCATCGTTTTCCCGGTACCGGTTTGCGCCTGCCCGGCTACGTCACGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTa                >  1:1604554/5‑150 (MQ=255)
                                                                                                                                   atcctGCGATTGCCGATCGCAGGATGAGAGAGAAGATAATGAAACGTTGACGTAAGAAACGCCATCGTTTTCCCGGTACCGGTTTGCGCCTGCCCGGCTACGTCACGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTa                >  1:1697770/5‑150 (MQ=255)
                                                                                                                                    ccTTGCGATCGGCAATCGCAGGATGAGAGAGAAGATAATGAAACGTTGACGTAAGAAACGCCATCGTTTTCCCGGTACCGGTTTGCGCCTGCCCGGCTACGTCACGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTAc               >  2:1651828/1‑150 (MQ=255)
                                                                                                                                    ccTTGCGATCGGCAATCGCAGGATGAGAGAGAAGATAATGAAACGTTGACGTAAGAAACGCCATAGTTTTCCCGGTACCGGTTTGCGCCTGCCCGGCTACGTCACGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTAc               <  1:1128497/150‑1 (MQ=255)
                                                                                                                                      ttGCGATCGGCAATCGCAGGATGAGAGAGAAGATAATGAAACGTTGACGTAAGAAACGCCATCGTTTTCCCGGTACCGGTTTGCGCCTGCCCGGCTACGTCACGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAg             <  2:1132212/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)
                                                                                                                                         cGATCGGCAATCGCAGGATGAGAGAGAAGATAATGAAACGTTGACGTAAGAAACGCCATCGTTTTCCCGGTACCGGTTTGCGCCTGCCCGGCTACGTCACGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTa          >  2:1044992/1‑150 (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)
                                                                                                                                              ggCAATCGCAGGATGAGAGAGAAGATAATGAAACGTTGACGTAAGAAACGCCATCGTTTTCCCGGTACCGGTTTGCGCCTGCCCGGCTACGTCACGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGaaa     >  2:1375315/1‑150 (MQ=255)
                                                                                                                                              ggCAATCGCAGGATGAGAGAGAAGATAATGAAACGTTGACGTAAGAAACGCCATCGTTTTCCCGGTACCGGTTTGCGCCTGCCCGGCTACGTCACGACCCGCCAGCGTCAGCGGAAGCGCCAGTGCCTGAATGGGCGTACAGTTATGaaa     >  1:73269/1‑150 (MQ=255)
                                                                                                                                               gCAATCGCAGGATGAGAGAGAAGCGACTCAAACGATGACGTAAGAAACGCGATCGGTTTCCCGGTACCGGTTTGCGCCTGCCCGGCTACGTCACGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAAc    <  1:290320/150‑1 (MQ=255)
                                                                                                                                               gCAATCGCAGGATGAGAGAGAAGATAATGAAACGTTGACGTAAGAAACGCCATCGTTTTCCCGGTACCGGTTTGCGCCTGCCCGGCTACGTCACGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAAc    >  1:1275077/1‑150 (MQ=255)
                                                                                                                                               gCAATCGCAGGATGAGAGAGAAGATAATGAAACGTTGACGTAAGAAACGCCATCGTTTTCCCGGTACCGGTATGCGCCTGACCGGCTACGTCACGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAAc    >  1:1530370/1‑150 (MQ=255)
                                                                                                                                                 |                                                                                                                                                     
TGTCGTAGCCATCACCACCGTAAGCCAGACCCAGCTTCAGGCCAGTAGCTTCCGCCAGCGGTTCTGCGTCGGCATGGATCTGCACGGCAAGTTCACGCGTCGGTGCCATAATTAAGGCACGCGGCTGATTCACCTTGCGATCGGCAATCGCAGGATGAGAGAGAAGATAATGAAACGTTGACGTAAGAAACGCCATCGTTTTCCCGGTACCGGTTTGCGCCTGCCCGGCTACGTCACGACCCGCCAGCGTCAGCGGAAGGGCCAGTGCCTGAATGGGCGTACAGTTATGAAACCC  >  NC_000913/3965255‑3965549

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