Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,932,096 C→G 6.6% intergenic (+37/+19) purT → / ← eda phosphoribosylglycinamide formyltransferase 2/KHG/KDPG aldolase; 2‑dehydro‑3‑deoxy‑phosphogluconate/4‑hydroxy‑2‑ oxoglutarate aldolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,932,0960CG 6.6% 528.6 / 20.9 151intergenic (+37/+19)purT/edaphosphoribosylglycinamide formyltransferase 2/KHG/KDPG aldolase; 2‑dehydro‑3‑deoxy‑phosphogluconate/4‑hydroxy‑2‑ oxoglutarate aldolase
Reads supporting (aligned to +/- strand):  ref base C (69/72);  new base G (6/4);  total (75/76)
Fisher's exact test for biased strand distribution p-value = 5.33e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.79e-01

Only 100 of 155 total aligned reads displayed.

GGAAATTGATGGCAGCCGTCGTCTGGGGGTGGCACTGGCTACTGCAGAGAGTGTTGTTGACGCCATTGAACGCGCGAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGTAATGCGGTCGTAATCGCCCGCTTCCAGCGCATCTGCCGGAACCAGCCAGGAACCACCGATGCACAGCACGCTTTTCAGCGCCAGGTA  >  NC_000913/1931948‑1932243
                                                                                                                                                    |                                                                                                                                                   
