Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,998,147 C→A 5.3% Y55* (TAC→TAA uvrD → DNA‑dependent ATPase I and helicase II

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,998,1470CA 5.3% 247.7 / 2.3 96Y55* (TAC→TAAuvrDDNA‑dependent ATPase I and helicase II
Reads supporting (aligned to +/- strand):  ref base C (34/56);  new base A (3/3);  total (37/59)
Fisher's exact test for biased strand distribution p-value = 6.73e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TGCTCGACAGCCTTAATGACAAACAGCGCGAAGCGGTGGCCGCGCCACGCAGCAACCTTCTGGTGCTGGCGGGCGCGGGCAGTGGTAAGACGCGCGTACTGGTGCATCGTATCGCCTGGTTGATGAGCGTGGAAAACTGCTCGCCATACTCGATTATGGCGGTGACGTTTACCAACAAAGCGGCGGCGGAGATGCGTCATCGTATCGGGCAACTGATGGGCACGAGCCAGGGCGGTATGTGGGTCGGCACCTTCCACGGGCTGGCGCACCGTTTGCTGCGTGCGCACCATATGGAC  >  NC_000913/3997999‑3998294
                                                                                                                                                    |                                                                                                                                                   
tGCTCGACAGCCTTAATGACAAACAGCGCGAAGCGGTGGCCGCGCCCCCCACCAACCATCTGGTGCTGGCGGGCGCGGGCAGTGGTAAGCCCCGCGTACTGGTGCATCGTATCGCCGGGTGGATGAGCGTGGAAAACTGCTAGCCATAc                                                                                                                                                     >  1:555809/1‑149 (MQ=255)
tGCTCGACAGCCTTAATGACAAACAGCGCGAAGCGGTGGCCGCGCCACGCAGCAACCTTCTGGTGCTGGCGGGCGCGGGCAGTGGTAAGACGCGCGTACTGGTGCATCGTATCGCCTGGTTGATGAGCGTGGAAAACTGCTCGCCATAc                                                                                                                                                     >  1:80381/1‑149 (MQ=255)
tGCTCGACAGCCTAAATGACAAACACCGCGAAGCGGTGGCCGCGCCACGCAGCAACCTTCTGGGGCTGGCGGGCGCGGGCAGTGGTAAGACGCGCGTACTGGTGCATCGTATCGCCTGGTTGATGAGCGTGGAAAACTGCTCGCCATAc                                                                                                                                                     <  1:887354/149‑1 (MQ=255)
 gCTCGACAGCCTTAATGACAAACAGCGCGAAGCGGTGGCCGCGCCACGCAGCAACCTTCTGGTGCTGGCGGGCGCGGGCAGTGGTAAGACGCGCGTTCTGGTGCATCGTATCGCCTGGTTGATGAGCGTGGAAAACTGCTCGCCCTAcg                                                                                                                                                    >  1:1353800/1‑148 (MQ=255)
  cTCGCCAGCCTAACTGACAAACCGAGCGAAGCGGGGGCCGCGCCACGCAGCAACCTGCTGGTGCTGGCGGGCGCGGGCAGTGGTAAGACGCGCGTACTGGTGCATCGTATCGCCTGGTTGATGAGCGTGGAAATCTGCTCGCCATACTc                                                                                                                                                   <  2:1455274/149‑1 (MQ=255)
   tCGACAGCCTTAATGACAAACAGCGCGAAGCGGTGGCCGCGCCACGCAGCAACCTTCTGGTGCTGGCGGGCGCGGGCAGTGGTAAGACGCGCGTACTGGTGCATCGTATCGCCTGGTTGATGAGCGTGGAAAACTGCTCGCCATACTCg                                                                                                                                                  >  2:20509/1‑149 (MQ=255)
     gACAGCCTTAATGACAAACAGCGCGAAGCGGTGGCCGCCCCCCGGACCAACCCTGTGGTGCTGGCGGGCGCGGGGAGTGGGAAGACGCGCGTACGGGTGGATCGTTTCGGCTGGTTGATGAGCGTGGAAAACTGCTCGCCAAACTCGAt                                                                                                                                                >  1:187993/1‑149 (MQ=255)
      ccctccTTAATGACAAACAGCGCGAAGCTGTGGCCGAGCCACGCAGCAACCTTCTGGTGCTGGCGGGCGCGGGCAGTGGTAAGACGCGCGTACTGGTGCATCGTATCGCCTGGTTGATGAGCGTGGAAAACTGCTCGCCATACTCGAtt                                                                                                                                               <  1:144142/145‑1 (MQ=255)
         cacTTAATGACAAACAGCCCGAAGCGGTGGCCGCGCCACGCAGCAACCTTCTGGTGCTGGCGGGCGCGGGCAGTGGTAAGACGCGCGTACTGGTGCATCGTATCGCCTGGTTGATGAGCGTGGAAAACTGCTCGCCATACTCGATTATg                                                                                                                                            <  1:1300838/147‑1 (MQ=255)
          cgttAAGGACAGACAGCGCGAAGCGGTGGCCGCGCCACGCATCAACCTTGTGGTGCTGGCGGGCGCGGTCAGTGGTAAGACGCGCGTACTGGTGGATCGTAGCGCCTGTTTGATGAGCGTGGAATACTGCTCGCCATACTCGATTATgg                                                                                                                                           <  1:1522564/147‑1 (MQ=255)
          ccTTAATGACAAACAGCGCGAAGCGGTGGCCGCGCCACGCAGCAACCTTCTGGTGGTGGCGGGCGCGGGCAGTGGTAAGACGCGCGTACTGGTGCATCGTAACGCCGGGTTGATGGGCGTGGAAAACTGCTCGCCATAATCGATTATgg                                                                                                                                           >  2:146108/1‑149 (MQ=255)
           cTTAATGACAAACAGCGCGAAGCGGTGGCCGCGCCACGCAGCAACCTTCTGGTGCTGGCGGGCGCGGGCAGTGGTAAGACGCGCGTACTGGTGCATCGTATCGCCTGGTTGATGAGCGTGGAAAACTGCTCGCCATACTCGATTATGGc                                                                                                                                          <  2:545301/149‑1 (MQ=255)
              aaggACAAACAGCCCGAAGCGGTGGCCGCGCCACGCAGCAACCTTCTGGTGCTGGCGGGCGCGGGCAGTGGTAAGACGCGCGTACTGGTGCATCGTATCGCCTGGTTGATGAGCGTGGAAAACTGCTCGCCATAATCGAGTATGGCgag                                                                                                                                       <  1:722051/146‑3 (MQ=255)
              aaTGACAAACAGCGCGAAGCGGTGGCCGCGCCACGCAGCAACCTTCTGGTGCTGGCGGGCGCGGGCAGTGGTAAGACGCGCGTACTGGTGCATCGTATCGCCTGGTTGATGAGCGTGGAAAACTGCTCGCCATACTCGATTATGGCGGt                                                                                                                                       <  1:1235835/149‑1 (MQ=255)
              aaTGACAAACAGCGCGAAGCGGTGGCCGCGCCACGCAGCAACCTTCTGGTGCTGGCGGGCGCGGGCAGTGGTAAGACGCGCGTACTGGTGCATCGTATCGCCTGGTTGATGAGCGTGGAAAACTGCTCGCCATAATCGAGTATGGCgag                                                                                                                                       <  2:1107829/149‑3 (MQ=255)
               aggacccaCAGCGCGCAGCGGCGGCCGCGCCACGCAGCAACCTTCTGGGGCTGGCGGGCGCGGGCAGTGGTAAGACGCGCGTACTGCTGCATCGTATCGCCTGGTTGATGAGCGTGGAAAACTGCTCGCCATACTCGATTATGGCGGTg                                                                                                                                      <  1:867245/142‑1 (MQ=255)
                 gccAAACAGCGCGCAGCGGTGGCCGCGCCACGCAGCAACCTTCTGGTGCTGGCGGGCGCGGGCAGTGGTAAGACGCGCGTACTGGTGCATCGTATCGCCTGGTTGATGAGCGGGGAAAACTGCTCGCCATACTCGATTATGGCGGTGAc                                                                                                                                    <  2:861411/147‑1 (MQ=255)
                 gACAAACAGCGCGAAGCGGTGGCCGCGCCACGCAGCAACCGTCTGGTGCTGGCGGGCGCGGGCAGTGGTAAGACGCGCGGACTGGTGCATCGTATCGCCGGGTTGATGAGCGTGGAAAACTGCTCGCCAAACTCGATTATTGCGGGGAc                                                                                                                                    >  2:625241/1‑149 (MQ=255)
                   cAAACAGCGCGAAGCGGTGGCCGCGCCACGCAGCAACCTTCTGGTGCTGGCGGGCGCGGGCAGTGGTAAGACGCGCGTACTGGTGCATCGTATCGCCTGGTTGATGAGCGTGGAAAACTGCTCGCCATACTCGATTATGGCGGTGACGt                                                                                                                                  <  1:763059/149‑1 (MQ=255)
                   cAAACAGCGCGAAGCGGTGGCCGCGCCACGCAGCAACCTTCTGGTGCTGGCGGGCGCGGGCAGTGGTAAGACGCGCGTACTGGTGCATCGTATCGCCTGGTTGATGAGCGTGGAAAACTGCTCGCCATACTCGATTATGGCGGTGACGt                                                                                                                                  >  2:975744/1‑149 (MQ=255)
                    aaaCAGCGCGAAGCGGTGGCCCCGCCACGCAGCAACCTTCTGGTGCTGGCGGGCGCGGGCAGTGGTAAGACGCGCGTACTGGTGCATCGTATAGCCGGGGTGATGAGCGTGGAAAACGGCTCGCCATACTCGATTATGGCGGGGACGtt                                                                                                                                 >  1:518460/1‑149 (MQ=255)
                       ccgcgcgcAGCGGTGGCCGCGGCACGCAGCAACCTTCTGTGGCGGGCGGGCGCGGGCAGTGGTAAGACGCGCGTACTGGTGCATCGTATCGCCTGGTTGATGAGCGTGGAAAACTGCTCGCCATACTCGATTATGGCGGTGACGTTTAc                                                                                                                              <  2:1017205/147‑1 (MQ=255)
                       catcgcgAAGCGGTGGCCGCGCCACGCAGCAACCTTCTGGTGCTGGCGGGCGCGGGCAGTGGTAAGACGCGCGTACTGGTGCATCGTATCGCCTGGTTGATGAGCGTGGAAAACTGCTCGCCATACTCGATTATGGCGGTGACGTTTAc                                                                                                                              <  2:174494/146‑1 (MQ=255)
                         gcgcGAAGCGGTGGCCGCGCCCCGCAGCCACCTGCGGGTGCTGGCGGGCGGGGGGAGTGGTAAGACGCGCGGACTGGTGCATAGTATCGCCTGGTGTATGAGGGTGGAAAACTGGACGCCATACTCGATTTTTGCGGTGACCTTTTCca                                                                                                                            >  2:1395933/1‑149 (MQ=255)
                          cgcgAAGCGGTGGCCGCGCCACGCAGCAACCTTCTGGGGCTGGCGGGCGCGGGCAGTGGTAAGACGCGCGTACTGGTGCATCGTATCGCCTGGTTGATGAGCGTGGAAAACTGCTCGCCATACTCGATTATGGCGGTGACGTTTACcaa                                                                                                                           <  1:478294/149‑1 (MQ=255)
                           gcgAAGCGGTGGCCGCGCCACGCAGCAACCTTGTGGTGCTGGCGGGCGCGGGCAGTGGTAAGACGCGCGTGATGGTGCCACGGATAGGCTGGTAGATTAGCGTGGTAAACGTCTAGGCGTTCTCGGTGATGGGGGTGACGTTTCCcgac                                                                                                                          >  2:228103/1‑149 (MQ=255)
                           gcgAAGCGGTGGCCGCGCCACGCAGCAACCTTCTGGTGCTGGCGGGCGCGGGCAGTGGTAAGACGCGCGTACTGGTGCATCGTCTCGCGTGGTTGATGGGAGTTGAAAACTGCTCGGGGTACTCGTTGATGGCGGTGACGTTTTCccac                                                                                                                          >  2:233915/1‑149 (MQ=255)
                             gAAGCGGTGGCCGCGCCACGCAGCAACCTTCTGGTGCTGGCGGGCGCGGGCAGTGGTAAGACGCGCGTACTGGTGCATCGTATCGCCTGGTTGATGAGCGTGGAAAACTGCTCGCCATACTCGATTATGGCGGTGACGTTTACcaacaa                                                                                                                        >  2:63270/1‑149 (MQ=255)
                                gCGGTGGCCGCGCCACGCAGCATCCTTCTGGTGCTGGCGTGCGCGGGCAGTGGTAAGACGCGCGTACTGGTGCATCGTTTCGCCTGGTTGATGAGCGTGGAAAACTGCTCGCCATACTCGATTATGGCGGTGACGTTTACGAACAAAgc                                                                                                                     <  2:5381/149‑1 (MQ=255)
                                gCGGTGGCCGCGCCACGCAGCAACCTTCTGGTGCGGGCGGGCGCGGGCAGTGGTTAGACGCGCGTACTGGGGCATCGGATAGCCTGGGTGATGCGCGTTGAAAAAGGCGCGCCATAATCGATTATGGCGGGGACGTTTACCAACAAAgc                                                                                                                     >  2:1499100/1‑149 (MQ=255)
                                     ggCCGCGCCACGCAGCAACCTTCTGGTGCTGGCGGGCGCGGGCAGTGGTAAGACGCGCGTACTGGTGCATCGTATCGCCTGGTTGATGAGCGTGGAAAACTGCTCGCCATACTCGATTATGGCGGTGACGGTTACCAACAAAgcggcgg                                                                                                                >  1:303981/1‑149 (MQ=255)
                                      gCCGCGCCACGCAGCAACCTTCTGGTGCTGGCGGGCGCGGGCAGTGGTAAGACGCGCGTACTGGTGCTTCGTATCGCCTGGTTGATGAGCGTGGAAAACTGCTCGCCATACTCGATTATGGCGGTGACGTTTACCAACAAAgcggcggc                                                                                                               >  1:9517/1‑149 (MQ=255)
                                      gCCGCGCCACGCAGCAACCTTCTGGTGCTGGCGGGCGCGGGCAGTGGTAAGACGCGCGTACTGGTGCATCGTATCGCCTGGTTGATGAGCGTGGAAAACTGCTCGCCATACTCGATTATGGCGGTGACGTTTACCAACAAAgcggcggc                                                                                                               <  1:1520650/149‑1 (MQ=255)
                                      gCCGCGCCACGCAGCAACCTTCTGGTGCTGGCGGGCGCGGGCAGTGGTAAGACGCGCGTACTGGTGCATCGTATCGCCTGGTTGATGAGCGTGGAAAACTGCTCGCCAAACTCGATTATGGCGGTGACGTTTACCAACAAAgcggcggc                                                                                                               <  2:1542478/149‑1 (MQ=255)
                                       ccGCGCCACGCAGCAACCTTCTGGTGCTGGCGGGCGCGGGCAGTGGTAAGACGCGCGTACTGGTGCATCGTATCGCCTGGTTGATGAGCGTGGAAAACTGCTCGCCATACTCGAGTATGGCGGTGACGTTTACCAACAAAgcggcggcg                                                                                                              >  2:1488445/1‑149 (MQ=255)
                                              aCGCAGCAACCTTCTGGTGCTGGCGGGCGCGGGCAGTGGTAAGACGCGCGTACTGGTGCATCGTATCGCCTGGTTGATGAGCGTGGAAAACTGCTCGCCATACTCGATTATGGCGGTGACGTTTACCAACAAAGCGGCGGCGGAGATGc                                                                                                       <  2:1208793/149‑1 (MQ=255)
                                                  agcaACCTTCTGGTGCTGGCGGGCGCGGGCAGTGGTACGACGCGCGTACTGGTGCATCGTATCGCCTGGTTGATGAGCGTGGAAAACTGCTCGCCATACTCGATTATGGCGGTGACGTTTACCAACAAAGCGGCGGCGGAGATGCGTCa                                                                                                   <  2:383394/149‑1 (MQ=255)
                                                    caACCTTCTGGTGCTGGCGGGCGCGGGCAGTGGTAAGACGCGCGTACTGGTGCATCGTATCGCCTGGTTGATGAGCGTGGAAAACTGCTCGCCATACTCGATTATGGCGGTGACGTTTACCAACAAAGCGGCGGCGGAGATGCGTCATc                                                                                                 <  2:781567/149‑1 (MQ=255)
                                                          tCTGGTGCTGGCGGGCGCGGGCAGTGGTAAGACGCGCGTACTGGTGCATCGTATCGCCTGGTTGATGAGCGTGGAAAACTGCTCGCCATACTCGATTATGGCGGTGACGTTTACCAACAAAGCGGCGGCGGAGATGCGTCATCGTATCg                                                                                           <  1:975744/149‑1 (MQ=255)
                                                            tgtcgCTGGGGTGCGCGTTCATTTGTAAGACGCGCGTACTGGTGCATCGTCTCGGGTGGTTGATTTGCGTGGATACCTGCTCGCCATAATCGATTATGGCGGTGACGTTTACCAACAAAGCGGCGGCGGAGATGCGTCATCGTATCggg                                                                                         <  1:1452240/145‑1 (MQ=255)
                                                                gCTGGCGGTCGCGGGCAGTTGTAAGACGCGCGTACTGGTGCATCGTATCGCCTGGTTGATGAGCGTGGAAAACTGCTCGCCATACTCGATTATGGCGGTGACGTTTACCAACAAAGCGGCGGCGGAGATGCGTCATCGTATCGGGCAAc                                                                                     <  2:1414165/149‑1 (MQ=255)
                                                                  tGGCGGGCGCGGGCAGTGGTAAGCCGCGCGTACTGGTGCATCGTATCGCCTGGTTGATGAGCGTGGAAAACTGCTCGCCATACTCGATTATGGCGGTGACGTTTACCAACAAAGCGGCGGCGGAGATGCGTCATCGTATCGGGCAACTg                                                                                   >  2:1416320/1‑149 (MQ=255)
                                                                  tGGCGGGCGCGGGCAGTGGTAAGACGCGCGTACTGGTGCATCGTATCGCCTGGTTGATGAGCGTGGAAAACTGCTCGCCATACTCGAGTATGGCGGTGACGTTTACCAACAAAGCGGCGGCGGAGATGCGTCATCGTATCGGGCAACTg                                                                                   <  1:584045/149‑1 (MQ=255)
                                                                     cgggcgcgggcAGTGGTAAGACGCGCGTACTGGTGCATCGGATCGCCTGGTTGATGAGCGTGGAAAACTGCTCGCCATACTCGATTATGGCGGTGACGTTTACCAACAAAGCGGCGGCGGAGATGCGTCATCGTATCGGGCAACTGATg                                                                                <  1:904186/149‑1 (MQ=255)
                                                                        gcgcgggcAGTGGTAAGACGCGCGTACTGGTGCATCGTATCGCCTGGTTGATGAGCGTGGAAAACTGCTCGCCATACTCGATTATGGCGGTGACGTTTACCAACAAAGCGGCGGCGGAGATGCGTCATCGTATCGGGCAACTGATGGGc                                                                             <  1:534290/149‑1 (MQ=255)
                                                                          gcgggcAGTGGTAAGACGCGCGTACTGGTGCATCGTATCGCCTGGTGGATGAGCGTGGAAAACTGCTCGCCATACTAGATTATGGCGGTGACGTTTACCAACAAAGCGGCGGCGGAGATGCGGCATCGTATCGGGCAACTGATGGGCAc                                                                           <  2:1137079/149‑1 (MQ=255)
                                                                            gggcAGTGGTAAGACGCGCGTACTGGTGCATCGTATCGCCTGGTTGAGGAGCGTGGAAAACTGCTCGCCATACTCGATTATGGCGGTGACGTTTACCAACAAAGCGGCGGCGGAGATGCGTCATCGTATCGGGCAACTGATGGGCACGa                                                                         <  2:1344227/149‑1 (MQ=255)
                                                                                aggggTAAGACGCGCGTACTGGGGCAACGTATCGCCGGGTTGCGGCGCGTGGAAACCTGCTCGCCAAACTCGAATATGGCGGTGACGTTTACCAACAACGCGGCGGCGGAGATGCGGCATCGTATCGGGCAACTGATGGGCACGAGCCa                                                                     <  2:1027845/146‑1 (MQ=255)
                                                                                aGTGGTAAGACTCGAGTAATGGGTCATCGTATCGGCTGGGTTATTTGAGTGGAAAACTCCTCGCCATACTCTATTATGGCGGTGACGTTTACCAACAAAGCGTCGGCGTAGATGCGTCATCGTATAGGTCAACTGATGGGCACGAGCCa                                                                     <  2:1023772/149‑1 (MQ=255)
                                                                                aGTGGTAAGACGCGCGTACTGGTGCATCGTATCGCCTGGTTGATGAGCGGGGAAAACTGCTCGCCATACTCGATTATGGCGGTGACGTTTACCAACAAAGCGGCGGCGGAGATGCGTCATCGTATCGGGCAACTGATGGGCACGAGCCa                                                                     <  1:1369944/149‑1 (MQ=255)
                                                                                    gTAAGACGCGCGTACTGGTGCATCGTATCGCCTGGTTGATGAGCGTGGAAAACTGCTCGCCATACTCGATTATGGCGGTGACGTTTACCAACAAAGCGGCGGCGGAGATGCGTCATCGTATCGGGCAACTGATGGGCACGAGCCAGGGc                                                                 >  1:535255/1‑149 (MQ=255)
                                                                                       cgACGCGCGTACTGGTGCAGCGTATCGCCTGGGTGATGAGCGTGGAAAACTGCTCGCCATACTCGATTATGGCGGTGACGTTTACCAACAAAGCGGCGGCGGAGATGCGTCATCGTATCGGGCAACTGATGGGCACGAGCCAGGGCGGt                                                              <  1:419026/148‑1 (MQ=255)
                                                                                        gACGCGCGTACTGGTGCATCGTATCGCCTGGTCGATGAGCGTGGAAAACTGCTCGCCATACTCGATTCTGGCGGTGACGTTTACCAACAAAGCGGCGGCGGAGATGCGTCATCGTATCGGGCAACTGATGGGCACGAGCCAGGGCGGTa                                                             <  1:1538504/149‑1 (MQ=255)
                                                                                           gcgcgTACTGGTGCAGCGTTTCGCCTGGTTGATGAGCGTGGAAAACTGCTCGCCATACTCGATTATGGCGGTGACGTTTACCAACAAAGCGGCGGCGGAGATGCGTCATCGTATCGGGCAACTGATGGGCACGAGCCAGGGCGGTAtgt                                                          <  1:1394036/149‑1 (MQ=255)
                                                                                               gTACTGGTGCATCGTATCGCCTGGTTGATGAGCGTGGAAAACTGCTCGCCATACTCGATTATGGCGGTGACGTTTACCAACAAAGCGGCGGCGGAGATGCGTCATCGTATCGGGCAACTGATGGGCACGAGCCAGGGCGGTATGTGGGt                                                      >  2:1531874/1‑149 (MQ=255)
                                                                                                 aCTGGTGCATCGTATCGCCTGTTTGATGAGCGTGGAAAACTGCTCGCCATACTCGATTATGGCGGTGACGTTTACCAACAAAGCGGCGGCGGAGATGCGTCATCGTATCGGGCAACTGATGGGCACGAGCCAGGGCGGTATGTGGGTCg                                                    <  1:146108/149‑1 (MQ=255)
                                                                                                  cTGGTGCATCGTATCGCCTGGTTGATGAGCGTGGAAAACTGCTCGCCATACTCGATTATGGCGGTGACGTTTACCAACAAAGCGGCGGCGGAGATGCGTCATCGTATCGGGCAACTGATGGGCACGAGCCAGGGCGGTATGTGGGTCgg                                                   <  1:1419706/149‑1 (MQ=255)
                                                                                                     gTGCATCGTATCGCCTGGTTGATGAGCGTGGAAAACTGCTCGCCATACTCGATTATGGCGGTGACGTTTACCAACAAAGCGGCGGCGGAGATGCGTCATCGTATCGGGCAACTGATGGGCACGAGCCAGGGCGGTATGTGGGTCGGCAc                                                >  1:1177991/1‑149 (MQ=255)
                                                                                                      tGCATCGTATCGCCTGGTTGCTGAGCGTGGAAAACTGCTCGCCATACTCGAGTATGGCGGTGACGTTTACCAACAAAGCGGCGGCGGAGATGCGTCATCGTATCGGGCAACTGATGGGCACGAGCCAGGGCGGTATGTGGGTCGGCAcc                                               <  2:518460/149‑1 (MQ=255)
                                                                                                       gCATCGTATCGCCTGGTTGATGAGCGTGGAAAACTGCTCGCCATACTCGATTATGGCGGTGACGTTTACCAACAAAGCGGCGGCGGAGATGCGTCATCGTATCGGGCAACTGATGGGCACGAGCCAGGGCGGTATGTGGGTCGGCACCt                                              <  1:983333/149‑1 (MQ=255)
                                                                                                       gCATCGTATCGCCTGGTTGATGAGCGTGGAAAACTGCTCGCCATACTCGATTATGGCGGTGACGTTTACCAACAAAGCGGCGGCGGAGATGCGTCATCGTATCGGGCAACTGATGGGCACGAGCCAGGGCGGTATGTGGGTCGGCACCt                                              <  2:80381/149‑1 (MQ=255)
                                                                                                       gCATCGTATCGCCTGGTTGATGAGCGTGGAAAACTGCTCGCCATACTCGATTATGGCGGTGACGTTTACCAACAAAGCGGCGGCGGAGATGCGTCATCGTATCGGGCAACTGATGGGCACGAGCCAGGGCGGTATGTGGGTCGGCACCt                                              >  1:969726/1‑149 (MQ=255)
                                                                                                       gCATCGTATCGCCTGGTTGATGAGCATGGAAAACTGCTCGCCTTACTCGATTCTGGCGGTGACGTTTACCAACAAAGCGGCGGCGGAGATGCGTCATCGTATCGGGCAACTGATGGGCACGAGCCAGGGCGGTATGTGGGTCGGCACCt                                              <  2:1045332/149‑1 (MQ=255)
                                                                                                         aTCGTATCACCTGGTTGATGAGCGTGGAAAACTGCTCGCCATACTCGATTATGGCGGTGACGTTTACCAACAAAGCGGCGGCGGAGATGCGTCATCGTATCGGGCAACTGATGGGCACGAGCCAGGGCGGTATGTGGGTCGGCACCTTc                                            <  2:242167/149‑1 (MQ=255)
                                                                                                           cGTATCGCCTGGTTGATGAGCGTGGAAAACTGCTCGCCATACTCGATTATGGCGGTGACGTTTACCAACAAAGCGGCGGCGGAGATGCGTCATCGTATCGGGCAACTGATGGGCACGAGCCAGGGCGGTATGTGGGTCGGCACCTTccc                                          >  1:1468674/1‑148 (MQ=255)
                                                                                                            gTATCGCCTGGTTGATGAGCGTGGAAAACTGCTCGCCATACTCGATTATGGCGGTGACGTTTACCAACAAAGCGGCGGCGGAGATGCGTCATCGTATCGGGCAACTGATGGGCACGAGCCAGGGCGGTATGTGGGTCGGCACCTTCCAc                                         <  2:1518952/149‑1 (MQ=255)
                                                                                                              aTCGCCTGGTTGATGAGCGTGGAAAACTGCTCGCCATACTCGATTATGGCGGTGACGTTTACCAACAAAGCGGCGGCGGAGATGCGTCATCGTATCGGGCAACTGATGGGCACGAGCCAGGGCGGTATGTGGGTCGGCACCTTCCACgg                                       >  2:580637/1‑149 (MQ=255)
                                                                                                              aTCGCCTGGTTGATGAGCGTGGAAAACTGCTCGCCATACTCGATTATGGCGGTGACGTTTACCAACAAAGCGGCGGCGGAGATGCGTCAGCGTATCGGGCAACTGATGGGCACGAGCCAGGGCGGTATGTGGGTCGGCACCTTCCACgg                                       >  2:456111/1‑149 (MQ=255)
                                                                                                                cGCCGGGTTGATGAGCGTGGAAAACTGCTCGCCATACTCGATTATTGCGGTGACGTTTACCAACAAAGCGGCGGCGGAGAGGCGTCATCGTATCGGGCAACTGATGGGCACGAGCCAGGGCGGTATGTGGGTCGGCACCTTCCACGGGc                                     <  1:270678/149‑1 (MQ=255)
                                                                                                                 gCCTGGTTGATGAGCGTGGAAAACTGCTCGCCATACTCGATTATGGCGGTGACGTTTACCAACAAAGCGGCGGCGGAGATGCGTCATCGTATCGGGCACCTGATGGGCACGAGCCAGGGCGGTATGTGGGTCGGCACCTTCCAAgggag                                    >  2:1458438/1‑147 (MQ=255)
                                                                                                                 gCCTGGTTGATGAGCGTGGAAAACTGCTCGCCATACTCGATTATGGCGGTGACGTTTACCAACAAAGCGGCGGCGGAGATGCGTCATCGTATCGGGCAACTGATGGGCCCGAGCCAGGGCGTTATGTGGGTCGGCACCTTCCACGGGCt                                    >  2:754104/1‑149 (MQ=255)
                                                                                                                  ccTGGTTGATGAGCGTGGAAAACTGCTCGCCATACTCGATTATGGCGGTGACGTTTACCAACAAAGCGGCGGCGGAGATGCGTCATCGTATCGGGCAACTGATGGGCACGAGCCAGGGCGGTATGTGGGTCGGCACCTTCCACGGGCTg                                   >  1:500400/1‑149 (MQ=255)
                                                                                                                     ggTTGATGAGCGGCGAAAACTGCTCGCCATACTCGATTATGGCGGCGACGATTACCAACAAAGCGGCGGCGGAGATGCGTCATCGTATCGGGCAACTGATGGGCACGAGCCAGGGCGGTATGTGGGTCGGCACCTTCCACGGGCTGgcg                                <  2:1352010/149‑1 (MQ=255)
                                                                                                                       ttGATGAGCGTGGAAAACTGCTCGCCATACTCGATTATGGCGGTGACGTTTACCAACAAAGCGGCGGCGGAGATGCGTCATCGTATCGGGCAACTGATGGGCACGAGCCAGGGCGGTATGTGGGTCGGCACCTTCCACGGGCTGGCGCa                              <  1:63270/149‑1 (MQ=255)
                                                                                                                           tgaGCGTGGAAAACTGCTCGCCATACTCGATTATGGCGGTGACGTTTACCAACAAAGCGGCGGCGGAGATGCGTCATCGTATCGGGCAACTGATGGGCACGAGCCAGGGCGGTATGTGGGTCGGCACCTTCCACGGGCTGGCGCACCGt                          >  1:451746/1‑149 (MQ=255)
                                                                                                                            gaGCGTGGAAACCTGCGCGCCATACTCGATTCTGGCGGTGACGTATACCAACAAAGCGGCGGCGGAGATGCGTCATCGTATCGGGCAACAGATGGGCACGAGCCAGGGCGGTATGTGGGTCGGCACCTTCCACGGGCTGGCGCACCGtt                         <  2:187993/149‑1 (MQ=255)
                                                                                                                              gcgcGGAAAACTGCTCGCCATACTCGATTATGGCGGTGACGTTTACCAACAAAGCGGCGGCGGAGATGCGTCATCGTAGCGGGCAACTGATGGGCACGAGCCAGGGCGGTATGTGGGTCGGCACCTTCCACGGGCTGGCGCACCGTTtg                       <  1:633815/149‑1 (MQ=255)
                                                                                                                              gCGTTGAAAACGGCTCGCCATACTCGATTATGGCGGTGACGTTTCCCAACAACGCGGCGGCGGAGATGCGTCTTCGTATCGGGCAACTGATGGGCACGAGCCAGGGCGGTATGTGGGTCGGCACCTTCCACGGGCTGGCGCACCGTTtg                       <  2:957874/149‑1 (MQ=255)
                                                                                                                              ccGTGGAAATCTGCTCGCCATACTCGCTTATGGCGGTGACGTTAACCAACAAAGCGGCGGCGGAGAAGCGTCATCGTATCGGGCAACTGATGGGCACGAGCCCGGGCGGGCTGTGGGTCGGCACCTTCCACGGGCGGGCGCACCGGTtg                       <  1:1019645/148‑1 (MQ=255)
                                                                                                                                gTGGAAAACTGCTCGCCATACTCGATTATGGCGGTGACGTTTACCAACAAAGCGGCGGCGGAGATGCGTCATCGTATCGGGCAACTGATGGGCACGAGCCAGGGCGGTATGTGGGTCGGCACCTTCCACGGGCTGGCGCACCGTTtgct                     <  2:1351717/149‑1 (MQ=255)
                                                                                                                                 tGGAAAACTGCTCGCCATACTCGATTATGGCGGTGACGTTTACCAACAAAGCGGCGGCGGAGATGCGTCATCGTATCGGGCAACTGATTGGCACGAGCCAGGGCGGTATGTGGGTCGGCACCTTCCACGGGCTGGCGCACCGTTtgctg                    >  1:1037113/1‑149 (MQ=255)
                                                                                                                                 tGGAAAACTGCTCGCCATACTCGATTATGGCGGTGACGTTTACCAACAAAGCGGCGGCGGAGATGCGTCATAGTATCGGGCAACTGATGGGCACGAGCCATGGCGGTATGTGGGTCGGCACCTTCCACGGGCAGGCGCAACGTTtgcag                    >  2:297400/1‑147 (MQ=255)
                                                                                                                                  ggAACACTGCGAGCCATACTCGATTATGCCGGTTCCGTGTAGCAACAATGCGGCGGCGGAGATGCGGCATCGTATCGGGCAACTGATGGGCTCGAGAGAGGGCGGTATGTGGGTCGGCACCTTCCACGGGCTGGCGCACCGTTtgctgc                   <  1:886785/149‑1 (MQ=255)
                                                                                                                                   gAAAACTGCTCGCCATACTCGATTATGGCGGTGACGTTTACCAACAAAGCGGCGGCGGAGATGCGTCATCGTATCGGGCAACTGATGGGCACGAGCCAGGGCGGTATGTGGGTCGGCACCTTCCACGGGCTGGCGCACCGTTTGCtgcg                  >  1:439664/1‑149 (MQ=255)
                                                                                                                                    aaaaCTGCTCGCCATACTCGATTATGGCGGTGACGTTTACCAACAAAGCGGCGGCGGAGATGCGTCATCGTATCGGGCAACTGATGGGCACGAGCCAGGGCGGTATGTGGGTCGGCACCTTCCACGGGCTGGCGCCCCGTTTGCtgcgt                 >  1:1292043/1‑149 (MQ=255)
                                                                                                                                    aaaaCTGCTCGCCATACTCGATTATGGCGGTGACGTTTACCAACAAAGCGGCGGCGGAGATGCGTCATCGTATCGGGCAACTGATGGGCACGAGCCAGGGCGGTATGTGGGTCGGCACCTTCCACGGGCTGGCGCCACGTTTGCtgcgt                 >  1:1487374/1‑149 (MQ=255)
                                                                                                                                    aaaaCTGCTCGCCATACTCGATTATGGCGGTGACGTTTACCAACAAAGCGGCGGCGGAGATGCGTCATCGTATCGGGCAACTGATGCGCACGAGCCAGGGCGGTATGTGGGTCGGCACCTTCCACGGGCTGGCGCACCGTTTGCtgcgt                 >  1:1060787/1‑149 (MQ=255)
                                                                                                                                     aaaCTGCTCGCCATACTCGATTATGGCGGTGACGTTTACCAACAAAGCGGCGGCGGAGATGCGTCATCGTATCGGGCAACTGATGGGCACGAGCCAGGGCGGTATGTGGGTCGGCACCTTCCACGGGCTGGCGCACCGTTTGCtgcgtg                <  2:1353800/149‑1 (MQ=255)
                                                                                                                                      aaCTGCTCGCCATACTCGATTATGGCGGTGACGTTTACCAACAAAGCGGCGGCGGAGATGCGTCATCGTATCGGGCAACTGATGGGCACGAGCCAGGGCGGTATGTGGGTCGGCACCTTCCACGGGCTGGCGCACCGTTTGCtgcgtgc               <  1:1416320/149‑1 (MQ=255)
                                                                                                                                       cgtcaccGCCATAATCGAGTATGGCGGTGACGTTTACCAACAAAGCGGCGGCGGAGATGCGTCATCGTATCGGGCAACTGATGGGCACGAGCCAGGGCGGTATGTGGGTCGGCACCTTCCCCGGGCTGGCGCCCCGTTTGCtgcgggcg              >  1:1062577/7‑149 (MQ=255)
                                                                                                                                       aCTGCTCGCCATACTCGCTTATGGCGGTGACGTTTACCAACAACGCGGCGGCGGAGATGCGGCATCGTATCGGGCACCTGATGGGCACGAGCCAGGGCGGTATGTGGGTCGGCACCTTCCACGGGCTGGCGCACCGTTTGCtgcgtgcg              <  1:1395933/149‑1 (MQ=255)
                                                                                                                                        cTGCTCGCCATACTCGATTCTGGCGGTGACGTTTACCCACAAAGCGGCGGCGGAGATGCGTCATCGTATCGGGCAACTGATGGGCACGAGCCAGGGCGGTATGTGGGTCGGCACCTTCCACGGGCTGGCGCACCGTTTGCTGCGTgcgc             <  2:303981/149‑1 (MQ=255)
                                                                                                                                         cgCTCGCCATACTCGCTTATGGCGGTGACGTTTACCAACAAAGCGGCGGCGGAGATGCGGCATCGTATCGGGCAACTGATGGGCACGAGCCAGGGCGGTATGTGGGGCGGCACCTTCCACGGGCTGGCGCACCGTTTGCTGCGTGCGCa            <  1:294007/148‑1 (MQ=255)
                                                                                                                                           cTCGCCATACTCGATTATGGCGGTGACGTTTACCAACAAAGCGGCGGCGGAGATGCGTCATCGTATCGGGCAACTGATGGGCACGAGCCAGGGCGGTATGTGGGTCGGCACCTTCCACGGGCTGGCGCACCGTTTGCTGCGTGCGCAcc          >  1:1230126/1‑149 (MQ=255)
                                                                                                                                             cGCCATACTCGATTTTGGCGGTGACGTATTCCAACAAAGCGGCGGCGGAGATGCGTCATCGTCTCGGGCAACTGATGGGCACGAGCCAGGGCGGTATGTGGGTCGGCACCTTCCACGGGCTGGCGCACCGTTTGTTGCGTGCGCACCat        <  1:174271/149‑1 (MQ=255)
                                                                                                                                             cGCCATACTCGATTATGGCGGTGACGTTTACCAACAAAGCGGCGGCGGAGATGCGTCATCGTATCGGGCAACTGATGGGCACGAGCCAGGGCGGTATGTGGGTCGGCACCTTCCACGGGCTGGCGCACCGTTTGCTGCGTGCGCACCat        <  2:1177991/149‑1 (MQ=255)
                                                                                                                                               ccATACTCGATTATGGCGGTGACGTTTACCAACAAAGCGGCGGCGGAGATGCGTCATCGTATCGGGCAACTGATGGGCACGAGCCAGGGCGGTATGTGGGTCGGCACCTTCCACGGGCTGGCGCACCGTTTGCTGCGTGCGCACCatat      >  1:835512/1‑149 (MQ=255)
                                                                                                                                                cATACTCGATTATGGCGGTGACGTTTACCAACAAAGCGGCGGCGGAGATGCGTCATCGTATCGGGCAACTGATGGGCACGAGCCAGGGCGGTATGTGGGTCGGCACCTTCCACGGGCTGGCGCACCGTTTGCTGCGTGCGCCCCATATg     >  1:141400/1‑149 (MQ=255)
                                                                                                                                                cATACTCGATTATGGCGGTGACGTTTACCAACAAAGCGGCGGCGGAGATGCGTCATCGTATCGGGCAACTGATGGGCACGAGCCAGGGCGGTATGTGGGTCGGCACCTTCCACGGGCTGGCGCACCGTTTGCTGCGTGCGCACCATATg     <  2:451746/149‑1 (MQ=255)
                                                                                                                                                   aCTCGATTATGGCGGTGACGTTTACCAACAAAGCGGCGGCGGAGATGCGTCATCGTATCGGGCAACTGATGGGCACGAGCCAGGGCGGTATGTGGGTCGGCACCTTCCACGGGCTGCCGCACCGTTTGCTGCGTGCGCACCATATGGAc  >  1:1467200/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
TGCTCGACAGCCTTAATGACAAACAGCGCGAAGCGGTGGCCGCGCCACGCAGCAACCTTCTGGTGCTGGCGGGCGCGGGCAGTGGTAAGACGCGCGTACTGGTGCATCGTATCGCCTGGTTGATGAGCGTGGAAAACTGCTCGCCATACTCGATTATGGCGGTGACGTTTACCAACAAAGCGGCGGCGGAGATGCGTCATCGTATCGGGCAACTGATGGGCACGAGCCAGGGCGGTATGTGGGTCGGCACCTTCCACGGGCTGGCGCACCGTTTGCTGCGTGCGCACCATATGGAC  >  NC_000913/3997999‑3998294

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: