Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,613,276 G→T 7.9% intergenic (‑25/+39) glsB ← / ← sad glutaminase 2/NAD(+)‑dependent succinate semialdehyde dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,613,2760GT 7.9% 377.7 / 16.0 114intergenic (‑25/+39)glsB/sadglutaminase 2/NAD(+)‑dependent succinate semialdehyde dehydrogenase
Reads supporting (aligned to +/- strand):  ref base G (53/52);  new base T (3/6);  total (56/58)
Fisher's exact test for biased strand distribution p-value = 4.90e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.62e-01

Only 100 of 115 total aligned reads displayed.

CCCAATCGGGAACCGTCTACTGTAGCCAGCGCCGGAATATAATCCGCGACTTTACCCTGACCAATGAGCGGCCGCACTTGCCGCAAGATGTTTTCTAAAATTGCATTATCCATGGCGACTGCCACTTTCTACTCCTGGACCGCAGGTCTGAAAAGACCTGCGAGTATATCAGAGCTGAATATGTCGCGTCAGATCCGGTCTTTCCACACCGTCTGGATATTACAGAATTCGTGTAAGCCGAAATGGGAAAGCTCACGACCAAAGCCACTCTTTTTCACGCCACCAAAGGCCACTCGCG  >  NC_000913/1613127‑1613424
                                                                                                                                                     |                                                                                                                                                    