ggAAATTGATGGCAGCCGTCGTCTGGGGGTGGCACTGGCTACTGCAGAGAGTGTTGTTGACGCCATTGAACGCGCGAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTc                                                                                                                                                     <  1:244304/149‑1 (MQ=255)
  aaaTTGATGGCAGCCGTCGTCTGGGGGTGGCACTGGCTACTGCAGAGAGTGTTGTTGACGCCATTGAACGCGCGAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTCGa                                                                                                                                                   >  4:1492456/1‑149 (MQ=255)
    aTTGATGGCAGCCGTCGTCTGGGGGTGGCACTGGCTACTGCAGAGAGTGTTGTTGACGCCATTGAACGCGCGAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCTACCAAAATGCCCGATCCTCGATc                                                                                                                                                 >  2:1534196/1‑149 (MQ=255)
      tGATGGCAGCCGTCGTCTGGGGGTGGCACTGGCTACTGCAGAGAGTGTTGTTGACGCCATTGAACGCGCGAGGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTCGATCgg                                                                                                                                               <  3:1562262/149‑1 (MQ=255)
        aTGGCAGCCGTCGTCTGGGGGTGGCACTGGCTACTGCAGAGAGTGTTGTTGACGCCATTGAACGCGCGAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGc                                                                                                                                             >  2:1298018/1‑149 (MQ=255)
             aGCCGTCGTCTGGGGGTGGCACTGGCTACTGCAGAGAGTGTTGTTGACGCCATTGAACGCGCGAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCAtttt                                                                                                                                        >  4:36520/1‑149 (MQ=255)
             aGCCGTCGTCTGGGGGTGGCACTGGCTACTGCAGAGAGTGTTGTTGACGCCATTGAACGCGCGAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCAtttt                                                                                                                                        >  3:486041/1‑149 (MQ=255)
             aGCCGTCGTCTGGGGGTGGCACTGGCTACTGCAGAGAGTGTTGTTGACGCCATTGAACGCGCGAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCAtttt                                                                                                                                        <  3:105435/149‑1 (MQ=255)
               ccGTCGTCTGGGGGTGGCACTGGCTACTGCAGAGAGTGTTGTTGACGCCATTGAACGCGCGAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCAtttt                                                                                                                                        <  3:1505179/147‑1 (MQ=255)
               ccGTCGTCTGGGGGTGGCACTGGCTACTGCAGAGAGTGTTGTTGACGCCATTGAACGCGCGAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCAtttt                                                                                                                                        >  4:1505179/1‑147 (MQ=255)
                cgtcttGTGGGGGTGGCACTGGCTGGTGCAGAGAGTGATGTTGGCCCCATTTAACGCGCGAGGCGCGCCGCAGGACTGGTAATAGTACATTGTTAAACCCGGCAAAATAACGCTACAAAAAAGCCCGATCCTCGATCGGGCATTTTGAc                                                                                                                                     <  3:194186/149‑1 (MQ=255)
                 gtcgtcTGGGGGTGGCACTGGCTACTGCAGAGAGTGTTGTTGACGCCATTGAACGCGCGAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACt                                                                                                                                    <  4:1036846/149‑1 (MQ=255)
                       tGGGGGTGGCACTGGCTACTGCAGAGAGTGTTGTTGACGCCATTGAACGCGCGAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTACa                                                                                                                              <  4:662693/149‑1 (MQ=255)
                       tGGGGGTGGCACTGGCTACTGCAGAGAGTGTTGTTGACGCCATTGAACGCGCGAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTACa                                                                                                                              <  2:1173539/149‑1 (MQ=255)
                        gggggTGGCACTGGCTACTGCAGAGAGTGTTGTTGACGCCATTGAACGCGCGAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTc                                                                                                                                                     >  3:748124/1‑125 (MQ=255)
                          gggTGGCACTGGCTACTGCAGAGAGTGTTGTTGACGCCATTGAACGCGCGAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCt                                                                                                                           <  2:75909/149‑1 (MQ=255)
                              ggCACTGGCTACTGCAGAGAGTGTTGTTGACGCCATTGAACGCGCGAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAACAACGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAgc                                                                                                                       <  1:1121474/149‑1 (MQ=255)
                              ggCACTGGCTACTGCAGAGAGTGTTGTTGACGCCATTGAACGCGCGAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAgc                                                                                                                       <  1:145435/149‑1 (MQ=255)
                                 aCTGGCTACTGCAGAGAGTGTTGTTGACGCCATTGAACGCGCGAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGcc                                                                                                                    >  4:167258/1‑149 (MQ=255)
                                  cTGGCTACTGCAGAGAGTGTTGTTGACGCCATTGAACGCGCGAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCt                                                                                                                   >  1:111252/1‑149 (MQ=255)
                                        aCTGCAGAGAGTTTTGTTGACGCCATTGAACGCGCGAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAGACGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACa                                                                                                             <  1:190892/149‑1 (MQ=255)
                                        aCTGCAGAGAGTGTTGTTGACGCCATTGAACGCGCGAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACa                                                                                                             <  3:878515/149‑1 (MQ=255)
                                            cAGAGAGTGTTGTTGACGCCATTGAACGCGCGAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACt                                                                                                                                    <  2:1394805/122‑1 (MQ=255)
                                             agagagTGTTGTTGACGCCATTGAACGCGCGAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCGAGGATCGGGCAtttttgta                                                                                                                                    <  4:1627238/121‑5 (MQ=39)
                                             agagagTGTTGTTGACGCCATTGAACGCGCGAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTc                                                                                                        <  4:265677/149‑1 (MQ=255)
                                               agagTGTTGTTGACGCCATTGAACGCGCGAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCAc                                                                                                      <  1:1274442/149‑1 (MQ=255)
                                                   tgttgtTGACGCCATTGAACGCGCGAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACtttt                                                                                                                                 <  3:1340959/118‑1 (MQ=255)
                                                   tgttgtTGACGCCATTGAACGCGCGAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACgcgc                                                                                                  >  1:487539/1‑149 (MQ=255)
                                                    gttgttGACGCCATTGAACGCGCGAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGcc                                                                                                 <  1:1603285/149‑1 (MQ=255)
                                                     ttgttgACGCCATTGAACGCGCGAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCa                                                                                                <  1:734522/149‑1 (MQ=255)
                                                     ttgttgACGCCATTGAACGCGCGAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCa                                                                                                <  2:1027545/149‑1 (MQ=255)
                                                        ttgACGCCATTGAACGCGCGAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGc                                                                                                           <  3:585770/135‑1 (MQ=255)
                                                        ttgACGCCATTGAACGCGCGAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGc                                                                                                           >  4:585770/1‑135 (MQ=255)
                                                        ttgACGCCATTGAACGCGCGAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCt                                                                                             >  4:1419430/1‑149 (MQ=255)
                                                          gACGCCATTGAACGCGCGAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTa                                                                                           >  2:743148/1‑149 (MQ=255)
                                                          gACGCCATTGAACGCGCGAAGCACGCCGCCGGACACGTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTa                                                                                           <  3:1154367/149‑1 (MQ=255)
                                                            cGCCATTGAACGCGCGAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGt                                                                                         >  2:560030/1‑149 (MQ=255)
                                                             gCCATTGAACGCGCGAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGTa                                                                                        >  4:1023454/1‑149 (MQ=255)
                                                             gCCATTGAACGCGCGAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGTa                                                                                        >  3:1181004/1‑149 (MQ=255)
                                                               cATTGAACGCGCGAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCAtttt                                                                                                                                        <  4:911217/99‑1 (MQ=255)
                                                                 ttGAACGCGCGAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGTAATGc                                                                                    >  4:1348113/1‑149 (MQ=255)
                                                                 ttGAACGCGCGAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGTAATGc                                                                                    <  2:1657079/149‑1 (MQ=255)
                                                                     aCGCGCGAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCGAGGATCGGGCa                                                                                                                                            >  1:576194/1‑89 (MQ=38)
                                                                     aCGCGCGAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGTAATGCGGTc                                                                                <  2:1258575/149‑1 (MQ=255)
                                                                         gcgAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGTAATGCGGTCGTaa                                                                            <  2:1690894/149‑1 (MQ=255)
                                                                          cgAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGTAATGCGGTCGTAAt                                                                           >  4:636887/1‑149 (MQ=255)
                                                                           gAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTCAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGTAATGCGGTCGTAATc                                                                          >  1:1650830/1‑149 (MQ=255)
                                                                             aGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCTACAAACATGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGTAATGCGGTCGTAATCGc                                                                        <  3:1023454/149‑1 (MQ=255)
                                                                               cACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGTAATGCGGTCGTAATCGccc                                                                      <  2:986117/149‑1 (MQ=255)
                                                                                aCGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGTAATGCGGTCGTAATCGCCCg                                                                     <  1:1497455/149‑1 (MQ=255)
                                                                                    cgccgGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCAc                                                                                                      <  4:480673/112‑1 (MQ=255)
                                                                                    cgccgGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGcc                                                                                                 <  4:1369711/117‑1 (MQ=255)
                                                                                       cgGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGTAATGCGGTCGTAATCGCCCGCTTCCAg                                                              <  3:1492456/149‑1 (MQ=255)
                                                                                            aGGTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGTAATGCGGTCGTAATCGCCCGCTTCCAGCGCAt                                                         >  1:76452/1‑149 (MQ=255)
                                                                                              gTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCtt                                                                                                                  >  4:1275682/1‑90 (MQ=255)
                                                                                              gTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCtt                                                                                                                  <  3:1275682/90‑1 (MQ=255)
                                                                                              gTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGTAATGCGGTCGTaa                                                                            <  3:1386568/128‑1 (MQ=255)
                                                                                               tAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGTAATGCGGTCGTAATCGCCCGCTTCCAGCGCATCt                                                       <  1:74626/148‑1 (MQ=255)
                                                                                               tAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGTAATGCGGTCGTAATCGCCCGCTTCCAGCGCATCTg                                                      >  2:629706/1‑149 (MQ=255)
                                                                                                   aGTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGTAATGCGGTCGTAATCGCCCGCTTCCAGCGCATCTGCCgg                                                  <  1:560030/149‑1 (MQ=255)
                                                                                                      aCAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGTAATGCGGTCGTAATCGCCCGCTTCCAGCGCATCTGCCGGAAc                                               <  2:1375629/149‑1 (MQ=255)
                                                                                                          ggTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGTAATGCGGTCGTAATCGCCCGCTTcc                                                                >  3:1233647/1‑128 (MQ=255)
                                                                                                            ttAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGTAATGCGGTCGTAATCGCCCGCTTCCAGCGCATCTGCCGGAAccagcca                                         >  4:31336/1‑149 (MQ=255)
                                                                                                            ttAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGTAATGCGGTCGTAATCGCCCGCTTCCAGCGCATCTGCCGGAAccagcca                                         <  4:15539/149‑1 (MQ=255)
                                                                                                            ttAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGTAATGCGGTCGTAATCGCCCGCTTCCAGCGCATCTGCCGGAAccagcca                                         >  3:946873/1‑149 (MQ=255)
                                                                                                            ttAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGTAATGCGGTCGTAATCGCCCGCTTCCAGCGCATCTGCCGGAAccagcca                                         <  3:850515/149‑1 (MQ=255)
                                                                                                            ttAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGTAATGCGGTCGTAATCGCCCGCTTCCAGCGCATCTGCCGGAAccag                                            <  1:320675/146‑1 (MQ=255)
                                                                                                            ttAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGTAATGCGGTCGTAATCGCCCGCTTCCAGCGCATCTGCCGGAAccag                                            >  2:320675/1‑146 (MQ=255)
                                                                                                               aaCCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGTAATGCGGTCGTAATCGCCCGCTTCCAGCGCATCTGCCGGAAcca                                             >  1:61330/1‑142 (MQ=255)
                                                                                                               aaCCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGTAATGCGGTCGTAATCGCCCGCTTCCAGCGCATCTGCCGGAAcca                                             <  2:61330/142‑1 (MQ=255)
                                                                                                                aCCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGTAATGCGGTCGTAATCGCCCGCTTCCAGCGCATCTGCCGGAACCAGCCAGGaa                                     <  2:76452/149‑1 (MQ=255)
                                                                                                                    ggCAAAAAAACGGTACAAAAATGCCCGATCCTCGATCGGGCATTGTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGTAATGCGGTCGTAATCGCCCGCTTCCAGCGCATCTGCCGGAACCAGCCAGGAaccac                                 <  3:1544118/149‑1 (MQ=255)
                                                                                                                    ggCAAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGTAATGCGGTCGTAATCGCCCGCTTCCAGCGCATCTGCCGGAACCAGCCAGGAaccac                                 <  3:36520/149‑1 (MQ=255)
                                                                                                                    ggCAAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGTAATGCGGTCGTAATCGCCCGCTTCCAGCGCATCTGCCGGAACCAGCCAGGAaccac                                 >  2:360859/1‑149 (MQ=255)
                                                                                                                        aaaaaaCGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGTAATGCGGTCGTAATCGCCCGCTTCCAGCGCATCTGCCGGAACCAGCCAGGAACCACCGAt                             <  4:599094/149‑1 (MQ=255)
                                                                                                                         aaaaaCGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGTAATGCGGTCGTAATCGCCCGCTTCCAGCGCATCTGCCGGAACCAGCCAGGAACCACCGATg                            >  2:942418/1‑149 (MQ=255)
                                                                                                                         aaaaaCGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGTAATGCGGTCGTAATCGCCCGCTTCCAGCGCATCTGCCGGAACCAGCCAGGAACCACCGATg                            >  3:595527/1‑149 (MQ=255)
                                                                                                                          aaaaCGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGTAATGCGGTCGTAATCGCCCGCTTCCAGCGCATCTGCCGGAACCAGCCAGGAACCACCGATGc                           >  3:1659371/1‑149 (MQ=255)
                                                                                                                            aaCGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGTAATGCGGTCGTAATCGCCCGCTTCCAGCGCATCTGCCGGAACCAGCCAGGAACCACCg                               <  2:381329/143‑1 (MQ=255)
                                                                                                                              cGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGTAATGCGGTGGTAATCGCCCGCTTCCAGCGCATCTGCCGGAACCAGCCAGGAACCACCGATGCACAg                       <  2:1400755/149‑1 (MQ=255)
                                                                                                                                 tACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCCACAGCTTCACGCGCCAGCTTAGTAATGCGGTCGTAATCGCCCGCTTCCAGCGCATCTGCCGGAACCAGCCAGGAACCACCGATGCACAGCAc                    >  1:128469/1‑149 (MQ=255)
                                                                                                                                 tACAAAAATGACCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCATGCGCCAGCTTAGTAATGCGGTCGTAATCGCCCGCTTCCACCTCATCTGCCGGAACCAGCCAGGAACCACCGATGCACAGCAc                    >  4:1096147/1‑149 (MQ=255)
                                                                                                                                    aaaaaTGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGTAATGCGGTCGTAATCGCCCGCTTCCAGCGCATCTGCCGGAACCAGCCAGGAACCACCGATGcac                         >  3:292248/1‑141 (MQ=255)
                                                                                                                                     aaaaTGCCCGATCGAGGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGTAATGCGGTCGTAATCGCCCGCTTCCAGCGCATCTGCCGGAACCAGCCAGGAACCACCGATGCACAGCACGCtt                >  3:1310348/1‑149 (MQ=255)
                                                                                                                                     aaaaTGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGTAATGCGGTCGTAATCGCCCGCTTCCAGCGCATCTGCCGGAACCAGCCAGGAACCACCGATGCACAGCACGCtt                <  3:384399/149‑1 (MQ=255)
                                                                                                                                        aTGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGTAATGCGGTCGTAATCGCCCGCTTCCAGCGCATCTGCCGGAACCAGCCAGGAACCACCGATGCACAGCACGCTTTTc             >  4:1313619/1‑149 (MQ=255)
                                                                                                                                        aTGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGTAATGCGGTCGTAATCGCCCGCTTCCAGCGCATCTGCCGGAACCAGCCAGGAACCACCGATGCACAGCACGCTTTTc             >  3:506479/1‑149 (MQ=255)
                                                                                                                                         tGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGTAATGCGGTCGTAATCGCCCGCTTCCAGCGCATCTGCCGGAACCAGCCAGGAACCACCGATGCACAGCACGCTTTTCa            <  3:167258/149‑1 (MQ=255)
                                                                                                                                          gCCCGATCGAGGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGTAATGCGGTCGTAATCGCCCGCTTCCAGCGCATCTGCCGGAACCAGCCAGGAACCACCGATGCACAGCACGCTTTTCAg           >  1:175647/1‑149 (MQ=255)
                                                                                                                                          gCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGTAATGCGGTCGTAATCGCCCGCTTCCAGCGCATCTGCCGGAACCAGCCAGGAACCACCGATGCACAGCACGCTTTTCAg           >  1:1019120/1‑149 (MQ=255)
                                                                                                                                          gCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGTAATGCGGTCGTAATCGCCCGCTTCCAGCGCATCTGCCGGAACCAGCCAGGAACCACCGATGCACAGCACGCTTTTCAg           <  4:506479/149‑1 (MQ=255)
                                                                                                                                           cccGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGTAATGCGGTCGTAATCGCCCGCTTCCAGCGCATCTGCCGGAACCAGCCAGGAACCACCGATGCACAGCACGCTTTTCAgc          <  3:1195724/149‑1 (MQ=255)
                                                                                                                                            ccGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGTAATGCgg                                                                                  <  1:1028052/76‑1 (MQ=255)
                                                                                                                                            ccGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGTAATGCgg                                                                                  >  2:1028052/1‑76 (MQ=255)
                                                                                                                                              gATCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGTAATGCGGTCGTAATCGCCCGCTTCCAGCGCATCTGCCGGAACCAGCCAGGAACCACCGATGCACAGCACGCTTTTCAGCGcc       >  3:439007/1‑149 (MQ=255)
                                                                                                                                              gATCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGTAATGCGGTCGTAATCGCCCGCTTCCAGCGCATCTGCCGGAACCAGCCAGGAACCACCGATGCACAGCACGCTTTGCAGCGcc       >  3:1450450/1‑149 (MQ=255)
                                                                                                                                               aTCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGTAATGCGGTCGTAATCGCCCGCTTCCAGCGCATCTGCCGGAACCAGCCAGGAACCACCGATGCACAGCACGCTTTTCAGCGCCa      <  4:1252873/149‑1 (MQ=255)
                                                                                                                                                  cTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGTAATGCGGTCGTAATCGCCCGCTTCCAGCGCATCTGCCGGAACCAGCCAGGAACCACCGATGCACAGCACGCTTTTCAGCGCCAGGt   >  4:1043515/1‑149 (MQ=255)
                                                                                                                                                   tCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGTAATGCGGTCGTAATCGCCCGCTTCCAGCGCATCTGCCGGAACCAGCCAGGAACCACCGATGCACAGCACGCTTTTCAGCGCCAGGTa  <  4:1181004/149‑1 (MQ=255)
                                                                                                                                                   tCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGTAATGCGGTCGTAATCGCCCGCTTCCAGCGCATCTGCCGGAACCAGCCAGGAACCACCGATGCACAGCACGCTTTTCAGCGCCAGGTa  >  1:1701671/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
GGAAATTGATGGCAGCCGTCGTCTGGGGGTGGCACTGGCTACTGCAGAGAGTGTTGTTGACGCCATTGAACGCGCGAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGTAATGCGGTCGTAATCGCCCGCTTCCAGCGCATCTGCCGGAACCAGCCAGGAACCACCGATGCACAGCACGCTTTTCAGCGCCAGGTA  >  NC_000913/1931948‑1932243

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