Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,277,228 A→T 7.3% D52E (GAT→GAA soxS ← DNA‑binding transcriptional dual regulator SoxS

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,277,2280AT 7.3% 611.6 / 21.2 194D52E (GAT→GAAsoxSDNA‑binding transcriptional dual regulator SoxS
Reads supporting (aligned to +/- strand):  ref base A (83/96);  new base T (6/9);  total (89/105)
Fisher's exact test for biased strand distribution p-value = 7.89e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 196 total aligned reads displayed.

GCTGGGAGTGCGATCAAACTGCCGACGGAAAACGCGGGAGAAGGTCTGCTGCGAGACATAACCCAGGTCCATTGCGATATCAAAAATCGGACGCTCGGTGGTGCGCAACTCAACGGCGGCCAGTAACAGGCGGCGTTGGCGAATGTAATCGCCAAGCGTCTGATGCGTCACCGTGCGGAACATTCGTTGCAAGTACCACTTTGAATAGCCTGATTTTTTTGCGACTACATCAATGTTAAGCGGCTGGTCAATATGCTCGTCAATCCATGCGATAAGATCCTGAATAATTTTCTGATG  >  NC_000913/4277081‑4277377
                                                                                                                                                   |                                                                                                                                                     
gCTGGGAGTGCGATCAAACTGCCGACGGAAAACGCGGGAGAAGGTCTGCTGCGAGACATAACCCAGGTCCATTGCGATATCAAAAATCGGACGCTCGGTGGTGCGCAACTCAACGGCGGCCAGTAACAGGCGGCGTTGGCGAATGTAATc                                                                                                                                                     >  2:91380/1‑150 (MQ=255)
 cTGGGAGTGCGATCAAACTGCCGACAGAAAACGCGGGAGAAGGTCTGCTGCGAGACATAACCCAGGTCCATTGCGATATCAAAAATCGGACGCTCGGTGGTGCGCAACTCAACGGCGGCCAGTAACAGGCGGCGTTGGCGAATGTAATCg                                                                                                                                                    <  2:1783785/150‑1 (MQ=255)
 cTGGGAGTGCCATCAAACTGCCGACGGAAAACGCGGTAGAAGGTCTGCTGCCAGACATAACCCAGGTCCATAGCGATATCAATAATCGGACGCTCGGTGGTGCGCAACTCAACGGCGGCCAGTAACAGGCGGCGTTGGCGAATGTAATCg                                                                                                                                                    >  2:1853942/1‑150 (MQ=255)
   gggAGTGCGATCAAACTGCCGACGGAAAACGCGGGAGAAGGTCTGCTGCGAGACATAACCCAGGTCCATTGCGATATCAAAAATCGGACGCTCGGTGGTGCGCAACTCAACGGCGGCCAGTAACAGGCGGCGTTGGCGAATGTAATCGcc                                                                                                                                                  >  2:612586/1‑150 (MQ=255)
    ggAGTGCGATCAAACTGCCGACGGAAAACGCGGGAGAAGGTCTGCTGCGAGACATAACCCAGGTCCATTGCGATATCAAAAATCGGACGCTCGGTGGTGCGCAACTCAACGGCGGCCAGTAACAGGCGGCGTTGGCGAATGTAATCGCCa                                                                                                                                                 <  2:1676254/150‑1 (MQ=255)
    ggAGTGCGATCAAACTGCCGACGGAAAACGCGGGAGAAGGTCTGCTGCGAGACATAACCCAGGTCCATTGCGATATCAAAAATCGGACGCTCGGTGGTGCGCAACTCAACGGCGGCCAGTAACAGGCGGCGTTGGCGAATGTAATCGCCa                                                                                                                                                 <  1:1998065/150‑1 (MQ=255)
    ggAGTGCGATCAAACTGCCGACGGAAAACGCGGGAGAAGGTCTGCTGCGAGACATAACCCAGGTCCATTGCGATATCAAAAATCGGACGCTCGGTGGTGCGCAACTCAACGGCGGCCAGTAACAGGCGGCGTTGGCGAATGTAATCGCCa                                                                                                                                                 <  2:1515473/150‑1 (MQ=255)
      aGTGCGATCAAACTGCCGACGGAAAACGCGGGAGAAGGTCTGCTGCGAGACATAACCCAGGTCCATTGCGATATCAAAAATCGGACGCTCGGTGGTGCGCAACTCAACGGCGGCCAGTAACAGGCGGCGTTGGCGATTACATTCGCCaac                                                                                                                                               <  1:15842/150‑2 (MQ=255)
      aGTGCGATCAAACTGCCGACGGAAAACGCGGGAGAAGGTCTGCTGCGAGACATAACCCAGGTCCATTGCGATATCAAAAATCGGACGCTCGGTGGTGCGCAACTCAACGGCGGCCAGTAACAGGCGGCGTTGGCGATTACATTCGCCaac                                                                                                                                               <  2:1763726/150‑2 (MQ=255)
      aGTGCGATCAAACTGCCGACGGAAAACGCGGGAGAAGGTCGGCTGCGAGACATAACCCAGGTCCATTGCGATATCAAAAATCGGACGCTCGGGGGTGCGCAACTCAACGGCGGCCAGTAACAGGCGGCGTTGGCGAATGTAATCGCCAAg                                                                                                                                               <  2:2380303/150‑1 (MQ=255)
       gTGCGATCAAACTGCCGACGGAAAACGCGGGAGAAGGTCTGCTGCGAGACATAACCCAGGTCCATTGCGATATCAAAAATCGGACGCTCGGTGGTGCGCAACTCAACGGCGGCCAGTAACAGGCGGCGTTGGCGAATGTAATCGCCAAGc                                                                                                                                              >  2:2124088/1‑150 (MQ=255)
        tGCGATCAAACTGCCGACGGAAAACGCGGGAGAAGGTCTGCTGCGAGACATAACCCAGGTCCATTGCGATATCAAAAATCGGACGCTCGGTGGTGCGCAACTCAACGGCGGCCAGTAACAGGCGGCGTTGGCGATTACATTCGCCaacgc                                                                                                                                             <  1:700626/150‑4 (MQ=255)
        tGCGATCAAACTGCCGACGGAAAACGCGGGAGAAGGTCTGCTGCGAGACATAACCCAGGTCCATTGCGATATCAAAAATCGGACGCTCGGTGGTGCGCAACTCAACGGCGGCCAGTAACAGGCGGCGTTGGCGAATGTAATCGCCAAGCg                                                                                                                                             <  2:594790/150‑1 (MQ=255)
        tGCGATCAAACTGCCGACGGAAAACGCGGGAGAAGGTCTGCTGCGAGACATAACCCAGGTCCATTGCGATATCAAAAATCGGACGCTCGGTGGTGCGCAACTCAACGGCGGCCAGTAACAGGCGGCGTTGGCGAATGTAATCGCCAAGCg                                                                                                                                             <  2:1129654/150‑1 (MQ=255)
         gCGATCAAACTGCCGACGGAAAACGCGGGAGAAGGTCTGCTGCGAGACATAACCCAGGTCCATTGCGATATCAAAAATCGGACGCTCGGTGGTGCGCAACTCAACGGCGGCCAGTAACAGGCGGCGTTGGCGATTACATTCGCCaacgcc                                                                                                                                            <  1:15481/150‑5 (MQ=255)
         gCGATCAAACTGCCGACGGAAAACGCGGGAGAAGGTCTGCTGCGAGACATAACCCAGGTCCATTGCGATATCAAAAATCGGACGCTCGGTGGTGCGCAACTCAACGGCGGCCAGTAACAGGCGGCGTTGGCGAATGTAATCGCCAAGCGt                                                                                                                                            >  1:1256621/1‑150 (MQ=255)
           gATCAAACTGCCGACGGAAAACGCGGGAGAAGGTCTGCTGCGAGACATAACCCAGGTCCATTGCGATATCAAAAATCGGACGCTCGGTGGTGCGCAACTCAACGGCGGCCAGTAACAGGCGGCGTTGGCGAATGTAATCGCCAAGCGTCt                                                                                                                                          <  1:1704819/150‑1 (MQ=255)
            aTCAAACTGCCGACGGAAAACGCGGGAGAAGGTCTGCTGCGAGACATAACCCAGGTCCATTGCGATATCAAAAATCGGACGCTCGGTGGTGCGCAACTCAACGGCGGCCAGTAACAGGCGGCGTTGGCGAATGTAATCGCCAAGCGTCTg                                                                                                                                         >  1:920084/1‑150 (MQ=255)
               aaaCTGCCGACGGAAAACGCGGGAGAAGGTCTGCTGCGAGACATAACCCAGGTCCATTGCGATATCAAAAATCGGACGCTCGGTGGTGCGCAACTCAACGGCGGCCAGTAACAGGCGGCGTTGGCGAATGTAATCGCCAAGCGTCTGATg                                                                                                                                      >  1:2358565/1‑150 (MQ=255)
               aaaCTGCCGACGGAAAACGCGGGAGAAGGTCTGCTGCGAGACATAACCCAGGTCCATTGCGATATCAAAAATCGGACGCTCGGTGGTGCGCAACTCAACGGCGGCCAGTAACAGGCGGCGTTGGCGAATGTAATCGCCAAGCGTCTGATg                                                                                                                                      <  2:323615/150‑1 (MQ=255)
                aaCTGCCGACGGAAAACGCGGGAGAAGGTCTGCTGCGAGACATAACCCAGGTCCATTGCGATATCAAAAATCGGACGCTCGGTGGTGCGCAACTCAACGGCGGCCAGTAACAGGCGGCGTTGGCGAATGTAATCGCCAAGCGTCTGATGc                                                                                                                                     <  1:110639/150‑1 (MQ=255)
                 aCTGCCGACGGAAAACGCGGGAGAAGGTCTGCTGCGAGACATAACCCAGGTCCATTGCGATATCAAAAATCGGACGCTCGGTGGTGCGCAACTCAACGGCGGCCAGTAACAGGCGGCGTTGGCGAATGTAATCGCCAAGCGTCTGATGCg                                                                                                                                    <  2:54118/150‑1 (MQ=255)
                    gCCGACGGAAAACGCGGGAGAAGGTCTGCTGCGAGACATAACCCAGGTCCATTGCGATATCAAAAATCGGACGCTCGGTGGTGCGCAACTCAACGGCGGCCAGTAACAGGCGGCGTTGGCGAATGTAATCGCCAAGCGTCTGATGCGTcc                                                                                                                                 >  2:2198155/1‑149 (MQ=255)
                     ccGACGGAAAACGCGGGAGAAGGTCTGCTGCGAGACATAACCCAGGTCCATTGCGATATCAAAAATCGGACGCTCGGTGGTGCGCAACTCAACGGCGGCAAGTAACAGGCGGCGTTGGCGAATGTAATCGCCAACCGTCTGATGcataac                                                                                                                                >  2:1048185/1‑145 (MQ=255)
                       gACGGAAAACGCGGGAGAAGGTCTGCTGCGAGACATAACCCAGGTCCATTGCGATATCAAAAATCGGACGCTCGGTGGTGCGCAACTCAACGGCGGCCAGTAACAGGCGGCGTTGGCGAATGTAATCGCCAAGCGTCTGATGCGTCACCg                                                                                                                              <  1:2368118/150‑1 (MQ=255)
                           gAAAACGCGGGAGAAGGTCTGCTGCGAGACATAACCCAGGTCCATTGCGATATCAAAAATCGGACGCTCGGTGGTGCGCAACTCAACGGCGGCCAGTAACAGGCGGCGTTGGCGAATGTAATCGCCAAGCGTCTGATGCGTCACCGTGc                                                                                                                           <  2:2007430/149‑1 (MQ=255)
                                cgcgGGAGAAGGTCTGCTGCGAGACATAACCCAGGTCCATTGCGATATCAAAAATCGGACGCTCGGTGGTGCGCAACTCAACGGCGGCCAGTAACAGGCGGCGTTGGCGAATGTAATCGCCAAGCGTCTGATGCGTCACCGTGCGGAACa                                                                                                                     >  2:1373997/1‑150 (MQ=255)
                                 gcgGGAGAAGGTCTGCTGCGAGACATAACCCAGGTCCATTGCGATATCAAAAATCGGACGCTCGGTGGTGCGCAACTCAACGGCGGCCAGTAACAGGCGGCGTTGGCGAATGTAATCGCCAAGCGTCTGATGCGTCACCGTGCGGAACAt                                                                                                                    <  1:909861/150‑1 (MQ=255)
                                 gcgGGAGAAGGTCTGCTGCGAGACATAACCCAGGTCCATTGCGATATCAAAAATCGGACGCTCGGTGGTGCGCAACTCAACGGCGGCCAGTAACAGGCGGCGTTGGCGAATGTAATCGCCAAGCGTCTGATGCGTCACCGTGCGGAACAt                                                                                                                    <  1:2076346/150‑1 (MQ=255)
                                  cgGGAGAAGGTCTGCTGCGAGACATAACCCAGGTCCATTGCGATATCAAAAATCGGACGCTCGGTGGTGCGCAACTCAACGGCGGCCAGTAACAGGCGGCGTTGGCGAATGTAATCGCCAAGCGTCTGATGCGTCACCGTGCGGAACAtt                                                                                                                   <  1:624693/150‑1 (MQ=255)
                                   gggAGAAGGTCTGCTGCGAGACATAACCCAGGTCCATTGCGATATCAAAAATCGGACGCTCGGTGGTGCGCAACTCAACGGCGGCCAGTAACAGGCGGCGTTGGCGAATGTAATCGCCAAGCGTCTGATGCGTCACCGTGCGGAACATTc                                                                                                                  >  2:1331873/1‑150 (MQ=255)
                                    ggAGAAGGTCTGCTGCGAGACATAACCCAGGTCCATTGCGATATCAAAAATCGGACGCTCGGTGGTGCGCAACTCAACGGCGGCCAGTAACAGGCGGCGTTGGCGAATGTAATCGCCAAGCGTCTGATGCGTCACCGTGCGGAACATTCg                                                                                                                 >  2:496344/1‑150 (MQ=255)
                                       gaAGGTCTGCTGCGAGACATAACCCAGGTCCATTGCGATATCAAAAATCGGACGCTCGGTGGTGCGCAACTCAACGGCGGCCAGTAACAGGCGGCGTTGGCGAATGTAATCGCCAAGCGTCTGATGCGTCACCGTGCGGAACATTCGTTg                                                                                                              <  2:68464/150‑1 (MQ=255)
                                       gaAGGTCTGCTGCGAGACATAACCCAGGTCCATTGCGATATCAAAAATCGGACGCTCGGTGGTGCGCAACTCAACGGCGGCCAGTAACAGGCGGCGTTGGCGAATGTAATCGCCAAGCGTCTGATGCGTCACCGTGCGGAACATTCGTTg                                                                                                              <  2:398354/150‑1 (MQ=255)
                                        aaGGTCTGCTGCGAGACATAACCCAGGTCCATTGCGATATCAAAAATCGGACGCTCGGTGGTGCGCAACTCAACGGCGGCCAGTAACAGGCGGCGTTGGCGAATGTAATCGCCAAGCGTCTGATGCGTCACCGTGCGGAACATTCGTTGc                                                                                                             <  1:2273649/150‑1 (MQ=255)
                                         aGGTCTGCTGCGAGACATAACCCAGGTCCATTGCGATATCAAAAATCGGACGCTCGGTGGTGCGCAACTCAACGGCGGCCAGTAACAGGCGGCGTTGGCGAATGTAATCGCCAAGCGTCTGATGCGTCACCGTGCGGAACATTCGTTGCa                                                                                                            <  1:1339161/150‑1 (MQ=255)
                                           gTCTGCTGCGAGACATAACCCAGGTCCATTGCGATATCAAAAATCGGACGCTCGGTGGTGCGCAACTCAACGGCGGCCAGTAACAGGCGGCGTTGGCGAATGTAATCGCCAAGCGTCTGATGCGTCACCGTGCGGAACATTCGTTGCAAg                                                                                                          <  2:940251/150‑1 (MQ=255)
                                           gTCTGCTGCGAGACATAACCCAGGTCCATTGCGATATCAAAAATCGGACGCTCGGTGGTGCGCAACTCAACGGCGGCCAGTAACAGGCGGCGTTGGCGAATGTAATCGCCAAGCGTCTGATGCGTCACCGTGCGGAACATTCGTTGCAAg                                                                                                          >  1:2157660/1‑150 (MQ=255)
                                               gctgcGAGACATAACCCAGGTCCATTGCGATATCAAAAATCGGACGCTCGGTGGTGCGCAACTCAACGGCGGCCAGTAACAGGCGGCGTTGGCGAATGTAATCGCCAAGCGTCTGATGCGTCACCGTGCGGAACATTCGTTGCAAGTAcc                                                                                                      <  1:737343/150‑1 (MQ=255)
                                                 tgcGAGACATAACCCAGGTCCATTGCGATATCAAAAATCGGACGCTCGGTGGTGCGCAACTCAACGGCGGCCAGTAACAGGCGGCGTTGGCGAATGTAATCGCCAAGCGTCTGATGCGTCACCGTGCGGAACATTCGTTGCAAGTACCAc                                                                                                    >  2:29177/1‑150 (MQ=255)
                                                  gcGAGACATAACCCAGGTCCATTGCGATATCAAAAATCGGACGCTCGGTGGTGCGCAACTCAACGGCGGCCAGTAACAGGCGGCGTTGGCGAATGTAATCGCCAAGCGTCTGATGCGTCACCGTGCGGAACATTCGTTGCAAGTACCACt                                                                                                   >  1:2032028/1‑150 (MQ=255)
                                                   cGAGTCATAACCCAGGTCCATTGCGATATCAAAAATCGGACGCTCGGTGGTGCGCAACTCAACGGCGGCCAGTAACAGGCGGCGTTGGCGAATGTAATCGCCAAGCGTCTGATGCGTCACCGTGCGGAACATTCGTTGCAAGTACCACtt                                                                                                  >  2:424818/1‑150 (MQ=255)
                                                     agaCATACCCCAGGTCCATTGCGATATCAAAAATCGGACGCTCGGTGGTGCGCAACTCAACGGCGGCCAGTAACAGGCGGCGTTGGCGTATGTAATCGACAAGCGTCTGATGCTTCAACGTGCTGTACATTGGTTGCAATTACCTCTTTg                                                                                                <  2:1650391/150‑1 (MQ=255)
                                                     agaCATAACCCAGGTCCATTGCGATATCAAAAATCGGACGCTCGGTGGTGCGCAACTCAACGGCGGCCAGTAACAGGCGGCGTTGGCGAATGTAATCGCCAAGCGTCTGATGCGTCACCGTGCGGAACATTCGTTGCAAGTACCACTTTg                                                                                                <  1:479311/150‑1 (MQ=255)
                                                        cATAACCCAGGTCCATTGCGATATCAAAAATCGGACGCTCGGTGGTGCGCAACTCAACGGCGGCCAGTAACAGGCGGCGTTGGCGAATGTAATCGCCAAGCGTCTGATGCGTCACCGTGCGGAACATTCGTTGCAAGTACCACTTTGAAt                                                                                             >  2:2387247/1‑150 (MQ=255)
                                                          tAACCCAGGTCCATTGCGATATCAAAAATCGGACGCTCGGTGGTGCGCAACTCAACGGCGGCCAGTAACAGGCGGCGTTGGCGAATGTAATCGCCAAGCGTCTGATGCGTCACCGTGCGGAACATTCGTTGCAAGTACCACTTTGAATAg                                                                                           >  1:191673/1‑150 (MQ=255)
                                                          tAACCCAGGTCCATTGCGATATCAAAAATCGGACGCTCGGTGGTGCGCAACTCAACGGCGGCCAGTAACAGGCGGCGTTGGCGAATGTAATCGCCAAGCGTCTGATGCGTCACCGTGCGGAACATTCGTTGCAAGTACCACTTTGAATAg                                                                                           <  2:1543906/150‑1 (MQ=255)
                                                           aaCCCAGGTCCATTGCGATATCAAAAATCGGACGCTCGGTGGTGCGCAACTCAACGGCGGCCAGTAACAGGCGGCGTTGGCGAATGTAATCGCCAAGCGTCTGATGCGTCACCGTGCGGAACATTCGTTGCAAGTACCACTTTGAATAGc                                                                                          >  1:2057568/1‑150 (MQ=255)
                                                                      aTTGCGATATCAAAAATCGGACGCTCGGTGGTGCGCAACTCAACGGCGGCCAGTAACAGGCGGCGTTGGCGAATGTAATCGCCAAGCGTCTGATGCGTCACCGTGCGGAACATTCGTTGCAAGTACCACTTTGAATAGCCTGAttttttt                                                                               <  2:2146392/150‑1 (MQ=255)
                                                                      aTTGCGATATCAAAAATCGGACGCTCGGTGGTGCGCAACTCAACGGCGGCCAGTAACAGGCGGCGTTGGCGAATGTAATCGCCAAGCGTCTGATGCGTCACCGTGCGGAACATTCGTTGCAAGTACCACTTTGAATAGCCTGAttttttt                                                                               <  2:342080/150‑1 (MQ=255)
                                                                        tGCGATATCAAAAATCGGACGCTCGGTGGTGCGCAACTCAACGGCGGCCAGTAACAGGCGGCGTTGGCGAATGTAATCGCCAAGCGTCTGATGCGTCACCGTGCGGAACATTCGTTGCAAGTACCACTTTGAATAGCCTGATTTTTTTGc                                                                             <  1:484860/150‑1 (MQ=255)
                                                                          cGATATCAAAAATCGGACGCTCGGTGGTGCGCAACTCAACGGCGGCCAGTAACAGGCGGCGTTGGCGAATGTAATCGCCAAGCGTCTGATGCGTCACCGTGCGGAACATTCGTTGCAAGTACCACTTTGAATAGCCTGATTTTTTTGCGa                                                                           >  2:236287/1‑150 (MQ=255)
                                                                           gATATCAAAAATCGGACGCTCGGTGGTGCGCAACTCAACGGCGGCCAGTAACAGGCGGCGTTGGCGAATGTAATCGCCAAGCGTCTGATGCGTCACCGTGCGGAACATTCGTTGCAAGTACCACTTTGAATAGCCTGATTTTTTTGCGAc                                                                          <  1:2290829/150‑1 (MQ=255)
                                                                              atCAAAAATCGGACGCTCGGTGGTGCGCAACTCAACGGCGGCCAGTAACAGGCGGCGTTGGCGAATGTAATCGCCAAGCGTCTGATGCGTCACCGTGCGGAACATTCGTTGCAAGTACCACTTTGAATAGCCTGATTTTTTTGCGACTAc                                                                       >  1:824363/1‑150 (MQ=255)
                                                                               tCAAAAATCGGACGCTCGGTGGTGCGCAACTCAACGGCGGCCAGTAACAGGCGGCGTTGGCGAATGTAATCGCCAAGCGTCTGATGCGTCACCGTGCGGAACATTCGTTGCAAGTACCACTTTGAATAGCCTGATTTTTTTGCGACTACa                                                                      <  1:1829914/150‑1 (MQ=255)
                                                                                cAAAAATCGGACGCTCGGTGGTGCGCAACTCAACGGCGGCCAGTAACAGGCGGCGTTGGCGAATGTAATCGCCAAGCGTCTGATGCGTCACCGTGCGGAACATTCGTTGCAAGTACCACTTTGAATAGCCTGATTTTTTTGCTACTACAt                                                                     >  1:187803/1‑150 (MQ=255)
                                                                                cAAAAATCGGACGCTCGGTGGTGCGCAACTCAACGGCGGCCAGTAACAGGCGGCGTTGGCGAATGTAATCGCCAAGCGTCTGATGCGTCACCGTGCGGAACATTCGTTGCAAGTACCACTTTGAATAGCCTGATTTTTTTGCGACTACAt                                                                     >  2:1420736/1‑150 (MQ=255)
                                                                                 aaaaaTCGGACGCTCGGTGGTGCGCAACTCAACGGCGGCCAGTAACAGGCGGCGTTGGCGAATGTAATCGCCAAGCGTCTGATGCGTCACCGTGCGGAACATTCGTTGCAAGTACCACTTTGAATAGCCTGATTTTTTTGCGACTACATc                                                                    >  2:2210162/1‑150 (MQ=255)
                                                                                 aaaaaTCGGACGCTCGGTGGTGCGCAACTCAACGGCGGCCAGTAACAGGCGGCGTTGGCGAATGTAATCGCCAAGCGTCTGATGCGTCACCGTGCGGAACATTCGTTGCAAGTACACTTTTAATTGCCTTGTTTTTTTTCGGCGACAtaa                                                                   >  1:1023429/1‑148 (MQ=255)
                                                                                     aTCGGACGCTCGGTGGTGCGCAACTCAACGGCGGCCAGTAACAGGCGGCGTTGGCGAATGTAATCGCCAAGCGTCTGATGCGTCACTGTGCGGAACATTCGTTGCAAGTACCACTTTGAATAGCCTGATTTTTTTGCGACTACATCAATg                                                                <  2:1256621/150‑1 (MQ=255)
                                                                                        ggACGCTCGGTGGTGCGCAACTCAACGGCGGCCAGTAACAGGCGGCGTTGGCGAATGTAATCTCCAAGCGTCTGATGCGTCACCGTGCGGAACATTCGTTGCAAGTACCACTTTGAATAGCCTGATTTTTTTGCGACTACATCAATGTTa                                                             <  2:260018/150‑1 (MQ=255)
                                                                                        ggACGCTCGGTGGTGCGCAACTCAACGGCGGCCAGTAACAGGCGGCGTTGGCGAATGTAATCGCCAAGCGTCTGATGCGTCACCGTGCGGAACATTCGTTGCAAGTACCACTTTGAATAGCCTGATTTTTTTGCGACTACATCAATGTTa                                                             >  1:1075371/1‑150 (MQ=255)
                                                                                           cGCTCGGTGGTGCGCAACTCAAAGGCGGCCAGTAACAGGCGGCGTTGTCGAATCTAATCGCCAAGCGTCTGATGCGTCACCGTCCGGAACATTCGTTGCAAGTACCACTTTGAATAGCCTGATTTTTTTGCGACTACATCAATCTTAAGc                                                          >  2:11030/1‑150 (MQ=255)
                                                                                            gCTCGGTGATGCGCAACTCAACGGCGGCCAGTAACAGGCGGCGTTGGCGAATGTAATCGCCAAGCGTCTGATGCGTCACCGTGCGGAACATTCGTTGCAAGTACCACTTTGAATAGCCTGATTTTTTTGCGACTACCTCAATGTTAAGCg                                                         >  2:2112302/1‑150 (MQ=255)
                                                                                             cTCGGTGGTGCGCAACTCAACGGCGGCCAGTAACAGGCGGCGTTGGCGAATGTAATCGCCAAGCGTCTGATGCGTCACCGTGCGGAACATTCGTTTCAAGTACCACTTTGAATAGCCTGATTTTTTTGCGACTACATCAATGTTAAGCgg                                                        >  1:1813390/1‑150 (MQ=255)
                                                                                             cTCGGTGGTGCGCAACTCAACGGCGGCCAGTAACAGGCGGCGTTGGCGAATGTAATCGCCAAGCGTCTGATGCGTCACCGTGCGGAACATTCGTTGCAAGTACCACTTTGAATAGCCTGATTTTTTTGCGACTACATCAATGTTAAGCgg                                                        >  2:157630/1‑150 (MQ=255)
                                                                                              tCGGTGGTGCGCAACTCAACGGCGGCCAGTAACAGGCGGCGTTGGCGAATGTAATCGCCAAGCGTCTGATGCGTCACCGTGCGGAACATTCGTTGCAAGTACCACTTTGAATAGCCTGATTTTTTTGCGACTACATCAATGTTAAGCGGc                                                       <  2:569376/150‑1 (MQ=255)
                                                                                               cGGTGGTGCGCAACTCAACGGCGGCCAGTAACAGGCGGCGTTGGCGAATGTAATCGCCAAGCGTCTGATGCGTCACCGTGCGGAACATTCGTTGCAAGTACCACTTTGAATAGCCTGATTTTTTTGCGACTACATCAATGTTAAGCGGCt                                                      <  2:926204/150‑1 (MQ=255)
                                                                                                ggtggtGCGCAACTCAACGGCGGCCAGTAACAGGCGGCGTTGGCGAATGTAATCGCCAAGCGTCTGATGCGTCACCGTGCGGAACATTCGTTGCAAGTACCACTTTGAATAGCCTGATTTTTTTGCGACTACATCAATGTTAAGCGGCTg                                                     <  2:2057568/150‑1 (MQ=255)
                                                                                                  tggtgCGCAACTCAACGGCGGCCAGTAACAGGCGGCGTTGGCGAATGTAATCGCCAAGCGTCTGATGCGTCACCGTGCGGAACATTCGTTGCAAGTACCACTTTGAATAGCCTGATTTTTTTGCGACTACATCAATGTTAAGCGGCTGGt                                                   <  1:2326613/150‑1 (MQ=255)
                                                                                                   ggtgCGCAACTCAACGGCGGCCAGTAACAGGCGGCGTTGGCGAATGTAATCGCCAAGCGTCTGATGCGTCACCGTGCGGAACATTCGTTGCAAGTACCACTTTGAATAGCCTGATTTTTTTGCGACTACATCAATGTTAAGCGGCTGGTc                                                  <  1:1260577/150‑1 (MQ=255)
                                                                                                   ggtgCGCAACTCAACGGCGGCCAGTAACAGGCGGCGTTGGCGAATGTAATCGCCAAGCGTCTGATGCGTCACCGTGCGGAACATTCGTTGCAAGTACCACTTTGAATAGCCTGATTTTTTTGCGACTACATCAATGTTAAGCGGCTGGTc                                                  >  1:376661/1‑150 (MQ=255)
                                                                                                      gcgcAACTCAACGGCGGCCAGTAACAGGCGGCGTTGGCGAATGTAATCGCCAAGCGTCTGATGCGTCACCGTGCGGAACATTCGTTGCAAGTACCACTTTGAATAGCCTGATTTTTTTGCGACTACATCAATGTTAAGCGGCTGGTCAat                                               >  1:1245560/1‑150 (MQ=255)
                                                                                                       ggcAACTCAACGGCGGCCAGTAACAGGCGGCGTTGGCGAATGTAATCGCCAAGCGTCTGATGCGTCACCGTGCGGAACATTCGTTGCAAGTACCACTTTGAATAGCCTGATTTTTTTGCGACTACATCAATGTTAAGCGGCTGGTCAata                                              <  2:2291659/149‑1 (MQ=255)
                                                                                                       cgcAACTCAACGGCGGCCAGTAACAGGCGGCGTTGGCGAATGTAATCGCCAAGCGTCTGATGCGTCACCGTGCGGAACATTCGTTGCAAGTACCACTTTGAATAGCCTGATTTTTTTGCGACTACATCAATGTTAAGCGGCTGGTCAata                                              <  2:1893890/150‑1 (MQ=255)
                                                                                                       cgcAACTCAACGGCGGCCAGTAACAGGCGGCGTTGGCGAATGTAATCGCCAAGCGTCTGATGCGTCACCGTGCGGAACATTCGTTGCAAGTACCACTTTGAATAGCCTGATTTTTTTGCGACTACATCAATGTTAAGCGGCTGGTCAata                                              >  2:1008795/1‑150 (MQ=255)
                                                                                                       cgcAACTCAACGGCGGCCAGTAACAGGCGGCGTTGGCGAATGTAATCGCCAAGCGTCTGATGCGTCACCGTGCGGAACATTCGTTGCAAGTACCACTTTGAATAGCCTGATTTTTTTGCGACTACATCAATGTTAAGCGGCTGGTCAata                                              <  2:1833853/150‑1 (MQ=255)
                                                                                                       cgcAACTCAACGGCGGCCAGTAACAGGCGGCGTTGGCGAATGTAATCGCCAAGCGTCTGATGCGTCACCGTGCGGAACATTCGTTGCAAGTACCACTTTGAATAGCCTGATTTTTTTGCGACTACATCAATGTTAAGCGGCTGGTCAata                                              <  1:1000600/150‑1 (MQ=255)
                                                                                                           aCTCAACGGCGGCCAGTAACAGGCGGCGTTGGCGAATGTAATCGCCAAGCGTCTGATGCGTCACCGTGCGGAACATTCGTTGCAAGTACCACTTTGAATAGCCTGATTTTTTTGCGACTACATCAATGTTAAGCGGCTGGTCAATATGCt                                          >  1:149666/1‑150 (MQ=255)
                                                                                                             tCAACGGCGGCCAGTAACAGGCGGCGTTGGCGAATGTAATCGCCAAGCGTCTGATGCGTCACCGTGCGGAACATTCGTTGCAAGTACCACTTTGAATAGCCTGATTTTTTTGCGACTACATCAATGTTAAGCGGCTGGTCAATATGCTCg                                        >  1:760835/1‑150 (MQ=255)
                                                                                                             tCAACGGCGGCCAGTAACAGGCGGCGTTGGCGAATGTAATCGCCAAGCGTCTGATGCGTCACCGTGCGGAACATTCGTTGCAAGTACCACTTTGAATAGCCTGATTTTTTTGCGACTACATCAATGTTAAGCGGCTGGTCAATATGCTCg                                        <  2:879462/150‑1 (MQ=255)
                                                                                                             tCAACGGCGGCCAGTAACAGGCGGCGTTGGCGAATGTAATCGCCAAGCGTCTGATGCGTCACCGTGCGGAACATTCGTTGCAAGTACCACTTTGAATAGCCTGATTTTTTTGCGACTACATCAATGTTAAGCGGCTGGTCAATATGCTCg                                        <  2:906639/150‑1 (MQ=255)
                                                                                                                 cggcggCCAGTAACAGGCGGCGTTGGCGAATGTAATCGCCAAGCGTCTGATGCGTCACCGTGCGGAACATTCGTTGCAAGTACCACTTTGAATAGCCTGATTTTTTTGCGACTACATCAATGTTAAGCGGCTGGTCAATATGCTCGTCaa                                    <  1:2124088/150‑1 (MQ=255)
                                                                                                                    cggcCAGTAACAGGCGGCGTTGGCGAATGTAATCGCCAAGCGTCTGATGCGTCACCGTGCGGAACATTCGTTGCAAGTACCACTTTGAATAGCCTGATTTTTTTGCGACTACATCAATGTTAAGCGGCTGGTCAATATGCTCGTCAATcc                                 >  2:47922/1‑150 (MQ=255)
                                                                                                                         aGTAACAGGCGGCGTTGGCGAATGTAATCGCCAAGCGTCTGATGCGTCACCGTGCGGAACATTCGTTGCAAGTACCACTTTGAATAGCCTGATTTTTTTGCGACTACATCAATGTTAAGCGGCTGGTCAATATGCTCGTCAATCCATGCg                            <  1:2210162/150‑1 (MQ=255)
                                                                                                                                ggcggcGTTGGCGAATGTAATCGCCAAGCGTCTGATGCGTCACCGTGCGGAACATTCGTTGCAAGTACCACTTTGAATAGCCTGATTTTTTTGCGACTACATCAATGTTAAGCGGCTGGTCAATATGCTCGTCAATCCATGCGATAAGAt                     >  2:1165594/1‑150 (MQ=255)
                                                                                                                                  cggcgTTTGCGAATGTAATCGCCAATCGGCTGATGCGTCACCGTGCGGAACATTCGTTGCAAGTACCACTTTGAATATCCTGATTTTTTTGCGACTACATCAATGTTAAGCGGCTGGGCAATATGCTCGTCAATCCATGCGATAAGATcc                   <  2:191673/150‑1 (MQ=255)
                                                                                                                                  cggcgTTGGCGAATGTAATCGCCAAGCGTCTGATGCGTCACCGTGCGGAACATTCGTTGCAAGTACCACTTTGAATAGCCTGATTTTTTTGCGACTACATCAATGTTAAGCGGCTGGTCAATATGCTCGTCAATCCATGCGATAAGATcc                   >  2:824065/1‑150 (MQ=255)
                                                                                                                                  cggcgTTGGCGAATGTAATCGCCAAGCGTCTGATGCGTCACCGTGCGGAACATTCGTTGCAAGTACCACTTTGAATAGCCGGATTTTTTTGCGACTACATCAATGTTAAGCGGCTGGTCAATATGCTCGTCAATCCATGCGATAAGATcc                   <  2:1372234/150‑1 (MQ=255)
                                                                                                                                    gcgTTGGCGAATGTAATCGCCAAGCGTCTGATGCGTCACCGTGCGGAACATTCGTTGCAAGTACCACTTTGAATAGCCTGATTTTTTTGCGACTACATCAATGTTAAGCGGCTGGTCAATATGCTCGTCAATCCATGCGATAAGATCCTg                 <  2:2464217/150‑1 (MQ=255)
                                                                                                                                    cgcttGGCGATTACATTCGCCAAGCGTCTGATGCGTCACCGTGCGGAACATTCGTTGCAAGTACCACTTTGAATAGCCTGATTTTTTTGCGACTACATCAATGTTAAGCGGCTGGTCAATATGCTCGTCAATCCATGCGATAAGATCCTg                 >  1:1147014/4‑150 (MQ=255)
                                                                                                                                      gTTGGCGAATGTAATCGCCAAGCGTCTGATGCGTCACCGTGCGGAACATTCGTTGCAAGTACCACTTTGAATAGCCTGATTTTTTTGCGACTACATCAATGTTAAGCGGCTGGTCAATATGCTCGTCAATCCATGCGATAAGATCCTGaa               >  2:1455793/1‑150 (MQ=255)
                                                                                                                                      cttGGCGATTTCATTCGCCAAGCGTCTGATGCGTCACCGTGCGGAACATTCGTTGCAAGTACCACTTTGAATAGCCTGATTTTTTTGCGACTACATCAATGTTAAGCGGCTGGTCAATATGCTCGTCAATCCATGCGATAAGATCCTGaa               >  1:945665/2‑150 (MQ=255)
                                                                                                                                      cttGGCGATTACATTCGCCAAGCGTCTGATGCGTCACCGTGCGGAACATTCGTTGCAAGTACCACTTTGAATAGCCTGATTTTTTTGCGACTACATCAATGTTAAGCGGCTGGTCAATATGGTCGTCAATCCATGCGATAAGATCCTGaa               >  2:655067/2‑150 (MQ=255)
                                                                                                                                      cttGGCGATTACATTCGCCAAGCGTCTGATGCGTCACCGTGCGGAACATTCGTTGCAAGTACCACTTTGAATAGCCTGATTTTTTTGCGACTACATCAATGTTAAGCGGCTGGTCAATATGCTCGTCAATCCATGCGATAAGATCCTgta               >  2:501988/2‑148 (MQ=255)
                                                                                                                                      cttGGCGATTACATTCGCCAAGCGTCTGATGCGTCACCGTGCGGAACATTCGTTGCAAGTACCACTTTGAATAGCCTGATTTTTTTGCGACTACATCAATGTTAAGCGGCTGGTCAATATGCTCGTCAATCCATGCGATAAGATCCTGaa               >  2:1287569/2‑150 (MQ=255)
                                                                                                                                         ggCGAATGTAATCGCCAAGCGTCTGATGCGTCACCGTGCGGAACATTCGTTGCAAGTACCACTTTGAATAGCCTGATTTTTTTGCGACTACATCAATGTTAAGCGGCTGGTCAATATGCTCGTCAATCCATGCGATAAGATCCTGaataa            >  1:1970346/1‑150 (MQ=255)
                                                                                                                                         ggCGAATGTAATCGCCAAGCGTCTGATGCGTCACCGTGCGGAACATTCGTTGCAAGTACCACTTTGAATAGCCTGATTTTTTTGCGACTACATCAATGTTAAGCGGCTGGTCAATATGCTCGTCAATCCATGCGATAAGATCCTGaataa            >  2:1790334/1‑150 (MQ=255)
                                                                                                                                                  aaTCGCCAAGCGTCTGATGCGTCACCGTGCGGAACATTCGTTGCAAGTACCACTTTGAATAGCCTGATTTTTTTGCGACTACATCAATGTTAAGCGGCTGGTCAATATGCTCGTCAATCCATGCGATAAGATCCTGAATAATTTTCTGAt   >  1:327361/1‑150 (MQ=255)
                                                                                                                                                   aTCGCCAAGCGTCTGATGCGTCACCGTGCGGAACATTCGTTGCAAGTACCACTTTGAATAGCCTGATTTTTTTGCGACTACATCAATGTTAAGCGGCTGGTCAATATGCTCGTCAATCCATGCGATAAGATCCTGAATAATTTTCTGATg  >  2:77575/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
GCTGGGAGTGCGATCAAACTGCCGACGGAAAACGCGGGAGAAGGTCTGCTGCGAGACATAACCCAGGTCCATTGCGATATCAAAAATCGGACGCTCGGTGGTGCGCAACTCAACGGCGGCCAGTAACAGGCGGCGTTGGCGAATGTAATCGCCAAGCGTCTGATGCGTCACCGTGCGGAACATTCGTTGCAAGTACCACTTTGAATAGCCTGATTTTTTTGCGACTACATCAATGTTAAGCGGCTGGTCAATATGCTCGTCAATCCATGCGATAAGATCCTGAATAATTTTCTGATG  >  NC_000913/4277081‑4277377

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