Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,762,221 G→T 5.9% I440I (ATC→ATA) ‡ yfjK ← radiation resistance protein; DEAD/H helicase‑like protein; CP4‑57 putative defective prophage

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,762,2210GT 5.9% 460.5 / 12.8 133I440I (ATC→ATA) ‡yfjKradiation resistance protein; DEAD/H helicase‑like protein; CP4‑57 putative defective prophage
Reads supporting (aligned to +/- strand):  ref base G (61/64);  new base T (2/6);  total (63/70)
Fisher's exact test for biased strand distribution p-value = 2.80e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.93e-01

Only 100 of 133 total aligned reads displayed.

CTCATCAAAGCGATCCTGTGAAAACTCGGACAAGTGGTCTGGTTGAACGCTGGCTAACAGGTTGATTGGGGTATCAATGCCCTGAATACCGAGCGGAATATCCACCTCTTGATTGAGATTGTCTTCCGGGATTTCCTCCAGGCAAAAGATTTTGCCGACGAAGTGAACGCCCATGCGCCCGGCGCGACCTTTGATGTTCCCATGGGTGAATTTGTCAATGCTGCGAGTTCCATCGCGGTTGTCGTAGATCACAACGTTCTTAGCGATAGTGTTGACTCCTTCAATGATCGTGGAG  >  NC_000913/2762074‑2762368
                                                                                                                                                   |                                                                                                                                                   
