Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,169,427 A→C 5.4% F85V (TTT→GTT)  PP_1922 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,169,4270AC 5.4% 279.5 / 4.6 92F85V (TTT→GTT) PP_1922hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (37/50);  new base C (3/2);  total (40/52)
Fisher's exact test for biased strand distribution p-value = 6.49e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.56e-01

GCGCAAGCTAGCTGCTACTGCCTCTGGCTTAGTATTTATAAAAGAAGCGATTATTCTTTGAGCCAAGGGGCGATGTTGGTAGGTCTGAGCGTAATTTAGCAAAGCGGAGTGCTGTGTTTCTTGCATTAGCTTAAGGCGCTCATCGACAAAGTTGTCGACGAATGCTTGGCTGAAACCGTGCTGCCCAAGCTCGATGATTAGGTTGGTGTGCTCGAATTGCTTCGCTCTTTCAGCCGGAGGTAGTGGGCGGGTGAACTCAAAGAACTCAGCTTCGAAGAATGTCGTAACTTG  >  NC_002947/2169278‑2169568
                                                                                                                                                     |                                                                                                                                             
gcgcAAGCTAGCTGCTACTGCCTCTGGCTTAGTATTTATAAAAGAAGCGATTATTCTTTGAGCCAAGGGGCGATGTTGGTAGGTCTGAGCGTAATTTAGCAAAGCGGAGTGCTGTGTTTCTTGCATTAGCTTAAGGCGCTCATCGACaaa                                                                                                                                               <  2:429935/150‑1 (MQ=255)
gcgcAAGCTAGCTGCTACTGCCTCTGGCTTAGTATTTATAAAAGAAGCGATTATTCTTTGAGCCAAGGGGCGATGTTGGTAGGTCTGAGCGTAATTTAGCAAAGCGAAGTGCTGTGTTTCTTGCATTAGCTTCAGGCGCTCATCGACaaa                                                                                                                                               >  2:963521/1‑150 (MQ=255)
      gctagctGCTACTGCCTCTGGCTTAGTATTTATAAAAGAAGCGATTATTCTTTGAGCCAAGGGGCGATGTTGGTAGGTCTGAGCGTAATTTAGCAAAGCGGAGTGCTGTGTTTCTTGCATTAGCTTAAGGCGCTCATCGACAAAGTTGTc                                                                                                                                         >  1:1718911/1‑150 (MQ=255)
           ctgctACTGCCTCTGGCTTAGTATTTATAAAAGAAGCGATTATTCTTTGAGCCAAGGGGCGATGTTGGTAGGTCTGAGCGTAATTTAGCAAAGCGGAGTGCTGTGTTTCTTGCATTAGCTTAAGGCGCTCATCGACAAAGTTGTcgacga                                                                                                                                    >  2:785489/1‑150 (MQ=255)
             gctACTGCCTCTGGCTTAGTATTTATAAAAGAAGCGATTATTCTTTGAGCCAAGAGGCGATGTTGGTAGGTCTGAGCGTAATTTAGCAAAGCGGAGTGCTGTGTTTCTTGCATTAGCTTAAGGCGCTCATCGACAAAGTTGTCGACGAAt                                                                                                                                  <  1:1519890/150‑1 (MQ=255)
               tACTGCCTCTGGCTTAGTATTTATAAAAGAAGCGATTATTCTTTGAGCCAAGGGGCGATGTTGGTAGGTCTGAGCGTAATTTAGCAAAGCGGAGTGCTGTGTTTCTTGCATTAGCTTAAGGCGCTCATCGACAAAg                                                                                                                                              >  2:1810320/1‑136 (MQ=255)
               tACTGCCTCTGGCTTAGTATTTATAAAAGAAGCGATTATTCTTTGAGCCAAGGGGCGATGTTGGTAGGTCTGAGCGTAATTTAGCAAAGCGGAGTGCTGTGTTTCTTGCATTAGCTTAAGGCGCTCATCGACAAAg                                                                                                                                              <  1:1810320/136‑1 (MQ=255)
               tACTGCCTCTGGCTTAGTATTTATAAAAGAAGCGATTATTCTTTGAGCCAAGGGGCGATGTTGGTAGGTCTGAGCGTAATTTAGCAAAGCGGAGTGCTGTGTTTCTTGCATTAGCTTAAGGCGCTCATCGACAAAGTTGTCGACGAATGc                                                                                                                                >  1:1081613/1‑150 (MQ=255)
               tACTGCCTCTGGCTTAGTATTTATAAAAGAAGCGATTATTCTTTGAGCCAAGGGGCGATGTTGGTAGGTCTGAGCGTAATTTAGCAAAGCGGAGTGCTGTGTTTCTTGCATTAGCTTAAGGCGCTCATCGACAAAGTTGTCGACGAATGc                                                                                                                                <  2:2266553/150‑1 (MQ=255)
                aCTGCCTCTGGCTTAGTATTTATAAAAGAAGCGATTATTCTTTGAGCCAAGGGGCGATGTTGGTAGGTCTGAGCGTAATTTAGCAAAGCGGAGTGCTGTGTGTCTTGCATTAGCTTAAGGCGCTCATCGACAAAGTTGTCGACGAATGCt                                                                                                                               <  1:663565/150‑1 (MQ=255)
                  tGCCTCTGGCTTAGTATTTATAAAAGAAGCGATTATTCTTTGAGCCAAGGGGCGATGTTGGTAGGTCTGAGCGTAATTTAGCAAAGCGGAGTGCTGTGTTTCTTGCATTAGCTTAAGGCGCTCATCGACAAAGTTGTCGACGAATGCTTg                                                                                                                             >  1:1616080/1‑150 (MQ=255)
                           cTTAGTATTTATAAAAGAAGCGATTATTCTTTGAGCCAAGGGGCGATGTTGGTAGGTCTGAGCGTAATTTAGCAAAGCGGAGTGCTGTGTTTCTTGCATTAGCTTAAGGCGCTCATCGACAAAGTTGTCGACGAATGCTTGGCTGAAAcc                                                                                                                    <  1:1663561/150‑1 (MQ=255)
                                tATTTATAAAAGAAGCGATTATTCTTTGAGCCAAGGGGCGATGTTGGTAGGTCTGAGCGTAATTTAGCAAAGCGGAGTGCTGTGTTTCTTGCATTAGCTTAAGGCGCTCATCGACAAAGTTGTCGACGAATGCTTGGCTGAAACCGtgct                                                                                                               >  1:291798/1‑150 (MQ=255)
                                 aTTTATAAAAGAAGCGATTATTCTTTGAGCCAAGGGGCGATGTTGGTAGGTCTGAGCGTAATTTAGCAAAGCGGAGTGCTGTGTTTCTTGCATTAGCTTAAGGCGCTCATCGACAAAGTTGTCGACGAATGCTTGGCTGAAACCGtgc                                                                                                                <  2:107281/148‑1 (MQ=255)
                                 aTTTATAAAAGAAGCGATTATTCTTTGAGCCAAGGGGCGATGTTGGTAGGTCTGAGCGTAATTTAGCAAAGCGGAGTGCTGTGTTTCTTGCATTAGCTTAAGGCGCTCATCGACAAAGTTGTCGACGAATGCTTGGCTGAAACCGtgc                                                                                                                >  1:107281/1‑148 (MQ=255)
                                  tttATAAAAGAAGCGATTATTCTTTGAGCCAAGGGGCGATGTTGGTAGGTCTGAGCGTAATTTAGCAAAGCGGAGTGCTGTGTTTCTTGCATTAGCTTAAGGCGCTCATCGACAAAGTTGTCGACGAATGCTTGGCTGAAACCGtgctgc                                                                                                             <  2:1690477/150‑1 (MQ=255)
                                      taAAAGAAGCGATTATTCTTTGAGTCAAGGGGCGATGTTGGTAGGTCTGAGCGTAATTTAGCAAAGCGGAGTGCTGTGTTTCTTGCATTAGCTTAAGGCGCTCATCGACAAAGTTGTCGACGAATGCTTGGCTGAAACCGTGCTGCCCaa                                                                                                         >  1:1843498/1‑150 (MQ=255)
                                        aaaGAAGCGATTATTCTTTGAGCCAAGGGGCGATGTTGGTAGGTCTGAGCGTAATTTAGCAAAGCGGAGTGCTGTGTTTCTTGCATTAGCTTAAGGCGCTCATCGACAAAGTTGTCGACGAATGCTTGGCTGAAACCGTGCTGCCCAAGc                                                                                                       <  2:165577/150‑1 (MQ=255)
                                             agCGATTATTCTTTGAGCCAAGGGGCGATGTTGGTAGGTCTGAGCGTAATTTAGCAAAGCGGAGTGCTGTGTTTCTTGCATTAGCTTAAGGCGCTCATCGACAAAGTTGTCGACGAATGCTTGGCTGAAACCGTGCTGCCCAAGCTCgct                                                                                                  >  1:1591144/1‑148 (MQ=255)
                                               cGATTATTCTTTGAGCCAAGGGGCGATGTTGGTAGGTCTGAGCGTAATTTAGCAAAGCGGAGTGCTGTGTTTCTTGCATTAGCTTAAGGCGCTCATCGACAAAGTTGTCGACGAATGCTTGGCTGAAACCGTGCTGCCCAAGCTCgatga                                                                                                <  1:1632358/150‑1 (MQ=255)
                                                 attattCTTTGAGCCAAGGGGCGATGTTGGTAGGTCTGAGCGTAATTTAGCAAAGCGGAGTGCTGTGTTTCTTGCATTAGCTTAAGGCGCTCATCGACAAAGTTGTCGACGAATGCTTGGCTGaaa                                                                                                                      <  2:1758398/126‑1 (MQ=255)
                                                 attattCTTTGAGCCAAGGGGCGATGTTGGTAGGTCTGAGCGTAATTTAGCAAAGCGGAGTGCTGTGTTTCTTGCATTAGCTTAAGGCGCTCATCGACAAAGTTGTCGACGAATGCTTGGCTGaaa                                                                                                                      >  1:1758398/1‑126 (MQ=255)
                                                 attattCTTTGAGCCAAGGGGCGATGTTGGTAGGTCTGAGCGTAATTTAGCAAAGCGGAGTGCTGTGTTTCTTGCATTAGCTTAAGGCGCTCATCGACAAAGTTGTCGACGAATGCTTGGCTGAAACCGTGCTGCCCAAGCTCGATGAtt                                                                                              >  2:175370/1‑150 (MQ=255)
                                                     ttCTTTGAGCCAAGGGGCGATGTTGGTAGGTCTGAGCGTAATTTAGCAAAGCGGAGTGCTGTGTTTCTTGCATTAGCTTAAGGCGCTCATCGACAAAGTTGTCGACGAATGCTTGGCTGAAACCGTGCTGCCCAAGCTCGATGATTAGGt                                                                                          <  2:2356033/150‑1 (MQ=255)
                                                       cTTTGAGCCAAGGGGCGATGTTGGTAGGTCTGAGCGTAATTTAGCAAAGCGGAGTGCTGTGTTTCTTGCATTAGCTTAAGGCGCTCATCGACAAAGTTGTCGACGAATGCTTGGCTGAAACCGTGCTGCCCAAGCTCGATGATTAGGTTg                                                                                        <  2:1846212/150‑1 (MQ=255)
                                                        tttGAGCCAAGGGGCGATGTTGGTAGGTCTGAGCGTAATTTAGCAAAGCGGAGTGCTGTGTTTCTTGCATTAGCTTAAGGCGCTCATCGACAAAGTTGTCGACGAATGCTTGGCTGAAACCGTGCTGCCCAAGCTCGATGATTAGGTTgg                                                                                       >  1:1344120/1‑150 (MQ=255)
                                                         ttGAGCCAAGGGGCGATGTTGGTAGGTCTGAGCGTAATTTAGCAAAGCGGAGTGCTGTGTTTCTTGCATTAGCTTAAGGCGCTCATCGACAAAGTTGTCGACGAATGCTTGGCTGAAACCGTGCTGCCCAAGCTCGATGATTAGGTTGgt                                                                                      <  2:1709821/150‑1 (MQ=255)
                                                           gAGCCAAGGGGCGATGTTGGTAGGTCTGAGCGTAATTTAGCAAAGCGGAGTGCTGTGTTTCTTGCATTAGCTTAAGGCGCTCATCGACAAAGTTGTCGACGAATGCTTGGCTGAAACCGTGCTGCCCAAGCTCGATGATTAGGTTGgtgt                                                                                    <  2:697335/150‑1 (MQ=255)
                                                             gCCAAGGGGCGATGTTGGTAGGTCTGAGCGTAATTTAGCAAAGCGGAGTGCTGTGTTTCTTGCATTAGCTTAAGGCGCTCATCGACAAAGTTGTCGACGAATGCTTGGCTGAAACCGTGCTGCCCAAGCTCGATGATTAGGTTGGTGTGc                                                                                  <  1:1606421/150‑1 (MQ=255)
                                                              ccAAGGGGCGATGTTGGTAGGTCTGAGCGTAATTTAGCAAAGCGGAGTGCTGTGTTTCTTGCATTAGCTTAAGGCGCTCATCGACAACGTTGTCGACGAATGCTTGGCTGAAACCGTGCTGCCCAAGCTCGATGATTAGGTTGGTGTGCt                                                                                 >  1:142198/1‑150 (MQ=255)
                                                              ccAAGGGGCGATGTTGGTAGGTCTGAGCGTAATTTAGCAAAGCGGAGTGCTGTGTTTCTTGCATTAGCTTAAGGCGCTCATCGACAAAGTTGTCGACGAATGCTTGGCTGAAACCGTGCTGCCCAAGCTCGATGATTAGGTTGGTGTGCt                                                                                 >  1:558612/1‑150 (MQ=255)
                                                              ccAAGGGGCGATGTTGGTAGGTCTGAGCGTAATTTAGCAAAGCGGAGTGCTGTGTTTCTTGCATTAGCTTAAGGCGCTCATCGACAAAGTTGTCGACGAATGCTTGGCTGAAACCGTGCTGCCCAAGCTCGATGATTAGGTTGGTGTGCt                                                                                 >  2:582038/1‑150 (MQ=255)
                                                               cAAGGGGCGATGTTGGTAGGTCTGAGCGTAATTTAGCAAAGCGGAGTGCTGTGTTTCTTGCATTAGCTTAAGGCGCTCATCGACAAAGTTGTCGACGAATGCTTGGCTGAAACCGTGCTGCCCAAGCTCGATGATTAGGTTGGTGTGCTc                                                                                >  2:1639549/1‑150 (MQ=255)
                                                               cAAGGGGCGATGTTGGTAGGTCTGAGCGTAATTTAGCAAAGCGGAGTGCTGTGTTTCTTGCATTAGCTTAAGGCGCTCATCGACAAAGTTGTCGACGAATGCTTGGCTGAAACCGTGCTGCCCAAGCTCGATGATTAGGTTGGTGTGCTc                                                                                <  1:2087231/150‑1 (MQ=255)
                                                                   gggCGATGTTGGTAGGTCTGAGCGTAATTTAGCAAAGCGGAGTGCTGTGTTTCTTGCATTAGCTTAAGGCGCTCATCGACAAAGTTGTCGACGAATGCTTGGCTGAAACCGTGCTGCCCAAGCTCGATGATTAGGTTGGTGTGCTCGAAt                                                                            <  1:1206695/150‑1 (MQ=255)
                                                                      cGATGTTGGTAGGTCTGAGCGTAATTTAGCAAAGCGGAGTGCTGTGTTTCTTGCATTAGCTTAAGGCGCTCATCGACAACGTTGTCGACGAATGCTTGGCTGAAACCGTGCTGCCCAAGCTCGATGATTAGGTTGGTGTGCTCGAATTGc                                                                         <  2:142198/150‑1 (MQ=255)
                                                                       gATGTTGGTAGGTCTGAGCGTAATTTAGCAAAGCGGAGTGCTGTGTTTCTTGCATTAGCTTAAGGCGCTCATCGACAACTTTGTcgatg                                                                                                                                     >  1:232566/1‑87 (MQ=255)
                                                                       gATGTTGGTAGGTCTGAGCGTAATTTAGCAAAGCGGAGTGCTGTGTTTCTTGCATTAGCTTAAGGCGCTCATCGACAACTTTGTcgatg                                                                                                                                     <  2:232566/89‑3 (MQ=255)
                                                                       gATGTTGGTAGGTCTGAGCGTAATTTAGCAAAGCGGAGTGCTGTGTTTCTTGCATTAGCTTAAGGCGCTCATCGACAAAGTTGTCGACGAATGCTTGGCTGAAACCGTGCTGCCCAAGCTCGATGATTa                                                                                             >  1:1208254/1‑129 (MQ=255)
                                                                       gATGTTGGTAGGTCTGAGCGTAATTTAGCAAAGCGGAGTGCTGTGTTTCTTGCATTAGCTTAAGGCGCTCATCGACAAAGTTGTCGACGAATGCTTGGCTGAAACCGTGCTGCCCAAGCTCGATGATTa                                                                                             <  2:1208254/129‑1 (MQ=255)
                                                                         tGTTGGTAGGTCTGAGCGTAATTTAACAAAGCGGAGTGCTGTGTTTCTTGCATTAGCTTAAGGCGCTCATCGACAAAGTTGTCGACGAATGCTTGGCTGAAACCGTGCTGCCCAAGCACGATGATTAGGTTGGTGTGCTCGAATTGCTTc                                                                      <  2:75459/150‑1 (MQ=255)
                                                                          gTTGGTAGGTCTGAGCGTAATTTAGCAAAGCGGAGTGCTGTGTTTCTTGCATTAGCTTAAGGCGCTCATCGACAAAGTTGTCGACGAATGCTTGGCTGAAACCGTGCTGCCCAAGCTCGATGATTAGGTTGGTGTGCTCGAATTGCTTCg                                                                     <  1:175370/150‑1 (MQ=255)
                                                                              gTAGGTCTGAGCGTAATTTAGCAAAGCGGAGTGCTGTGTTTCTTGCATTAGCTTAAGGCGCTCATCGACAAAGTTGTCGACGAATGCTTGGCTGAAACCGTGCTGCCCAAGCTCGATGATTAGGTTGGTGTGCTCGAATTGCTTCGctct                                                                 <  1:785489/150‑1 (MQ=255)
                                                                              gTAGGTCTGAGCGTAATTTAGCAAAGCGGAGTGCTGTGTTTCTTGCATTAGCTTAAGGCGCTCATCGACAAAGTTGTCGACGAATGCTTGGCTGAAACCGTGCTGCCCAAGCTCGATGATTAGGTTGGTGTGCTCGAATTGCTTCGctct                                                                 <  2:291798/150‑1 (MQ=255)
                                                                                aGGTCTGAGCGTAATTTAGCAAAGCGGAGTGCTGTGTTTCTTGCATTAGCTTAAGGCGCTCATCGACAAAGTTGTCGACGAATGCTTGGCTGAAACCGTGCTGCCCAAGCTCGATGATTAGGTTGGTGTGCTCGAATTGCTTCGCTCttt                                                               <  2:1081613/150‑1 (MQ=255)
                                                                                          gTAATTTAGCAAAGCGGAGTGCTGTGTTTCTTGCATTAGCTTAAGGCGCTCATCGACAAAGTTGTCGACGAATGCTTGGCTGAAACCGTGCTGCGCAAGCTCGATGATTATGTTGGTGTGCTCGAATTGCTTCGCTCTTTCAGCCGGAgg                                                     <  2:696249/150‑1 (MQ=255)
                                                                                          gTAATTTAGCAAAGCGGAGTGCTGTGTTTCTTGCATTAGCTTAAGGCGCTCATCGACAAAGTTGTCGACGAATGCTTGGCTGAAACCGTGCTGCCCAAGCTCgat                                                                                                  >  1:69015/1‑105 (MQ=255)
                                                                                          gTAATTTAGCAAAGCGGAGTGCTGTGTTTCTTGCATTAGCTTAAGGCGCTCATCGACAAAGTTGTCGACGAATGCTTGGCTGAAACCGTGCTGCCCAAGCTCgat                                                                                                  <  2:69015/105‑1 (MQ=255)
                                                                                                tAGCAAAGCGGAGTGCTGTGTTTCTTGCATTAGCTTAAGGCGCTCATCGACAAAGTTGTCGACGAATGCTTGGCTGAAACCGTGCTGCCCAAGCTCGATGATTAGGTTGGTGTGCTCGAATTGCTTCGCTCTTTCAGCCGGAGGTAGTgg                                               <  1:1639549/150‑1 (MQ=255)
                                                                                                   cAAAGCGGAGTGCTGTGTTTCTTGCATTAGCTTAAGGCGCTCATCGACAAAGTTGTCGACGAATGCTTGGCTGAAACCGTGCTGCCCAAGCTCGATGATTAGGTTGGTGTGCTCGAATTGCTTCGCTCTTTCAGCCGGAGGTAGTGGGCg                                            <  1:1256082/150‑1 (MQ=255)
                                                                                                   cAAAGCGGAGTGCTGTGTTTCTTGCATTAGCTTAAGGCGCTCATCGACAAAGTTGTCGACGAATGCTTGGCTGAAACCGTGCTGCCCAAGCTCGATGATTAGGTTGGTGTGCTCGAATTGCTTCGCTCTTTCAGCCGGAGGTAGTGGGCg                                            >  2:702580/1‑150 (MQ=255)
                                                                                                   cAAAGCGGAGTGCTGTGTTTCTTGCATTAGCTTAAGGCGCTCATCGACAAAGTTGCCGACGAATGCTTGGCTGAAACCGTGCTGCCCAAGCTCGATGATTAGGTTGGTGTGCTCGAATTGCTTCGc                                                                    <  1:1490658/126‑1 (MQ=255)
                                                                                                   cAAAGCGGAGTGCTGTGTTTCTTGCATTAGCTTAAGGCGCTCATCGACAAAGTTGCCGACGAATGCTTGGCTGAAACCGTGCTGCCCAAGCTCGATGATTAGGTTGGTGTGCTCGAATTGCTTCGc                                                                    >  2:1490658/1‑126 (MQ=255)
                                                                                                    aaaGCGGAGTGCTGTGTTTCTTGCATTAGCTTAAGGCGCTCATCGACAAAGTTGTCGACGAATGCTTGGCTGAAACCGTGCTGCCCAAGCTCGATGATTAGGTTGGTGTGCTCGAATTGCTTCGctc                                                                  >  1:1336898/1‑127 (MQ=255)
                                                                                                    aaaGCGGAGTGCTGTGTTTCTTGCATTAGCTTAAGGCGCTCATCGACAAAGTTGTCGACGAATGCTTGGCTGAAACCGTGCTGCCCAAGCTCGATGATTAGGTTGGTGTGCTCGAATTGCTTCGctc                                                                  <  2:1336898/127‑1 (MQ=255)
                                                                                                              gCTGTGTTTCTTGCATTAGCTTAAGGCGCTCATCGACAAAGTTGTCGACGAATGCTTGGCTGAAACCGTGCTGCCCAAGCTCGATGATTAGGTTGGTGTGCTCGAATTGCTTCGCTCTTTCAGCCGGAGGTAGTGGGCGGGTGAACTCaa                                 >  2:2029979/1‑150 (MQ=255)
                                                                                                               cTGTGTTTCTTGCATTAGCTTAAGGCGCTCATCGACAAAGTTGTCGACGAATGCTTGGCTGAAACCGTGCTGCCCAAGCTCGATGATTAGGTTGGTGTGCTCGAATTGCTTCGCTCTTTCAGCCGGAg                                                      <  1:1447461/128‑1 (MQ=255)
                                                                                                               cTGTGTTTCTTGCATTAGCTTAAGGCGCTCATCGACAAAGTTGTCGACGAATGCTTGGCTGAAACCGTGCTGCCCAAGCTCGATGATTAGGTTGGTGTGCTCGAATTGCTTCGCTCTTTCAGCCGGAg                                                      >  2:1447176/1‑128 (MQ=255)
                                                                                                               cTGTGTTTCTTGCATTAGCTTAAGGCGCTCATCGACAAAGTTGTCGACGAATGCTTGGCTGAAACCGTGCTGCCCAAGCTCGATGATTAGGTTGGTGTGCTCGAATTGCTTCGCTCTTTCAGCCGGAg                                                      >  2:1447461/1‑128 (MQ=255)
                                                                                                               cTGTGTTTCTTGCATTAGCTTAAGGCGCTCATCGACAAAGTTGTCGACGAATGCTTGGCTGAAACCGTGCTGCCCAAGCTCGATGATTAGGTTGGTGTGCTCGAATTGCTTCGCTCTTTCAGCCGGAg                                                      <  1:1447176/128‑1 (MQ=255)
                                                                                                               cTGTGTTTCTTGCATTAGCTTAAGGCGCTCATCGACAAAGTTGTCGACGAATGCTTGGCTGAAACCGTGCTGCCCAAGCTCGATGATTAGGTTGGTGTGCTCGAATTGCTTCGCTCTTTCAGCCGGAGGTAGTGGGCGGGTGATCTCaat                                <  2:1843498/150‑2 (MQ=255)
                                                                                                                tgtgTTTCTTGCATTAGCTTAAGGCGCTCATCGACAAAGTTGTCGACGAATGCTTGGCTGAAACCGTGCTGCCCAAGCTCGATGAtt                                                                                              <  1:1505567/87‑1 (MQ=255)
                                                                                                                tgtgTTTCTTGCATTAGCTTAAGGCGCTCATCGACAAAGTTGTCGACGAATGCTTGGCTGAAACCGTGCTGCCCAAGCTCGATGAtt                                                                                              >  2:1505567/1‑87 (MQ=255)
                                                                                                                tgtgTTTCTTGCATTAGCTTAAGGCGCTCATCGACAAAGTTGTCGACGAATGCTTGGCTGAAACCGTGCTGCCCAAGCTCGATGATTAGGTTGGTGTGCTCGAATTGCTTCGCTCTTTCAGCCGGAGGTAGTGGGCGGGTGAACTCAAAg                               <  1:1121547/150‑1 (MQ=255)
                                                                                                                 gtgtTTCTTGCATTAGCTTAAGGCGCTCATCGACAAAGTTGTCGACGAATGCTTGGCTGAAACCGTGCTGCCCAAGCTCGATGATTAGGTTGGTGTGCTCGAATTGCTTCGCTCTTTCAGCCGGAGGTAGTGGGCGGGTGAACTCAAAGa                              >  2:2393992/1‑150 (MQ=255)
                                                                                                                   gtTTCTTGCATTAGCTTAAGGCGCTCATCGACAAAGTTGTCGACGAATGCTTGGCTGAAACCGTGCTGCCCAAGCTCGATGATTAgg                                                                                           >  1:1514839/1‑87 (MQ=255)
                                                                                                                   gtTTCTTGCATTAGCTTAAGGCGCTCATCGACAAAGTTGTCGACGAATGCTTGGCTGAAACCGTGCTGCCCAAGCTCGATGATTAgg                                                                                           <  2:1514839/87‑1 (MQ=255)
                                                                                                                   gtTTCTTGCATTAGCTTAAGGCGCTCATCGACAAAGTTGTCGACGAATGCTTGGCTGAAACCGTGCTGCCCAAGCTCGATGATTAGGTTGGTGTGCTCGAATTGCTTCGCTCTTTCAGCCGGAGGTAGTGGGCGGGTGAACTCAAAGAAc                            <  2:558612/150‑1 (MQ=255)
                                                                                                                    tttCTTGCATTAGCTTAAGGCGCTCATCGACAAAGTTGTCGACGAATGCTTGGCTGAAACCGTGCTGCCCAAGCTc                                                                                                     >  2:477900/1‑76 (MQ=255)
                                                                                                                    tttCTTGCATTAGCTTAAGGCGCTCATCGACAAAGTTGTCGACGAATGCTTGGCTGAAACCGTGCTGCCCAAGCTc                                                                                                     <  1:477900/76‑1 (MQ=255)
                                                                                                                    tttCTTGCATTAGCTTAAGGCGCTCATCGACAAAGTTGTCGACGAATGCTTGGCTGAAACCGTGCTGCCCAAGCTCGATGATTAGGTTGGTGTGCTCGAAt                                                                            <  1:1116653/101‑1 (MQ=255)
                                                                                                                    tttCTTGCATTAGCTTAAGGCGCTCATCGACAAAGTTGTCGACGAATGCTTGGCTGAAACCGTGCTGCCCAAGCTCGATGATTAGGTTGGTGTGCTCGAAt                                                                            >  2:1116653/1‑101 (MQ=255)
                                                                                                                     ttgttGCCTTAGCTTAAGGCGCTCAGCGACAAAGTTGTCGACGAATGCTTGGCTGAAACCGTGCTGCCCAAGCTCGATGATTAGGTTGGTGTGCTCGAATTGCTTCGCTCTTTCAGCCGGAGGTAGTGGGCGGGTGAACTCAAAGAACTc                          <  2:1344120/147‑1 (MQ=255)
                                                                                                                        ttGCATTAGCTTAAGGCGCTCATCGACAAAGTTGTCGACGAATGCTTGGCTGAAACCGTGCTGCCCAAGCTCGATGATTAGGTTGGTGTGCTCGAATTGCTTCGCTCTTTCAGCCGGAGGTAGTGGGCGGGTGAACTCAAAGAACTCAGc                       >  2:2039606/1‑150 (MQ=255)
                                                                                                                          gCATTAGCTTAAGGCGCTCATCGACAAAGTTGTCGACGAATGCTTGGCTGAAACCGTGCTGCCCAAGCTCGATGATTAGGTTGGTGTGCTCGAATTGCTTCGCTCTTTCAGCCGGAGGTAGTGGGCGGGTGAACTCAAAGAACTCAGCtt                     <  1:2039606/150‑1 (MQ=255)
                                                                                                                           cATTAGCTTAAGGCGCTCATCGACAAAGTTGTCGACGAATGCTTGGCTGaaa                                                                                                                      >  1:637763/1‑52 (MQ=255)
                                                                                                                           cATTAGCTTAAGGCGCTCATCGACAAAGTTGTCGACGAATGCTTGGCTGaaa                                                                                                                      <  2:637763/52‑1 (MQ=255)
                                                                                                                            aTTAGCTTAAGGCGCTCATCGACAAAGTTGTCGACGAATGCTTGGCTGAAACCGTGCTGCCCAAGc                                                                                                       <  2:1466752/66‑1 (MQ=255)
                                                                                                                            aTTAGCTTAAGGCGCTCATCGACAAAGTTGTCGACGAATGCTTGGCTGAAACCGTGCTGCCCAAGc                                                                                                       >  1:1466752/1‑66 (MQ=255)
                                                                                                                            aTTAGCTTAAGGCGCTCATCGACAAAGTTGTCGACGAATGCTTGGCTGAAACCGTGCTGCCCAAGCTCGATGATTAGGTTGGTGTGCTCGAATTGCTTCGCTCTTTCAGCCGGAGGTAGTggg                                              >  1:1383861/1‑123 (MQ=255)
                                                                                                                            aTTAGCTTAAGGCGCTCATCGACAAAGTTGTCGACGAATGCTTGGCTGAAACCGTGCTGCCCAAGCTCGATGATTAGGTTGGTGTGCTCGAATTGCTTCGCTCTTTCAGCCGGAGGTAGTggg                                              <  2:1383861/123‑1 (MQ=255)
                                                                                                                            aTTAGCTTAAGGCGCTCATCGACAAAGTTGTCGACGAATGCTTGGCTGAAACCGTGCTGCCCAAGCTCGATGATTAGGTTGGTGTGCTCGAATTGCTTCGCTCTTTCAGCCGGAGGTAGTGGGCGGGTGAAc                                     >  2:789152/1‑132 (MQ=255)
                                                                                                                            aTTAGCTTAAGGCGCTCATCGACAAAGTTGTCGACGAATGCTTGGCTGAAACCGTGCTGCCCAAGCTCGATGATTAGGTTGGTGTGCTCGAATTGCTTCGCTCTTTCAGCCGGAGGTAGTGGGCGGGTGAAc                                     <  1:789152/132‑1 (MQ=255)
                                                                                                                              tAGCTTAAGGCGCTCATCGACAAAGTTGTCGACGAATGCTTGGCTGAAACCGTGCTGCCCAAGCTCGATGATTAGGTTGGTGTGCTCGAATTGCTTCGCTCTTTCAGCCGGAGGTAGTGGGCGGGTGAACTCAAAGAACTCAGCTTCgaa                 <  1:582038/150‑1 (MQ=255)
                                                                                                                                gCTTAAGGCGCTCATCGACAAAGTTGTCGACGAATGCTTGGCTGAAACCGTGCTGCCCAAGCTCGATGATTAGGt                                                                                          >  1:512109/1‑75 (MQ=255)
                                                                                                                                gCTTAAGGCGCTCATCGACAAAGTTGTCGACGAATGCTTGGCTGAAACCGTGCTGCCCAAGCTCGATGATTAGGt                                                                                          <  2:512109/75‑1 (MQ=255)
                                                                                                                                 cTTAAGGCGCTCATCGACAAAGTTGTCGACGAATGCTTGGCTGAAACCGTGCTGCCCAAGCTCGATGATTAGGTTGGTGTGCTCGAATTGCTTCGCTCTTTCAGCCGGAGGTAGTGGGCGGGTGAACTCAAAGAACTCAGCTTCgaagaa              <  1:702580/150‑1 (MQ=255)
                                                                                                                                     cggCTCTCATCGACAAAGTTGTCGACGAATGCTTGGCTGAAACCGTGCTGCCCAAGCTCGATGATTAGGTTGGTGTGCTCGAATTGCTTCGCTCTTTCAGCCGGAGGTAGTGGGCGGGTGAACTCAAAGAACTCAGCTTCGAAGAATGTc          <  1:963521/149‑1 (MQ=255)
                                                                                                                                     aGGCGCTCATCGACAAAGTTGTCGACGAATGCTTGGCTGAAACCGTGCTGCCCAAGCTCGATGATTAGGTTGgt                                                                                      <  1:97802/74‑1 (MQ=255)
                                                                                                                                     aGGCGCTCATCGACAAAGTTGTCGACGAATGCTTGGCTGAAACCGTGCTGCCCAAGCTCGATGATTAGGTTGgt                                                                                      >  2:97802/1‑74 (MQ=255)
                                                                                                                                         gcTCATCGACAAAGTTGTCGACGAATGCTTGGCTGAAACCGTGCTGCCCAAGCTCGATGATTAGGTTGGTGTGCTCGAATTGCTTCGCTCTTTCAGCCGGAGGTAGTGGGCGGGTGAACTCAAAGAACTCAGCTTCGAAGAATGTCGTaa      >  2:503017/1‑150 (MQ=255)
                                                                                                                                             gtCGACAACTTTGTCGATGAATGCTTGGCTGAAACCGTGCTGCCCAAGCTCGATGATTAGGTTGGTGTGCTCGAATTGCTTCGCTCTTTCAGCCGGAGGTAGTGGGCGGGTGAACTCAAAGAACTCAGCTTCGAAGAATGTCGTAACTTg  >  1:1465724/2‑150 (MQ=255)
                                                                                                                                             aTCGACAAAGTTGTCGACGAATGCTTGGCTGAAACCGTGCTGCCCAAGCTCGATGATTAg                                                                                            >  2:902448/1‑60 (MQ=255)
                                                                                                                                             aTCGACAAAGTTGTCGACGAATGCTTGGCTGAAACCGTGCTGCCCAAGCTCGATGATTAg                                                                                            <  1:902448/60‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                             
GCGCAAGCTAGCTGCTACTGCCTCTGGCTTAGTATTTATAAAAGAAGCGATTATTCTTTGAGCCAAGGGGCGATGTTGGTAGGTCTGAGCGTAATTTAGCAAAGCGGAGTGCTGTGTTTCTTGCATTAGCTTAAGGCGCTCATCGACAAAGTTGTCGACGAATGCTTGGCTGAAACCGTGCTGCCCAAGCTCGATGATTAGGTTGGTGTGCTCGAATTGCTTCGCTCTTTCAGCCGGAGGTAGTGGGCGGGTGAACTCAAAGAACTCAGCTTCGAAGAATGTCGTAACTTG  >  NC_002947/2169278‑2169568

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