Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,518,197 C→G 13.2% noncoding (35/76 nt) alaW ← tRNA‑Ala

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,518,1970CG13.2% 236.9 / 25.6 89noncoding (35/76 nt)alaWtRNA‑Ala
Reads supporting (aligned to +/- strand):  ref base C (52/25);  new base G (7/5);  total (59/30)
Fisher's exact test for biased strand distribution p-value = 5.29e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.77e-01

CGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGCCCCACGATGCGTTTACGTACCAAGTTTGCTGGGTGCAAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGCCCCGTCACGTAAAGCTTGTCGAGTTGACGGGCGGCA‑TCATATGAATTCCGCCCGAATGTGTCAACGGCAAATTGCAACGTGTAGTTTCAATCGCTGAAAAATCAGGCAAATGAACAAT  >  NC_000913/2518054‑2518346
                                                                                                                                               |                                                                                                                                                      
ctggATCGAACCTCTGACCTCTTTCATGGCATTAAAGCGCTCTACCATCTGAGCTATATCCTCACTATTATTTTACGTTCGTTGTTTGCTGGGTTGTAATTTTTGTGTAGCTAAGCGGGATCGAACCGCTGACCTCTTTCATGGCATGCa                                                                                                                                                  <  2:1760773/148‑1 (MQ=12)
   gATCGAACCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGCCCCACGATGCGTTTACGTACCAAGTTTGCTGGGTGCAAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAgc                                                                                                                                               >  1:928158/1‑150 (MQ=255)
    ctctaaCCTCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATATCCACACGATGATTTTACGTACCAAGTTTGCTGGGTTGAAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGGCATGCAAgag                                                                                                                                              <  2:1171598/146‑3 (MQ=25)
     tCGAACCGCTGACCTCTTTCATGCCATGCTAGCGCATTACCAGAGCAGCTATATCCCCAGTATTAGTTTACGTACCAAGTTTGCTGGGTGCAAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGGCATGCAAgagg                                                                                                                                             <  2:841570/150‑4 (MQ=21)
     tCGAACCGCTGACCTCTTGCATGCCATGCTAGCGCTCTCCCAGCTGAGCTATAGCCCCACGATGCGTATACGTACCAAGTTTGCTGGGTGCAAAATTTGGTGGAGCTAAGCGGGATAGAACCGCTGACCTCTTGCATGGCATGCAAgagg                                                                                                                                             <  1:543457/150‑4 (MQ=38)
     tCGAACCGCTGACCTCTTGCATGCAATGCAAGAGCTGTCCAATCTGAGCTATATGCAGACGATCCTTTTACGTACCAAGTTTGCTGGGTGCTAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGGCATGCAAgagg                                                                                                                                             <  2:261447/150‑4 (MQ=21)
            ccggACCTCTTTCATGCACTGCTACCGCTCTCCAATCTGCGCTATATCCCCACTCTCCTTTCACGTACCAAGTTTGCTGGGTGCAAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCa                                                                                                                                      <  1:976965/147‑1 (MQ=18)
              tGACCTCTTGCATGCGATGCAAGCTGTCTCCCATCTGCGCTATAGCCCGAGTATTCTTTTACGTACCAAGTTTGGTGGGTGGTAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGc                                                                                                                                    <  2:69074/150‑1 (MQ=18)
                      tGCATGGCATGCAAGCGCTCTCCCAGCTGAGCTATAGCCCCACGATGCGTTTACGTACCAAGTTTGCTGGGTGCAAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCtat                                                                                                                            >  1:1944177/1‑150 (MQ=255)
                       gCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGCCCCACGAAGAGTATACGTACCAATTTTCCTGCGTGCAAAATTTGGTGGAGCTAAGCTGGATCGAACCGCTGACCTCTTGCATGCCATGCAACCGATCTCACAGCTGAGCtata                                                                                                                           >  2:911720/1‑150 (MQ=21)
                          tgccatgcAAGCGCTCTCCCAGCTGAGCTATAGCCCCACGATGCGTTTACGTACCAAGTTTGCTGGGTGCAAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCACAGCTGAGCTATAGcc                                                                                                                        >  1:2281032/1‑150 (MQ=255)
                                gcAAGCGCTCTCCCAGCTGAGCTATAGCCCCAACAAGAAAATAATTTCCAATTTACCTGCGTAAAAAATATGATGCAGCTAATCAGAAACGAACCGCTTACATATTGCATACAATGAAAACACTATAACATCTGAGCTATAGcccattca                                                                                                                  >  2:866199/1‑145 (MQ=255)
                                 ccaGCGCTCTCCCAGCGGAGCTATAGCCCCACGATCCGTTTACGTCCCAAGTTTGCTGGGTTGAAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGCCCCGTCAc                                                                                                                 <  1:1119601/148‑1 (MQ=21)
                                                cTGAGCTATAGCCCCACGATGCGTTTACGTACCAAGTTTGCTGGGTGCAAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGCCCCGTCACGTAAAGCTTGTCGAg                                                                                                  >  1:461262/1‑150 (MQ=255)
                                                     cTATAGCCCCACGATGCGTTTACGTACCAAGTTTGCTGGGTGCAAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGCCCCGTCACGTAAAGCTTGTCGAGTTGAc                                                                                             >  1:348895/1‑150 (MQ=255)
                                                      tataGCCCCACGATGCGTTTACGTACCAAGTTTGCTGGGTGCAAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGCCCCGTCACGTAAAGCTTGTCGAGTTGACg                                                                                            >  1:1470293/1‑150 (MQ=255)
                                                          gCCCCACGATGCGTTTACGTACCAAGTTTGCTGGGTGCAAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAGCGCTGTCCCAGCTGCGCTATAGCCCCGTCAGGTAAAGCTTGTCGAGTTGACGggcg                                                                                        >  1:2420980/1‑150 (MQ=255)
                                                            cccACGATGCGTTTACGTACCAAGTTTGCTGGGTGCAAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGCCCCGTCACGTAAAGCTTGTCGAGTTGACGggcggc                                                                                      >  2:1525974/1‑150 (MQ=255)
                                                            cccACGATGCGTTTACGTACCAAGTTTGCTGGGTGCAAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGCCCCGTCACGTAAAGCTTGTCGAGTTGACGggcggc                                                                                      >  2:1218620/1‑150 (MQ=255)
                                                            cccACGATGCGTTTACGTACCAAGTTTGCTGGGTGCAAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGCCCCGTCACGTAAAGCTTGTCGAGTTGACGggcggc                                                                                      >  2:1209514/1‑150 (MQ=255)
                                                            cccACGATGCGTTTACGTACCAAGTTTGCTGGGTGCAAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGCCCCGTCACGTAAAGCTTGTCGAGTTGACGggcggc                                                                                      <  2:1944177/150‑1 (MQ=255)
                                                            cccACGATGCGTTTACGTACCAAGTTTGCTGGGTGCAAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGCCCCGTCACGTAAAGCTTGTCGAGTTGACGggcggc                                                                                      >  2:958342/1‑150 (MQ=255)
                                                             ccACGATGCGTTTACGTACCAAGTTTGCTGGGTGCAAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGCCCCGTCACGTAAAGGGTGTCGAGTTGGCGggcggcg                                                                                     >  1:1070972/1‑149 (MQ=255)
                                                                  atccctatACGCACCAAGCTTGCTGGGTGGTAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGCCCCGTCACGTAAAGCTTGTCGAGTTGACGGGCGGCA‑TCata                                                                               <  1:1716734/143‑1 (MQ=21)
                                                                     cGTTTACGTACCAAGTTTGCTGGGTGCAAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGCCCCGTCACGTAAAGCTTGTCGAGTTGACGGGCGGCA‑TCATATGa                                                                            >  1:769651/1‑150 (MQ=255)
                                                                            gtccgaaGTTTGCTGGGTTCTAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGCCCCGTCACGTAAAGCTTGTCGAGTTGACGGGCGGCA‑TCATATGAATTCCGc                                                                     <  1:128722/145‑1 (MQ=255)
                                                                             tACCAAGTTTGCTGGGTGCAAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGCCCCGTCACGTAAAGCTTGTCGAGTTGACGGGCGGCA‑TCATATGAATTCCGcc                                                                    >  1:1413451/1‑150 (MQ=255)
                                                                             tACCAAGTTTGCTGGGTGCAAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGCCCCGTCACGTAAAGCTTGTCGAGTTGACGGGCGGCA‑TCATATGAATTCCGcc                                                                    >  2:2311759/1‑150 (MQ=255)
                                                                              aCCAAGTTTGCTGGGTGCAAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGCCCCGTCACGTAAAGCTTGTCGAGTTGACGGGCGGCA‑TCATATGAATTCCGccc                                                                   >  1:228446/1‑150 (MQ=255)
                                                                              aCCAAGTTTGCTGGGTGCAAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGCCCCGTCACGTAAAGCTTGTCGAGTTGACGGGCGGCA‑TCATATGAATTCCGccc                                                                   >  2:266407/1‑150 (MQ=255)
                                                                               ccAAGTTTGCTGGGTGCAAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGCCCCGTCACGTAAAGCTTGTCGAGTTGACGGGCGGCA‑TCATATGAATTCCGCCCg                                                                  >  1:2078977/1‑150 (MQ=255)
                                                                               ccAAGGATGCTGGGTGCAAAATTTGGTGGAGCTAAGCTGGATCGAACCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGCCCCGTCACGTAAAGCTTGTCGAGTTGACGGGCGGAATTCATATG‑ATGCCGCCCg                                                                  <  1:181984/150‑1 (MQ=21)
                                                                                cAAGTTTGCTGGGTGCAAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGCCCCGTCACGTAAAGCTTGTCGAGTTGACGGGCGGCA‑TCATATGAATTACGCCCGa                                                                 >  2:2128036/1‑150 (MQ=255)
                                                                                 aaGTTTGCTCGGTGCAAAATTTGGTGCAGCGACGCGGGCGCGAACCGCTGACCTCTTGCACGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGCCCCGTCACGTAAAGCTTGTCGAGTTGACGGGCGGAATTCATATG‑ATGCCGCCCgtc                                                                <  1:866199/150‑3 (MQ=21)
                                                                                     ttGCTGGGTTCAAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGCCCCGTCACGTAAAGCTTGTCGAGTTGACGGGCGGCA‑TCATATGAATTCCGCCCGAAtgtg                                                            <  2:401202/150‑1 (MQ=255)
                                                                                     ttGCTGGGTGCAAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGCCCCGTCACGTAAAGCTTGTCGAGTTGACGGGCGGCA‑TCATATGAATTCCGCCCGAAtgtg                                                            <  1:1186851/150‑1 (MQ=255)
                                                                                        cTGGGTGCAAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGCCCCGTCACGTAAAGCTTGTCGAGTTGACGGGCGGCA‑TCATATGAATTCCGCCCGAATGTGTCa                                                         >  2:2146571/1‑150 (MQ=255)
                                                                                            gTGCAAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGCCCAGTCACGTAAAGCTTGTCGAGTTGACGGGCGGCA‑TCATATGAATTCCGCCCGAATGTGTCAAAgg                                                     >  1:605779/1‑150 (MQ=255)
                                                                                             tGCAAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGCCCCGTCACGTAAAGCTTGTCGAGTTGACGGGCGGCA‑TCATATGAATTCCGCCCGAATGTGTCAACGGc                                                    >  1:133248/1‑150 (MQ=255)
                                                                                             tGCAAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGCCCCGTCACGTAAAGCTTGTCGAGTTGACGGGCGGCA‑TCATATGAATTCCGCCCGAATGTGTCAACGGc                                                    <  1:2146571/150‑1 (MQ=255)
                                                                                              gCAAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGCCCCGTCACGTAAAGCTTGTCGAGTTGACGGGCGGCA‑TCATATGAATTCCGCCCGAATGTGTCAACGGCa                                                   >  1:27867/1‑150 (MQ=255)
                                                                                               ggaaaTTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGCCCCGTCAGGAAAAGAATGTCGAGTTGACGTGCGGCA‑TCATATTAATTCCGCACGAATGTGTCAACGGCaa                                                  >  1:432206/3‑150 (MQ=18)
                                                                                               ggaaaTTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGCCCCGGCAAGTAAAGAATTTCGAGTTGACGGGCGGCA‑TCATATTAATTCCGCAAAGATGGGTCAACCGCaa                                                  >  2:2155311/3‑150 (MQ=17)
                                                                                               cAAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGCCCCGTCACGTAAAGCTTGTCGAGTTGACGGGCGGCA‑TCATATGAATTCCGCCCGAATGTGTCACCGGCaa                                                  >  1:1019240/1‑150 (MQ=255)
                                                                                               cAAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGCCCCGTCACGTAAAGCTTGTCGAGTTGACGGGCGGCA‑TCATATGAATTCCGCCCGAATGTGTCAACGGCaa                                                  >  2:442646/1‑150 (MQ=255)
                                                                                               cAAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGCCCCGTCACGTAAAGCTTGTCGAGTTGACGGGCGGCA‑TCATATGAATTCCGCCCGAATGTGTCAACGGCaa                                                  >  1:791112/1‑150 (MQ=255)
                                                                                               cAAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGCCCCGTCACGTAAAGCTTGTCGAGTTGACGGGCGGCA‑TCATATGAATTCCGCCCGAATGTGTCAACGGCaa                                                  >  1:2306210/1‑150 (MQ=255)
                                                                                               cAAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGCCCCGTCACGTAAAGCTTGTCGAGTTGACGGGCGGCA‑TCATATGAATTCCGCCCGAATGTGTCAACGGCaa                                                  >  2:675966/1‑150 (MQ=255)
                                                                                                 aaaTTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGCCCCGTCACGTAAAGCTTGTCGAGTTGACGGGCGGCA‑TCATATGAATTCCGCCCGAATGTGTCAACGGCAAAt                                                >  2:487870/1‑150 (MQ=255)
                                                                                                   aTTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGCCCCGTCATGTAAATTTTGTCGAGTTGACGGGAGTCA‑TAATATGAATTCCTCGCGAATGTGTCAACTGGAAATTg                                              <  2:1122788/150‑1 (MQ=18)
                                                                                                   aTTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGCCCCGTCACGTAAAGCTTGTCGAGTTGACGGGCGGCA‑TCATATGAATTCCGCCCGAATGTGTCAACGGCAAATTg                                              <  1:602499/150‑1 (MQ=255)
                                                                                                    tttGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGCCCCGTCACGTAAAGCTTGTCGAGTTGACGGGCGGCA‑TCATATGAATTCCGCCCGAATGTGTCAACGGCAAATTGc                                             <  1:262054/150‑1 (MQ=255)
                                                                                                      tggtggAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGCCCCGTCACGTAAAGCTTGTCGAGTTGACGGGCGGCA‑TCATATGAATTCCGCCCGAATGTGTCAACGGCAAATTGCaa                                           >  2:1355550/1‑150 (MQ=255)
                                                                                                       ggtggAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGCCCCGTCACGTAAAGCTTGTCGAGTTGACGGGCGGCA‑TCATATGAATTCCGCCCGAATGTGTCAACGGCAAATTGCAAc                                          <  1:958342/150‑1 (MQ=255)
                                                                                                           gAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGCCCCGTCACGTAAAGCTTGTCGAGTTGACGGGCGGCA‑TCATATGAATTCCGCCCGAATGTGTCAACGGCAAATTGCAACgtgt                                      <  1:1265146/150‑1 (MQ=255)
                                                                                                            aGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGCCCCGTCACGTAAAGCTTGTCGAGTTGACGGGCGGCA‑TCATATGAATTCCGCCCGAATGTGTCAACGGCAAATTGCAACGTGTa                                     >  2:1219260/1‑150 (MQ=255)
                                                                                                            aGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGCCCCGTCACGTAAAGCTTGTCGAGTTGACGGGCGGCA‑TCATATGAATTCCGCCCGAATGTGTCAACGGCAAATTGCAACGTGTa                                     >  2:1254266/1‑150 (MQ=255)
                                                                                                             gCTAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGCCCCGTCACGTAAAGCTTGTCGAGTTGACGGGGGGCA‑TCATATGAATTCCGCCCGAATGTGTCAACGGCAAATTGCAACGTGTAg                                    >  1:2187142/1‑150 (MQ=255)
                                                                                                             gCTAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGCCCCGTCACGTAAAGCTTGTCGAGTTGACGGGCGGCA‑TCATATGAATTCCGCCCGAATGTGTCAACGGCAAAATGCAACGTGTAg                                    >  2:1001936/1‑150 (MQ=255)
                                                                                                                 aGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGCCCCGTCACGTAAAGCTTGTCGAGTTGACGGGCGGCA‑TCATATGAATTCCGCCCGAATGTGTCAACGGCAAATTGCAACGTGTAGTTTc                                <  1:2348733/150‑1 (MQ=255)
                                                                                                                 aGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGCCCCGTCACGTAAAGCTTGTCGAGTTGACGGGCGGCA‑TCATATGAATTCCGCCCGAATGTGTCAACGGCAAATTGCAACGTGTAGTTTc                                >  2:2055900/1‑150 (MQ=255)
                                                                                                                  gCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGCCCCGTCACGTAAAGCTTGTCGAGTTGACGGGCGGCA‑TCATATGAATTCCGCCCGAATGTGTCAACGGCAAATTGCAACg                                         <  2:1614318/140‑1 (MQ=255)
                                                                                                                  gCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGCCCCGTCACGTAAAGCTTGTCGAGTTGACGGGCGGCA‑TCATATGAATTCCGCCCGAATGTGTCAACGGCAAATTGCAACg                                         >  1:1614318/1‑140 (MQ=255)
                                                                                                                  gCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGCCCCGTCACGTAAAGCTTGTCGAGTTGACGGGCGGCA‑TCATATGAATTCCGCCCGAATGTGTCAACGGCAAATTGCAACGTGTAGTTTCa                               >  1:2203939/1‑150 (MQ=255)
                                                                                                                  gCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGCCCCGTCACGTAAAGCTTGTCGAGTTGACGGGCGGCA‑TCATATGAATTCCGCCCGAATGTGTCAACGGCAAATTGCAACGTGTAGTTTCa                               <  1:1219260/150‑1 (MQ=255)
                                                                                                                   cGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGCCCCGTCACGTAAAGCTTGTCGAGTTGACGGGCGGCA‑TCATATGAATTCCGCCCGAATGTGTCAACGGCAAATTGCAACGTGTAGTTTCaa                              <  2:285350/150‑1 (MQ=255)
                                                                                                                   cGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGCCCAGTCACGTAAAGCTTGTCGAGTTGACGGGCGGCA‑TCATATGAATTACGCCCGAATGTGTCAACGGCAAATTGCAACGTGTAGTTTCaa                              >  2:420172/1‑150 (MQ=255)
                                                                                                                     ggATCGAACCGCTGACCTCTTGCATGCCAAGCAAGCGCTCTCCCAGCTGAGCTATAGCCCCGTCACGTAAAGCTTGTCGAGTTGACGGGCGGCA‑TCATATGAATTCCGCCCGAATGTGTCAACGGCAAATTGCAACGTGTAGTTTCAATc                            >  1:249305/1‑150 (MQ=255)
                                                                                                                         cGAACCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGCCCCGTCACGTAAAGCTTGTCGAGTTGACGGGCGGCA‑TCATATGAATTCCGCCCGAATGTGTCAACGGCAAATTGCAACGTGTAGTTTCAATCGCTg                        >  1:2302719/1‑150 (MQ=255)
                                                                                                                          gAACCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGCCCCGTCACGTAAAGCTTGTCGAGTTGACGGGCGGCA‑TCATATGAATTCCGCCCGAATGTGTCAACGGCAAATTGCAACGTGTAGTTTCAATCGCTGa                       >  1:2272180/1‑150 (MQ=255)
                                                                                                                           aaCCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGCCCCGTCACGTAAAGCTTGTCGAGTTGACGGGCGGCA‑TCATATGAATTCCGCCCGAATGTGTCAACGGCAAATTGCAACGTGTAGTTTCAATCGCTGaa                      >  1:2340856/1‑150 (MQ=255)
                                                                                                                                 tGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGCCCCGTCACGTAAAGCTTGTCGAGTTGACGGGCGGCA‑TCATATGAATTCCGCCCGAATGTGTCAACGGCAAATTGCAACGTGTAGTTTCAATCGCTGAAAAATCa                >  2:2203278/1‑150 (MQ=255)
                                                                                                                                 tGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGCCCCGTCACGTAAAGCTTGTCGAGTTGACGGGCGGCA‑TCATATGAATTCCGCCCGAATGTGTCAACGGCAAATTGCAACGTGTAGTTTCAATCGCTGAAAAATCa                <  2:2340856/150‑1 (MQ=255)
                                                                                                                                  gACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGCCCCGTCACGTAAAGCTTGTCGAGTTGACGGGCGGCA‑TCATATGAATTCCGCCCGAATGTGTCAACGGCAAATTGCAACGTGTAGTTTCAATCGCTGAAAAATCAg               >  2:1913531/1‑150 (MQ=255)
                                                                                                                                  gACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGCCCCGTCACGTAAAGCTTGTCGAGTTGACGGGCGGCA‑TCATATGAATTCCGCCCGAATGTGTCAACGGCAAATTGCAACGTGTAGTTTCAATCGCTGAAAAATCAg               <  1:2318891/150‑1 (MQ=255)
                                                                                                                                   aCCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGCCCCGTCACGTAAAGCTTGTCGAGTTGACGGGCGGCA‑TCATATGAATTCCGCCCGAATGTGTCAACGGCAAATTGCAACGTGTAGTTTCAATCGCTGAAAAATCAgg              <  2:2119493/150‑1 (MQ=255)
                                                                                                                                    gcgctTGCATGGCATGCAAGCGCTCTCCCAGCTGAGCTATAGCCCCGTCACGTAAAGCTTGTCGAGTTGACGGGCGGCA‑TCATATGAATTCCGCCCGAATGTGTCAACGGCAAATTGCAACGTGTAGTTTCAATCGCTGAAAAATCAGGc             >  2:978199/4‑150 (MQ=255)
                                                                                                                                    gcgctTGCATGGCATGCAAGCGCTCTCCCAGCTGAGCTATAGCCCCGTCACGTAAAGCTTGTCGAGTTGACGGGCGGCA‑TCATATGAATTCCGCCCGAATGTGTCAACGGCAAATTGCAACGTGTAGTTTCAATCGCTGAAAAATCAGGc             >  1:1987222/4‑150 (MQ=255)
                                                                                                                                     ctctTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGCCCCGTCACGTAAAGCTTGTCGAGTTGACGGGCGGCA‑TCATATGAATTCCGCCCGAATGTGTCAACGGCAAATTGCAACGTGTAGTTTCAATCGCTGAAAAATCAGGCa            <  1:1797583/150‑1 (MQ=255)
                                                                                                                                     cgctTGCATGGCATGCAAGCGCTCTCCCAGCTGAGCTATAGCCCCGTCACGTAAAGCTTGTCGAGTTGACGGGCGGCA‑TCATATGAATTCCGCCCGAATGTGTCAACGGCAAATTGCAACGTGTAGTTTCAATCGCTGAAAAATCAGGCa            >  1:63791/3‑150 (MQ=255)
                                                                                                                                     cgctTGCATGGCATGCAAGCGCTCTCCCAGCTGAGCTATAGCCCCGTCACGTAAAGCTTGTCGAGTTGACGGGCGGCA‑TCATATGAATTCCGCCCGAATGTGTCAACGGCAAATTGCAACGTGTAGTTTCAATCGCTGAAAAATCAGGCa            >  1:550611/3‑150 (MQ=255)
                                                                                                                                     cgctTGCATGGCATGCAAGCGCTCTCCCAGCTGAGCTATAGCCCCGTCACGTAAAGCTTGTCGAGTTGACGGGCGGCA‑TCATATGAATTCCGCCCGAATGTGTCAACGGCAAATTGCAACGTGTAGTTTCAATCGCTGAAAAATCAGGCa            >  1:292279/3‑150 (MQ=255)
                                                                                                                                     cgctTGCATGGCATGCAAGCGCTCTCCCAGCTGAGCTATAGCCCCGTCACGTAAAGCTTGTCGAGTTGACGGGCGGCA‑TCATATGAATTCCGCCCGAATGTGTCAACGGCAAATTGCAACGTGTAGTTTCAATCGCTGAAAAATCAGGCa            >  1:1243980/3‑150 (MQ=255)
                                                                                                                                     cgctTGCATGGCATGCAAGCGCTCTCCCAGCTGAGCTATAGCCCCGTCACGTAAAGCTTGTCGAGTTGACGGGCGGCA‑TCATATGAATTCCGCCCGAATGTGTCAACGGCAAATTGCAACGTGTAGTTTCAATCGCTGAAAAATCAGGCa            >  1:1790654/3‑150 (MQ=255)
                                                                                                                                      tctTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGCCCCGTCACGTAAAGCTTGTCGAGTTGACGGGCGGCA‑TCATATGAATTCCGCCCGAATGTGTCAACGGCAAATTGCAACGTGTAGTTTCAATCGCTGAAAAATCAGGCaa           >  1:95716/1‑150 (MQ=255)
                                                                                                                                      tctTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGCCCCGTCACGTAAAGCTTGTCGAGTTGACGGGCGGCA‑TCATATGAATTCCGCCCGAATGTGTCAACGGCAAATTGCAACGTGTAGTTTCAATCGCTGAAAAATCAGGCaa           >  1:883907/1‑150 (MQ=255)
                                                                                                                                      tctTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGCCCCGTCACGTAAAGCTTGTCGAGTTGACGGGCGGCA‑TCATATGAATTCCGCCCGAATGTGTCAACGGCAAATTGCAACGTGTAGTTTCAATCGCTGAAAAATCAGGCaa           >  1:25003/1‑150 (MQ=255)
                                                                                                                                      tctTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGCCCCGTCACGTAAAGCTTGTCGAGTTGACGGGCGGCA‑TCATATGAATTCCGCCCGAATGTGTCAACGGCAAATTGCAACGTGTAGTTTCAATCGCTGAAAAATCAGGCaa           >  1:1742912/1‑150 (MQ=255)
                                                                                                                                         tGCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGCCCCGTCACGTAAAGCTTGTCGAGTTGACGGGCGGCA‑TCATATGAATTCCGCCCGAATGTGTCAACGGCAAATTGCAACGTGTAGTTTCAATCGCTGAAAAATCAGGCAAATg        <  1:911720/150‑1 (MQ=255)
                                                                                                                                            atggcatgcAAGCGCTCTCCCAGCTGAGCTATAGCCCCGTCACGTAAAGCTTGTCGAGTTGACGGGCGGCA‑TCATATGAATTCCGCCCGAATGTGTCAACGGCAAATTGCAACGTGTAGTTTCAATCGCTGAAAAATCAGGCAAATGAAc     <  2:292279/150‑1 (MQ=255)
                                                                                                                                            atgccatgcAAGCGCTCTCCCAGCTGAGCTATAGCCCCGTCACGTAAAGCTTGTCGAGTTGACGGGCGGCA‑TCATATGAATTCCGCCCGAATGTGTCAACGGCAAATTGCAACGTGTAGTTTCAATCGCTGAAAAATCAGGCAAATGAAc     <  1:992428/150‑1 (MQ=255)
                                                                                                                                               ccatgcAAGCGCTCTCCCAGCTGAGCTATAGCCCCGTCTCGTAAAGCTTGTCGAGCTGACGGGCGGCA‑TCATATGAATTCCGCCCGAATGTGTCAACGGCAAATTGCAACGTGTAGTTTCAATCGCTGAAAAATCAGGCAAATGAACAAt  >  1:2453285/1‑150 (MQ=255)
                                                                                                                                               |                                                                                                                                                      
CGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGCCCCACGATGCGTTTACGTACCAAGTTTGCTGGGTGCAAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGCCCCGTCACGTAAAGCTTGTCGAGTTGACGGGCGGCA‑TCATATGAATTCCGCCCGAATGTGTCAACGGCAAATTGCAACGTGTAGTTTCAATCGCTGAAAAATCAGGCAAATGAACAAT  >  NC_000913/2518054‑2518346

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