Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,603,724 A→C 7.8% T236P (ACC→CCC)  lsrC → autoinducer 2 import system permease protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,603,7240AC 7.8% 273.4 / 6.7 103T236P (ACC→CCC) lsrCautoinducer 2 import system permease protein
Reads supporting (aligned to +/- strand):  ref base A (50/44);  new base C (2/7);  total (52/51)
Fisher's exact test for biased strand distribution p-value = 9.23e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 103 total aligned reads displayed.

ATGCCACGGGCGATAATTTACAGGGCGCTCGTCAACTGGGCGTTCGTACTGAAGCCATTCGCATTGTGGCATTTTCGTTGAACGGCTGCATGGCGGCACTGGCGGGAATTGTGTTTGCTTCGCAGATTGGTTTTATCCCCAACCAGACCGGTACCGGGCTGGAGATGAAAGCAATTGCAGCCTGCGTGCTGGGCGGCATTAGTTTGCTCGGTGGTTCCGGTGCGATCATTGGTGCGGTACTCGGCGCATGGTTCCTGACGCAGATCGATAGCGTACTGGTGCTGTTGCGCATTCC  >  NC_000913/1603578‑1603872
                                                                                                                                                  |                                                                                                                                                    
aTGCCACGGGCGATAATTTACAGGGCGCTCGTCAACTGGGCGTTCGTACTGAAGCCATTCGCATTGTGGCATTTTCGTTGAACGGCTGCATGGCGGCACTGGCGGGAATTGTGTTTGCTTCGCAGATTGGTTTTATCCCCAACCAGACCg                                                                                                                                                   >  1:21072/1‑150 (MQ=255)
 tGCCACGGGCGATAATTTACAGGGCGCTCGTCAACTGGGCGTTCGTACTGAAGCCATTCGCATTGTGGCATTTTCGTTGAACGGCTGCATGGCGGCACTGGCGGGAATTGTGTTTGCTTCGCAGATTGGTTTTATCCCCAACCAGACCgg                                                                                                                                                  <  1:51399/150‑1 (MQ=255)
 tGCCACGGGCGATAATTTACAGGGCGCTCGTCAACTGGGCGTTCGTACTGAAGCCATTCGCATTGTGGCATTTTCGTTGAACGGCTGCATGGCGGCACTGGCGGGAATTGTGTTTGCTTCGCAGATTGGTTTTATCCCCAACCAGACCgg                                                                                                                                                  >  2:1810709/1‑150 (MQ=255)
  gCCACGGGCGATAATTTACAGGGCGCTCGTCAACTGGGCGTTCGTACTGAAGCCATTCGCATTGTGGCATTTTCGTTGAACGGCTGCATGGCGGCACTGGCGGGAATTGTGTTTGCTTCGCAGATTGGTTTTATCCCCAACCAGACCGGt                                                                                                                                                 <  2:716238/150‑1 (MQ=255)
  gCCACGGGCGATAATTTACAGGGCGCTCGTCAACTGGGCGTTCGTACTGAAGCCATTCGCATTGTGGCATTTTCGTTGAACGGCTGCATGGCGGCACTGGCGGGAATTGTGTTTGCTTCGCAGATTGGTTTTATCCCCAACCAGACCGGt                                                                                                                                                 >  2:217844/1‑150 (MQ=255)
  gCCACGGGCGAGAATTTACAGGGCGCTCGTCAACTGGGCGTTCGTACTGAATCAATTCTCATTGTGGAATTGTCGTTGAAAGGCTGCATGGCGGCACTGTCGGGAATGGTGTTTGCTTCGCAGATTTGTTTTATCCCCAACCAGCCCGGt                                                                                                                                                 >  2:1058752/1‑150 (MQ=255)
        ggCGATAATTTACAGGGCGCTCGTCAACTGGGCGTTCGTACTGAAGCCATTCGCATTGTGGCATTTTCGTTGAACGGCTGCATGGCGGCACTGGCGGGAATTGTGTTTGCTTCGCAGATTGGTTTTATCCCCAACCAGACCGGTACCggg                                                                                                                                           >  1:231832/1‑150 (MQ=255)
         gCGATAATTTACAGGGCGCTCGTCAACTGGGCGTTCGTACTGAAGCCATTCGCATTGTGGCATTTTCGTTGAACGGCTGCATGGCGGCACTGGCGGGAATTGTGTTTGCTTCGCAGATTGGTTTTATCCCCAACCAGACCGGTACCGGGc                                                                                                                                          >  2:866317/1‑150 (MQ=255)
         gCGATAATTTACAGGGCGCTCGTCAACTGGGCGTGCGTACTGAAGCCATTCGCATTGTGGCATTTTCGTGGAACGGCTGCATGGCGGCACTGGCGGGAATTGTGTTTGCTTCGCAGATTGGTTTTATCCCCAACCAGACCGGTACCGGGc                                                                                                                                          <  2:1631923/150‑1 (MQ=255)
          cGATAATTTACAGGGCGCTCGTCAACTGGGCGTTCGTACTGCAGCCATTCGCATTGTGGCATTTTCGTTGAACGGCTGCATGGCGGCACTGGCGGGAATTGTGTTTGCTTCGCAGATTGGTTTTATCCCCAACCAGACCGGTACCGGGCt                                                                                                                                         >  2:1242253/1‑150 (MQ=255)
           gATAATTTACAGGGCGCTCGTCAACTGGGCGTTCGTACTGAAGCCATTCGCATTGTGGCATTTTCGTTGAACGGCTGCATGGCGGCACTGGCGGGAATTGTGTTTGCTTCGCAGATTGGTTTTATCCCCAACCAGACCGGTACCGGGCTg                                                                                                                                        <  1:101878/150‑1 (MQ=255)
             aaaTTTACCGGGCGCTCGTCTACTGGGCGTTCGAACTGAAGCAATTCGCATTGTGGCATTTTCGTTGAACGGCTGCATGGCGGCACTGGCGGGAATTGTGTTTGCTTCGCAGATTGGTTTTATCCCCAACCAGACCGGTACCGGGCTGga                                                                                                                                      <  1:1349720/149‑1 (MQ=255)
              aaTTTCCCGGGCGCTCGTCAACTGGGCGTTCGTACTGAAGCCATTCGCATTGTGGCATTTTCGTTGAACGGCTGCATGGCGGCACTGGCGGGAATTGTGTTTGCATCGCAGATTGGTTATATCCCCAACCAGACCGGTACCGGGCTGgag                                                                                                                                     <  2:1913156/150‑1 (MQ=255)
              aaTTTACAGGGCGCTCGTCAACTGGGCGTTCGTACTGAAGCCATTCGCATTGTGGCATTTTCGTTGAACGGCTGCATGGCGGCACTGGCGGGAATTGTGTTTGCTTCGCAGATTGGTTTTATCCCCAACCAGACCGGTACCGGGGGGgag                                                                                                                                     >  1:888536/1‑150 (MQ=255)
              aaTTTACAGGGCGCTCGTCAACTGGGCGTTCGTACTGAAGCCATTCGCATTGTGGCATTTTCGTTGAACGGCTGCATGGCGGCACTGGCGGGAATTGTGTTTGCTTCGCAGATTGGTTTTATCCCCAACCAGACCGGTACCGGGCTGgag                                                                                                                                     >  2:1340175/1‑150 (MQ=255)
              aaTTTACAGGGCGCTCGTCAACTGGGCGTTCGTACTGAAGCCATTCGCATTGTGGCATTTTCGTTGAACGGCTGCATGGCGGCACTGGCGGGAATTGTGTTTGCTTCGCAGATTGGTTTTATCCCCAACCAGACCGGTACCGGGCTGgag                                                                                                                                     <  2:1852906/150‑1 (MQ=255)
                 ttACAGGGCGCTGCTCAACTGGGCGTTCGTACTGAAGCCATTCGCATTGTGGCATTTTCGTTGAACGGCTGCATGGCGGCACTGGCGGGAATTGTGTTTGCTTCGCAGATTGGTTTTATCCGCATCCAGCCCGGTACCGGTCTggttggg                                                                                                                                  <  1:365123/150‑6 (MQ=255)
                 ttACAGGGCGCTCGTCAACTGGGCGTTCGTACTGAAGCCATTCGCATTGTGGCATTTTCGTTGAACGGCTGCATGGCGGCACTGGCGGGAATTGTGTTTGCTTCGCAGATTGGTTTTATCCCCAACCAGCCCGGTACCGGTCTggttggg                                                                                                                                  <  2:550687/150‑6 (MQ=255)
                 ttACAGGGCGCTCGTCAACTGGGCGTTCGTACTGAAGCCATTCGCATTGTGGCATTTTCGTTGAACGGCTGCATGGCGGCACTGGCGGGAATTGTGTTTGCTTCGCAGATTGGTTTTATCCCCAACCAGCCCGGTACCGGTCTggttggg                                                                                                                                  <  1:1800307/150‑6 (MQ=255)
                 ttACAGGGCGCTCGTCAACTGGGCGTTCGTACTGAAGCCATTCGCATTGTGGCATTTTCGTTGAACGGCTGCATGGCGGCACTGGCGGGAATTGTGTTTGCTTCGCAGATTGGTTTTATCCCCAACCAGACCGGTACCGGGCTGGAGATg                                                                                                                                  >  2:899192/1‑150 (MQ=255)
                 ttACAGGGCGCTCGTCAACTGGGCGTTCGTACTGAAGCCATTCGCATTGTGGCATTTTCGTTGAACGGCTGCATGGCGGCACTGGCGGGAATTGTGTTTGCTTCGCAGATTGGTTTTATCCCCAACCAGACCGGTACCGGGCTGGAGATg                                                                                                                                  <  1:217844/150‑1 (MQ=255)
                 ttACAGGGCGCTCGTCAACTGGGCGTTCGTACTGAAGCCATTCGCATTGTGGCATTATCGTTGAACGGCTGCATGGCGGCACTGGCGGGAATTGTGTTTGCTTCGCAGATTGGTTTTATCCCCAACCAGCCCGGTACCGGTCTggttggg                                                                                                                                  <  2:1339820/150‑6 (MQ=255)
                     acggggcTCGTAAACTGGGCGTTAGTACTGATCCCATTCGCATTGTGGCATTTTCGTTGAACGGCTGCATGGCGGCACTGGCCGGAATTGTGTTTGCATCGCAGATTGGTTTTATCCCGAACCTGACCGGGACCGGGCTGGAGCTGAAAg                                                                                                                              <  2:204297/145‑1 (MQ=255)
                      gggCGCTCGTCAACTGGGCGTTCGTACTGAAGCCATTCGCATTGTGGCATTTTCGTTGAACGGCTGCATGGCGGCACTGGCGGGAATTGTGTTTGCTTCGCAGATTGGTTTTATCCCCAACCAGACCGGTACCGGGCTGGAGATGAAAGc                                                                                                                             >  2:663426/1‑150 (MQ=255)
                      gggCGCTCGTCAACTGGGCGTTCGTACTGAAGCCATTCGCATTGTGGCATTTTCGTTGAACGGCTGCATGGCGGCACTGGCGGGAATTGAGTTTGCTTCGCAGATTGGTTTTATCCCCAACCAGACCGGTACCGGGCTGGAGATGAAAGc                                                                                                                             >  2:1407363/1‑150 (MQ=255)
                         cgcTCGTCAACTGGGCGTTCGTACTGAAGCCATTCGCATTGTGGCATTTTCGTTGAACGGCTGCATGGCGGCACTGGCGGGAATTGTGTTTGCTTCGCAGATTGGTTTTATCCCCAACCAGACCGGTACCGGGCTGGAGATGAAAGCAAt                                                                                                                          >  2:2155929/1‑150 (MQ=255)
                          gcTCGTCAACTGGGCGTTCGTACTGAAGCCATTCGCATTGTGGCATTTTCGTTGAACGGCTGCATGGCGGCACTGGCGGGAATTGTGTTTGCTTCGCAGATTGGTTTTATCCCCAACCAGACCGGTACCGGGCTGGAGATGAAAGCAAtt                                                                                                                         <  1:1733771/150‑1 (MQ=255)
                           cTCGTCAACTGGGCGTTCGTACTGAAGCCATTCGCATTGTGGCATTTTCGTTGAACGGCTGCATGGCGGCACTGGCGGGAATTGTGTTTGCTTCGGAGATTGGTTTTATCCCCAACCAGACCGGTACCGGGCTGGAGATGAAAGCAATTg                                                                                                                        >  2:1981215/1‑150 (MQ=255)
                            tCGTCAACTGGGCGTTCGTACTGAAGCCATTCTCATTGTGGCATTTTCGTTGAACGGCTGCATGGCGGCACTGGCGGGAATTGTGTTTGCTTCGCAGATTGGTTTTATCCCCAACCAGACCGGTACCGGGCTGGAGATGAAAGCAATTGc                                                                                                                       <  1:2207609/150‑1 (MQ=255)
                            gcGTCAACTGGGCGTGAGTACTGAAGCCATACGCATTTTGGCATTTTCGTTTAACGGCTGCATGTCGGCACTGGCGGGATTTGTGTTTGCTTCGCAGATTGGTTTTATCCCCAACCAGACCGGTACCGGGCTGGAGATGAAAGCAATTGc                                                                                                                       <  2:1399736/149‑1 (MQ=255)
                             cGTCAACTGGGCGTTCGTACTGAAGCCATTCGCATTGTGGCATTTTCGTTGAACGGCTGCATGGCGGCACTGGCGGGAATTGTGTTTGCTTCGCAGATTGGTTTTATCCCCAACCAGACCGGTACCGGGCTGGAGATGAAAGCAATTGCa                                                                                                                      <  2:2107458/150‑1 (MQ=255)
                                caccTGGGCGTTCGTACTGAAGCCATTCGCATTGTGGCATTTTCGTTGAACGGCTGCATGGCGGCACTGGCGGGAATTGTGTTTGCTTCGCAGATTGGTTTTATCCCCAACCAGACCGGTACCGGGCTGGAGATGAAAGCAATTGCAGcc                                                                                                                   >  1:1663895/4‑150 (MQ=255)
                                cAACTGGGCGTTCGTACTGAAGCCATTCGCATTGTGGCATTTTCGTTGAACGGCTGCATGGCGGCACTGGCGGGAATTGTGTTTGCTTCGCAGATTGGTTTTATCCCCAACCAGACCGGTACCGGGCTGGAGATGAAAGCAATTGCAGcc                                                                                                                   >  2:1550306/1‑150 (MQ=255)
                                cAACTGGGCGTTCGTACTGAAGCCATTCGCATTGTGGCATTTTCGTTGAACGGCTGCATGGCGGCACTGGCGGGAATTGTGTTTGCTTCGCAGATTGGTTTTATCCCCAACCAGACCGGTACCGGGCTGGAGATGAAAGCAATTGCAGcc                                                                                                                   >  1:1663889/1‑150 (MQ=255)
                                  aCTGGGCGTTCGCACTCAACCCATTCGCATTGTGGCATTTTCGTTGAACGGCTGCATGGCGGCACTGGCTGGAATTGTGTTTGCTTCGCAGATTGGTTTTATCCCCAACCAGACCGGTACCGGGCTGGAGATGAAAGCAATTGCAGACTg                                                                                                                 <  2:1497723/150‑1 (MQ=255)
                                    tGGGCGTTCGTACTGAAGCCATTCGCATTGTGGCATTTTCGTTGAACGGCTGCATGGCGGCACTGGCGGGAATTGTGTTTGCTTCGCAGATTGGTTTTATCCCCAACCAGACCGGTACCGGGCTGGAGATGAAAGCAATTGCAGCCTGCg                                                                                                               <  2:1593391/150‑1 (MQ=255)
                                    tGGGCGTTCGTACTGAAGCCATTCGCATTGTGGCATTTTCGTTGAACGGCTGCATGGCGGCACTGGCGGGAATTGTGTTTGCTTCGCAGATTGGTTTTATCCCCAACCAGACCGGTACCGGGCTGGAGATGAAAGCAATTGCAGCCTGCg                                                                                                               >  1:788067/1‑150 (MQ=255)
                                          ttCGTACTGAAGCCATTCGCATTGTGGCATTTTCGTTGAACGGCTGCATGGCGGCACTGGCGGGAATTGTGTTTGCTTCGCAGATTGGTTTTATCCCCAACCAGACCGGTACCGGGCTGGAGATGAAAGCAATTGCAGCCTGCGTGCTgg                                                                                                         >  2:400766/1‑150 (MQ=255)
                                             gTACTGCAGCCATTCGCATTGTGGCATTTTCGTTGAACGGCTGCATGGCGGCACTGGCGGGAATTGTGTTTGCTTCGCAGATTGGTTTTATCCCCAACCAGACCGGTACCGGGCTGGAGATGAAAGCAATTGCAGCCTGCGTGCTGggcg                                                                                                      <  1:1563871/150‑1 (MQ=255)
                                             gTACTGAAGCCATTCGCATTGTGGCATTTTCGTTGAACGGCTGCATGGCGGCACTGGCGGGAATTGTGTTTGCTTCGCAGATTGGTTTTATCCCCAACCAGACCGGTACCGGGCTGGAGATGAAAGCAATTGCAGCCTGCGTGCTGggcg                                                                                                      >  1:751293/1‑150 (MQ=255)
                                               aCTGAAGCCATTCGCATTGTGGCATTTTCGTTGAACGGCTGCATGGCGGCACTGGCGGGAATTGTGTTTGCTTCGCAGATTGGTTTTATCCCCAACCAGACCGGTACCGGGCTGGAGATGAAAGCAATTGCAGCCTGCGTGCTGggcggc                                                                                                    <  2:698402/150‑1 (MQ=255)
                                                cTGAAGCCATTCGCATTGTGGCATTTTCGTTGAACGGCTGCATGGCGGCACTGGCGGGAATTGTGTTTGCTTCGCAGATTGGTTTTATCCCCAACCAGACCGGTACCGGGCTGGAGATGAAAGCAATTGCAGCCTGCGTGCTGGGCGGCa                                                                                                   >  2:1870548/1‑150 (MQ=255)
                                                cTGAAGCCATTCGCATTGTGGCATTTTCGTTGAACGGCTGCATGGCGGCACTGGCGGGAATTGTGTTTGCTTCGCAGATTGGTTTTATCCCCAACCAGACCGGTACCGGGCTGGAGATGAAAGCAATTGCAGCCTGCGTGCTGGGCGGCa                                                                                                   >  2:743909/1‑150 (MQ=255)
                                                  gAAGCCATTCGCATTGTGGCATTTTCGTTGAACGGCTGCATGGCGGCACTGGCGGGAATTGTGTTTGCTTCGCAGATTGGTTTTATCCCCAACCAGACCGGTACCGGGCTGGAGATGAAAGCAATTGCAGCCTGCGTGCTGGGCGGCAtt                                                                                                 <  1:1594736/150‑1 (MQ=255)
                                                  caaGCCATTCGCATTGTGGCATTTTCGTTGAACGGCTGCATGGCGGCACTGGCGGGAATTGTGTTTGCTTCGCAGATTGGTTTTATCCCCAACCAGCCCGGTCACGGGCAGGAGAGGAAAGCAATTGCAGCCTGCGTGCTGGTCGGCagt                                                                                                 >  1:1349906/2‑148 (MQ=255)
                                                   aaGCCATTCGCATTGTGGCATTTTCGTTGAACGGCTGCATGGCGGCACTGGCGGGAATTGTGTTTGCTTCGCAGATTGGTTTTATCCCCAACCAGACCGGTACCGGGCTGGAGATGAAAGCAATTGCAGCCTGCGTGCTGGGCGGCATTa                                                                                                <  2:1795169/150‑1 (MQ=255)
                                                   aaGCCATTCGCATTGTGGCATTTTCGTTGAACGGCGGCATGGCGGCACTGGCGGGAATTGTGTTTGCTTCGCAGATTGGTTTTATCCCCAACCAGACCGGTACCGGGCTGGAGATGAAAGCAATTGCAGCCTGCGTGCTGGGCGGCATTa                                                                                                <  2:1623508/150‑1 (MQ=255)
                                                    aGCCATTCGCATTGTGGCATTTTCGTTGAACTGCTGCATGGCGGCACTGGCGGGAATTGTGTTTGCTTCGCAGATTGGTTTTATCCCCAACCAGACCGGTACCGGGCTGGAGATGAAAGCAATTGCAGCCTGCGTGCTGGGCGGCATTAg                                                                                               >  1:1122642/1‑150 (MQ=255)
                                                      ccATTCGCATTGTGGCATTTTCGTTGAACGGCTGCATGGCGGCACTGGCGGGAATTGTGTTTGCTTCGCAGATTGGTTTTATCCCCAACCAGACCGGTACCGGGCTGGAGATGAAAGCAATTGCAGCCTGCGTGCTGGGCGGCATTAGtt                                                                                             <  1:899192/150‑1 (MQ=255)
                                                       cATTCGCATTGTGGCATTTTCGTTGTCTGGCTGCATGGCGTCACTGGCGGGAATTGTGTTTGATTCGCAGATTGGTTTTATCCCCAACCAGACCGGTACCGGGCTGGAGATGAAAGCAATTGCAGCCTTCGTGCTGGGCGGCATTAGttt                                                                                            <  1:284223/150‑1 (MQ=255)
                                                         ttCGCATTGTGGCATTTTCGTTGAACGGCTGCATGGCGGCACTGGCGGTAATTGTGTTTGCTTCGCAGATTGGTTTTATCCCCAACCAGACCGGGACCGGGCTGGAGATGAAAGCAATTGCAGCCTGCGTGCTGGGCGGCATTAGTTTGc                                                                                          <  2:412367/150‑1 (MQ=255)
                                                             cATTGTGGCATTTTCGTTGAACGGCTGCCTGGCGGCACTGGCGGGAATTGTGTTTGCTTCGCAGATTGGTTTTATCCCCAACCAGACCGGTACCGGGCTGGAGATGAAAGCAATTGCAGCCTGCGTGCTGGGCGGCATTAGTTTGCTCgg                                                                                      <  1:2132268/150‑1 (MQ=255)
                                                             cATTGTGGCATTTTCGTTGAACGGCTGCATGGCGGCACTGGCGGGAATTGTGTTTGCTTCGCAGATTGGTTTTATCCCCAACCAGACCGGTACCGGGCTGGAGATGAAAGCAATTGCAGCCTGCGTGCTGGGCGGCATTAGTTTGCTCgg                                                                                      <  2:927467/150‑1 (MQ=255)
                                                                  tgGCATTTTCGTTGAACGGCTGCATGGCGGCACTGGCGGGAATTGTGTTTGCTTCGCAGATTGGTATTATCACCAACCAGACAGGTACCGGGCTGGAGATGAAAGCAATTGCAGCCTGCGTGCTGGGCGGCATTAGTTTGCTCGGTGGtt                                                                                 >  1:1828619/1‑150 (MQ=255)
                                                                     cATTTTCGTTGAACGGCTGCATGGCGGCACTGGCGGGAATTGTGTTTGCTTCGCAGATTGGTTTTATCCCCAACCAGACCGGTACCGGGCTGGAGATGAAAGCAATTGCAGCCTGCGTGCTGGGCGGCATTAGTTTGCTCGGTGGTTCCg                                                                              >  2:1396764/1‑150 (MQ=255)
                                                                      aTTTTCGTTGAACGGCTGCATGGCGGCACTGGCGGGAATTGTGTTTGCTTCGCAGATTGGTTTTATCCCCAACCAGACCGGTACCGGGCTGGAGATGAAAGCAATTGCAGCCTGCGTGCTGGGCGGCATTAGTTTGCTCGGTGGTTCCgg                                                                             <  1:2113119/150‑1 (MQ=255)
                                                                      aTTTTCGTTGAACGGCTGCATGGCGGCACTGGCGGGAATTGTGTTTGCTTCGCAGATTGGTTTTATCCCCAACCAGACCGGTACCGGGCTGGAGATGAAAGCAATTGCAGCCTGCGTGCTGGGCGGCATTAGTTTGCTCGGTGGTTCCgg                                                                             <  1:772072/150‑1 (MQ=255)
                                                                        tttCGTTGAACGGCTGCATGGCGGCACTGGCGGGAATTGTGTTTGCTTCGCAGATTGGTTTTATCCCCAACCAGACCGGTACCGGGCTGGAGATGAAAGCAATTGCAGCCTGCGTGCTGGGCGGCATTAGTTTGCTCGGTGGTTCCGGTg                                                                           <  2:662097/150‑1 (MQ=255)
                                                                          tCGTTCCACGGCTGCATGGAGGCCAGGGCGGGATTTGTCTTTGCTGCGCAGATTGGTGTTATCCCCAACCCGACGGGTACCGTGGTGGAGATGAAAGCAATTGCAGCCTGCGTGGTGGGCGGCATTAGTTTGCTCGGTGGTTCCGGGGCg                                                                         <  2:640878/150‑1 (MQ=255)
                                                                            gTTGAACGGCTGCAGGGCGGCACTGGCGGGAATTGTGTTTGCTGCGCAGATTGGTTTTATCCCCAACCAGACCGGTACCGGGCTGGAGATGAAAGCAATTGCAGCCTGCGTGCTGGGCGGCATTAGTTTGCTCGGTGGTTCCGGTGCGAt                                                                       <  1:1271763/150‑1 (MQ=255)
                                                                             ttGAACGGCTGCATGGCGGCACTGGCGGGAATTGTGTTTGCTTCGCAGATTGGTTTTATCCCCAACCAGACCGGTACCGGGCTGGAGATGAAAGCAATTGCAGCCTGCGTGCTGGGCGGCATTAGTTTGCTCGGTGGTTCCGGTGCGATc                                                                      >  1:1777214/1‑150 (MQ=255)
                                                                                aaCGGCTGCATGGCGGCACTGGCGGGAATTGTGTTTGCTTCGCAGATTGGTTTTATCCCCAACCAGACCGGTACCGGGCTGGAGATGAAAGCAATTGCAGCCTGCGTGCTGGGCGGCATTAGTTTGCTCGGTGGTTCCGGTGCGATCAtt                                                                   >  1:994649/1‑150 (MQ=255)
                                                                                    gCTGCATGGCGGCACTGGCGGGAATTGTGTTTGCTTCGCAGATTGGTTTTATCCCCAACCAGACCGGTACCGGGCTGGAGATGAAAGCAATTGCAGCCTGCGTGCTGGGCGGCATTAGTTTGCTCGGTGGTTCCGGTGCGATCATTGGTg                                                               <  2:1145912/150‑1 (MQ=255)
                                                                                     cTGCATGGCGGCACTGGCGGGAATTGTGTTTGCTTCGCAGATTGGTTTTATCCCCAACCAGACCGGTACCGGGCTGGAGATGAAAGCAATTGCAGCCTGCGTGCTGGGCGGCATTAGTTTGCTCGGTGGTTCCGGTGCGATCATTGGTGc                                                              <  2:852161/150‑1 (MQ=255)
                                                                                     cTGCATGGCGGCACTGGCGGGAATTGTGTTTGCTTCGCAGATTGGTTTTATCCCCAACCAGACCGGTACCGGGCTGGAGATGAAAGCAATTGCAGCCTGCGTGCTGGGCGGCATTAGTTTGCTCGGTGGTTCCGGTGCGATCATTGGTGc                                                              <  1:1810709/150‑1 (MQ=255)
                                                                                       gCATGGCGGCATTGGGTGGAATTTTGTTAGCTTCGAAGATTGGTTTTATCCCCAATCAGACCGGTACCGGGCTGGAGATTAAAGCAATTGCTGCCTGCGTGCTGGGCGGCATTAGTTTGCTCGGTGGTTCCGGTGCGATCATTGGTGCgg                                                            <  2:1074731/150‑1 (MQ=255)
                                                                                        cATGGCGGCACTGGCGGGAATTGTGTTTGCTTCGCAGATTGGTTTTATCCCCAACCAGACCGGTACCGGGCTGGAGATGAAAGCAATTGCAGCCTGCGTGCTGGGCGGCATTAGTTTGCTCGGTGGTTCCGGTGCGATCATTGGTGCGGt                                                           <  1:1242253/150‑1 (MQ=255)
                                                                                               gcACTGGCGGGAATTGTGTTTGCTTCGCAGATTGGTTTTATCCCCAACCAGACCGGTACCGGGCTGGAGATGAAAGCAATTGCAGCCTGCGTGCTGGGCGGCATAAGTTTGCTCGGTGGTTCAGGTGCGATCATAGGTGCGGTACTCGgc                                                    >  1:760690/1‑150 (MQ=255)
                                                                                                cACTGGCGGGAATTGTGTTTGCTTCGCAGATTGGTTTTATCCCCAACCAGACCGGTACCGGGCTGGAGATGAAAGCAATTGCAGCCTGCGTGCTGGGCGGCATTAGTTTGCTCGGTGGTTCCGGTGCGATCATTGGTGCGGTACTCGgcg                                                   >  1:667440/1‑150 (MQ=255)
                                                                                                 aCTGGCGGGAATTGTGTTTGCTTCGCAGATTGGTTTTATCCCCAACCAGACCGGTACCGGGCTGGAGATGAAAGCAATTGCAGCCTGCGTGCTGGGCGGCATTAGTTTGCTCGGTGGTTCCGGTGCGATCATTGGTGCGGTACTCGgcgc                                                  <  1:940996/150‑1 (MQ=255)
                                                                                                  cTGGCGGGAATTGTGTTTGCTTCGCAGATTGGTTTTATCCCCAACCAGACCGGTACCGGGCTGGAGATGAAAGCAATTGCAGCCTGCGTGCTGGGCGGCATTAGTTTGCTCGGTGGTTCCGGTGCGATCATTGGTGCGGTACTCGGCGCa                                                 >  1:532784/1‑150 (MQ=255)
                                                                                                      cGGGAATTGTGTTTGCTTCGCAGATTGGTTTTATCCCCAACCAGACCGGTACCGGGCTGGAGATGAAAGCAATTGCAGCCTGCGTGCTGGGCGGCATTAGTTTGCTCGGTGGTTCCGGTGCGATCATTGGTGCGGTACTCGGCGCATGGt                                             >  2:1225965/1‑150 (MQ=255)
                                                                                                         tatgtgtgtgCGCTTCGCAGATAGGTTTTAACCCCAACCAGAACGGTACCGGGCTGGAGATTACAGCAATTGCAGCCGGCGTGCTGGGCGGGATTAGTTTGCTCGGTTGTTCCGGTGCGATCTTTGGTGCGGTACTCGGCGCATGGTTcc                                          <  2:1515963/146‑1 (MQ=255)
                                                                                                         gAATTGTGTTTGCTTCGCAGATTGGTTTTATCCCCAACCAGACCGGTACCGGGCTGGAGATGAAAGCAATTGCAGCCTGCGTGCTGGGCGGCATTAGTTTGCTCGGTGGTTCCGGTGCGATCATTGGTGCGGTACTCGGCGCATGGTTcc                                          <  2:1302848/150‑1 (MQ=255)
                                                                                                              gtgtTTGCTTCGCAGATTGGTTTTATCCCCAACCAGACCGGTACCGGGCTGGAGATGAAAGCAATTGCAGCCTGCGTGCTGGGCGGCATTAGTTTGCTCGGTGGTTCCGGTGCGATCATTGGTGCGGTACTCGGCGCATGGTTCCTGACg                                     >  2:2209425/1‑150 (MQ=255)
                                                                                                                   tGCTTCGCAGATTGGTTTTATCCCCAACCAGACCGGTACCGGGCTGGAGATGAAAGCAATTGCAGCCTGCGTGCTGGGCGGCATTAGTTTGCTCGGTGGTTCCGGTGCGATCATTGGTGCGGTACTCGGCGCATGGTTCCTGACGCAGAt                                >  1:45820/1‑150 (MQ=255)
                                                                                                                     cTTCGCAGATTGGTTTTATCCCCAACCAGACCGGTACCGGGCTGGAGATGAAAGCAATTGCAGCCTGCGTGCTGGGCGGCATTAGTTTGCTCGGTGGTTCCGGTGCGATCATTGGTGCGGTACTCGGCGCATGGTTCCTGACGCAGATCg                              >  1:1593840/1‑150 (MQ=255)
                                                                                                                        cGCAGATTGGTTTTATCCCCAACCAGACCGGTACCGGGCTGGAGATGAAAGCAATTGCAGCCTGCGTGCTGGGCGGCATTAGTTTGCTCGGTGGTTCCGGTGCGATCATTGGTGCGGTACTCGGCGCATGGTTCCTGCCGCAGATCGATa                           >  1:1786314/1‑150 (MQ=255)
                                                                                                                            gATTGGTTTTATCCCCAACCAGACCGGTACCGGGCTGGAGATGAAATCAATTGCAGCCTGCGTGCTGGGCGGCATTAGTTTGCTCGGTGGTTCCGGTGCGATCATTGGTGCGGTACTCGGCGCATGGTTCCTGACGCAGATCGATAGCGt                       >  2:262528/1‑150 (MQ=255)
                                                                                                                            gATTGGTTTTATCCCCAACCAGACCGGTACCGGGCTGGAGATGAAAGCAATTGCAGCCTGCGTGCTGGGCGGCATTAGTTTGCTCGGTGGTTCCGGTGCGATCATTGGTGCGGTACTCGGCGCATGGTTCCTGACGCAGATCGATAGCGt                       <  1:1523554/150‑1 (MQ=255)
                                                                                                                               tGGTTTTATCCCCAACCAGACCGGTACCGGGCTGGAGATGAAAGCAATTGCAGCCTGCGTGCTGGGCGGCATTAGTTTGCTCGGTGGTTCCGGTGCGATCATTGGTGCGGTACTCGGCGCATGGTTCCTGACGCAGATCGATAGCGTACt                    >  2:1848425/1‑150 (MQ=255)
                                                                                                                                ggTTTTATCCCCAACCCGCCCGGTACCGGGCTGGAGATGAACGCAATTGCAGCCTGCGTGCTGGGCGGCATTAGTTTGCTCGGTGGTTCCGGTGCGATCATTGGTGCGGTACTCGGCGCATGGTTCCTGACGCAGATCGATAGCGTACTg                   <  1:177208/150‑1 (MQ=255)
                                                                                                                                ggTTTTATCCCCAACCAGACCGGTACCGGGCTGGAGATGAAAGCAATTGCAGCCTGCGTGCTGGGCGGCATTAGTTTGCTCGGTGGTTCCGGTGCGATCATTGGTGCGGTAGTCGGAGCATGGTTCCTGACGCAGAGCGATAGCGTACTg                   >  2:326361/1‑150 (MQ=255)
                                                                                                                                ggTTTTATCCCCAACCAGACCGGTACCGGGCTGGAGATGAAAGCAATTGCAGCCTGCGTGCTGGGCGGCATTAGTTTGCTCGGTGGTTCCGGTGCGATCATTGGTGCGGTACTCGGCGCATGGTTCCTGACGCAGATCGATGGCGTACTg                   >  2:1310603/1‑150 (MQ=255)
                                                                                                                                 gTTTTATCCCCAACCAGACCGGTACCGGGCTGGAGATGAAAGCAATTGCAGCCTGCGTGCTGGGCGGCATTAGTTTGCTCGGTGGTTCCGGTGCGATCATTGGTGCGGTACTCGGCGCATGGTTCCTGACGCAGATCGATAGCGTACTgg                  >  1:1811492/1‑150 (MQ=255)
                                                                                                                                  ttttATCACCAACCAGACCGGTACCGGGCTGGAGATGAAATCAATTGCAGCCTGCGTGCTGGGCGGCATTAGTTTGCTCGGTGGTTCCGGTGCGATCATTGGTGCGGTACTCGGCGCATGGTTCCTGACGCAGATCGATAGCGTACTGGt                 >  2:69434/1‑150 (MQ=255)
                                                                                                                                   tttATCCCCAACCAGACCGGTACCGGGCTGGAGATGAAAGCAATTGCAGCCTGCGTGCTGGGCGGCATTAGTTTGCTCGGTGGTTCCGGTGCGATCATTGGTGCGGTACTCGGCGCATGGTTCCTGACGCAGATCGATAGCGTACTGGTg                >  2:1664114/1‑150 (MQ=255)
                                                                                                                                    ttATCCCCACCCAGACCGGTACCGGGCTGGAGATGAAAGCAATTGCAGCCTGCGGGCTGGGCGGCATTAGTTTGCTCGGTGGGTCCGGTGCGATCATTGGTGCGGTACTCGGCGCATGGGTCCTGACGCAGATCGATAGCGTACTGGTGc               <  1:1340175/150‑1 (MQ=255)
                                                                                                                                      aTCCCCAACCAGACCGGTACCGGGCTGGAGATGAAAGCAATTGCAGCCTGCGTGCTGGGCGGCATTAGTTTGCTCGGTGGTTCCGGTGCGATCATTGGTGCGGTACTCGGCGCATGGTTCCTGACGCAGATCGATAGCGTACTGGTGCTg             >  1:1036389/1‑150 (MQ=255)
                                                                                                                                       tacccAACCAGACCGGTACCGGGCTGGAGATGAAAGCAATTGTAGCCTGCGTGCTGGGCGGCATTAGTTTGCTCGGTGGTTCCGGTGCGATCATTGGTGCGGTACTCGGCGCATGGTTCCTGACGCAGATCGATTGCGTACTGGTGCTGt            >  2:607650/3‑150 (MQ=255)
                                                                                                                                       tCCCCATCCAGACCTGTACCGGGCTGGAGATGAAAGCAATTGCAGCCTGCGTGCTGGGCGGCATTAGTTTGCTCGGTGGTTCCGGTGCGATCATTGGTGCGGTACTCGGCGCATGGTTCCTGACGCAGATCGATAGCGTACTGGTGCTGt            <  1:295243/150‑1 (MQ=255)
                                                                                                                                        cccctcccAGCCGGGTACGGGGCTGGTGATGAAAGGTATTGCAGCCTGCGTGCTGGGCGGCATTAGTTTGCTCGGTGGTTCCGGTGCGATCATTGGTGCGGTACTCGGCGCATGGTTCCTGACGCAGATCGATAGCGTACTGGTGCTGtt           <  1:400766/144‑1 (MQ=255)
                                                                                                                                         cccAACCAGACCGGTACCGGGCTGGAGATGAAAGCAATTGCAGCCTGCGTGCTGGGCGGCATTAGTTTGCTCGGTGGTTCCGGTGCGATCATTGGTGCGGTACTCGGCGCATGGTTCCTGACGCAGATCGCTTGCGTACTGGTGCTGTTg          >  1:48497/1‑150 (MQ=255)
                                                                                                                                         cccAACCAGACCGGTACCGGGCTGGAGATGAAAGCAATTGCAGCCTGCGTGCTGGGCGGCATTAGTTTGCTCGGTGGTTCCGGTGCGATCATTGGTGCGGTACTCGGCGCATGGTTCCTGACGCAGATCGATAGCGTACTGGTGCTGTTg          >  1:790317/1‑150 (MQ=255)
                                                                                                                                           cacccAGACCGGTACCGGGCGGGAGATCAAAGCAATTGCAGCCTGCGTGCTGGGCGGCATAAGTTTGCGCGGTGGTTCCGGTGCGATCATTGGTGCGGTACTCGGCGCATGGTTCCTGACGCAGATCGATAGCGTACTGGTGCTGTTgcg        <  2:231832/147‑1 (MQ=255)
                                                                                                                                            aaCCAGACCGGTACCGGGCTGGAGATGAAAGCAATTGCAGCCTGCGTGCTGGGCGGCATTAGTTTGCTCGGTGGTTCCGGTGCGATCATTGTTGCGGTACTCTGAGCATGGTTCCTGACGCAGATCGATAGCGTACTGGTGCTGTTgcgc       >  2:1049786/1‑150 (MQ=255)
                                                                                                                                             cccAGACCGGTACCGGGCGGGAGAGGAAAGCAATTGCAGCCTGCGTGCTGGGCGGCATTAGTTTGCTCGGTGGTTCCGGTGCGATCATTGGTGCGGTACTCGGCGCATGGTTCCTGACGCAGATCGATAGCGTACTGGTGCTGTTGCGCa      <  1:1808998/149‑1 (MQ=255)
                                                                                                                                                aGACCGGTACCGGGCTGGAGATGAAAGCAATTGCAGCCTGCGTGCTGGGCGGCATTAGTTTGCTCGGTGGTTCCGGTGCGATCATTGGTGCGGTACTCGGCGCATGGTTCCTGACGCAGATCGATAGCGTACTGGTGCTGTTGCGCATTc   >  1:2044979/1‑150 (MQ=255)
                                                                                                                                                 gACCGGTACCGGGCTGGAGATGAAAGCAATTGCAGCCTGCGTGCTGGGCGGCATTAGTTTGCTCGGTGGTTCCGGTGCGATCATTGGTGCGGTACTCGGCGCATGGTTCCTGACGCAGATCGATAGCGTACTGGTGCTGTTGCGCATTcc  >  1:1942229/1‑150 (MQ=255)
                                                                                                                                                 gACCGGTACCGGGCTGGAGATGAAAGCAATTGCAGCCTGCGTGCTGGGCGGCATTAGTTTGCTCGGTGGTTCCGGTGCGATCATTGGTGCGGTACTCGGCGCATGGTTCCTGACGCAGATCGATAGCGTACTGGTGCTGTTGCGCATTcc  >  1:1631152/1‑150 (MQ=255)
                                                                                                                                                  |                                                                                                                                                    
ATGCCACGGGCGATAATTTACAGGGCGCTCGTCAACTGGGCGTTCGTACTGAAGCCATTCGCATTGTGGCATTTTCGTTGAACGGCTGCATGGCGGCACTGGCGGGAATTGTGTTTGCTTCGCAGATTGGTTTTATCCCCAACCAGACCGGTACCGGGCTGGAGATGAAAGCAATTGCAGCCTGCGTGCTGGGCGGCATTAGTTTGCTCGGTGGTTCCGGTGCGATCATTGGTGCGGTACTCGGCGCATGGTTCCTGACGCAGATCGATAGCGTACTGGTGCTGTTGCGCATTCC  >  NC_000913/1603578‑1603872

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: ‑