Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,967,531 C→T 8.8% T372M (ACG→ATG)  rho → transcription termination factor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,967,5310CT 8.8% 290.0 / 13.8 91T372M (ACG→ATG) rhotranscription termination factor
Reads supporting (aligned to +/- strand):  ref base C (34/49);  new base T (4/4);  total (38/53)
Fisher's exact test for biased strand distribution p-value = 7.15e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TACCGGTTCTAAAATGGACGAAGTTATCTACGAAGAGTTTAAAGGTACAGGCAACATGGAACTGCACCTCTCTCGTAAGATCGCTGAAAAACGCGTCTTCCCGGCTATCGACTACAACCGTTCTGGTACCCGTAAAGAAGAGCTGCTCACGACTCAGGAAGAACTGCAGAAAATGTGGATCCTGCGCAAAATCATTCACCCGATGGGCGAAATCGATGCAATGGAATTCCTCATTAATAAACTGGCAATGACCAAGACCAATGACGATTTCTTCGAAATGATGAAACGCTCATAAATTT  >  NC_000913/3967382‑3967680
                                                                                                                                                     |                                                                                                                                                     
tACCGGTTCTAAAATGGACGAAGTTATCTACGAAGAGTTTAAAGGTACAGGCAACATGGAACTGCACCTCTCTCGTAAGATCGCTGAAAAACGCGTCTTCCCGGCTATCGACTACAACCGTTCTGGTACCCGTAAAGAAGAGCTGCTCAc                                                                                                                                                       <  2:522106/150‑1 (MQ=255)
 aCCGGTTCTAAAATGGACGAAGTTATCTACGAAGAGTTTAAAGGTACAGGCACCATGGAACTGCACCTCTCTCGTAAGATCGCTGAAAAACGCGTCTTCCCGGCTATCGACTACAACCGTTCTGGTACCCGTAAAGAAGAGCTGCTCACg                                                                                                                                                      >  2:623414/1‑150 (MQ=255)
  ccGGTTCTAAAATGGACGAAGTTATCTACGAAGAGTTTAAAGGTACAGGCAACATGGAACTGCACCTCTCTCGTAAGATCGCTGAAAAACGCGTCTTCCCGGCTATCGACTACAACCGTTCTGGTACCCGTAAAGAAGAGCTGCTCACGa                                                                                                                                                     <  2:1020300/150‑1 (MQ=255)
  ccGGTTCTAAAATGGACGAAGTTATCTACGAAGAGTTTAAAGGTACAGGCAACATGGAACTGCACCTCTCTCGTAAGATCGCTGAAAAACGCGTCTTCCCGGCTATCGACTACAACCGTTCTGGTACCCGTAAAGAAGAGCTGCTCACGa                                                                                                                                                     >  2:308183/1‑150 (MQ=255)
     gTTCTAAAATGGACGAAGTTATCTACGAAGAGTTTAAAGGTACAGGCAACATGGAACTGCACCTCTCTCGTAAGATCGCTGAAAAACGCGTCTTCCCGGCTATCGACTACAACCGTTCTGGTACCCGTAAAGAAGAGCTGCTCACGACTc                                                                                                                                                  >  1:139517/1‑150 (MQ=255)
      ttCTAAAATGGACGAAGTTATCTACGAAGAGTTTAAAGGTACAGGCAACATGGAACTGCACCTCTCTCGTAAGATCGCTGAAAAACGCGTCTTCCCGGCTATCGACTACAACCGTTCTGGTACCCGTAAAGAAGAGCTGCTCACGACTCa                                                                                                                                                 <  2:420752/150‑1 (MQ=255)
       tCTAAAATGGACGAAGTTATCTACGAAGAGTTTAAAGGTACAGGCAACATGGAACTGCACCTCTCTCGTAAGATCGCTGAAAAACGCGTCTTCCCGGCTATCGACTACAACCGTTCTGGTACCCGTAAAGAAGAGCTGCTCACGACTCAg                                                                                                                                                <  2:240861/150‑1 (MQ=255)
              tGGACGAAGTTATCTACGAAGAGTTTAAAGGTACAGGCAACATGGAACTGCACCTCTCTCGTAAGATCGCTGAAAAACGCGTCTTCCCGGCTATCGACTACAACCGTTCTGGTACCCGTAAAGAAGAGCTGCTCATGACTCAGGAAGAAc                                                                                                                                         <  1:83500/150‑1 (MQ=255)
                gACGAAGTTATATACGAAGATTTTAAAGGTACTGGCAACATGGAACTGCACCTCTCTCGTATGATCGCTGAAAAACGCGTCTTCCCGGCTATCGACTACAACCGTTCTGGTACCCGTAAAGAAGAGCTGCTCACGACTCAGGAAGAACTg                                                                                                                                       <  1:485776/150‑1 (MQ=255)
                     aGTTATCTACGAAGAGTTTAAAGGTACAGGCAACATGGAACTGCACCTCTCTCGTAAGATCGCTGAAAAACGCGTCTTCCCGGCTATCGACTACAACCGTTCTGGTACCCGTAAAGAAGAGCTGCTCACGACTCAGGAAGAACTGCAGaa                                                                                                                                  <  2:166903/150‑1 (MQ=255)
                       ttATCTACGAAGAGTTTAAAGGTACAGGCAACATGGAACTGCACTTCTCTCGTAAGATCGCTGAAAAACGCGTCTTCCCGGCTATCGACTACAACCGTTCTGGTACCCGTAAAGAAGAGCTGCTCACGACTCAGGAAGAACTGCAGaaaa                                                                                                                                <  2:894407/150‑1 (MQ=255)
                        tATCTACGAAGAGTTTAACGGTACAGGCAACATGGAACTGCACCTCTCTCGTAAGATCGCTGAAAAACGCGTCTTCCCGGCTATCGACTACAACCGTTCTGGTACCCGTAAAGAAGAGCTGCTCACGACTCAGGAAGAACTGCAGAAAAt                                                                                                                               <  1:433129/150‑1 (MQ=255)
                          tCTACGAAGAGTTTAAAGGTACAGGCAACATGGAACTGCACCTCTCTCGTAAGATCGCTGAAAAACGCGTCTTCCCGGCTATCGACTACAACCGTTCTGGTACCCGTAAAGAAGAGCTGCTCACGACTCAGGAAGAACTGCAGAAAAtgt                                                                                                                             <  1:1370257/150‑1 (MQ=255)
                           cTACGAAGAGTTTAAAGGTACAGGCAACATGGAACTGCACCTCTCTCGTAAGATCGCTGAAAAACGCGTCTTCCCGGCTATCGACTACAACCGTTCTGGTACCCGTAAAGAAGAGCTGCTCACGACTCAGGAAGAACTGCAGAAAAtgtg                                                                                                                            <  1:461813/150‑1 (MQ=255)
                             aCGAAGAGTTTAAAGGTACAGGCAACATGGAACTGCACCTCTCTCGTAAGATCGCTGAAAAACGCGTCTTCCCGGCTATCGACTACAACCGTTCTGGTACCCGTAAAGAAGAGCTGCTCACGACTCAGGAAGAACTGCAGAAAATGTGGa                                                                                                                          >  1:178908/1‑150 (MQ=255)
                               gAAGAGTTTAAAGGTACAGGCAACATGGAACTGCACCTCTCTCGTAAGATCGCTGAAAAACGCGTCTTCCCGGCTATCGACTACAACCGTTCTGGTACCCGTAAAGAAGAGCTGCTCACGACTCAGGAAGAACTGCAGAAAATGTGGATc                                                                                                                        >  2:936229/1‑150 (MQ=255)
                               gAAGAGTTTAAAGGTACAGGCAACATGGAACTGCACCTCTCTCGTAAGATCGCTGAAAAACGCGTCTTCCCGGCTATCGACTACAACCGTTCTGGTACCCGTAAAGAAGAGCTGCTCACGACTCAGGAAGAACTGCAGAAAATGTGGATc                                                                                                                        >  2:936092/1‑150 (MQ=255)
                                aaGAGTTTAAAGGTACAGGCAACATGGAACTGCACCTCTCTCGTAAGATCGCTGAAAAACGCGTCTTCCCGGCTATCGACTACAACCGTTCTGGTACCCGTAAAGTAGAGCTGCTCACGACTCAGGAAGAACTGCAGAAAATGTGGATcc                                                                                                                       >  2:110539/1‑150 (MQ=255)
                                aaGAGTTTAAAGGTACAGGCAACATGGAACTGCACCTCTCTCGTAAGATCGCTGAAAAACGCGTCTTCCCGGCTATCGACTACAACCGTTCTGGTACCCGTAAAGAAGAGCTGCTCACGACTCAGGAAGAACTGCAGAAAATGTGGATcc                                                                                                                       <  1:898316/150‑1 (MQ=255)
                                aaGAGTTTAAAGGTACAGGCAACATGGAACTGCACCTCTCTCGTAAGATCGCTGAAAAACGCGTCTTCCCGGCTATCGACTACAACCGTTCTGGTACCCGTAAAGAAGAGCTGCTCACGACTCAGGAAGAACTGCAGAAAATGTGGATcc                                                                                                                       <  1:109665/150‑1 (MQ=255)
                                aaGAGTTTAAAGGTACAGGCAACATGGAACTGCACCTCTCTCGTAAGATCGCTGAAAAACGCGTCTTCCCGGCTATCGACTACAACCGTTCTGGTACCCGTAAAGAAGAGCTGCTCACGACTCAGGAAGAACTGCAGAAAATGTGGATcc                                                                                                                       <  1:109669/150‑1 (MQ=255)
                                aaGAGTTTAAAGGTACAGGCAACATGGAACTGCACCTCTCTCGTAAGATCGCTGAAAAACGCGTCTTCCCGGCTATCGACTACAACCGTTCTGGTACCCGTAAAGAAGAGCTGCTCACGACTCAGGAAGAACTGCAGAAAATGTGGATcc                                                                                                                       <  1:109741/150‑1 (MQ=255)
                                 agagTTTAAAGGTACAGGCAACATGGAACTGCGCCTCTCTCGTAAGATCGCTGAAAAACGCGTCTTCCCGGCTTTCGACTACAACCGTTATGGTACCCGTAAAGAATTGCTGCTCACGACTCAGGAAGAACTGCAGAAAATGTGGATCCt                                                                                                                      <  1:1309857/150‑1 (MQ=255)
                                 agagTTTAAAGGTACAGGCAACATGGAACTGCACCTCTCTCGTAAGATCGCTGAAAAACGCGTCTTCCCGGCTATCGACTACAACCGTTCTGGTACCCGTAAAGAAGAGCTGCTCACGACTCAGGAAGAACTGCAGAAAATGTGGATCCt                                                                                                                      <  1:609665/150‑1 (MQ=255)
                                 agagTTTAAAGGTACAGGCAACATGGAACTGCACCTCTCTCGTAAGATCGCTGAAAAACGCGTCTTCCCGGCTATCGACTACAACCGTTCTGGTACCCGTAAAGAAGAGCTGCTCACGACTCAGGAAGAACTGCAGAAAATGTGGATCCt                                                                                                                      >  1:615070/1‑150 (MQ=255)
                                  gagTTTAAAGGTACAGGCAACATGGAACTGCACCTCTCTCGTACGATCGCTGAAAAACGCGTCTTCCCGGCTATCGACTACAACCGTTCTGGTACCCGTAAAGAAGAGCTGCTCACGACTCAGGAAGACCTGCAGACAATGTGGATCCTg                                                                                                                     >  1:401114/1‑150 (MQ=255)
                                   agTTTAAAGGTACAGGCAACATGGAACTGCACCTCTCTCGTAAGATCGCTGAAAAACGCGTCTACCCGGCTATCGACTACAACCGTTCTGGTACCCGTAAAGAAGAGCTGCTCACGACTCAGGAAGAACTGCAGAAAATGTGGATCCTgc                                                                                                                    <  1:269559/150‑1 (MQ=255)
                                     tttAAAGGTACAGGCAACATGGAACTGCACCTCTCTCGTAAGATCGCTGAAAAACGCGTCTTCCCGGCTATCGACTACAACCGTTCTGGTACCCGTAAAGACGAGCTGCTCACGACTCAGGAAGAACTGCAGAAAATGTGGATCCTgcgc                                                                                                                  >  2:1353249/1‑150 (MQ=255)
                                     tttAAAGGTACAGGCAACATGGAACTGCACCTCTCTCGTAAGATCGCTGAAAAACGCGTCTTCCCGGCTATCGACTACAACCGTTCTGGTACCCGTAAAGAAGAGCTGCTCACGACTCAGGAAGAACTGCAGAAAATGTGGATCCTgcgc                                                                                                                  >  2:1353244/1‑150 (MQ=255)
                                      ttAAAGGTACAGGCAACATGGAACTGCACCTCTCTCGTAAGATCGCTGAAAAACGCGTCTTCCCGGCTATCGACTACAACCGTTCTGGTACCCGTAAAGAAGAGCTGCTCATGACTCAGGAAGAACTGCAGAAAATGTGGATCCTGCGCa                                                                                                                 >  2:299561/1‑150 (MQ=255)
                                      ttAAAGGTACAGGCAACATGGAACTGCACCTCTCTCGTAAGATCGCTGAAAAACGCGTCTTCCCGGCTATCGACTACAACCGTTCTGGTACCCGTAAAGAAGAGCTGCTCATGACTCAGGAAGAACTGCAGAAAATGTGGATCCTGCGCa                                                                                                                 >  2:1214892/1‑150 (MQ=255)
                                      ttAAAGGTACAGGCAACATGGAACTGCACCTCTCTCGTAAGATCGCTGAAAAACGCGTCTTCCCGGCTATCGACTACAACCGTTCTGGTACCCGTAAGAAGAGCTGCTCATGACTCAGGAAGAACTGCAGAAAATGTGGATCCTGCGCaa                                                                                                                >  2:299568/1‑150 (MQ=255)
                                           ggTACAGGCAACATGGAACTGCACCTCTCTCGTAAGATCGCTGAAAAACGCGTCTTCCCGGCTATCGACTACAACCGTTCTGGTACCCGTAAAGAAGAGCTGCTCACGACTCAGGAAGAACTGCAGAAAATGTGGATCCTGCGCAAAATc                                                                                                            >  1:576329/1‑150 (MQ=255)
                                            gTACAGGCAACATGGAACTGCACCTCTCTCGTAAGATCGCTGAAAAACGCGTCTTCCCGGCTATCGACTACAACCGTTCTGGTACCCGTAAAGAAGAGCTGCTCATGACTCAGGAAGAACTGCAGAAAATGTGGATCCTGCGCAAAATCa                                                                                                           >  1:531491/1‑150 (MQ=255)
                                            gTACAGGCAACATGGAACTGCACCTCTCTCGTAAGATCGCTGAAAAACGCGTCTTCCCGGCTATCGACTACAACCGTTCTGGTACCCGTAAAGAAGAGCTGCTCACGACTCAGGAAGAACTGCAGAAAATGTGGATCCTGCGCAAAATCa                                                                                                           >  2:337774/1‑150 (MQ=255)
                                                aGGCAACATGGAACTGCACCTCTCTCGTAAGATCGCTGAAAAACGCGTCTTCCCGGCTATCGACTACAACCGTTCTGGTACCCGTAAAGAAGAGCTGCTCACGACTCAGGAAGAACTGCAGAAAATGTGGATCCTGCGCAAAATCATTCa                                                                                                       <  2:903821/150‑1 (MQ=255)
                                                 ggCAACATGGAACTGCACCTCTCTCGTAAGATCGCTGAAAAACGCGTCTTCCCGGCTATCGACTACAACCGTTCTGGTACCCGTAAAGAAGAGCTGCTCACGACTCAGGAAGAACTGCAGAAAATGTGGATCCTGCGCAAAATCATTCAc                                                                                                      <  1:190146/150‑1 (MQ=255)
                                                 ggCAACATGGAACTGCACCTCTCTCGTAAGATCGCTGAAAAACGCGTCTTCCCGGCTATCGACTACAACCGTTCTGGTACCCGTAAAGAAGAGCTGCTCACGACTCAGGAAGAACTGCAGAAAATGTGGATCCTGCGCAAAATCATTCAc                                                                                                      <  1:1355291/150‑1 (MQ=255)
                                                  gCAACATGGAACTGCACCTCTCTCGTAAGATCGCTGAAAAACGCGTCTTCCCGGCTATCGACTACAACCGTTCTGGTACCCGTAAAGAAGAGCTGCTCACGACTCAGGAAGAACTGCAGAAAATGTGGATCCTGCGCAAAATCATTCAcc                                                                                                     <  2:504137/150‑1 (MQ=255)
                                                     aCATGGAACTGCACCTCTCTCGTAAGATCGCTGAAAAACGCGTCTTCCCGGCTATCGACTACAACCGTTCTGGTACCCGTAAAGAAGAGCTGCTCACGACTCAGGAAGAACTGCAGAAAATGTGGATCCTGCGCAAAATCATTCACCCGa                                                                                                  >  1:959518/1‑150 (MQ=255)
                                                     aCATGGAACTGCACCTCTCTCGTAAGATCGCTGAAAAACGCGTCTTCCCGGCTATCGACTACAACCGTTCTGGTACCCGTAAAGAAGAGCTGCTCACGACTCAGGAAGAACTGCAGAAAATGTGGATCCTGCGCAAAATCATTCACCCGa                                                                                                  <  1:745241/150‑1 (MQ=255)
                                                       aTGGAACTGCACCTCTCTCGTAAGATCGCTGAAAAACGCGTCTTCCCGGCTATCGACTACAACCGTTCTGGTACCCGTAAAGAAGAGCTGCTCACGACTCAGGAAGAACTGCAGAAAATGTGGATCCTGCGCAAAATCATTCACCCGATg                                                                                                >  1:303480/1‑150 (MQ=255)
                                                       aTGGAACAGCACCTCTCTCGTAAGATCGCTGAAAAACGCGTCTTCCCGGCTATCGACTACAACCGTTCTGGTACCCGTAAAGAAGAGCTGATCACGACTCAGGAAGAACTGCAGAAAATGTGGATCCTGCGCAAAATCATTCACCCGATg                                                                                                >  1:88892/1‑150 (MQ=255)
                                                          gAACTGCACCTCTCTCGTAAGATCGCTGAAAAACGCGTCTTCCCGGCTATCGACTACAACCGTTCTGGTACCCGTAAAGAAGAGCTGCTCACGACTCAGGAAGAACTGCAGAAAATGTGGATCCTGCGCAAAATCATTCACCCGATGGGc                                                                                             <  1:863700/150‑1 (MQ=255)
                                                                  ccTCTCTCGTAAGATCGCTGAAAAACGCGTCGTCCCGGCTATCGACTACAACCGTTCTGGTACCCGTAAAGAAGAGCTGCTCACGACTCAGGAAGAACTGCAGAAAATGTGGATCCTGCGCAAAATCATTCACCCGATGGGCGAAATCGa                                                                                     <  2:474374/150‑1 (MQ=255)
                                                                    tctctcGTAAGATCGCTGAAAAACGCGTCTTCCCGGCTATCGACTACAACCGTTCTGGTACCCGTAAAGAAGAGCTGCTCACGACTCAGGAAGAACTGCAGAAAATGTGGATCCTGCGCAAAATCATTCACCCGATGGGCGAAATCGATg                                                                                   <  1:789690/150‑1 (MQ=255)
                                                                    tctctcGTAAGATCGCTGAAAAACGCGTCTTCCCGGCTATCGACTACAACCGTTCTGGTACCCGTAAAGAAGAGCTGCTCACGACTCAGGAAGAACTGCAGAAAATGTGGATCCTGCGCAAAATCATTCACCCGATGGGCGAAATCGATg                                                                                   <  2:1050103/150‑1 (MQ=255)
                                                                        tcGTAAGATCGCTGAAAAACGCGTCTTCCCGGCTATCGACTACAACCGTTCTGGTACCCGTAAAGAAGAGCTGCTCACGACTCAGGAAGAACTGCAGAAAATGTGGATCCTGCGCAAAATCATTCACCCGATGGGCGAAATCGATGCact                                                                               >  2:1280616/1‑148 (MQ=255)
                                                                          gTAAGATCGCTGAAAAACGCGTCTTCCCGGCTATCGACTACAACCGTTCTGGTACCCGTAAAGAAGAGCTGCTCACGACTCAGGAAGAACTGCAGAAAATGTGGATCCTGCGCAAAATCATTCACCCGATGGGCGAAATCGATGCAATgg                                                                             <  1:894371/150‑1 (MQ=255)
                                                                               aTCGCTGAAAAACGCGTCTTCCCGGCTATCGACTACAACCGTTCTGGTACCCGTAAAGAAGAGCTGCTCACGACTCAGGAAGAACTGCAGAAAATGTGGATCCTGCGCAAAATCATTCACCCGATGGGCGAAAACGATGCCATGGAATTc                                                                        >  1:1309428/1‑150 (MQ=255)
                                                                                    tGAAAAACGCGTCTTCCCGGCTATCGACTACAACCGTTCTGGTACCCGTAAAGAAGAGCTGCTCACGACTCAGGAAGAACTGCAGAAAATGTGGATCCTGCGCAAAATCATTCACCCGATGGGCGAAATCGATGCAATGGACTTCCTCAt                                                                   >  2:635743/1‑150 (MQ=255)
                                                                                    tGAAAAACGCGTCTTCCCGGCTATCGACTACAACCGTTCTGGTACCCGTAAAGAAGAGCTGCTCACGACTCAGGAAGAACTGCAGAAAATGTGGATCCTGCGCAAAATCATTCACCCGATGGGCGAAATCGATGCAATGGAATTCCTCAt                                                                   <  1:590165/150‑1 (MQ=255)
                                                                                          aCGCGTCTTCCCGGCTATCGACTACAACCGTTCTGGTACCCGTAAAGAAGAGCTGCTCACGACTCAGGAAGAACTGCAGAAAATGTGGATCCTGCGCAAAATCATTCACCCGATGGGCGAAATCGATGCAATGGAATTCCTCATtaataa                                                             <  1:456255/150‑1 (MQ=255)
                                                                                            gcgTCTTCCCGGCTATCGACTACAACCGTTCTGGTACCCGTAAAGAAGAGCTGCTCACGACTCAGGAAGAACTGCAGAAAATGTGGATCCTGCGCAAAATCATTCACCCGATGGGCGAAATCGATGCAATGGAATTCCTCATTAATAAAc                                                           <  1:799837/150‑1 (MQ=255)
                                                                                                 ttCCCGGCTATCGACTACAACCGTTCTGGTACCCGTAAAGAAGAGCTGCTCACGACTCAGGAAGAACTGCAGAAAATGTGGATCCTGCGCAAAATCATTCACCCGATGGGCGAAATCGATGCAATGGAATTCCTCATTAATAAACTGGCa                                                      >  2:531601/1‑150 (MQ=255)
                                                                                                  tCCCGGCTATCGACTACAACCGTTCTGGTACCCGTAAAGAAGAGCTGCTCACGACTCAGGAAGAACTGCAGAATATGTGGATCCTGCGCAAAATCATTCACCCGATGGGCGAAATCGATGCAATGGAATTCCTCATTAATAAACTGGCaa                                                     <  1:1060369/150‑1 (MQ=255)
                                                                                                  tCCCGGCTATCGACTACAACCGTTCTGGTACCCGTAAAGAAGAGCTGCTCACGACTCAGGAAGAACTGCAGAATATGTGGATCCTGCGCAAAATCATTCACCCGATGGGCGAAATCGATGCAATGGAATTCCTCATTAATAAACTGGCaa                                                     <  1:1060340/150‑1 (MQ=255)
                                                                                                   cccGGCTATCGACTACAACCGTTCTGGTACCCGTAAAGAAGAGCTGCTCACGACTCAGGAAGAACTGCAGAAAATGTGGATCCTGCGCAAAATCATTCACCCGATGGGCGAAATCGATGCAATGGAATTCCTCATTAATAAACTGGCAAt                                                    <  2:541220/150‑1 (MQ=255)
                                                                                                   cccGGCTATCGACTACAACCGTTCTGGTACCCGTAAAGAAGAGCTGCTCACGACTCAGGAAGAACTGCAGAAAATGTGGATCCTGCGCAAAATCATTCACCCGATGGGCGAAATCGATGCAATGGAATTCCTCATTAATAAACTGGCAAt                                                    <  2:541216/150‑1 (MQ=255)
                                                                                                   cccGGCTATCGACTACAACCGTTCTGGTACCCGTAAAGAAGAGCTGCTCACGACTCAGGAAGAACTGCAGAAAATGTGGATCCTGCGCAAAATCATTCACCCGATGGGCGAAATCGATGCAATGGAATTCCTCATTAATAAACTGGCAAt                                                    <  2:541186/150‑1 (MQ=255)
                                                                                                       gCTATCGACTACAACCGTTCTGGTACCCGTAAAGAAGAGCTGCTCACGACTCAGGAAGAACTGCAGAAAATGTGGATCCTGCGCAAAATCATTCACCCGATGGGCGAAATCGATGCAATGGAATTCCTCATTAATAAACTGGCAATgacc                                                <  1:623414/150‑1 (MQ=255)
                                                                                                         tATCGACTACAACCGTTCTGGTACCCGTAAAGAAGAGCTGCTCACGACTCAGGAAGAACTGCAGAAAATGTGGATCCTGCGCAAAATCATTCACCCGATGGGCGAAATCGATGCAATGGAATTCCTCATTAATCAACTGGCAATgaccaa                                              <  2:43256/150‑1 (MQ=255)
                                                                                                         tATCGACTACAACCGTTCTGGTACCCGTAAAGAAGAGCTGCTCACGACTCAGGAAGAACTGCAGAAAATGTGGATCCTGCGCAAAATCATTCACCCGATGGGCGAAATCGATGCAATGGAATTCCTCATTAATAAACTGGCAATgaccaa                                              <  1:491659/150‑1 (MQ=255)
                                                                                                           tCGACTACAACCGTTCTGGTACCCGTAAAGAAGAGCTGCTCACGACTCAGGAAGAACTGCAGAAAATGTGGATCCTGCGCAAAATCATTCACCCGATGGGCGAAATCGATGCCATGGAATTCCTCATTAATAAACTGGCAATgaccaaga                                            >  1:121634/1‑150 (MQ=255)
                                                                                                           tCGACTACAACCGTTCTGGTACCCGTAAAGAAGAGCTGCTCACGACTCAGGAAGAACTGCAGAAAATGTGGATCCTGCGCAAAATCATTCACCCGATGGGCGAAATCGATGCAATGGAATTCCTCATTAATAAACTGGCAATgaccaaga                                            >  1:121783/1‑150 (MQ=255)
                                                                                                            cGACTACAACCGTTCTGGTACCCGTAAAGAAGAGCTGCTCACGACTCAGGAAGAACTGCAGAAAATGTGGATCCTGCGCAAAATCATTCACCCGATGGGCGAAATCGATGCAATGGAATTCCTCATTAATAAACTGGCAATgaccaagac                                           <  1:785542/150‑1 (MQ=255)
                                                                                                              aCTACAACCGTTCTGGTACCCGTAAAGAAGAGCTGCTCACGACTCAGGAAGAACTGCAGAAAATGTGGATCCTGCGCAAAATCATTCACCCGATGGGCGAAATCGATGCAATGGAATTCCTCATTAATAAACTGGCAATgaccaagacca                                         >  1:928336/1‑150 (MQ=255)
                                                                                                                tACAACCGTTCTGGTACCCGTAAAGAAGAGCTGCTCACGACTCAGGAAGAACTGCAGAAAATGTGGATCCTGCGCAAAATCATTCACCCGATGGGCGAAATCGATGCAATGGAATTCCTCATTAATAAACTGGCAATGACCAAGACCAAt                                       <  1:263073/150‑1 (MQ=255)
                                                                                                                tACAACCGTTCTGGTACCCGTAAAGAAGAGCTGCTCACGACTCAGGAAGAACTGCAGAAAATGTGGATCCTGCGCAAAATCATTCACCCGATGGGCGAAATCGATGCAATGGAATTCCTCATTAATAAACTGGCAATGACCAAGACCAAt                                       <  1:770940/150‑1 (MQ=255)
                                                                                                                         tCTGGTACCCGTAAAGAAGAGCTGCTCATGACTCAGGAAGAACTGCAGAAAATGTGGATCCTGCGCAAAATCATTCACCCGATGGGCGAAATCGATGCTATGGAATTCCTCATTAATAAACTGGCAATGACCAAGACCAATGACGATttc                              <  2:773185/150‑1 (MQ=255)
                                                                                                                            ggTACCCGTAAAGAAGAGCTGCTCACGACTCAGGAAGAACTGCAGAAAATGTGGATCCTGCGCAAAATCATTCACCCGATGGGCGAAATCGATGCAATGGAATTCCTCATTAATAAACTGGCAATGACCAAGACCAATGACGATttcttc                           <  1:173879/150‑1 (MQ=255)
                                                                                                                            ggTACCCGTAAAGAAGAGCTGCTCACGACTCAGGAAGAACTGCAGAAAATGTGGATCCTGCGCAAAATCATTCACCCGATGGGCGAAATCGATGCAATGGAATTCCTCATTAATAAACTGGCAATGACCAAGACCAATGACGATttcttc                           <  1:1289760/150‑1 (MQ=255)
                                                                                                                             gTACCCGTAAAGAAGAGCTGCTCACGACTCAGGAAGAACTGCAGAAAATGTGGATCCTGCGCAAAATCATTCACCCGATGGGCGAAATCGATGCAATGGAATTCCTCATTAATAAACTGGCAATGACCAAGACCAATGACTATTTCTTCg                          >  2:512685/1‑150 (MQ=255)
                                                                                                                               aCCCGTAAAGAAGAGCTGCTCACGACTCAGGAAGAACTGCAGAAAATGTGGATCCTGCGCAAAATCATTCACCCGATGGGCGAAATCGATGCAATGGAATTCCTCATTAATAAACTGGCAATGTCCAAGACCAATGACGATTTCTTCTaa                        >  2:281154/1‑150 (MQ=255)
                                                                                                                                 ccGTAAAGAAGAGCTGCTCACGACTCAGGAAGAACTGCAGAAAATGTGGATCCTGCGCACAATCATTCACCCGATGGGCGAAATCGATGCACTGGAATTCCTCATTAATAAACTGGCAATGACCAAGACCAATGACGATTTCTTCGAAat                      >  2:641478/1‑150 (MQ=255)
                                                                                                                                 ccGTAAAGAAGAGCTGCTCACGACTCAGGAAGAACTGCAGAAAATGTGGATCCTGCGCAAAATCATTCACCCGATGGGCGAAATCGATGCAATGGAATTCCTCATTAATAAACTGGCAATGACCAAGACCAATGACGATTTCTTCGAAat                      >  2:641502/1‑150 (MQ=255)
                                                                                                                                   gTAAAGAAGAGCTGCTCACGACTCAGGAAGAACTGCAGAAAATGTGGATCCTGCGCAAAATCATTCACCCGATGGGCGAAATCGATGCAATGGAATTCCTCATTAATAAACTGGCAATGACCAAGACCAATGACGATTTCTTCGAAatga                    >  1:329933/1‑150 (MQ=255)
                                                                                                                                   gTAAAGAAGAGCTGCTCACGACTCAGGAAGAACTGCAGAAAATGTGGATCCTGCGCAAAATCATTCACCCGATGGGCGAAATCGATGCAATGGAATTCCTCATTAATAAACTGGCAATGACCAAGACCAATGACGATTTCTTCGAAatga                    <  1:308183/150‑1 (MQ=255)
                                                                                                                                      aagaagAGCTGCTCACGACTCAGGAAGAACTGCAGAAAATGTGGATCCTGCGCAAAATCATTCACCCGATGGGCGAAATCGATGCAATGGAATTCCTCATTAATAAACTGGCAATGACCAAGACCAATGACGATTTCTTCGAAAtgatga                 <  1:494891/150‑1 (MQ=255)
                                                                                                                                       agaagaGCTGCTCACGACTCAGGAAGAACTGCAGAAAATGTGGATCCTGCGCAAAATCATTCACCCGATGGGCGAAATCGATGCAATGGAATTCCTCATTAATAAACTGGCAATGACCAAGACCAATGACGATTTCTTCGAAATGATGaa                <  1:94119/150‑1 (MQ=255)
                                                                                                                                        gaagaGCTGCTCACGACTCAGGAAGAACTGCAGAAAATGTGGATCCTGCGCAAAATCATTCACCCGATGGGCGAAATCGATGCAATGGAATTCCTCATTAATAAACTGGCAATGACCAAGACCAATGACGATTTCTTCGAAATGATGaaa               <  2:95041/150‑1 (MQ=255)
                                                                                                                                        gaagaGCTGCTCACGACTCAGGAAGAACTGCAGAAAATGTGGATCCTGCGCAAAATCATTCACCCGATGGGCGAAATCGATGCAATGGAATTCCTCATTAATAAACTGGCAATGACCAAGACCAATGACGATTTCTTCGAAATGATGaaa               <  2:839754/150‑1 (MQ=255)
                                                                                                                                            agCTGCTCATGACTCAGGAAGAACTGCAGAAAATGTGGATCCTGCGCAAAATCATTCACCCGATGGGCGAAATCGATGCAATGGAATTCCTCATTAATAAACTGGCAATGACCAAGACCAATGACGATTTCTTCGAAATGATGAAACGCt           <  1:1367716/150‑1 (MQ=255)
                                                                                                                                            agCTGCTCATGACTCAGGAAGAACTGCAGAAAATGTGGATCCGGCGCAAAATCATTCACCCGATGGGCGAAATCGATGCAATGGAATTCCTCATTAATAAACTGGCAATGACCAAGACCAATGACGATTTCTTCGAAATGATGAAACGCt           <  1:692217/150‑1 (MQ=255)
                                                                                                                                            agCTGCTCACGACTCAGGAAGAACTGCAGAATATGTGGATCCTGCGCAAAATCATTCACCCGATGGGCGAAATCGATGCAATGGAATTCCTCATTAATAAACTGGCAATGACCAAGACCAATGACGATTTCTTCGAAATGATGAAACGCt           <  2:1208389/150‑1 (MQ=255)
                                                                                                                                              ctgctCACGACTCAGGAAGAACTGCAGAAAATGTGGATCCTGCGCAAAATCATTCACCCGATGGGCGAAATCGATGCAATGGAATTCCTCATTAATAAACTGGCAATGACCAAGACCAATGACGATTTCTTCGAAATGATGAAACGCTCa         >  1:243078/1‑150 (MQ=255)
                                                                                                                                               tgctCACGACTCAGGAAGAACTGCAGAAAATGTGGATCCTGCGCAAAATCATTCACCCGATGGGCGAAATCGATGCAATGGAATTCCTCATTAATAAACTGGCAATGACCAAGACCAATGACGATTTCTTCGAAATGATGAAACGCTCAt        >  1:318083/1‑150 (MQ=255)
                                                                                                                                               tgctCACGACTCAGGAAGAACTGCAGAAAATGTGGATCCTGCGCAAAATCATTCACCCGATGGGCGAAATCGATGCAATGGAATTCCTCATTAATAAACTGGCAATGACCAAGACCAATGACGATTTCTTCGAAATGATGAAACGCTCAt        <  1:463294/150‑1 (MQ=255)
                                                                                                                                               tgctCACGACTCAGGAAGAACTGCAGAAAATGTGGATCCTGCGCAAAATCATTAACCAGATGGGCGAAATCGATGCAATGGAATTCCTCATTAATAAACTGGCAATGACCAAGACCAAATACGATTTCTTCGAAATGATGAAACGCTCAt        >  1:1180375/1‑150 (MQ=255)
                                                                                                                                                gctCACGACTCAGGAAGAACTACAGAAAATGTGGATCCTGCGCAAAATCATTCACCCGATGGGCGAAATCGATGCAATGGAATTCCTCATTAATAAACTGGCAATGACCAAGACCAATGACGATTTCTTCGAAATGATGAAACGCTCATa       >  1:191756/1‑150 (MQ=255)
                                                                                                                                                  tCACGACTCAGGAAGAACTGCAGAAAATGTGGATCCTGCGCAAAATCATTCACCCGATGGGCGAAATCGATGCAATGGAATTCCTCATTAATAAACTGGCAATGACCAAGACCAATGACGATTTCTTCGAAATGATGAAACGCTCATaaa     >  2:433531/1‑150 (MQ=255)
                                                                                                                                                  tCACGAATCAGGAAGAACTGCAGAAAATGTGGATCCTCCGCAAAATCATACACAAGATGTGACAAATCGATGCAATGGAATTCCTCAAAAATAAACAGGCAATGACCAAGACCAAAGACGATTTCTACGAAATGATCAAACGCTCAaaaa     >  2:139188/1‑150 (MQ=255)
                                                                                                                                                     cGACTCAGGAAGAACTGCAGAAAATGTGGATCCTGCGCAAAATCATTCACCCGATGGGCGAAATCGATGCAATGGAATTCCTCATTAATAAACTGGCAATGACCAAGACCAATGACGATTTCTTCGAAATGATGAAACGCTCATAAAttt  <  1:5645/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
TACCGGTTCTAAAATGGACGAAGTTATCTACGAAGAGTTTAAAGGTACAGGCAACATGGAACTGCACCTCTCTCGTAAGATCGCTGAAAAACGCGTCTTCCCGGCTATCGACTACAACCGTTCTGGTACCCGTAAAGAAGAGCTGCTCACGACTCAGGAAGAACTGCAGAAAATGTGGATCCTGCGCAAAATCATTCACCCGATGGGCGAAATCGATGCAATGGAATTCCTCATTAATAAACTGGCAATGACCAAGACCAATGACGATTTCTTCGAAATGATGAAACGCTCATAAATTT  >  NC_000913/3967382‑3967680

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