cccAATCGGGAACCGTCTACTGTAGCCAGCGCCGGAATATAATCCGCGACTTTACCCTGACCAATGAGCGGCCGCACTTGCCGCAAGATGTTTTCTAAAATTGCATTATCCATGGCGACTGCCACTTTCTACTCCTGGACCGCAGGTCTg                                                                                                                                                      <  1:1656132/150‑1 (MQ=255)
  cAATCGGGAACCGTCTACTGTAGCCAGCGCCGGAATATAATCCGCGACTTTACCCTGACCAATGAGCGGCCGCACTTGCCGCAAGATGTTTTCTAAAATTGCATTATCCATGGCGACTGCCACTTTCTACTCCTGGACCGCAGGTCTGaa                                                                                                                                                    >  1:1031103/1‑150 (MQ=255)
  cAATCGGGAACCGTCTACTGTAGCCAGCGCCGGAATATAATCCGCGACTTTACCCTGACCAATGAGCGGCCGCACTTGCCGCAAGATGTTTTCTAAAATTGCATTATCCATGGCGACTGCCACTTTCTACTCCTGGACCGCAGGTCTGaa                                                                                                                                                    >  1:1207609/1‑150 (MQ=255)
   aaTCGGGAACCGTCTACTGTAGCCAGCGCCGGAATATAATCCGCGACTTTACCCTGACCAATGAGCGGCCGCACTTGCCGCAAGATGTTTTCTAAAATTGCATTATCCATGGCGACTGCCACTTTCTACTCCTGGACCGCAGGTCTGaaa                                                                                                                                                   <  1:1365705/150‑1 (MQ=255)
     tCGGGAACCGTCTACTGTAGCCAGCGCCGGAATATAATCCGCGACTTTACCCTGACCAATGAGCGGCCGCACTTGCCGCAAGATGTTTTCTAAAATTGCATTATCCATGGCGACTGCCACTTTCTACTCCTGGACCGCAGGTCTGAAAAg                                                                                                                                                 <  1:469262/150‑1 (MQ=255)
      cGGGAACCGTCTACTGTAGCCAGCGCCGGAATATAATCCGCGACTTTACCCTGACCAATGAGCGGCCGCACTTGCCGCAAGATGTTTTCTAAAATTGCATTATCCATGGCGACTGCCACTTTCTACTCCTGGACCGCAGGTCTGAAAAGa                                                                                                                                                <  2:168925/150‑1 (MQ=255)
        ggAACCGTCTACTGTAGCCAGCGCCGGAATATAATCCGCGACTTTACCCTGACCACTGAGCGGCCGCAGTTGCCGCAAGATGTTTTCTAAATTTGCATTAACCATGGCGCCGTCCACTTTCTACTCCTGGACCGCAGGTCTGAAAAGAcc                                                                                                                                              <  2:1432731/150‑1 (MQ=255)
           aCCGTCTACTGTAGCCAGCGCCGGAATATAATCCGCGACTTTACCCTGACCAATGAGCGGCCGCACTTGCCGCAAGATGTTTTCTAAAATTGCATTATCCATGGCGACTGCCACTTTCTACTCCTGGACCGCAGGTCTTTTCAGACCTGc                                                                                                                                           <  1:192126/150‑1 (MQ=255)
           aCCGTCTACTGTAGCCAGCGCCGGAATATAATCCGCGACTTTACCCTGACCAATGAGCGGCCGCACTTGCCGCAAGATGTTTTCTAAAATTGCATTATCCATGGCGACTGCCACTTTCTACTCCTGGACCGCAGGTCTGAAAAGACCTGc                                                                                                                                           <  1:1532200/150‑1 (MQ=255)
            ccGTCTACTGTAGCCAGCGCCGGAATATAATCCGCGACTTTACCCTGACCAATGAGCGGCCGCACTTGCCGCAAGATGTTTTCTAAAATTGCATTATCCATGGCGACTGCCACTTTCTACTCCTGGACCGCAGGTCTTTTCAGACCTGCg                                                                                                                                          <  1:661664/150‑1 (MQ=255)
             cGTCTACTGTAGCCAGCGCCGGAATATAATCCGCGACTTTACCCTGACCAATGAGCGGCCGCACTTGCCGCAAGATGTTTTCTAAAATTGCATTATCCATGGCGACTGCCACTTTCTACTCCTGGACCGCAGGTCTTTTCAGACCTGCgg                                                                                                                                         <  1:373147/150‑2 (MQ=255)
             cGTCTACTGTAGCCAGCGCCGGAATATAATCCGCGACTTTACCCTGACCAATGAGCGGCCGCACTTGCCGCAAGATGTTTTCTAAAATTGCATTATCCATGGCGACTGCCACTTTCTACTCCTGGACCGCAGGTCTTTTCAGACCTGCgg                                                                                                                                         <  1:988112/150‑2 (MQ=255)
              gTCTACTGTAGCCAGCGCCGGAATATAATCCGCGACTTTACGCTGACCTATGAGCGGCCGCACTTGCCGCAAGATGTTTTCTAAAATTGCATTATCCATGGCGACTGCCACTTTCTACTCCTGGACCGCAGGTCTTTTCAGACCTGCggt                                                                                                                                        <  2:247843/150‑3 (MQ=255)
              gTCTACTGTAGCCAGCGCCGGAATATAATCCGCGACTTTACCCTGACCAATGAGCGGCCGCACTTGCCGCAAGATGTTTTCTAAAATTGCATTATCCATGGCGACTGCCACTTTCTACTCCTGGACCGCAGGTCTGAAAAGACCTGCGAg                                                                                                                                        <  1:369316/150‑1 (MQ=255)
               tcttctGTCGCGAGCGCCGGAATATAATCCGCGACTTTACCCTGACCAATTAGCGGCCGCACTTGCCGAAAGATGTTTTCTAAAATTGCATTATCCATGGCGACTGCCACTTTCTACTCCTGGACCGCAGGTCTGAAAAGACCTGCGAGt                                                                                                                                       <  2:648954/150‑1 (MQ=255)
                 tACTGTAGCCAGCGCCGGAATATAATCCGCGACTTTACCCTGACCAATGAGCGGCCGCACTTGCCGCAAGATGTTTTCTAAAATTGCATTATCCATGGCGACTGCCACTTTCTACTCCTGGACCGCAGGTCTGAAAAGACCTGCGAGtat                                                                                                                                     >  1:1001733/1‑150 (MQ=255)
                  aCTGTAGCCAGCGCCGGAATATAATCCGCGACTTTACCCTGACCAATGAGCGGCCGCACTTGCCGCAAGATGTTTTCTAAAATTGCATTATCCATGGCGACTGCCACTTTCTACTCCTGGACCGCAGGTCTGAAAAGACCTGCGAGtata                                                                                                                                    <  1:223720/150‑1 (MQ=255)
                   cTGTAGCCAGCGCCGGAATATAATCCGCGTCTTTACCCTGACCAATGAGCGGCCGCACTTGCCGCAAGATGTTTTCTAAAATTGCATTATCCATGGCGACTGCCACTTTCTACTCCTGGACCGCAGGTCTGAAAAGACCTGCGAGTatat                                                                                                                                   <  2:1128591/150‑1 (MQ=255)
                   cTGTAGCCAGCGCCGGAATATAATCCGCGACTTTACCCTGACCAATGAGCGGCCGCACTTGCCGCAAGATGTTTTCTAAAATTGCATTATCCATGGCGACTGCCACTTTCTACTCCTGGACCGCAGGTCTTTTCAGACCTGCggt                                                                                                                                        <  1:1260440/145‑3 (MQ=255)
                   cTGTAGCCAGCGCCGGAATATAATCCGCGACTTTACCCTGACCAATGAGCGGCCGCACTTGCCGCAAGATGTTTTCTAAAATTGCATTATCCATGGCGACTGCCACTTTCTACTCCTGGACCGCAGGTCTTTTCAGACCTGCggt                                                                                                                                        >  2:1260440/1‑143 (MQ=255)
                   cTGTAGCCAGCGCCGGAATATAATCCGCGACTTTACCCTGACCAATGAGCGGCCGCACTTGCCGCAAGATGTTTTCTAAAATTGCATTATCCATGGCGACTGCCACTTTCTACTCCTGGACCGCAGGTCTGAAAAGACCTGCGAGTatat                                                                                                                                   >  1:325551/1‑150 (MQ=255)
                    tGTAGCCAGCGCCGGAATATAATCCGCGACTTTACCCTGACCAATGAGCGGCCGCACTTGCCGCAAGATGTTTTCTAAAATTGCATTATCCATGGCGACTGCCACTTTCTACTCCTGGACCGCAGGTCTGAAAAGACCTGCGAGTATATc                                                                                                                                  <  1:126339/150‑1 (MQ=255)
                    tGTAGCCAGCGCCGGAATATAATCCGCGACTTTACCCTGACCAATGAGCGGCCGCACTTGCCGCAAGATGTTTTCTAAAATTGCATTATCCATGGCGACTGCCACTTTCTACTCCTGGACCGCAGGTCTGAAAAGACCTGCGAGTATATc                                                                                                                                  >  1:557912/1‑150 (MQ=255)
                    tGTAGCCAGCGCCGGAATATAATCCGCGACTTTACCCTGACCAATGAGCGGCCGCACTTGCCGCAAGATGTTTTCTAAAATTGCATTATCCATGGCGACTGCCACTTTCTACTCCTGGACCGCAGGTCTGAAAAGACCTGCGAGTATATc                                                                                                                                  >  2:149645/1‑150 (MQ=255)
                    tGTAGCCAGCGCAGGAATATAATCCGCGACTTTACCCTGACCAATGAGCGGCCGCACTTGCCGCAAGATGTTTTCTAAAATTGCATTATCCATGGCGACTGCCACTTTCTACTCCTGGACCGCAGGTCTGAAAAGACCTGCGAGTATATc                                                                                                                                  <  1:461191/150‑1 (MQ=255)
                      tAGCCAGCGCCGGAATATAATCCGCGACTTTACCCTGACCAATGAGCGGCCGCACTTGCCGCAAGATGTTTTCTAAAATTGCATTATCCATGGCGACTGCCACTTTCTACTCCTGGACCGCAGGTCTGAAAAGACCTGCGAGTATATCag                                                                                                                                <  2:1623006/150‑1 (MQ=255)
                        gCCAGCGCCGGAATATAATCCGCGACTTTACCCTGACCAATGAGCGGCCGCACTTGCCGCAAGATGTTTTCTAAAATTGCATTATCCATGGCGACTGCCACTTTCTACTCCTGGACCGCAGGTCTGAAAAGACCTGCGAGTATATCagag                                                                                                                              >  1:252719/1‑150 (MQ=255)
                        gCCAGCGCCGGAATATAATCCGCGACTTTACCCTGACCAATGAGCGGCCGCACTTGCCGCAAGATGTTTTCTAAAATTGCATTATCCATGGCGACTGCCACTTTCTACTCCTGGACCGCAGGTCTGAAAAGACCTGCGAGTATATCagag                                                                                                                              <  2:557912/150‑1 (MQ=255)
                          cAGCTCCGGAATATAATCCGCGACTTTACCCTGACCAATGAGCGGCCGCACTTGCCGCAAGATGTTTTCTAAAATTGCATTATCCATGGCGACTGCCACTTTCTACTCCTGGACCGCAGGTCTGAAAAGACCTGCGAGTATATCAGAGCt                                                                                                                            <  2:1207609/150‑1 (MQ=255)
                           aGCGCCGGAATATAATCCGCGACTTTACCCTGACCAATGAGCGGCCGCACTTGCCGCAAGATGTTTTCTAAAATTGCATTATCCATGGCGACTGCCACTTTCTACTCCTGGACCGCAGGTCTGAAAAGACCTGCGAGTATATCAGAGCTg                                                                                                                           <  2:451359/150‑1 (MQ=255)
                           aGCGCCGGAATATAATCCGCGACTTTACCCTGACCAATGAGCGGCCGCACTTGCCGCAAGATGTTTTCTAAAATTGCATTATCCATGGCGACTGCCACTTTCTACTCCTGGACCGCAGGTCTGAAAAGACCTGCGAGTATATCAGAGCTg                                                                                                                           <  2:1259382/150‑1 (MQ=255)
                              gcCGGAATATAATCCGCGACTTTACCCTGACCAATGAGCGGCCGCACTTGCCGCAAGATGTTTTCTAAAATTGCATTATCCATGGCGACTGCCACTTTCTACTCCTGGACCGCAGGTCTGAAAAGACCTGCGAGTATATCAGAGCTGAat                                                                                                                        <  2:1212300/150‑1 (MQ=255)
                                     tataATCCGCGACTTTACCCTGACCAATGAGCGGCCGCACTTGCCGCAAGATGTTTTCTAAAATTGCATTATCCATGGCGACTGCCACTTTCTACTCCTGGACCGCAGGTCTGAAAAGACCTGCGAGTATATCAGAGCTGAATATGTcgc                                                                                                                 >  1:1245182/1‑150 (MQ=255)
                                        aaTCCGCGACTTTACCCTGACCAATGAGCGGCCGCACTTGCCGCAAGATGTTTTCTAAAATTGCATTATCCATGGCGACTGCCACTTTCTACTCCTGGACCGCAGGTCTGAAAAGACCTGCGAGTATATCAGAGCTGAATATGTCGCGTc                                                                                                              <  1:243437/150‑1 (MQ=255)
                                         aTCCGCGACTTTACCCTGACCAATGAGCGGCCGCACTTGCCGCAAGATGTTTTCTAAAATTGCATTATCCATGGCGACTGCCACTTTCTACTCCTGGACCGCAGGTCTGAAAAGACCTGCGAGTATATCAGAGCTGAATATGTCGCGTCa                                                                                                             <  1:382885/150‑1 (MQ=255)
                                         aTCCGCGACTTTACCCTGACCAATGAGCGGCCGCACTTGCCGCAAGATGTTTTCTAAAATTGCATTATCCATGGCGACTGCCACTTTCTACTCCTGGACCGCAGGTCTGAAAAGACCTGCGAGTATATCAGAGCTGAATATGTCGCGTCa                                                                                                             >  2:135191/1‑150 (MQ=255)
                                          tCCGCGACTTTACCCTGACCAATGAGCGGCCGCACTTGCCGCAAGATGTTTTCTAAAATTGCATTATCCATGGCGACTGCCACTTTCTACTCCTGGACCGCAGGTCTGAAAAGACCTGCGAGTATATCAGAGCTGAat                                                                                                                        >  1:507701/1‑138 (MQ=255)
                                          tCCGCGACTTTACCCTGACCAATGAGCGGCCGCACTTGCCGCAAGATGTTTTCTAAAATTGCATTATCCATGGCGACTGCCACTTTCTACTCCTGGACCGCAGGTCTGAAAAGACCTGCGAGTATATCAGAGCTGAat                                                                                                                        <  2:507701/138‑1 (MQ=255)
                                           ccGCGACTTTACCCTGACCAATGAGCGGCCGCACTTGCCGCAAGATGTTTTCTAAAATTGCATTATCCATGGCGACTGCCACTTTCTACTCCTGGACCGCAGGTCTGAAAAGACCTGCGAGTATATCAGAGCTGAATATGTCGCGTCAGa                                                                                                           >  2:1571955/1‑150 (MQ=255)
                                            cgcgACTTTACCCTGACCAATGAGCGGCCGCACTTGCCGCAAGATGTTTTCTAAAATTGCATTATCCATGGCGACTGCCACTTTCTACTCCTGGACCGCAGGTCTGAAAAGAcc                                                                                                                                              <  1:993692/114‑1 (MQ=255)
                                              cgACTTTACCCTGACCAATGAGCGGCCGCACTTGCCGCAAGATGTTTTCTAAAATTGCATTATCCATGGCGACTGCCACTTTCTACTCCTGGACCGCAGGTCTGAAAAGACCTGCGAGTATATCAGAGCTGAATATGTCGCGTCAGATcc                                                                                                        <  1:1167094/150‑1 (MQ=255)
                                                aCTTTACCCTGACCAATGAGCGGCCGCACTTGCCGCAAGATGTTTTCTAAAATTGCATTATCCATGGCGACTGCCACTTTCTACTCCTGGACCGCAGGTCTGAAAAGACCTGCGAGta                                                                                                                                      >  1:1234531/1‑118 (MQ=255)
                                                aCTTTACCCTGACCAATGAGCGGCCGCACTTGCCGCAAGATGTTTTCTAAAATTGCATTATCCATGGCGACTGCCACTTTCTACTCCTGGACCGCAGGTCTGAAAAGACCTGCGAGta                                                                                                                                      <  2:1234531/118‑1 (MQ=255)
                                                aCTTTACCCTGACCAATGAGCGGCCGCACTTGCCGCAAGATGTTTTCTAAAATTGCATTATCCATGGCGACTGCCACTTTCTACTCCTGGACCGCAGGTCTGAAAAGACCTGCGAGTATATCAGAGCTGAATATGTCGCGTCAGATCCgg                                                                                                      >  2:513543/1‑150 (MQ=255)
                                                    tACCCTGACCAATGAGCGGCCGCACTTGCCGCAAGATGTTTTCTAAAATTGCATTATCCATGGCGACTGCCACTTTCTACTCCTGGACCGCAGGTCTGAAAAGACCTGCGAGTATATCAGAGCTGAATATGTCGCGTCAGATCCGGTCtt                                                                                                  <  2:274402/150‑1 (MQ=255)
                                                    tACCCTGACCAATGAGCGGCCGCACTTGCCGCAAGATGTTTTCTAAAATTGCATTATCCATGGCGACTGCCACTTTCTACTCCTGGACCGCAGGTCTGAAAAGACCTGCGAGTATATCAGAGCTGAATATGTCGCGTCAGATCCGGTCtt                                                                                                  >  1:289747/1‑150 (MQ=255)
                                                     aCCCTGACCAATGAGCGGCCGCACTTGCCGCAAGATGTTTTCTAAAATTGCATTATCCATGGCGACTGCCACTTTCTACTCCTGGACCGCAGGTCTGAAAAGACCTGCGAGTATATCAGAGCTGAATATGTCGCGTCAGATCCGGTCttt                                                                                                 <  2:1245182/150‑1 (MQ=255)
                                                      cccTGACCAATGAGCGGCCGCACTTGCCGCAAGATGTTTTCTAAAATTGCATTATCCATGGCGCCTGCCACTTTCTACTCCTGGACCGCAGGTCTGAAAAGACCTGCGAGTATATCAGAGCTGAATATGTCGCGTCAGATCCGGTCTTTc                                                                                                >  1:800001/1‑150 (MQ=255)
                                                      cccTGACCAATGAGCGGCCGCACTTGCCGCAAGATGTTTTCTAAAATTGCATTATCCATGGCGACTGCCACTTTCTACTCCTGGACCGCAGGTCTGAAAAGACCTGCGAGTATATCAGAGCTGAATATGTCGCGTCAGATCCGGTCTTTc                                                                                                >  2:538398/1‑150 (MQ=255)
                                                       ccTGACCAATGAGCGGCCGCACTTGCCGCAAGATGTTTTCTAAAATTGCATTATCCATGGCGACTGCCACTTTCTACTCCTGGACCGCAGGTCTGAAAAGACCTGCGAGTATATCAGAGCTGAATATGTCGCGTCAGATCCGGTCTTTcc                                                                                               >  2:605632/1‑150 (MQ=255)
                                                        cTGACCAATGAGCGGCCGCACTTGCCGCAAGATGTTTTCTAAAATTGCATTATCCATGGCGACTGCCACTTTCTACTCCTGGACCGCAGGTCTGAAAAGACCTGCGAGTATATCAGAGCTGAATATGTCGCGTCAGATCCGGTCTTTCca                                                                                              <  2:1001733/150‑1 (MQ=255)
                                                              aaTGAGCGGCCGCACTTGCCGCAAGATGTTTTCTAAAATTGCATTATCCATGGCGACTGTCACTTTCTACTCCTGGACCGCAGGTCTGAAAAGACCTGCGAGTATATCAGAGCTGAATATGTCGCGTCAGATCCGGTCTTTCCACACCGt                                                                                        >  1:37359/1‑150 (MQ=255)
                                                              aaTGAGCGGCCGCACTTGCCGCAAGATGTTTTCTAAAATTGCATTATCCATGGCGACTGCCACTTTCTACTCCTGGACCGCAGGTCTGAAAAGACCTGCGAGTATATCAGAGCTGAATATGTCGCGTCAGATCCGGTCTTTCCACACCGt                                                                                        >  1:1312724/1‑150 (MQ=255)
                                                                   gCGGCCGCACTTGCCGCAAGATGTTTTCTAAAATTGCATTATCCATGGCGACTGCCACTTTCTACTCCTGGACCGCAGGTCTGAAAAGACCTGCGAGTATATCAGAGCTGAATATGTCGCGTCAGATCCGGTCTTTCCACACCGTCTGGa                                                                                   >  1:99329/1‑150 (MQ=255)
                                                                   gCGGCCGCACTTGCCGCAAGATGTTTTCTAAAATTGCATTATCCATGGCGACTGCCACTTTCTACTCCTGGACCGCAGGTCTGAAAAGACCTGCGAGTATATCAGAGCTGAATATGTCGCGTCAGATCCGGTCTTTCCACACCGTCTGGa                                                                                   >  1:93276/1‑150 (MQ=255)
                                                                     ggCCGCACTTTCCGCAAGATGTTTTCTAAAATTGCATTATCCAAGGCGACTGCCACTTTCTACTCCTGGACCGCAGGTCCGAAAAGACCTCCAAGTATATAAGATCAGAATATGTCGCGACAGAACCGATCTTTACCCACCGTCTgcaaa                                                                                 >  2:704006/1‑146 (MQ=255)
                                                                     ggCCGCACTTGCCGCAAGATGTTTTCTAAAATTGCATTATCCATGGCGACTGCCACTTTCTACTCCTGGACCGCAGGTCTGAAAAGACCTGCGAGTATATCAGAGCTGAATATGTCACGTCAGATCCGGTCTTTCCACACCGTCTGGata                                                                                 >  1:343875/1‑150 (MQ=255)
                                                                        cGCACTTGCCGCAAGATGTTTTCTAAAATTGCATTATCCATGGCGACTGCCACTTTCTACTCCTGGACCGCAGGTCTGAAAAGACCTGCGAGTATATCAGAGCTGAATATGTCGCGTCAGATCCGGTCTTTCCACACCGTCTGGATATTa                                                                              <  2:1259618/150‑1 (MQ=255)
                                                                          cACTTGCCGCAAGATGTTTTCTAAAATTGCATTATCCATGGCGACTGCCACTTTCTACTCCTGGACCGCAGGTCTGAAAAGACCTGCGAGTATATCAGAGCTGAATATGTCGCGTCAGATCCGGTCTTTCCACACCGTCTGGATATTACa                                                                            <  1:90130/150‑1 (MQ=255)
                                                                           aCTTGCCGCAAGATGTTTTCTAAAATTGCATTATCCATGGCGACTGCCACTTTCTACTCCTGGACCGCAGGTCTGAAAAGACCTGCGAGTATATCAGAGCTGAATATGTCGCGTCAGATCCGGTCTTTCCACACCGTCTGGATATTACAg                                                                           >  2:1211095/1‑150 (MQ=255)
                                                                            cTTGCCGCAAGATGTTTTCTAAAATTGCATTATCCATGGCGACTGCCACTTTCTACTCCTGGACCGCAGGTCTGAAAAGACCTGCGAGTATATCAGAGCTGAATATGTCGCGTCAGATCCGGTCTTTCCACACCGTCTGGATATTACAGa                                                                          >  1:649946/1‑150 (MQ=255)
                                                                                 cgaaaGATGTTTTCTAAAATTGCATTATCCATGGCTACTTCCACTTTCTACTCCTTGACCGCAGGTCTGAAAAGACCTGCGAGTATATCAGAGCTGAATATGTCGCGTCAGATCCGGTCTTTCCACACCGTCTGGATATTACAGAATTCg                                                                     <  2:845410/147‑1 (MQ=255)
                                                                                 cGCAAGATGTTTTCTAAAATTGCATTATCCATGGCGACTGCCACTTTCTACTCCTGGACCGCAGGTCTGAAAAGACCTGCGAGTATATCAGAGCTGAATATGTCGCGTCAGATCCGGTCTTTCCACACCGTCTGGATATTACAGAATTCg                                                                     >  2:1003662/1‑150 (MQ=255)
                                                                                 cGCAAGATGTTTTCTAAAATTGCATTATCCATGGCGACTGCCACTTTCTACTCCTGGACCGCAGGTCTGAAAAGACCTGCGAGTATATCAGAGCTGAATATGTCGCGTCAGATCCGGTCTTTCCACACCGTCTGGATATTACAGAATTCg                                                                     <  1:538398/150‑1 (MQ=255)
                                                                                   cAAGATGTTTTCTAAAATTGCATTATCCATGGCGACTGCCACTTTCTACTCCTGGACCGCAGGTCTGAAAAGACCTGCGAGTATATCAGAGCTGAATATGTCGCGTCAGATCCGGTCTTTCCACACCGTCTGGATATTACAGAATTCgtg                                                                   <  2:1551200/150‑1 (MQ=255)
                                                                                   cAAGATGTTTTCTAAAATTGCATTATCCATGGCGACTGCCACTTTCTACTCCTGGAACGCAGGTCTGAAAAGACCTGCGAGTATATCAGAGCTGAATATGTCGCGTCAGATCCGGTCTTTCCACACCGTCTGGATTTTACAGAATTGgtg                                                                   <  2:1072362/150‑1 (MQ=255)
                                                                                    aaGATGTTTTCTAAAATTGCATTATCCATGGCGACTGCCACTTTCTACTCCTGGACCGCAGGTCTGAAAAGACCTGCGAGTATATCAGAGCTGAATATGTCGCGTCAGATCCGGTCTTTCCACACCGTCTGGATATTACAGAATTCgtgt                                                                  <  1:1231611/150‑1 (MQ=255)
                                                                                         gTTTTCTAAAATTGCATTATCCATGGCGACTGCCACTTTCTACTCCTGGACCGCAGGTCTGAAAAGACCTGCGAGTATATCAGAGCTGAATATGTCGCGTCAGATCCGGTCTTTCCACACCGTCTGGATATTACAGAATTCGTGTAAGcc                                                             <  1:971508/150‑1 (MQ=255)
                                                                                           tttattaaaTTTTATTATTCATTGCGACTGACACTTTCTACTCCTGGACCGCAGGTCGGAAAAGACCTGCGAGTATATCAGAGCTGAATATGTCGCGTCAGATCCGGTCTTTCCACACCGTCTGGATATTACAGAATTCGTGTAAGCCGa                                                           <  1:704006/144‑1 (MQ=255)
                                                                                             tCTAAAATTGCATTATCCATGGCGACTGCCACTTTCTACTCCTGGACCGCAGGTCTGAAAAGACCTGCGAGTATATCAGAGCTGAATATGTCGCGTCAGATCCGGTCTTTCCACACCGTCTGGATATTACAGAATTCGTGTAAGCCGaaa                                                         >  1:345493/1‑150 (MQ=255)
                                                                                              cTAAAATTGCATTATCCATGGCGACTGCCACTTTCTACTCCTGGACCGCAGGTCTGAAAAGACCTGCGAGTATATCAGAGCTGAATATGTCGCGTCAGATCCGGTCTTTCCACACCGTCTGGATATTACAGAATTCGTGTAAGCCGAAAt                                                        >  1:1629629/1‑150 (MQ=255)
                                                                                               tAAAATTGCATTATCCATGGCGACTGCCACTTTCTACTCCTGGACCGCAGGTCTGAAAAGACCTGCGAGTATATCAGAGCTGAATATGTCGCGTCAGATCCGGTCTTTCCACACCGTCTGGATATTACAGAATTCGTGTAAGCCGAAATg                                                       <  2:1526039/150‑1 (MQ=255)
                                                                                                      gCATTATCCATGGCGACTGCCACTTTCTACTCCTGGACCGCAGGTCTGAAAAGACCTGCGAGTATATCAGAGCTGAATATGTCGCGTCAGATCCGGTCTTTCCACACCGTCTGGATATTACAGAATTCGTGTAAGCCGAAATGGGAAAGc                                                >  2:28163/1‑150 (MQ=255)
                                                                                                             ccATGGCGACTGCCACTTTCTACTCCTGGACCGCAGGTCTGAAAAGACCTGCGAGTATATCAGAGCTGAATATGTCGCGTCAGATCCGGTCTTTCCACACCGTCTGGATATTACAGAATTCGTGTAAGCCGAAATGGGAAAGCTCACGAc                                         <  2:93276/150‑1 (MQ=255)
                                                                                                             ccATGGCGACTGCCACTTTCTACTCCTGGACCGCAGGTCTGAAAAGACCTGCGAGTATATCAGAGCTGAATATGTCGCGTCAGATCCGGTCTTTCCACACCGTCTGGATATTACAGAATTCGTGTAAGCCGAAATGGGAAAGCTCACGAc                                         <  2:99329/150‑1 (MQ=255)
                                                                                                                  gCGACTGCCACTTTCTACTCCTGGACCGCAGGTCTGAAAAGACCTGCGAGTATATCAGAGCTGAATATGTCGCGTCAGATCCGGTCTTTCCACACCGTCTGGATATTACAGAATTCGTGTAAGCCGAAATGGGAAAGCTCACGACCAAAg                                    <  2:1531043/150‑1 (MQ=255)
                                                                                                                  gCGACTGCCACTTTCTACTCCTGGACCGCAGGTCTGAAAAGACCTGCGAGTATATCAGAGCTGAATATGTCGCGTCAGATCCGGTCTTTCCACACCGTCTGGATATTACAGAATTCGTGTAAGCCGAAATGGGAAAGCTCACGACCAAAg                                    >  1:1024462/1‑150 (MQ=255)
                                                                                                                  gCGACTGCCACTTTCTACTCCTGGACCGCAGGTCTGAAAAGACCTGCGAGTATATCAGAGCTGAATATGTCGCGTCAGATCCGGTCTTTCCACACCGTCTGGATATTACAGAATTCGTGTAAGCCGAAATGGGAAAGCTCACGACCAAAg                                    <  2:817609/150‑1 (MQ=255)
                                                                                                                    gACTGCCACTTTCTACTCCTGGACCGCAGGTCTGAAAAGACCTGCGAGTATATCAGAGCTGAATATGTCGCGTCAGATCCGGGCTTTCCACACCGTCTGGATATTACAGAATTCGTGTAAGCCGAAATGGGAAAGCGCACGACCAAAGcc                                  <  2:660645/150‑1 (MQ=255)
                                                                                                                        gCCACTTTCTACTCCTGGACCGCAGGTCTGAAAAGACCTGCGAGTATATCAGAGCTGAATATGTCGCGTCAGATCCGGTCTTTCCACACCGTCTGGATATTACAGAATTCGTGTAAGCCGAAATGGGAAAGCTCACGACCAAAGCCActc                              >  1:1665748/1‑150 (MQ=255)
                                                                                                                         ccACTTTCTACTCCTGGACCGCAGGTCTGAAAAGACCTGCGAGTATATCAGAGCTGAATATGTCGCGTCAGATCCGGTCTTTCCACACCGTCTGGATATTACAGAATTCGTGTAAGCCGAAATGGGAAAGCTCACGACCAAAGCCActct                             >  1:805371/1‑150 (MQ=255)
                                                                                                                          cACTTTCTACTCCTGGACCGCAGGTCTGAAAAGACCTGCGAGTATATCAGAGCTGAATATGTCGCGTCAGATCCGGTCTTTCCACACCGTCTGGATATTACAGAATTCGTGTAAGCCGAAATGGGAAAGCTCACGACCAAAGCCACTCtt                            >  2:1265831/1‑150 (MQ=255)
                                                                                                                           aCTTTCTACTCCTGGACCGCAGGTCTGAAAAGACCTGCGAGTATATCAGAGCTGAATATGTCGCGTCAGATCCGGTCTTTCCACACCGTCTGGATATTACAGAATTCGTGTAAGCCGAAATGGGAAAGCTCACGACCAAAGCCACTCttt                           >  2:1589763/1‑150 (MQ=255)
                                                                                                                           aCTTTCTACTCCTGGACCGCAGGTCTGAAAAGACCTGCGAGTATATCAGAGCTGAATATGTCGCGTCAGATCCGGTCTTTCCACACCGTCTGGATATTACAGAATTCGTGTAAGCCGAAATGGGAAAGCTCACGACCAAAGCCACTCttt                           >  1:239943/1‑150 (MQ=255)
                                                                                                                                cTACTCCTGGACCGCAGGTCTGAAAAGACCTGCGAGTATATCAGAGCTGAATATGTCGCGTCAGATCCGGTCTTTCCACACCGTCTGGATATTACAGAATTCGTGTAAGCCGAAATGGGAAAGCTCACGACCAAAGCCACTCTTTTTCAc                      >  2:267755/1‑150 (MQ=255)
                                                                                                                                cTACTCCTGGACCGCAGGTCTGAAAAGACCTGCGAGTATATCAGAGCTGAATATGTCGCGTCAGATCCGGTCTTTCCACACCGTCTGGATATTACAGAATTCGTGTAAGCCGAAATGGGAAAGCTCACGACCAAAGCCACTCTTTTTCAc                      >  1:969206/1‑150 (MQ=255)
                                                                                                                                 tACTCCTGGACCGCAGGTCTGAAAAGACCTGCGAGTATATCAGAGCTGAATATGTCGCGTCAGATCCGGTCTTTCCACACCGTCTGGATATTACAGAATTCGTGTAAGCCGAAATGGGAAAGCTCACGACCAAAGCCACTCTTTTTCACg                     <  2:1632337/150‑1 (MQ=255)
                                                                                                                                 tACTCCTGGACCGCAGGTCTGAAAAGACCTGCGAGTATATCAGAGCTGAATATGTCGCGTCAGATCCGGTCTTTCCACACCGTCTGGATATTACAGAATTCGTGTAAGCCGAAATGGGAAAGCTCACGACCAAAGCCACTCTTTTTCACg                     <  2:289747/150‑1 (MQ=255)
                                                                                                                                   cTCCTGGACCGCAGGTCTGAAAAGACCTGCGAGTATATCAGAGCTGAATATGTCGCGTCAGATCCGGTCTTTCCACACCGTCTGGATATTACAGAATTCGTGTAAGCCGAAATGGGAAAGCTCACGACCAAAGCCACTCTTTTTCACGcc                   >  1:1240365/1‑150 (MQ=255)
                                                                                                                                     ccTGGACCGCAGGTCTGAAAAGACCTGCGAGTATATCAGAGCTGAATATGTCGCGTCAGATCCGGTCTTTCCACACCGTCTGGATATTACAGAATTCGTGTAAGCCGAAATGGGAAAGCTCACGACCAAAGCCACTCTTTTTCACGccac                 >  2:818478/1‑150 (MQ=255)
                                                                                                                                      cTGGACCGCAGGTCTGAAAAGACCTGCGAGTATATCAGAGCTGAATATGTCGCGTCAGATCCGGTCTTTCCACACCGTCTGGATATTACAGAATTCGTGTAAGCCGAAATGGGAAAGCTCACGACCAAAGCCACTCTTTTTCACGccacc                >  2:1112646/1‑150 (MQ=255)
                                                                                                                                        ggACCGCAGGTCTGAAAAGACCTGCGAGTATATCAGAGCTGAATATGTCGCGTCAGATCCGGTCTTTCCACACCGTCTGGATATTACAGAATTCGTGTAAGCCGAAATGGGAAAGCTCACGACCAAAGCCACTCTTTTTCACGCCACCaa              >  2:210787/1‑150 (MQ=255)
                                                                                                                                          aCCGCAGGTCTGAAAAGACCTGCGAGTATATCAGAGCTGAATATGTCGCGTCAGATCCGGTCTTTCCACACCGTCTGGATATTACAGAATTCGTGTAAGCCGAAATGGGAAAGCTCACGACCAAAGCCACTCTTTTTCACGCCACCAAAg            <  2:37359/150‑1 (MQ=255)
                                                                                                                                          aCCGCAGGTCTGAAAAGACCTGCGAGTATATCAGAGCTGAATATGTCGCGTCAGATCCGGTCTTTCCACACCGTCTGGATATTACAGAATTCGTGTAAGCCGAAATGGGAAAGCTCACGACCAAAGCCACTCTTTTTCACGCCACCAAAg            >  1:222080/1‑150 (MQ=255)
                                                                                                                                            cGCAGGTCTTTTCAGACCTGCGAGTATATCAGATCTGAATATGTCGCGTCAGATCCGGTCTTTCCACACCGTCTGGATATTACAGAATTCGTGTAAGCCGAAATGGGAAAGCTCACGACCAAAGCCACTCTTTTTCACGCCACCAAAGGc          >  1:578313/1‑150 (MQ=255)
                                                                                                                                            cGCAGGTCTTTTCAGACCTGCGAGTATATCAGAGCTGAATATGTCGCGTCAGATCCGGTCTTTCCACACCGTCTGGATATTACAGAATTCGTGTAAGCCGAAATGGGAAAGCTCACGACCAAAGCCACTCTTTTTCACGCCACCAAAGGc          >  1:101768/1‑150 (MQ=255)
                                                                                                                                             gCAGGTCTGAAAAGACCTGCGAGTATATCAGAGCTGAATATGTCGCGTCAGATCCGGTCTTTCCACACCGTCTGGATATTACAGAATTCGTGTAAGCCGAAATGGGAAAGCTCACGACCAAAGCCACTCTTTTTCACGCCACCAAAGGcc         <  1:1374089/150‑1 (MQ=255)
                                                                                                                                             gCAGGTCTGAAAAGACCTGCGAGTATATCAGAGCTGAATATGTCGCGTCAGATCCGGTCTTTCCACACCGTCTGGATATTACAGAATTCGTGTAAGCCGAAATGGGAAAGCTCACGACCAAAGCCACTCTTTTTCACGCCACCAAAGGcc         <  1:1042985/150‑1 (MQ=255)
                                                                                                                                                ggTCTGAAAAGACATGCGAGTATATCAGAGCTGAATATGTCGCGTCAGATCCGGTCTTTCCACACCGTCTGGATATTACAGAATTCGTGTAAGCCGAAATGGGAAAGCTCACGACCAAAGCCACTCTTTTTCACGCCACCAAAGGCCACt      >  2:25737/1‑150 (MQ=255)
                                                                                                                                                    tGAAAAGACCTGCGAGTATATCAGAGCTGAATATGTCGCGTCAGATCCGGTCTTTCCACACCGTCTGGATATTACAGAATTCGTGTAAGCCGAAATGGGAAAGCTCACGACCAAAGCCACTCTTTTTCACGCCACCAAAGGCCACTcgcg  <  2:1470538/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
CCCAATCGGGAACCGTCTACTGTAGCCAGCGCCGGAATATAATCCGCGACTTTACCCTGACCAATGAGCGGCCGCACTTGCCGCAAGATGTTTTCTAAAATTGCATTATCCATGGCGACTGCCACTTTCTACTCCTGGACCGCAGGTCTGAAAAGACCTGCGAGTATATCAGAGCTGAATATGTCGCGTCAGATCCGGTCTTTCCACACCGTCTGGATATTACAGAATTCGTGTAAGCCGAAATGGGAAAGCTCACGACCAAAGCCACTCTTTTTCACGCCACCAAAGGCCACTCGCG  >  NC_000913/1613127‑1613424

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