cTCATCAAAGCGATCCTGTGAAAACTCGGACAAGTGGTCTGGTTGAACGCTGGCTAACAGGTTGATTGGGGTATCAATGCCCTGAATACCGAGCGGAATATCCACCTCTTGATTGAGATTGTCTTCCGGGATTTCCTCCAGGCAAAAGa                                                                                                                                                    >  3:904375/1‑149 (MQ=255)
 tcatcaAAGCGATCCTGTGAAAACTCGGACAAGTGGTCTGGTTGAACGCTGGCTAACAGGTTGATTGGGGTATCAATGCCCTGAATACCGAGCGGAATATCCACCTCTTGATTGAGATTGTCTTCCGGGATTTCCTCCAGGCAAAAGAt                                                                                                                                                   >  4:875562/1‑149 (MQ=255)
  catcaAAGCGATCCTGTGAAAACTCGGACAAGTGGTCTGGTTGAACGCTGGCTAACAGGTTGATTGGGGTATCAATGCCCTGAATACCGAGCGGAATATCCACCTCTTGATTGAGATTGTGTTCCGGGATTTCCTCCAGGCAAAAGAtt                                                                                                                                                  <  2:729394/149‑1 (MQ=255)
    tcaAAGCGATCCTGTGAAAACTCGGACAAGTGGTCTGGTTGAACGCTGGCTAACAGGTTGATTGGGGTATCAATGCCCTGAATACCGAGCGGAATATCCACCTCTTGATTGAGATTGTCTTCCGGGATTTCCTCCAGGCAAAAGAtttt                                                                                                                                                >  3:337839/1‑149 (MQ=255)
    tcaAAGCGATCCTGTGAAAACTCGGACAAGTGGTCTGGTTGAACGCTGGCTAACAGGTTGATTGGGGTATCAATGCCCTGAATACCGAGCGGAATATCCACCTCTTGATTGAGATTGTCTTCCGGGATTTCCTCCAGGCAAAAGAtttt                                                                                                                                                <  4:710361/149‑1 (MQ=255)
    tcaAAGCGATCCTGTGAAAACTCGGACAAGTGGTCTGGTTGAACGCTGGCTAACAGGTTGATTGGGGTATCAATGCCCTGAATACCGAGCGGAATATCCACCTCTTGATTGAGATTGTCTTCCGGGATTTCCTCCAGGCAAAAGAtttt                                                                                                                                                <  4:715539/149‑1 (MQ=255)
    tcaAAGCGATCCTGTGAAAACTCGGACAAGTGGTCTGGTTGAACGCTGGCTAACAGGTTGATTGGGGTATCAATGCCCTGAATACCGAGCGGAATATCCACCTCTTGATTGAGATTGTCTTCCGGGATTTCCTCCAGGCAAAAGAtttt                                                                                                                                                <  3:827476/149‑1 (MQ=255)
      aaaGCGATCCTGTGAAAACTCGGACAAGTGGTCTGGTTGAACGCTGGCTAACAGGTTGATTGGGGTATCAATGCCCTGAATACCGAGCGGAATATCCACCTCTTGATTGAGATTGTCTTCCGGGATTTCCTCCAGGCAAAATCTTTTGc                                                                                                                                              >  1:640205/1‑149 (MQ=255)
      aaaGCGATCCTGTGAAAACTCGGACAAGTGGTCTGGTTGAACGCTGGCTAACAGGTTGATTGGGGTATCAATGCCCTGAATACCGAGCGGAATATCCACCTCTTGATTGAGATTGTCTTCCGGGATTTCCTCCAGGCAAAAGATTTTGc                                                                                                                                              >  3:183886/1‑149 (MQ=255)
      aaaGCGATCCTGTGAAAACTCGGACAAGTGGTCTGGTTGAACGCTGGCTAACAGGTTGATTGGGGTATCAATGCCCTGAATACCGAGCGGAATATCCACCTCTTGATTGAGATTGTCTTCCGGGATTTCCTCCAGGCAAAAGATTTTGc                                                                                                                                              <  4:486895/149‑1 (MQ=255)
      aaaGCGATCCTGTGAAAACTCGGACAAGTGGTCTGGTTCAACGCTGGCTAACAGGTTGATTGGGGTATCAATGCCCTGAATACCGAGCGGAATATCCACCTCTTGATTGAGATTGTCTTCCGGGATTTCCTCCAGGCAAAAGATTTTGc                                                                                                                                              >  4:891848/1‑149 (MQ=255)
        aGCGATCCTGTGAAAACTCGGACAAGTGGTCTGGTTGAACGCTGGCTAACAGGTTGATTGGGGTATCAATGCCCTGAATACCGAGCGGAATATCCACCTCTTGATTGAGATTGTCTTCCGGGATTTCCTCCAGGCAAAATCTTTTGcct                                                                                                                                            <  1:250470/149‑2 (MQ=255)
        aGCGATCCTGTGAAAACTCGGACAAGTGGTCTGGTTGAACGCTGGCTAACAGGTTGATTGGGGTATCAATGCCCTGAATACCGAGCGGAATATCCACCTCTTGATTGAGATTGTCTTCCGGGATTTCCTCCAGGCAAAAGATTTTGCcg                                                                                                                                            >  2:1010890/1‑149 (MQ=255)
         gCGATCCTGTGAAAACTCGGACAAGTGGTCTGGTTGAACGCTGGCTAACAGGTTGATTGGGGTATCAATGCCCTGAATACCGAGCGGAATATCCACCTCTTGATTGAGATTGTCTTCCGGGATTTCCTCCAGGCAAAATCTTTTGcctg                                                                                                                                           <  4:900403/149‑3 (MQ=255)
          cGATCCTGTGAAAACTCGGACAAGTGGTCTGGTTGAACGCTGGCTAACAGGTTGATTGGGGTATCAATGCCCTGAATACCGAGCGGAATATCCACCTCTTGATTGAGATTGTCTTCCGGGATTTCCTCCAGGCAAAATCTTTTGcctgg                                                                                                                                          <  4:350001/149‑4 (MQ=255)
           gATCCTGTGAAAACTCGGACAAGTGGTCTGGTTGAACGCTGGCTAACAGGTTGATTGGGGTATCAATGCCCTGAATACCGAGCGGAATATCCACCTCTTGATTGAGATTGTCTTCCGGGATTTCCTCCAGGCAAAATCTTTTGcctgga                                                                                                                                         <  2:640205/149‑5 (MQ=255)
            aTCCTGTGAAAACTCGGACAAGTGGTCTGGTTGAACGCTGGCTAACAGGTTGATTGGGGTATCAATGCCCTGAATACCGAGCGGAATATCCACCTCTTGATTGAGATTGTCTTCCGGGATTTCCTCCAGGCAAAAGATTTTGCcgacga                                                                                                                                        <  1:196749/149‑1 (MQ=255)
             tCCTGTGAAAACTCGGACAAGTGGTCTGGTTGAACGCTGGCTAACAGGTTGATTGGGGTATCAATGCCCTGAATACCGAGCGGAATATCCACCTCTTGATTGAGATTGTCTTCCGGGATTTCCTCCAGGCAAAAGATTTTGCCGACGaa                                                                                                                                       <  2:959750/149‑1 (MQ=255)
             tCCTGTGAAAACTCGGACAAGTGGTCTGGTTGAACGCTGGCTAACAGGTTGATTGGGGTATCAATGCCCTGAATACCGAGCGGAATATCCACCTCTTGATTGAGATTGTCTTCCGGGATTTCCTCCAGGCAAAAGATTTTGCCGACGaa                                                                                                                                       <  4:885339/149‑1 (MQ=255)
               cTGTGAAAACTCGGACAAGTGGTCTGGTTGAACGCTGGCTAACAGGTTGATTGGGGTATCAATGCCCTGAATACCGAGCGGAATATCCACCTCTTGATTGAGATTGTCTTCCGGGATTTCCTCCAGGCAAAAGATTTTGCCGACGAAGt                                                                                                                                     <  3:544512/149‑1 (MQ=255)
               cTGTGAAAACTCGGACAAGTGGTCTGGTTGAACGCTGGCTAACAGGTTGATTGGGGTATCAATGCCCTGAATACCGAGCGGAATATCCACCTCTTGATTGAGATTGTCTTCCGGGATTTCCTCCAGGCAAAAGATTTTGCCGACGAAGt                                                                                                                                     >  2:782037/1‑149 (MQ=255)
                       aCTCGGACAAGTGGTCTGGTTGAACGCTGGCTAACAGGTTGATTGGGGTATCAATGCCCTGAATACCGAGCGGAATATCCAACTCTTGATTGAGATTGTCTTCCGGGATTTCCTCCAGGCAAAAGATTTTGCCGACGAAGTGAACGccc                                                                                                                             <  2:1034305/149‑1 (MQ=255)
                         tCGGACAAGTGGTCTGGTTGAACGCTGGCTAACAGGTTGATTGGGGTATCAATGCCCTGAATACCGAGCGGAATATCCACCTCTTGATTGAGATTGTCTTCCGGGATTTCCTCCAGGCAAAAGATTTTGCCGACGAAGTGAACGCCCAt                                                                                                                           <  2:415322/149‑1 (MQ=255)
                         tCGGACAAGTGGTCTGGTTGAACGCTGGCTAACAGGTTGATTGGGGTATCAATGCCCTGAATACCGAGCGGAATATCCACCTCTTGATTGAGATTGTCTTCCGGGATTTCCTCCAGGCAAAAGATTTTGCCGACGAAGTGAACGCCCAt                                                                                                                           <  1:172440/149‑1 (MQ=255)
                         tCGGACAAGTGGTCTGGTTGAACGCTGGCTAACAGGTTGATTGGGGTATCAATGCCCTGAATACCGAGCGGAATATCCACCTCTTGATTGAGATTGTCTTCCGGGATTTCCTCCAGGCAAAAGATTTTGCCGACGAAGTGAACGCCCAt                                                                                                                           >  3:149152/1‑149 (MQ=255)
                               aaGTGGTCTGGTTGAACGCTGGCTAACAGGTTGATTGGGGTATCAATGCCCTGAATACCGAGCGGAATATCCACCTCTTGATTGAGATTGTCTTCCGGGATTTCCTCCAGGCAAAAGATTTTGCCGACGAAGTGAACGCCCATGCGccc                                                                                                                     <  4:173992/149‑1 (MQ=255)
                                  tGGTCTGGTTGAACGCTGGCTAACAGGTTGATTGGGGTATCAATGCCCTGAATACCGAGCGGAATATCCACCTCTTGCTTGAGATTGTCTTCCGGGATTTCCTCCAGGCAAAAGATTTTGCCGACGACGTGAACGCCCATGCGCCCGgc                                                                                                                  >  4:642581/1‑149 (MQ=255)
                                    gTCTGGTTGAACGCTGGCTAACAGGTTGATTGGGGTATCAATGCCCTGAATACCGAGCGGAATATCCACCTCTTGATTGAGATTGTCTTCCGGGATTTCCTCCAGGCAAAAGATTTTGCCGACGAAGTGAACGCCCATGCGCCCGgcgc                                                                                                                >  3:663309/1‑149 (MQ=255)
                                       tGGTTGAACGCTGGCTAACAGGTTGATTGGGGTATCAATGCCCTGAATACCGAGCGGAATATCCACCTCTTGATTGAGATTGTCTTCCGGGATTTCCTCCAGGCAAAAGATTTTGCCGACGAAGTGAACGCCCATGCGCCCGGCGCGAc                                                                                                             >  1:88802/1‑149 (MQ=255)
                                       tGGTTGAACGCTGGCTAACAGGTTGATTGGGGTATCAATGCCCTGAATACCGAGCGGAATATCCACCTCTTGATTGAGATTGTCTTCCGGGATTTCCTCCAGGCAAAAGATTTTGCCGACGAAGTGAACGCCCATGCGCCCGGCGCGAc                                                                                                             <  2:318688/149‑1 (MQ=255)
                                         gTTGAACGCTGGCTAACAGGTTGATTGGGGTATCAATGCCCTGAATACCGAGCGGAATATCCACCTCTTGATTGAGATTGTCTTCCGGGATTTCCTCCAGGCAAAAGATTTTGCCGACGAAGTGAACGCCCATGCGCCCGGCGCGACCt                                                                                                           >  3:312453/1‑149 (MQ=255)
                                          ttGAACGCTGGCTAACAGGTTGATTGGGGTATCAATGCCCTGAATACCGAGCGGAATATCCACCTCTTGATTGAGATTGTCTTCCGGGATTTCCTCCAGGCAAAAGATTTTGCCGACGAAGTGAACGCCCATGCGCCCGGCGCGACCtt                                                                                                          >  3:24748/1‑149 (MQ=255)
                                             aaCGCTGGCTAACAGGTTGATTGGGGTATCAATGCCCTGAATACCGAGCGGAATATCCACCTCTTGATTGAGATTGTCTTCCGGGATTTCCTCCAGGCAAAAGATTTTGCCGACGAAGTGAACGCCCATGCGCCCGGCGCGACCTTTGa                                                                                                       <  4:321270/149‑1 (MQ=255)
                                                 cTGGCTAACAGGTTGATTGGGGTATCAATGCCCTGAATACCGAGCGGAATATCCACCTCTTGATTGAGATTGTCTTCCGGGATTTCCTCCAGGCAAAAGATTTTGCCGACGAAGTGAACGCCCATGCGCCCGGCGCGACCTTTGATGtt                                                                                                   <  3:803140/149‑1 (MQ=255)
                                                   ggCTAACAGGTTGATTGGGGTATCAATGCCCTGAATACCGAGCGGAATATCCACCTCTTGATTGAGATTGTCTTCCGGGATTTCCTCCAGGCAAAAGATTTTGCCGACGAAGTGAACGCCCATGCGCCCGGCGCGACCTTTGATGTTcc                                                                                                 <  1:41623/149‑1 (MQ=255)
                                                    gCTAACAGGTTGATTGGGGTATCAATGCCCTGAATACCGAGCGGAATATCCACCTCTTGATTGAGATTGTCTTCCGGGATTTCCTCCAGGCAAAAGATTTTGCCGACGAAGTGAACGCCCATGCGCCCGGCGCGACCTTTGATGTTccc                                                                                                <  3:191500/149‑1 (MQ=255)
                                                      tAACAGGTTGATTGGGGTATCAATGCCCTGAATACCGAGCGGAATATCCACCTCTTGATTGAGATTGTCTTCCGGGATTTCCTCCAGGCAAAAGATTTTGCCGACGAAGTGAACGCCCATGCGCCCGGCGCGACCTTTGATGTTCCCa                                                                                               <  2:1049772/148‑1 (MQ=255)
                                                      tAACAGGTTGATTGGGGTATCAATGCCCTGAATACCGAGCGGAATATCCACCTCTTGATTGAGATTGTCTTCCGGGATTTCCTCCAGGCAAAAGATTTTGCCGACGAAGTGAACGCCCATGCGCCCGGCGCGACCTTTGATGTTCCCa                                                                                               >  1:1049772/1‑148 (MQ=255)
                                                      tAACAGGTTGATTGGGGTATCAATGCCCTGAATACCGAGCGGAATATCCACCTCTTGATTGAGATTGTCTTCCGGGATTTCCTCCAGGCAAAAGATTTTGCCGACGAAGTGAACGCCCATGCGCCCGGCGCGACCTTTGATGTTCCCAt                                                                                              >  2:81052/1‑149 (MQ=255)
                                                       aaCAGGTTGATTGGGGTATCAATGCCCTGAATACCGAGCGGAATATCCACCTCTTGATTGAGATTGTCTTCCGGGATTTCCTCCAGGCAAAAGATTTTGCCGACGAAGTGa                                                                                                                                   >  3:470828/1‑111 (MQ=255)
                                                       aaCAGGTTGATTGGGGTATCAATGCCCTGAATACCGAGCGGAATATCCACCTCTTGATTGAGATTGTCTTCCGGGATTTCCTCCAGGCAAAAGATTTTGCCGACGAAGTGAACGCCCATGCGCCCGGCGCGACCTTTGATGTTCCCATg                                                                                             <  2:88802/149‑1 (MQ=255)
                                                       aaCAGGTTGATTGGGGTATCAATGCCCTGAATACCGAGCGGAATATCCACCTCTTGATTGAGATTGTCTTCCGGGATTTCCTCCAGGCAAAAGATTTTGCCGACGAAGTGAACGCCCATGCGCCCGGCGCGACCTTTGATGTTCCCATg                                                                                             >  1:130444/1‑149 (MQ=255)
                                                        aCAGGTTGATTGGGGTATCAATGCCCTGAATACCGAGCGGAATATCCACCTCTTGATTGAGATTGTCTTCCGGGATTTCCTCCAGGCAAAAGATTTTGcc                                                                                                                                             >  4:252607/1‑100 (MQ=255)
                                                        aCAGGTTGATTGGGGTATCAATGCCCTGAATACCGAGCGGAATATCCACCTCTTGATTGAGATTGTCTTCCGGGATTTCCTCCAGGCAAAAGATTTTGcc                                                                                                                                             <  3:252607/100‑1 (MQ=255)
                                                        aCAGGTTGATTGGGGTATCAATGCCCTGAATACCGAGCGGAATATCCACCTCTTGATTGAGATTGTCTTCCGGGATTTCCTCCAGGCAAAAGATTTTGCCGACGAAGTGAACGCCCATGCGCCCGGCGCGACCTTTGATGTTCCCATgg                                                                                            <  2:808150/149‑1 (MQ=255)
                                                        aCAGGTTGATTGGGGTATCAATGCCCTGAATACCGAGCGGAATATCCACCTCGTGATTGAGATTGTCTTCCGGGATTTCCTCCACGCAAAAGACTTTGCCGACGAAGTGAACGCCCATGCGCCCGGCGCGACCTTTGATGTTCCCATgg                                                                                            >  2:1034868/1‑149 (MQ=255)
                                                            gTTGATTGGGGTATCAATGCCCTGAATACCGAGCGGAATATCCACCTCTTGATTGAGATTGTCTTCCGGGATTTCCTCCAGGCAAAAGATTTTGCCGACGAAGTGAACGCCCATGCGCCCGGCGCGACCTTTGATGTTCCCATGGGTGa                                                                                        >  1:700354/1‑149 (MQ=255)
                                                                aTTGGGGTATCAATGCCCTGAATACCGAGCGGAATATCCACCTCTTGATTGAGATTGTCTTCCGGGATTTCCTCCAGGCAAAAGATTTTGCCGACGAAGTGAACGCCCATGCGCCCGGCGCGACCTTTGATGTTCCCATGGGTGAAttt                                                                                    <  1:914108/149‑1 (MQ=255)
                                                                aTTGGGGAATCAATGCCCTGAATACCGAGCGGAATATCCACCTCTTGATTGAGATTGTCTTCCGGGATTTCCTCCAGGCAAAAGATTTTGCCGACGAAGTGAACGCCCATGCGCCCGGCGCGACCTTTGATGTTCCCATGGGTGAAttt                                                                                    >  2:10831/1‑149 (MQ=255)
                                                                 ttGGGGTATCAATGCCCTGAATACCGAGCGGAATATCCACCTCTTGATTGAGATTGTCTTCCGGGATTTCCTCCAGGCAAAAGATTTTGCCGACGAAGTGAACGCCCATGCGCCCGGCGCGACCTTTGATGTTCCCATGGGTGAATTTg                                                                                   >  4:767042/1‑149 (MQ=255)
                                                                           aaTGCCCTGAATACCGAGCGGAATATCCACCTCTTGATTGAGATTGTCTTCCGGGATTTCCTCCAGGCAAAAGATTTTGCCGACGAAGTGAACGCCCATGCGCCCGGCGCGACCTTTGATGTTCCCATGGGTGAATTTGTCAAtgctgc                                                                         <  4:183886/149‑1 (MQ=255)
                                                                           aaTGCCCTGAATACCGAGCGGAATATCCACCTCTTGATTGAGATTGTCTTCCGGGATTTCCTCCAGGCAAAAGATTTTGCCGACGAAGTGAACGCCCATGCGCCCGGCGCGACCTTTGATGTTCCCATGGGTGAATTTGGCAAtgctgc                                                                         >  1:440190/1‑149 (MQ=255)
                                                                            aTGCCCTGAATACCGAGCGGAATATCCACCTCTTGATTGAGATTGTCTTCCGGGATTTCCTCCAGGCAAAAGATTTTGCCGACGAAGTGAACGCCCATgt                                                                                                                         >  3:927112/1‑99 (MQ=255)
                                                                             tGCCCTGAATACCGAGCGGAATATCCACCTCTTGATTGAGATTGTCTTCCGGGATTTCCTCCAGGCAAAAGATTTTGCCGACGAAGTGAACGCCCATGCGCCCGGCGCGACCTTTGATGTTCCCATGGGTGAATTTGTCAATGCTGCGa                                                                       <  2:960601/149‑1 (MQ=255)
                                                                               cccTGAATACCGAGCGGAATATCCACCTCTTGATTGAGATTGTCTTCCGGGATTTCCTCCAGGCAAAAGATTTTGCCGACGAAGTGAACGCCCATGCGCCCGGCGCGACCTTTGATGTTCCCATGGGTGAATTTGTCAATGCTGCGAGt                                                                     >  3:904936/1‑149 (MQ=255)
                                                                                ccTGAATACCGAGCGGAATATCCACCTCTTGATTGAGATTGTCTTCCGGGATTTCCTCCAGGCAAAATCTTTTGcctgga                                                                                                                                         >  1:1062655/1‑76 (MQ=255)
                                                                                ccTGAATACCGAGCGGAATATCCACCTCTTGATTGAGATTGTCTTCCGGGATTTCCTCCAGGCAAAATCTTTTGcctgga                                                                                                                                         <  2:1062655/80‑5 (MQ=255)
                                                                                  tGAATACCGAGCGGAATATCCACCTCTTGATTGAGATTGTCTTCCGGGATTTCCTCCAGGCAAAAGATTTTGCCGACGAAGTg                                                                                                                                    >  2:159299/1‑83 (MQ=255)
                                                                                           aGCGGAATATCCACCTCTTGATTGAGATTGTCTTCCGGGATTTCCTCCAGGCAAAAGATTTTGCCGACGAAGTGAACGCCCATGCGCCCGGCGCGACCTTTGATGTTCCCATGGGTGAATTTGTCAATGCTGCGAGTTCCATCGCGGtt                                                         <  3:1039923/149‑1 (MQ=255)
                                                                                           aGCGGAATATCCACCTCTTGATTGAGATTGTCTTCCGGGATTTCCTCCAGGCAAAAGATTTTGCCGACGAAGTGAACGCCCATGCGCCCGGCGCGACCTTTGATGTTCCCATGGGTGAATTTGTCAATGCTGCGAGTTCCATCGCGGtt                                                         <  1:136177/149‑1 (MQ=255)
                                                                                           aGCGGAATATCCACCTCTTGATTGAGATTGTCTTCCGGGATTTCCTCCAGGCAAAAGATTTTGCCGACGAAGTGAACGCCCATGCGCCCGGCGCGACCTTTGATGTTCCCATGGGTGAATTTGTCAATGCTGCGAGTTCCATCGCGGtt                                                         <  3:222413/149‑1 (MQ=255)
                                                                                            gCGGGATATCCACCTCTTGATTGAGATTGTCTACCATGATTTCCTCCAGGCAAGAGATTTTGCCGACGAAGTGAACGCCCATGCGCCTGGCCCGACCTTTGATGTTCCCATGGGTGAATTTGTCAATGCTGCGCGTTCCATCGCggacg                                                        >  4:581997/1‑146 (MQ=255)
                                                                                               gAATATCCACCTCTTGATTGAGATTGTCTTCCGGGATTTCCTCCAGGCAAAAGATTTTGCCGACGAAGTGAACGCCCATGCGCCCGGCGCGACCTTTGATGTTCCCATGGGTGAATTTGTCAATGCTGCGAGTTCCATCGCGGTTGTCg                                                     <  3:636273/149‑1 (MQ=255)
                                                                                                 atatCCACCTCTTGATTGAGATTGTCTTCCGGGATTTCCTCCAGGCAAAAGATTTTGCCGACGAAGTGAACGCCCATGCGCCCGGCGCGACCTTTGATGTTCCCATGGGTGAATTTGTCAATGCTGCGAGTTCCATCGCGGTTGTCGTa                                                   <  1:782037/149‑1 (MQ=255)
                                                                                                    cccACCTCTTGATTGAGATTGTCTTCCGGGATTTCCTCAAGTCAAAAGATTTTGCGGACGAAGTGAACGACCATGCGCCCGGAGCGACCTTTGAGGTTCCCATTAGTGAATTTGTCAATGCTGCGAGTTCCATCGCGGTTGTCGTAGAt                                                <  2:941067/148‑1 (MQ=255)
                                                                                                       aCCTCTTGATTGAGATTGTCTTCCGGGATTTCCTCCAGGCAAAAGATTTTGCCGACGAAGTGAAc                                                                                                                                 >  3:236716/1‑65 (MQ=255)
                                                                                                       aCCTCTTGATTGAGATTGTCTTCCGGGATTTCCTCCAGGCAAAAGATTTTGCCGACGAAGTGAAc                                                                                                                                 <  4:236716/65‑1 (MQ=255)
                                                                                                       aCCTCTTGATTGAGATTGTCTTCCGGGATTTCCTCCAGGCAAAAGATTTTGCCGACGAAGTGAACGCCCATGCGCCCGGCGCGACCTTTGATGTTCCCATGGGTGAAttt                                                                                    >  1:1049422/1‑110 (MQ=255)
                                                                                                        ccTCTTGATTGAGATTGTCTTCCGGGATTTCCTCCAGGCAAAAGATTTTGCCGACGAAGTGAACGCCCATGCGCCCGGCGCGACCTTTGATGTTCCCATGGGTGAATTTGTCAATGCTGCGAGTTCCATCGCGGTTGTCGTAGATcaca                                            >  2:11121/1‑149 (MQ=255)
                                                                                                         ctctTGATTGAGATTGTCTTCCGGGATTTCCTCCAGGCAAAAGATTTTGCCGACGAAGTGAACGCCCATGCGCCCGGCGCGACCTTTGATGTTCCCATGGGTGAATTTGTCAATGCTGCGAGTTCCATCGCGGTTGTCGTAGATCACaa                                           <  1:619349/149‑1 (MQ=255)
                                                                                                          tctTGATTGAGATTGTCTTCCGGGATTTCCTCCAGGCAAAAGATTTTGCCGACGAAGTGAACGCCCATGCGCCCGGCGCGACCTTTg                                                                                                        >  3:1006764/1‑87 (MQ=255)
                                                                                                          tctTGATTGAGATTGTCTTCCGGGATTTCCTCCAGGCAAAAGATTTTGCCGACGAAGTGAACGCCCATGCGCCCGGCGCGACCTTTg                                                                                                        <  4:1006764/87‑1 (MQ=255)
                                                                                                            ttgattgaGATTGTCTTCCGGGATTTCCTCCAGGCAAAAGATTTTGCCGACGAAGTGAACGCCCATGCGCCCGGCGCGACCTTTGATGTTCCCATGGGTGAATTTGTCAATGCTGCGAGTTCCATCGCGGTTGTCGTAGATCACAACGt                                        <  4:372559/149‑1 (MQ=255)
                                                                                                              gattgaGATTGTCTTCCGGGATTTCCTCCAGGCAAAAGATTTTGCCGACGAAGTGAACGCCCATGCGCCCGGCGCGACCTTTGATGTTCCCATGGGTGAATTTGTCAATGCTGCGAGTTCCATCGCGGTTGTCGTAGATCACAACGTTc                                      <  4:312453/149‑1 (MQ=255)
                                                                                                               attgaGATTGTCTTCCGGGATTTCCTCCAGGCAAAAGATTTTGCCGACGAAGTGAACGCCCATGCGCCCGGCGCGACCTTTGATGTTCCCATGGGTGAATTTGTCAATGCTGCGAGTTCCATcgcg                                                            <  1:830211/126‑1 (MQ=255)
                                                                                                               attgaGATTGTCTTCCGGGATTTCCTCCAGGCAAAAGATTTTGCCGACGAAGTGAACGCCCATGCGCCCGGCGCGACCTTTGATGTTCCCATGGGTGAATTTGTCAATGCTGCGAGTTCCATcgcg                                                            >  2:830211/1‑126 (MQ=255)
                                                                                                               attgaGATTGTCTTCCGGGATTTCCTCCAGGCAAAAGATTTTGCCGACGAAGTGAACGCCCATGCGCCCGGCGCGACCTTTGATGTTCCCATGGGTGAATTTGTCAATGCTGCGAGTTCCATCGCGGTTGTCGTAGATCACAACGTTCt                                     >  4:466984/1‑149 (MQ=255)
                                                                                                               attgaGATTGTCTTCCGGGATTTCCTCCAGGCAAAAGATTTTGCCGACGAAGTGAACGCCCATGCGCCCGGCGCGACCTTTGATGTTCCCATGGGTGAATTTGTCAATGCTGCGAGTTCCATCGCGGTTGTCGTAGATCACAACGTTCt                                     <  2:652744/149‑1 (MQ=255)
                                                                                                               attgaGATTGTCTTCCGGGATTTCCTCCAGGCAAAAGATTTTGCCGACGAAGTGAACGCCCATGCGCCCGGCGCGACCTTTGATGTTCCCATGGGTGAATTTGTCAATGCTGCGAGTTCCATCGCGGTTGTCGTAGATCACAACGTTCt                                     >  3:325417/1‑149 (MQ=255)
                                                                                                                ttgaGATTGTCTTCCGGGATTTCCTCCAGGCAAAAGATTTTGCCGACGAAGTGAACGCCCATGCGCCCGGCGCGACCTTTGATGTTCCCATGGGTGAATTTGTCAATGCTGCGAGTTCCATCGCGGTTGTCGTAGATCACAACGTTCtt                                    >  3:435136/1‑149 (MQ=255)
                                                                                                                 tgaGATTGTCTTCCGGGATTTCCTCCAGGCAAAAGATTTTGCCGACGAAGTGAACGCCCATGCGCCCGGCGCGACCTTTGATGTTCCCATGGGTGAATTTGTCAATGCTGCGAGTTCCATCGCGGTTGTCGTAGATCACAACGTTCTTa                                   <  3:891848/149‑1 (MQ=255)
                                                                                                                           ttCCGGGATTTCCTCCAGGCAAAAGATTTTGCCGACGAAGTGAACGCCCATGCGTCCGGCGCGACCTTTGATGTTCCCATGGGTGAATTTGTCAATGCTGCGAGTTCCATCGCGGTTGTCGTAGATCACAACGTTCTTAGCGATAgtgt                         >  1:710209/1‑149 (MQ=255)
                                                                                                                            tCCGGGATTTCCTCCAGGCAAAAGATTTTGCCGACGAAGTGAACGCCCATGCGCCCGGCGCGACCTTTGATGTTCCCATGGGTGAATTTGTCAATGCTGCGa                                                                       >  4:162229/1‑102 (MQ=255)
                                                                                                                            tCCGGGATTTCCTCCAGGCAAAAGATTTTGCCGACGAAGTGAACGCCCATGCGCCCGGCGCGACCTTTGATGTTCCCATGGGTGAATTTGTCAATGCTGCGAGTTCCATCGCGGTTGTCGTAGATCACAACGTTCTTAGCGATAGTGtt                        >  4:1026943/1‑149 (MQ=255)
                                                                                                                            tCCGGGATTTCCTCCAGGCAAAAGATTTTGCCGACGAAGTGAACGCCCATGCGCCCGGCGCGACCTTTGATGTTCCCATGGGTGAATTTGTCAATGCTGCGAGTTCCATCGCGGTTGTCGTAGATCACAACGTTCTTAGCGATAGTGtt                        <  1:923850/149‑1 (MQ=255)
                                                                                                                               gggATTTCCTCCAGGCAAAAGATTTTGCCGACGAAGTGAACGCCCATGCGCCCGGCGCGACCTTTGATGTTCCCATGGGTGAATTTGTCAATGCTGCGAGTTCCATCGCGGTTGTCGTAGATCACAACGTTCTTAGCGATAGTGtt                        >  4:797118/1‑146 (MQ=255)
                                                                                                                               gggATTTCCTCCAGGCAAAAGATTTTGCCGACGAAGTGAACGCCCATGCGCCCGGCGCGACCTTTGATGTTCCCATGGGTGAATTTGTCAATGCTGCGAGTTCCATCGCGGTTGTCGTAGATCACAACGTTCTTAGCGATAGTGtt                        <  3:797118/146‑1 (MQ=255)
                                                                                                                               gggATTTCCTCCAGGCAAAAGATTTTGCCGACGAAGTGAACGCCCATGCGCCCGGCGCGACCTTTGATGTTCCCATGGGTGAATTTGTCAATGCTGCGAGTTCCATCGCGGTTGTCGTAGATCACAACGTTCTTAGCGATAGTGTTGAc                     <  1:839577/149‑1 (MQ=255)
                                                                                                                                   tttCCTCCAGGCAAAAGATTTTGCCGACGAAGTGAACGCCCATGCGCCCGGCGCGACCTTTGATGTTCCCATGGGTGAATTTGTCAATGCTGCGAGTTCCATCGCGGTTGTCGTAGATCACAACGTTCTTAGCGATAGTGTTGACTCCt                 >  1:332679/1‑149 (MQ=255)
                                                                                                                                   tttCCTCCAGGCAAAAGATTTTGCCGACGAAGTGAACGCCCATGCGCCCGGCGCGACCTTTGATGTTCCCATGGGTGAATTTGTCAATGCTGCGAGTTCCATCGCGGTTGTCGTAGATCACAACGTTCTTAGCGATAGTGTTGACTCCt                 >  1:115127/1‑149 (MQ=255)
                                                                                                                                       ctccAGGCAAAAGATTTTGCCGACGAAGTGAACGCCCATGCGCCCGGCGCGACCTTTGATGTTCCCATGGGTGAATTTGTCAATGCTGCGAGTTCCATCGCGGTTGTCGTAGATCACAACGTTCTTAGCGATAGTGTTGACTCCTTCaa             <  2:566102/149‑1 (MQ=255)
                                                                                                                                          cAGGCAAAAGATTTTGCCGACGAAGTGAACGCCCATGCGCCCGGCGCGACCTTTGATGTTCCCATGGGTGAATTTGTCAATGCTGCGAGTTCCATCGCGGTTGTCGTAGATCACAACGTTCTTAGCGATAGTGTTGACTCCTTCAATGa          >  4:226855/1‑149 (MQ=255)
                                                                                                                                             gCAAAAGATTTTGCCGACGAAGTGAACGCCCATGCGCCCGGCGCGACCTTTGATGTTCCCATGGGTGAATTTGTCAATGCTGCGAGTTCCATCGCGGTTGTCGTAGATCACAACGTTCTTAGCGATAGTGTTGACTCCTTCAATGATCg       >  4:356165/1‑149 (MQ=255)
                                                                                                                                               aaaaGATTTTGCCGACGAAGTGAACGCCCATGCGCCCGGCGCGACCTTTGATGTTCCCATGGGTGAATTTGTCAATGCTGCGAGTTCCATCGCGGTTGTCGTAGATCACAACGTTCTTAGCGATAGTGTTGACTCCTTCAATGATCGTg     >  3:442442/1‑149 (MQ=255)
                                                                                                                                                aaaGATTTTGCCGACGAAGTGAACGCCCATGCGCCCGGCGCGACCTTTGATGTTCCCATGGGTGAATTTGTCAATGCTGCGAGTTCCATCGCGGTTGTCGTAGATCACAACGTTCTTAGCGATAGTGTTGACTCCTTCAATGATCGt      >  1:915553/1‑147 (MQ=255)
                                                                                                                                                aaaGATTTTGCCGACGAAGTGAACGCCCATGCGCCCGGCGCGACCTTTGATGTTCCCATGGGTGAATTTGTCAATGCTGCGAGTTCCATCGCGGTTGTCGTAGATCACAACGTTCTTAGCGATAGTGTTGACTCCTTCAATGATCGt      <  2:915553/147‑1 (MQ=255)
                                                                                                                                                aaaGATTTTGCCGACGAAGTGAACGCCCATGCGCCCGGCGCGACCTTTGATGTTCCCATGGGTGAATTTGTCAATGCTGCGAGTTCCATCGCGGTTGTCGTAGATCACAACGTTCTTAGCGATAGTGTTGACTCCTTCAATGATCGTgg    <  1:1010890/149‑1 (MQ=255)
                                                                                                                                                aaaGATTTTGCCGACGAAGTGAACGCCCATGCGCCCGGCGCGACCTTTGATGTTCCCATGGGTGAATTTGTCAATGCTGCGAGTTCCATCGCGGTTGTCGTAGATCACAACGTTCTTAGCGATAGTGTTGACTCCTTCAATGATCGTgg    >  2:147783/1‑149 (MQ=255)
                                                                                                                                                 aaGATTTTGCCGACGAAGTGAACGCCCATGCGCCCGGCGCGACCTTTGATGTTCCCATGGGTGAATTTGTCAATGCTGCGAGTTCCATCGCGGTTGTCGTAGATCACAACGTTCTTAGCGATAGTGTTGACTCCTTCAATGATCGTGGa   >  1:751830/1‑149 (MQ=255)
                                                                                                                                                  aGATTTTGCCGACGAAGTGAACGCCCATGCGCCCGGCGCGACCTTTGATGTTCCCATGGGTGAATTTGTCAATGCTGCGAGTTCCATCGCGGTTGTCGTAGATCACAACGTTCTTAGCGATAGTGTTGACTCCTTCAATGATCGTGGAg  <  4:749223/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
CTCATCAAAGCGATCCTGTGAAAACTCGGACAAGTGGTCTGGTTGAACGCTGGCTAACAGGTTGATTGGGGTATCAATGCCCTGAATACCGAGCGGAATATCCACCTCTTGATTGAGATTGTCTTCCGGGATTTCCTCCAGGCAAAAGATTTTGCCGACGAAGTGAACGCCCATGCGCCCGGCGCGACCTTTGATGTTCCCATGGGTGAATTTGTCAATGCTGCGAGTTCCATCGCGGTTGTCGTAGATCACAACGTTCTTAGCGATAGTGTTGACTCCTTCAATGATCGTGGAG  >  NC_000913/2762074‑2762368

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