Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,090,780 C→T Q197* (CAA→TAA)  hisG → ATP phosphoribosyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,090,7800CT84.4% 116.3 / 12.5 45Q197* (CAA→TAA) hisGATP phosphoribosyltransferase
Reads supporting (aligned to +/- strand):  ref base C (5/2);  new base T (21/17);  total (26/19)
Fisher's exact test for biased strand distribution p-value = 6.81e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGGACGGTCCGCTCTCCTTAAACGGTAAACGTATCGCCACCTCTTATCCTCACCTGCTCAAGCGTTATCTCGACCAGAAAGGCATCTCTTTTAAATCCTGCTTACTGAACGGTTCTGTTGAAGTCGCCCCGCGTGCCGGACTGGCGGATGCGATTTGCGATCTGGTTTCCACCGGTGCCACGCTGGAAGCTAACGGCCTGCGCGAAGTCGAAGTTATCTATCGCTCGAAAGCCTGCCTGATTCAACGCGATGGCGAAATGGAAGAATCCAAACAGCAACTGATCGACAAACTGCTGACCCGTATTCAGGGTGTGATCCAGGCGCGCGAATCAAAATACATCATGATGCACGCACCGACCGAACGTCTGGATGAAGTCATCGCCCTGCTGCCAGGTGCCGAACGCCCAACTATTCTGCCGCTGGCGGGTGACCAACAGCGCGTAGCGATGCACATGGTCAGCAGCGAAACCCTGTTCTGGGAAACCAT  >  NC_000913/2090538‑2091024
                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                      
gggACGGTCCGCTCTCCTTAAACGGTAAACGTATCGCCACCTCTTATCCTCACCTGCTCAAGCGTTATCTCGACCAGAAAGGCATCTCTTTTAAATCCTGCTTACTGAACGGTTCTGTTGAAGTCGCCCCGCGTGCCGGACTGGCGGATGCGATTTGCGATCTGGTTTCCACCGGTGCCACGCTGGAAGCTAACGGCCTGCGCGAAGTCGAAGTTATCTATCGCTCGAAAGCCTGCCTGATTTAAcgcg                                                                                                                                                                                                                                                  <  1:250092/249‑1 (MQ=255)
    cGGTCCGCTCTCCTTAAACGGTAAACGTATCGCCTCCTCTTATCCTCACCTGCTCAAGCGTTATCTCGACCAGAAAGGCATCTCTTTTAAATCCTGCTTACTGAACGGTTCTGTTGAAGTCGCCCCGCGTGCCGGACTGGCGGATGCGATTTGCGATCTGGTTTCCACCGGTGCCACGCTGGAAGCTAACGGCCTGCGCGAAGTCGAAGTTATCTATCGCTCGAAAGCCTGCCTGATTTAACGCTATGGCg                                                                                                                                                                                                                                            <  1:132776/251‑1 (MQ=255)
        ccGCTCTCCTTAAACGGTAAACGTATCGCCACCTCTTATCCTCACCTGCTCAAGCGTTATCTCGACCAGAAAGGCATCTCTTTTAAATCCTGCTTACTGAACGGTTCTGTTGAAGTCGCCCCGCGTGCCGGACTGGCGGATGCGATTTGCGATCTGGTTTCCACCGGTGCCACGCTGGAAGCTAACGGCCTGCGCGAAGTCGAAGTTATCTATCGCTCGAAAGCCTGCCTGATTTAACGCGATGGCGaaa                                                                                                                                                                                                                                         <  1:148545/250‑1 (MQ=255)
                 ttAAACGGTAAACGTATCGCCACCTCTTATCCTCACCTGCTCAAGCGTTATCTCGACCAGAAAGGCATCTCTTTTAAATCCTGCTTACTGAACGGTTCTGTTGAAGTCGCCCCGCGTGCCGGACTGGCGGATGCGATTTGCGATCTGGTTTCCACCGGTGCCACGCTGGAAGCTAACGGCCTGCGCGAAGTCGAAGTTATCTATCGCTCAAAAGCCTGCCTGATTTAACGCGATGGCGAAATGGAAGAATc                                                                                                                                                                                                                               >  1:183670/1‑251 (MQ=255)
                        gTAAACGTATCGCCACCTCTTATCCTCACCTGCTCAAGCGTTATCTCGACCAGAAAGGCATCTCTTTTAAATCCTGCTTACTGAACGGTTCTGTTGAAGTCGCCCCGCGTGCCGGACTGGCGGATGCGATTTGCGATCTGGTTTCCACCGGTGCCACGCTGGAAGCTAACGGCCTGCGCGAAGTCGAAGTTATCTATCGCTCGAAAGCCTGCGTGATTTAAGGCGATGGCGAAATGGAAGAATCCAAACAg                                                                                                                                                                                                                        <  1:378594/251‑1 (MQ=255)
                               tATCGCCACCTCTTATCCTCACCTGCTCAAGCGTTATCTCGACCAGAAAGGCATCTCTTTTAAATCCTGCTTACTGAACGGTTCTGTTGAAGTCGCCCCGCGTGCCGGACTGGCGGATGCGATTTGCGCTCTGGTTTCCCCCGGTGCCACGCTGGAAGCTAACGGCCTGCGCGCAGTCGAAGTTATCTCTCGCTCGAAAGCCTGCCTGATTCAACGCGATGGCGAAATTGAATAATCCAAACAGCCACTg                                                                                                                                                                                                                  >  2:354116/1‑250 (MQ=255)
                                  cGCCACCTCTTATCCTCACCTGCTCAAGCGTTATCTCGACCAGAAAGGCATCTCTTTTAAATCCTGCTTACTGAACGGTTCTGTTGAAGTCGCCCCGCGTGCCGGACTGGCGGATGCGATTTGCGATCTGGTTTCCACCGGTGCCACGCTGGAAGCTAACGGCCTGCGCGAAGTCGAAGTTATCTATCGCTCGAAAGCCTGCCTGATTTAACGCGATGGCGAAATGGAAGAATCCAAACAGCAACTGATc                                                                                                                                                                                                               <  2:143264/250‑1 (MQ=255)
                                       ccTCTTATCCTCACCTGCTCAAGCGTTATCTCGACCAGAAAGGCATCTCTTTTAAATCCTGCTTACTGAACGGTTCTGTTGAAGTCGCCCCGCGTGCCGGACTGGCGGATGCGATTTGCGATCTGGTTTCCACCGGTGCCACGCTGGAAGCTAACGGCCTGCGCGAAGTCGAAGTTATCTATCGCTCGAAAGCCTGCCTGATTTAACGCGATGGCGAAATGGAAGAATCCAAACAGCAACTGATCGACaaa                                                                                                                                                                                                         <  1:420439/251‑1 (MQ=255)
                                                  cACCTGCTCAAGCGTTATCTCGACCAGAAAGGCATCTCTTTTAAATCCTGCTTACTGAACGGTTCTGTTGAAGTCGCCCCGCGTGCCGGACTGGCGGATGCGATTTGCGATCTGGTTTCCACCGGTGCCACGCTGGAAGCTAACGGCCTGCGCGAAGTCGAAGTTATCTATCGCTCAAAAGCCTGCCTGATTTAACGCGATGGCGAAATGGAAGAATCCAAACAGCAACTGATCGACAAACTGCTGAcc                                                                                                                                                                                                <  2:183670/249‑1 (MQ=255)
                                                           aaGCGTTATCTCGACCAGAAAGGCATCTCTTTTAAATCCTGCTTACTGAACGGTTCTGTTGAAGTCGCCCCGCGTGCCGGACTGGCGGATGCGATTTGCGATCTGGTTTCCACCGGTGCCACGCTGGAAGCTAACGGCCTGCGCGAAGTCGAAGTTATCTATCGCTCGAAAGCCTGCCTGATTTAACGCGATGGCGAAATGGAAGAATCCAAACAGCAACTGATCGACAAACTGCTGACCCGTATTCAggg                                                                                                                                                                                     >  1:372992/1‑251 (MQ=255)
                                                              cGTTATCTCGACCAGAAAGGCATCTCTTTTAAATCCTGCTTACTGAACGGTTCTGTTGAAGTCGCCCCGCGTGCCGGACTGGCGGATGCGATTTGCGATCTGGTTTCCACCGGTGCCACGCTGGAAGCTAACGGCCTGCGCGAAGTCGAAGTTATCTATCGCTCGAAAGCCTGCCTGATTTAACGCGATGGCGAAATGGAAGAATCCAAACAGCAACTGATCGACAAACTACTGACCCGTATTCAGGgtgt                                                                                                                                                                                  >  1:89576/1‑251 (MQ=255)
                                                                      cGACCAGAAAGGCATCTCTTTTAAATCCTGCTTACTGAACGGTTCTGTTGAAGTCGCCCCGCGTGCCGGACTGGCGGATGCGATTTGCGATCTGGTTTCCACCGGGGCCACGCTGGAAGCTAACGGCCTGCGGGAAGTCGAAGTTTTCTATCGCTCGAAAGCCTGCCTGCTTTAACGCGATGGCGAAATGGAAGAATCCAAACAGAAACTGATCGAAAAACTGCTGACCCGTATTCAGGGTGTGATCCAgg                                                                                                                                                                          >  2:125163/1‑251 (MQ=255)
                                                                       gACCAGAAAGGCATCTCTTTTAAATCCTGCTTACTGAACGGTGCTGTTGAAGTCGCCCCGCGTGCCGGACTGGCGGATGCGATTTGCGATCTGGTTTCCACCGGTGCCACGCTGGAAGCTAACGGCCTGCGCGAAGTCGAAGTTATCTATCGCTCGAAAGCCTGCCTGATTTAACGCGATGGCGAAATGGAAGAATCCAAGCAGCAACTGATCGACAAACTGCTGACCCGTATTCAGGGTGTGATCCAGgc                                                                                                                                                                         <  1:456420/251‑1 (MQ=255)
                                                                                ggCATCTCTTTTAAATCCTGCTTACTGAACGGTTCTGTTGAAGTCGCCCCGCGTGCCGGACTGGCGGATGCGATTTGCGATCTGGTTTCCACCGGTGCCACGCTGGAAGCTAACGGCCTGCGCGAAGTCGAAGTTATCTATCGCTCGAAAGCCTGCCTGATTTAACGCGATGGCGAAATGGAAGAATCCAAACAGCAACTGATCGACAAACTGCTGACCCGTATTCAGGGTGTGATCCAGGCGCGCGAATc                                                                                                                                                                <  2:7386/251‑1 (MQ=255)
                                                                                  cATCTCTTTTAAATCCTGCTTACTGAACGGTTCTGTTGAAGTCGCCCCGCGTGCCGGACTGGCGGATGCGATTTGCGATCTGGTTTCCACCGGTGCCACGCTGGAAGCTAACGGCCTGCGCGAAGTCGAAGTTATCTATCGCTCGAAAGCCTGCCTGATTTAACGCGATGGCGAAATGGAAGAATCCAAACAGCAACTGATCGACAAACTGCTGACCCGTATTCAGGGTGTGATCCAGGCGCGCGAATCa                                                                                                                                                               >  1:23290/1‑250 (MQ=255)
                                                                                            aaaTCCTGCTTACTGAACGGTTCTGTTGAAGTCGCCCCGCGTGCCGGACTGGCGGATGCGATTTGCGATCTGGTTTCCACCGGTGCCAAGCTGGAAGCTAACGGCCTGCGCGAAGTCGAAGTTATCTATCGCTCGAAAGCCTGCCTGATTTAACGCGATGGCGAAATGGAAGAATCCAAACAGCAACTGATCGACAAACTGCTGACCCGTATTCAGGGTGTGATCCAGGCGCGCGAATCAAAATAcatca                                                                                                                                                     >  2:460960/1‑250 (MQ=255)
                                                                                                ccTGCTTACTGAACGGTTCTGTTGAAGTCGCCCCGCGTGCCGGACTGGCGGATGCGATTTGCGATCTGGTTTCCACCGGTGCCACGCTGGAAGCTAACGGCCTGCGCGAAGTCGAAGTTATCTATCGCTCGAAAGCCTGCCTGATTTAACGCGATGGCGAAATGGAAGAATCCAAACAGCAACTGATCGACAAACTGCTGACCCGTATTCAGGGTGTGATCCAGGCGCGCGAATCAAAATACATCatgatg                                                                                                                                                <  1:356922/251‑1 (MQ=255)
                                                                                                 cTGCTTACTGAACGGTTCTGTTGAAGTCGCCCCGCGTGCCGGACTGGCGGATGCGATTTGCGATCTGGTTTCCACCGGTGCCACGCTGGAAGCTAACGGCCTGCGCGAAGTCGAAGTTATCTATCGCTCGAAAGCCTGCCTGATTCAACGCGATGGCGAAATGGAAGAATCCAAACAGCAACTGATCGACAAACTGCTGACCCGTATTCAGGGTGTGATCCAGGCGCGCGAATCAAAATACATCATGATgc                                                                                                                                               >  1:260834/1‑251 (MQ=255)
                                                                                                    cTTACTGAACGGTTCTGTTGAAGTCGCCCCGCGTGCCGGACTGGCGGATGCGATTTGCGATCTGGTTTCCACCGGTGCCACGCTGGAAGCTAACGGCCTGCGCGAAGTCGAAGTTATCTATCGCTCGAAAGCCTGCCTGATTTAACGCGATGGCGAAATGGAAGAATCCAAACAGCAACTGATCGACAAACTGCTGACCCGTATTCAGGGTGTGATCCAGGCGCGCGAATCAAAATACATAATGATgcacg                                                                                                                                            >  2:175415/1‑251 (MQ=255)
                                                                                                       aCTGAACGGTTCTGTTGAAGTCGCCCCGCGTGCCGGACTGGCGGATGCGATTTGCGATCTGGTTTCCACCGGTGCCACGCTGGAAGCTAACGGCCTGCGCGAAGTCGAAGTTATCTATCGCTCGAAAGCCTGCCTGATTTAACGCGATGGCGAAATGGAAGAATCCAAACAGCAACTGATCGACAAACTGCTGACCCGTATTCAGGGTGTGATCCAGGCGCGCGAATCAAAATACATCATGATgcacgca                                                                                                                                          <  1:175415/250‑1 (MQ=255)
                                                                                                          gAACGGTTCTGTTGAAGTCGCCCCGCGTGCCGGACTGGCGGATGCGATTTGCGATCTGGTTTCCACCGGTGCCACGCTGGAAGCTAACGGCCTGCGCGAAGTCGAAGTTATCTATCGCTCGAAAGCCTGCCTGATTTAACGCGATGGCGAAATGGAAGAATCCAAACAGCAACTGATCGACAAACTGCTGACCCGTATTCAGGGTGTGATCCAGGCGCGCGAATCAAAATACATCATGATGCACGCaccga                                                                                                                                      <  2:126445/251‑1 (MQ=255)
                                                                                                          gAACGGTTCTGTTGAAGTCGCCCCGCGTGCCGGACTGGCGGATGCGAGTTGCGATCTGGTTTCCACCGGTGCCACGCTGGAAGCTCACGGCCTGCGCGCAGGCGCAGCTATCTATCGCTCGCAAGCCTGCCTGATTTAACGCGCTGGCGAAATGGGAGAATCCCAACAGCAACTGATCGACAAACTGCTGACCCGCATTCAGGGTGTGATCCAGGCGCGCGAATCAAAATACATCATGATGCACGCaccgc                                                                                                                                      >  1:238247/1‑250 (MQ=255)
                                                                                                           aaCGGTTCTGTTGAAGTCGCCCCGCGTGCCGGACTGGCGGATGCGATTTGCGATCTGGTTTCCACCGGTGCCACGCTGGAAGCTAACGGCCTGCGCGAAGTCGAAGTTATCTATCGCTCGAAAGCCTGCCTGATTTAACGCGATGGCGAAATGGAAGAATCCAAACAGCAACTGATCGACAAACTGCTGACCCGTATTCAGGGTGTGATCCAGGCGCGCGAATCAAAATACATCATGATGCACGCaccgac                                                                                                                                     >  1:366277/1‑251 (MQ=255)
                                                                                                              ggTTCTGTTGAAGTCGCCCCGCGTGCCGGACTGGCGGATGCGATTTGCGATCTGGTTTCCACCGGTGCCACGCTGGAAGCTAACGGCCTGCGCGAAGTCGAAGTTATCTATCGCTCGAAAGCCTGCCTGATTTAACGCGATGGCGAAATGGAAGAATCCAAACAGCAACTGATCGACAAACTGCTGACCCGTATTCAGGGTGTGATCCAGGCGCGCGAATCAAAATACATCATGATGCACGCaccgaccg                                                                                                                                   >  1:122039/1‑250 (MQ=255)
                                                                                                                               cccGCGTGCCGGACTGGCGGATGCGATTTGCGATCTGGTTTCCACCGGTGCCACGCTGGAAGCTAACGGCCTGCGCGAAGTCGAAGTTATCTATCGCTCGAAAGCCTGCCTGATTCAACGCGATGGCGAAATGGAAGAATCCAAACAGCAACTGATCGACAAACTGCTGACCCGTATTCAGGGTGTGATCCAGGCGCGCGAATCAAAATACATCATGATGCACGCACCGACCGAACGTCTGGATGAAGTc                                                                                                                  >  1:167558/1‑250 (MQ=255)
                                                                                                                                 cgcgTGCCGGACTGGCGGATGCGATTTGCGATCTGGTTTCCACCGGTGCCACGCTGGAAGCTAACGGCCTGCGCGAAGTCGAAGTTATCTATCGCTCGAAAGCCTGCCTGATTTAACGCGATGGCGAAATGGAAGAATCCAAACAGCAACTGATCGAAAAACTGCTGACCCGGATTCAGGGGGTGATCCAGGCGCGCGAATCAAAATACATCATGATGCACGCACCGACCGAACGTCTGGATGAAGTCATc                                                                                                               >  2:52479/1‑251 (MQ=255)
                                                                                                                                     tGCCGGACTGGCGGATGCGATTTGCGATCTGGTTTCCAACGGTGCCACGCTGGAAGCTAACGGCCTGCGCGAAGTCGAAGTTATCTATCGCTCGAAAGCCTGCCTGATTCAACGCGATGGCGAAATGGAAGAATCCAAACAGCAACTGATCGACAAACTGCTGACCCGTATTCAGGGTGTGATCCAGGCGCGCGAATCAAAATACATCATGATGCACGCACCGACCGAACGTCTGGATGAAGTCATCGccc                                                                                                           >  1:395843/1‑251 (MQ=255)
                                                                                                                                          gACTGGCGGATGCGATTTGCGATCTGGTTTCCACCGGTGCCACGCTGGAAGCTAACGGCCTGCGCGAAGTCGAAGTTATCTATCGCTCGAAAGCCTGCCTGATTTAACGCGATGGCGAAATGGAAGAATCCAAACAGCAACTGATCGACAAACTGCTGACCCGTATTCAGGGTGTGATCCAGGCGCGCGAATCAAAATACATCATGATGCACGCACCGACCGAACGTCTGGATGAAGTCATCGCCctgctg                                                                                                      >  1:25378/1‑251 (MQ=255)
                                                                                                                                              ggCGGATGCGATTTGCGATCTGGTTTCCACCGGTGCCACGCTGGAAGCTAACGGCCTGCGCGAAGTCGAAGTTATCTATCGCTCGAAAGCCTGCCTGATTCAACGCGATGGCGAAATGGAAGAATCCAAACAGCAACTGATCGACAAACTGCTGACCCGTATTCAGGGTGTGATCCAGGCGCGCGAATCAAAATACATCATGATGCACGCACCGACCGAACGTCTGGATGAAGTCATCGCCCTGCTGcc                                                                                                    <  1:354116/249‑1 (MQ=255)
                                                                                                                                                       gATTTGCGATCTGGTTTCCACCGGTGCCACGCTGGAAGCTAACGGCCTGCGCGAAGTCGAAGTTATCTATCGCTCGAAAGCCTGCCTGATTTAACGCGATGGCGAAATGGAAGAATCCAAACAGCAACTGATCGACAAACTGCTGACCCGTATTCAGGGTGTGATCCAGGCGCGCGAATCAAAATACATCATGATGCACGCACCGACCGAACGTCTGGATGAAGTCATCGCCCTGCTGCCAGGTGCCGAAc                                                                                         <  1:150459/251‑1 (MQ=255)
                                                                                                                                                                        ccACCGGTGCCACGCTGGAAGCTAACGGCCTGCGCGAAGTCGAAGTTATCTATCGCTCGAAAGCCTGCCTGATTTAACGCGATGGCGAAATGGAAGAATCCAAACAGCAACTGATCGACAAACTGCTGACCCGTATTCAGGGTGTGATCCAGGCGCGCGAATCAAAATACATCATGATGCACGCACCGACCGAACGTCTGGATGAAGTCATCGCCCTGCTGCCAGGTGCCGAACGCCCAACTATTCTGCCg                                                                        >  1:226891/1‑251 (MQ=255)
                                                                                                                                                                            cGGTGCCACGCTGGAAGCTAACGGCCTGCGCGAAGTCGAAGTTATCTATCGCTCGAAAGCCTGCCTGATTTAACGCGATGGCGAAATGGAAGAATCCAAACAGCAACTGATCGACAAACTGCTGACCCGTATTCAGGGTGTGATCCAGGCGCGCGAATCAAAATACATCATGATGCACGCACCGACCGAACGTCTGGATGAAGTCATCGCCCTGCTGCCAGGTGCCGAACGCCCAACTATTCTGCCGCTgg                                                                    <  2:238247/251‑1 (MQ=255)
                                                                                                                                                                              gTGCCACGCTGGAAGCTAACGGCCTGCGCGAAGTCGAAGTTATCTATCGCTCGAAAGCCTGCCTGATTTAACGCGATGGCGAAATGGAAGAATCCAAACAGCAACTGATCGACAAACTGCTGACCCGTATTCAGGGTGTGATCCAGGCGCGCGAATCAAAATACATCATGATGCACGCACCGACCGAACGTCTGGATGAAGTCATCGCCCTGCTGCCAGGTGCCGAACGCCCAACTATTCTGCCGCTGGCg                                                                  <  1:341661/251‑1 (MQ=255)
                                                                                                                                                                                               aaCGGCCTGCGCGAAGTCGAAGTTATCTATCGCTCGAAAGCCTGCCTGATTTAACGCGATGGCGAAATGGAAGAATCCAAACAGCAACTGATCGACAAACTGCTGACCCGTATTCAGGGTGTGATCCAGGCGCGCGAATCAAAATACATCATGATGCACGCACCGACCGAACGTCTGGATGAAGTCATCGCCCTGCTGCCAGGTGCCGAACGCCCAACTATTCTGCCGCTGGCGGGTGACCAACAGCGCGt                                                 >  2:123002/1‑251 (MQ=255)
                                                                                                                                                                                                        cgcgAAGTCTAAGTTATCTATCGCTCGAAAGCCTGCCTGATTTAACGCGATGGCGAAATGGAAGAATCCAAACAGCAACTGATCGAAAAACTGCTGACCCGTATTCAGGGTGTGATCCAGGCGCGCGAATCAAAATACATCATGATGCACGCACCGACCGAACGTCTGGATGAAGTCATCGCCCTGCTGCCAGGTGCCGAACGCCCAACTCTTCTGCCGCTGGCGGCTGACCAACAGCGCCAAGCGATGc                                         >  2:137436/1‑250 (MQ=255)
                                                                                                                                                                                                        cgcgAAGTCGAAGTTATCTATCGCTCGAAAGCCTGCCTGATTTAACGCGATGGCGAAATGGAAGAATCCAAACAGCAACTGATCGACAAACTGCTGACCCGTATTCAGGGTGTGATCCAGGCGCGCGAATCAAAATACATCATGATGCACGCACCGACCGAACGTCTGGATGAAGTCATCGCCCTGCTGCCAGGTGCCGAACGCCCAACTATTCTGCCGCTGGCGGGTGACCAACAGCGCGTAGCGATGc                                         >  1:303607/1‑250 (MQ=255)
                                                                                                                                                                                                            aagtcgaagtTATCTATCGCTCGAAAGCCTGCCTGATTTAACGCGATGGCGAAATGGAAGAATCCAAACAGCAACTGATCGACAAACTGCTGACCCGTATTCAGGGTGTGATCCAGGCGCGCGAATCAAAATACATCATGATGCACGCACCGACCGAACGTCTGGATGAAGTCATCGCCCTGCTGCCAGGTg                                                                                               <  2:165116/192‑1 (MQ=255)
                                                                                                                                                                                                            aagtcgaagtTATCTATCGCTCGAAAGCCTGCCTGATTTAACGCGATGGCGAAATGGAAGAATCCAAACAGCAACTGATCGACAAACTGCTGACCCGTATTCAGGGTGTGATCCAGGCGCGCGAATCAAAATACATCATGATGCACGCACCGACCGAACGTCTGGATGAAGTCATCGCCCTGCTGCCAGGTg                                                                                               >  1:165116/1‑192 (MQ=255)
                                                                                                                                                                                                                      tatctatcGCTCGAAAGCCTGCCTGATTTAACGCGATGGCGAAATGGAAGAATCCAAACAGCAACTGATCGACAAACTGCTGACCCGTATTCAGGGTGTGATCCAGGCGCGCGAATCAAAATACATCATGATGCACGCACCGACCGAACGTCTGGATGAAGTCATCGCCCTGCTGCCAGGTGCCGAACGCCCAACTATTCTGCCGCTGGCGGGTGACCCACAGCGCGTGGCGATGCACATggcccggagcg                        >  2:360638/1‑242 (MQ=255)
                                                                                                                                                                                                                          tatcGCTCGAAAGCCTGCCTGATTTAACGCGATGGCGAAATGGAAGAATCCAAACAGCAACTGATCGACAAACTGCTGACCCGTATTCAGGGTGTGATCCAGGCGCGCGAATCAAAATACATCATGATGCACGCACCGACCGAACGTCTGGATGAAGTCATCGCCCTGCTGCCAGGTGCCGAACGCCCAACTATTCTGCCGCTGGCGGGTGACCAACAGCGCGTAGCGATGCACATGGTCAGCAGCGAAAc                    >  1:67229/1‑251 (MQ=255)
                                                                                                                                                                                                                           atcGCTCGAAAGCCTGCCTGATTCAACGCGATGGCGAAATGGAAGAATCCAAACAGCAACTGATCGACAAACTGCTGACCCGTATTCAGGGTGTGATCCAGGCGCGCGAATCAAAATACATCATGATGCACGCACCGACCGaa                                                                                                                                 <  1:160690/143‑1 (MQ=255)
                                                                                                                                                                                                                           atcGCTCGAAAGCCTGCCTGATTCAACGCGATGGCGAAATGGAAGAATCCAAACAGCAACTGATCGACAAACTGCTGACCCGTATTCAGGGTGTGATCCAGGCGCGCGAATCAAAATACATCATGATGCACGCACCGACCGaa                                                                                                                                 >  2:160690/1‑143 (MQ=255)
                                                                                                                                                                                                                                tCGAAAGCCTGCCTGATTTAACGCGATGGCGAAATGGAAGAATCCAAACAGCAACTGATCGACAAACTGCTGACCCGTATTCAGGGTGTGATCCAGGCGCGCGAATCAAAATACATCATGATGCACGCACCGACCGAACGTCTGGATGAAGTCATCGCCCTGCTGCCAGGTGCCGAACGCCCAACTATTCTTCCGCTGGCGGGTGACACACAGCGCGAAGCGCTGCGACTGGTCAGCAGCGaacctctggt              >  2:354201/1‑243 (MQ=255)
                                                                                                                                                                                                                                      gcctgcctGATTTAACGCGATGGCGAAATGGAAGAATCCAAACAGCAACTGATCGACAAACTGCTGACCCGTATTCAGGGTGTGATCCAGGCGCGCGAATCAAAATACATCATGATGCACGCACCGACCGAACGTCTGGATGAAGTCATCGCCCTGCTGCCAGGTGCCGAACGCCCAACTATTCTGCCGCTGGCGGGTGACCAACAGCGCGTAGCGATGCACATGGTCAGCAGCGAAACCCTGTTCTGGGa        <  2:366277/251‑1 (MQ=255)
                                                                                                                                                                                                                                            ctgATTTAACGCGATGGCGAAATGGAAGAATCCAAACAGCAACTGATCGACAAACTGCTGACCCGTATTCAGGGTGTGATCCAGGCGCGCGAATCAAAATACATCATGATGCACGCACCGACCGAACGTCTGGATGAAGTCATCGCCCTGCTGCCAGGTGCCGAACGCCCAACTATTCTGCCGCTGGCGGGTGACCAACAGCGCGGAGCGATGACAATGGTCAGCAGCGAAACCCTGTTCTGGGAAACCAt  >  2:410303/1‑251 (MQ=255)
                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                      
GGGACGGTCCGCTCTCCTTAAACGGTAAACGTATCGCCACCTCTTATCCTCACCTGCTCAAGCGTTATCTCGACCAGAAAGGCATCTCTTTTAAATCCTGCTTACTGAACGGTTCTGTTGAAGTCGCCCCGCGTGCCGGACTGGCGGATGCGATTTGCGATCTGGTTTCCACCGGTGCCACGCTGGAAGCTAACGGCCTGCGCGAAGTCGAAGTTATCTATCGCTCGAAAGCCTGCCTGATTCAACGCGATGGCGAAATGGAAGAATCCAAACAGCAACTGATCGACAAACTGCTGACCCGTATTCAGGGTGTGATCCAGGCGCGCGAATCAAAATACATCATGATGCACGCACCGACCGAACGTCTGGATGAAGTCATCGCCCTGCTGCCAGGTGCCGAACGCCCAACTATTCTGCCGCTGGCGGGTGACCAACAGCGCGTAGCGATGCACATGGTCAGCAGCGAAACCCTGTTCTGGGAAACCAT  >  NC_000913/2090538‑2091024

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