Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,932,095 T→A 5.3% intergenic (+36/+20) 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,0950TA 5.3% 285.4 / 4.0 108intergenic (+36/+20)purT/edaphosphoribosylglycinamide formyltransferase 2/KHG/KDPG aldolase; 2‑dehydro‑3‑deoxy‑phosphogluconate/4‑hydroxy‑2‑ oxoglutarate aldolase
Reads supporting (aligned to +/- strand):  ref base T (51/48);  new base A (2/4);  total (56/52)
Fisher's exact test for biased strand distribution p-value = 4.37e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.01e-01

Only 100 of 110 total aligned reads displayed.

GAAATTGATGGCAGCCGTCGTCTGGGGGTGGCACTGGCTACTGCAGAGAGTGTTGTTGACGCCATTGAACGCGCGAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGTAATGCGGTCGTAATCGCCCGCTTCCAGCGCATCTGCCGGAACCAGCCAGGAACCACCGATGCACAGCACGCTTTTCAGCGCCAGGTA  >  NC_000913/1931949‑1932243
                                                                                                                                                  |                                                                                                                                                    
gAAATTGATGGCAGCCGTCGTCTGGGGGTGGCACTGGCTACTGCAGAGAGTGTTGTTGACGCCATTGAACGCGCGAAGCACTCCGCCGGACAGGTAAAAGTGCAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTCg                                                                                                                                                    >  1:1031924/1‑149 (MQ=255)
 taaTTGATGGCAGCAGGCGTCTGGGGGTGGCACTTGCTACTGCAGAGAGTTTTGTTGACGCCATTGAACGCGCGAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAATCGCTACAAAAATGCCCGATCCTCGa                                                                                                                                                   <  2:418857/148‑1 (MQ=255)
  aaTTGCAGGCAGCCGTCGTCGTGTGGGGGCACTGGCTACTGCAGAGAGTGTTGTTGACGCCATTGAACGCGCGAAGCTCGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGTTCCTCGAt                                                                                                                                                  <  2:474962/149‑1 (MQ=255)
     tGATGGCAGCCGGCGTCTGGGGGTGGCACTGGCTACTGCAGAGAGTGTTGTTGACGCCATTGAACGCGCGAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTCGATCgg                                                                                                                                               <  1:250422/149‑1 (MQ=255)
       aTGGCAGCCGTCGTCTGGGGGTGGCACTGGCTACTGCAGAGAGTGTTGTTGACGCCATTGAACGCGCGAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCGAGGATCGGGc                                                                                                                                             <  2:586735/149‑1 (MQ=255)
        tGGCAGCGGTCGTCGGGGGGTGGCACTGGCTACTGCAGAGAGTGTTGTTGCCGCCATTGAACGCGCTAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAATAATCGCTACAAAAATGCCCGATCGAGGATCGGGCa                                                                                                                                            <  1:169875/149‑1 (MQ=255)
          gCAGCCGTCGTCTGGGGGTGGCACTGGCTACTGCAGAGAGTGTTGTTGACGCCATTGAACGCGCGAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCAtt                                                                                                                                          >  1:1118289/1‑149 (MQ=255)
          gCAGCCGGCGGCTGGGGGTGGCACTGGCTACTGCAGAGAGTGTTGTTGACGCCATTGAACGCGCGAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCGAGGATCGGGCAtt                                                                                                                                          <  2:1351304/149‑1 (MQ=255)
           cAGCCGTCGTCTGGGGGTGGCACTGGCTCCTGCAGAGAGGGTTGTTGACGCCATTGAACGCGCGAAGCACGCCGCCGGACAGGTAAAAGTACAGGGATAAAACCGGCAAAAAAACGCTACAAAAAAGCCCGATCCTCGATCGGGCAttt                                                                                                                                         <  1:1535534/149‑1 (MQ=255)
           cAGCCGTCGTCTGGGGGTGGCACTGGCTACTGCAGATAGTGTTGTTGACGCCATTGAACGCGCGAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCAttt                                                                                                                                         <  2:882826/149‑1 (MQ=255)
           cAGCCGTCGTCTGGGGGTGGCACTGGCTACTGCAGAGAGTGTTGTTGACGCCATTGAACGCGCGAAGTACGCCGCCGGACAGGAAAAAGTACAGGGATAAACCCGGCAAAAACCCGCTACAAAAATGCCCGCTCCGCGATCGGGcctat                                                                                                                                         >  2:1015772/1‑145 (MQ=255)
           cAGCCGTCGTCTGGGGGGTGCACTGGCTACTGCAGCGAGTGTTGTTGACGCCATTGAACGCGCGAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCAttt                                                                                                                                         <  1:1541641/149‑1 (MQ=255)
            aGCCGTCGTCTGGGGGTGGCACTGGCGACTGCAGAGAGTGTAGTTGACGCCATTGAACGCGCGAAGCACTCCGCCCGACAGGTAAAAGTACAGCGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCAtttt                                                                                                                                        <  1:471144/149‑1 (MQ=255)
              ccGGCGGCTGGGGGGGGCACTGGCTACTGCAGAGAGTGTTGTTGACGCCATTGAACGCGCGAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAAACCCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGa                                                                                                                                      <  1:318141/149‑1 (MQ=255)
                    tcTGGGGGTGGCACTGGCTACTGCAGAGAGTGTTGTTGACGCCATTGAACGCGCGAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTa                                                                                                                                <  1:258586/149‑1 (MQ=255)
                    gcggggggTGGCACTGGCTACTGCAGAGAGTGTTGTTGACGCCATTGAACGCGCGAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCTACAACAATGCCCGATCCTCGATCGGGCATTTTGACTTTTa                                                                                                                                <  1:350016/146‑1 (MQ=255)
                       gggggTGGCACTGGCTACTGCAGAGAGTGTTGTTGACGCCATTGAACGCGCGAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTACAg                                                                                                                             >  1:1204302/1‑149 (MQ=255)
                             ggCACTGGCTACTGCAGAGAGTGTTGTTGACGCCATTGAACGCGCGAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGAGCCTCGATCGGGCATTTTGACTTTTACAGCTTAgc                                                                                                                       <  1:1311434/149‑1 (MQ=255)
                                      tACTGCAGTGAGTGTTGTTGACGCCATTGAACGCGCGAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTAc                                                                                                              <  1:25025/149‑1 (MQ=255)
                                      tACTGCAGAGAGTGTTGTTGACGGAATTGTACGCGCGAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTAc                                                                                                              <  2:280020/149‑1 (MQ=255)
                                      tACTGCAGAGAGTGTTGTTGACGCCATTGACCGCGCGAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGAACCTCGCTCGGGCATTTTGACTTATACAGCTTAGCGCCTTCTAc                                                                                                              >  2:286506/1‑149 (MQ=255)
                                      tACTGCAGAGAGTGTTGTTGACGCCATTGAACGCGCGAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTAc                                                                                                              <  1:1283463/149‑1 (MQ=255)
                                      tACTGCAGAGAGTGTTGTTGACGCCATTGAACGCGCGAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTAc                                                                                                              >  2:626951/1‑149 (MQ=255)
                                      tACTGCAGAGAGTGTTGTTGACGCCATTGAACGCGCGAAGCACGCCGCAGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTAc                                                                                                              >  1:959375/1‑149 (MQ=255)
                                         tGCAGAGTGTGGTGTTGACGCCATTGAACGGGCGAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACGCGGCAAAAAAACGCTACAAAAATGCCCGAACCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGc                                                                                                           <  2:174513/149‑1 (MQ=255)
                                         tGCAGAGAGTGTTGTTGACGCCATTGAACGCGCGAAGCACGCCGCCGGACAGGTAAAAGAACAGGGTTACACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGAATTTTGCCTTTTACAGCATAGCGCCGTCTACAGc                                                                                                           >  2:370555/1‑149 (MQ=255)
                                          gCAGAGAGTGTTGTTGACGCCATTGAACGCGCGAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCt                                                                                                          <  2:245287/149‑1 (MQ=255)
                                           cAGAGAGTGTTGTTGACGCCATTGAACGCGCGAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTCGATCGAGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCtt                                                                                                         <  2:1343797/149‑1 (MQ=255)
                                            agagagTGTTGTTGACGCCATTGAACGCGCGAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTc                                                                                                        <  1:391533/149‑1 (MQ=255)
                                            agagagTGTTGTTGACGCCATTGAACGCGCGAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTc                                                                                                        <  1:1138147/149‑1 (MQ=255)
                                             gagagTGTTGTTGACGCCATTGAACGCGCGAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCa                                                                                                       <  2:591043/149‑1 (MQ=255)
                                                  tgttgtTGACGCCATTGAACGCGCGAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACgcgc                                                                                                  >  2:132023/1‑149 (MQ=255)
                                                  tgttgtTGACGCCATTGAACGCGCGAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACgcgc                                                                                                  >  2:457607/1‑149 (MQ=255)
                                                     tgttgACGCCATTGAACGCGCGAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAg                                                                                               >  1:891301/1‑149 (MQ=255)
                                                         gACGCCATTGAACGCGCGAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTa                                                                                           >  1:301722/1‑149 (MQ=255)
                                                         gACGCCATTGAACGCGCGAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTa                                                                                           >  1:780006/1‑149 (MQ=255)
                                                         gACGCCATTGAACGCGCGAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACCCGCCAGCTTa                                                                                           >  1:950902/1‑149 (MQ=255)
                                                            gCCATTGAACTCGCGAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAACACGCTACAAAAATGCCCGAGCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGACTTCTACATCTTCACGTGCCAGCttatta                                                                                        >  1:1143710/1‑149 (MQ=255)
                                                            gCCATTGAACGCGCGAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTACCGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGTa                                                                                        >  2:320339/1‑149 (MQ=255)
                                                            gCCATTGAACGCGCGAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGTa                                                                                        >  2:401860/1‑149 (MQ=255)
                                                             ccATTGAACGCGCGAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCATCGATCGGGCATTATGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTAAGTaa                                                                                       >  1:213643/1‑149 (MQ=255)
                                                              cATTGAACGCGCGAAGCACGCCGCCGGACAGGTAAAAGTAGAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTCGCTCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTCCAGCTTCACGCGCCAGCTTAGTAAt                                                                                      >  2:1233654/1‑149 (MQ=255)
                                                              cATTGAACGCGCGAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGTAAt                                                                                      >  1:212269/1‑149 (MQ=255)
                                                                    aCGCGCGAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGTAATGCgggc                                                                                <  2:455901/149‑3 (MQ=255)
                                                                      gcgcgAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTGCTACAGCTTCCCGCGCCAGCGTAGTAATGCGGTCGt                                                                              >  2:222515/1‑149 (MQ=255)
                                                                          gAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGTAATGCGGTCGTAATc                                                                          >  1:506852/1‑149 (MQ=255)
                                                                          gAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGTAATGCGGTCGTAATc                                                                          >  1:954869/1‑149 (MQ=255)
                                                                            aGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGTAATGCGGTCGTAATCGc                                                                        >  1:203173/1‑149 (MQ=255)
                                                                            aGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGTAATGCGGTCGTAATCGc                                                                        <  2:891301/149‑1 (MQ=255)
                                                                            aGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCGCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGTAATGCGGTCGTAATCGc                                                                        >  2:376963/1‑149 (MQ=255)
                                                                               aCGCCGCCGGACATGTAAAAGTACATGGTTAAACCCGTCTAAAAAACGCTACAAAAATGCCCTATCCTCGATCGGGCATTTTGACTTTTACATCTTAGCGCCTTCTTCAGCTTCACGCGCCAGCTTAGTAATGCGGTCGTAATCGCCCg                                                                     <  2:1239365/149‑1 (MQ=255)
                                                                                 gccgccGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATTCACGATCCGCGATCGGGCAGTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGGGCCAGCTTAGTAGTGCGGTTGTAATCGCCCGCt                                                                   >  1:1183444/1‑149 (MQ=255)
                                                                                    gccgGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGTAATGCGGTCGAAATCGCCCGCTTcc                                                                >  1:882645/1‑149 (MQ=255)
                                                                                          cAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGTAATGCGGTCGTAATCGCCCGCTTCCAGCGCa                                                          <  1:565047/149‑1 (MQ=255)
                                                                                             gTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGTAATGCGGTCGTAATCGCCCGCTTCCAGCGCATCt                                                       >  2:1228364/1‑149 (MQ=255)
                                                                                             gTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGTAATGCGGTCGTAATCGCCCGCTTCCAGCGCATCt                                                       >  1:1081239/1‑149 (MQ=255)
                                                                                               aaaaGTCCAGGTTTAACCCCGGCATAAAAACGCAACAAAAATGCCCGCTCCTCGATCGTGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGTAATGCTGTCGTAATCGCCCGCTTCCAGCGCATCTGc                                                     <  2:1031924/149‑1 (MQ=255)
                                                                                                aaaGTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGTAATGCGGTCGTAATCGCCCGCTTCCAGCGCATCTGcc                                                    >  1:436428/1‑149 (MQ=255)
                                                                                                 aaGTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGTAATGCGGTCGTAATCGCCCGCTTCCAGCGCATCTGCCg                                                   <  1:141277/149‑1 (MQ=255)
                                                                                                      cAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGTAATGCGGTCGTAATCGCCCGCTTCCAGCGCATCTGCCGGAAcc                                              >  1:667582/1‑149 (MQ=255)
                                                                                                      cAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGTAAAGCGGTCGTAATCGCCCGCTTCCAGCGCATCTGCCGGAAcc                                              >  1:780284/1‑149 (MQ=255)
                                                                                                        gggTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGTAATGCGGGCGTAATCGCCCGCTTCCAGCGCATCTGCCGGAAccag                                            <  1:728285/149‑1 (MQ=255)
                                                                                                         ggTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGTAATGCGGTCGTAATCGCCCGCTTCCAGCGCATCTGCCGGAAccagc                                           <  2:212269/149‑1 (MQ=255)
                                                                                                           ttAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGTAATGCGGTCGTAATCGCCCGCTTCCAGCGCATCTGCCGGAAccagcca                                         <  2:1118289/149‑1 (MQ=255)
                                                                                                               aCCCGGAAAAAAAACGCTACAAAAATGCCCGATACTCGATCGGGCATTTTGACTTTTACAGCTTAGCTCCTTCTACCGCTTCACGCGCCAGCTTAGTAATGAGGTCGTAATCGCCCGCTACCAGCCCATCTGACGGAACCAGCCAggca                                     >  1:479179/1‑147 (MQ=255)
                                                                                                                cccGGCAAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGTAATGCGGGCGTAATCGCCCGCTTCCAGCGCATCTGCCGGAACCAGCCAGGAac                                    <  2:854566/149‑1 (MQ=255)
                                                                                                                 ccGGCAAAAAACCGCTACAAAAATGCCCGCTCCTCGATCGGGCATTTTGCCTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGTAATTCGGTCGTAATCGCCCGCTTCCCGCGCATCTGCCGGAACCAGCCCGGAacc                                   >  1:619863/1‑149 (MQ=255)
                                                                                                                 ccGGCAAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGTAATGCGGTCGTAATCGCCCGCTTCCAGCGCATCTGCCGGAACCAGCCAGCAacc                                   >  2:665189/1‑149 (MQ=255)
                                                                                                                    gCAAAAAAACGCTACAAAAATGCCCGATCCTCTATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGTAATGCGGTCGTAATCGCCCGCTTCCAGCGCATCTGCCGGAACCAGCCAGGAaccacc                                <  1:132023/149‑1 (MQ=255)
                                                                                                                    gCAAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGTAATGCGGTCGTAATCGCCCGCTTCCAGCGCATCTGCCGGAACCAGCCAGGAaccacc                                <  2:438949/149‑1 (MQ=255)
                                                                                                                      aaaaaaaCGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCCCGCGCCAGCTTAGTAATGCGGTCGTAATCGAACGCTTCCAGCGCATCTTCCGGAACAAGCCAGGAACCACCGa                              >  1:245212/1‑149 (MQ=255)
                                                                                                                      aaaaaaaCGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGTAATGCGGTCGTAATCGCCCGCTTCCAGCGCATCTGCCGGAACCAGCCATGAACCACCGa                              <  2:1490617/149‑1 (MQ=255)
                                                                                                                      aaaaaaaCGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGTAATGCGGTCGTAATCGCCCGCTTCCAGCGCATCTGCCGGAACCAGCCAGGAACCACCGa                              >  2:728630/1‑149 (MQ=255)
                                                                                                                      aaaaaaaCGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGTAATGCGGTCGTAATCGCCCGCTTCCAGCGCATCTGCCGGAACCAGCCAGGAACCACCGa                              >  1:1108819/1‑149 (MQ=255)
                                                                                                                      aaaaaaaCGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGTAATGCGGTCGTAATCGCCCGCTTCCAGCGCATCTGCCGGAACCAGCCAGGAACCACCGa                              <  1:457607/149‑1 (MQ=255)
                                                                                                                          aaccGCTACAATAATGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGTAATGCGGTCGTAATCGCCCGCTTCCAGCGCATCTGCCGGAACCAGCCAGGAACCACCGATGca                          <  1:1461831/146‑1 (MQ=255)
                                                                                                                           aaCGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGTAATGCGGTCGTAATCGCCCGCTTCCAGCGCATCTGCCGGAACCAGCCAGGAACCACCGCTGcac                         >  2:510364/1‑149 (MQ=255)
                                                                                                                           aaCGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGTAATGCGGTCGTAATCGCCCGCTTCCAGCGCATCTGCCGGAACCAGCCAGGAACCACCGATGcac                         <  2:301722/149‑1 (MQ=255)
                                                                                                                           aaCGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGTAATGCGGTCGTAATCGCCCGCTTCCAGCGCATCTGCCGGAACCAGCCAGGAACCACCGATGcac                         <  1:159949/149‑1 (MQ=255)
                                                                                                                                tACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGTAATGCGGTCGTAATCGCCCGCTTCCAGCGCATCTGCCGGAACCAGCCAGGAACCACCGATGCACAGCAc                    >  1:1148935/1‑149 (MQ=255)
                                                                                                                                tACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGTAATGCGGTCGTAATCGCCCGCTTCCAGCGCATCTGCCGGAACCAGCCAGGAACCACCGATGCACAGCAc                    >  2:1334387/1‑149 (MQ=255)
                                                                                                                                tACAAAAAAGCCCGATCCTCGCGCGGGCATTTTGACATTTACAGCTTAGCGCCTTCTACAGCATCAAGCGCGAGCTTTGTAATGCGGTCGTAATCGACCGCTTCCAGCGCATCTGCCGGAACCAGACAGGAACCACCGATGCACAGCAc                    >  2:1349063/1‑149 (MQ=255)
                                                                                                                                  cAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGTAATGCGGTCGTAATCGCCCGCTTCCAGCGCATCTGCCGGACCCAGCCAGGAACCACCGATGCACAGCACGc                  >  1:1515910/1‑149 (MQ=255)
                                                                                                                                  cAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGTAATGCGGTCGTAATCGCCCGCTTCCAGCGCATCTGCCGGAACCAGCCAGGAACCACCGATGCACAGCACGc                  >  2:661292/1‑149 (MQ=255)
                                                                                                                                  cAAAAATGCCCGATCCTCGATCGGGCAATTTGACTTTTACAGCTTCGCGCCTTCTACAGCTTCACGCGCCAGCTTACTAATGCGGTCTTAATCGCCCGCTTCCAGCGCAACTGCCGGAACCAGCCAGGAACCACCGAAGCCCAGCAccc                  >  1:676435/1‑147 (MQ=255)
                                                                                                                                   aaaaaTGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGTAATGCGGTCGTAATCGCCCGCTTCCAGCGCATCTGCCGGACCCAGCCAGGAACCACCGATGCACAGCACGCt                 >  2:266951/1‑149 (MQ=255)
                                                                                                                                   aaaaaTGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGTAATGCGGTCGTAATCGCCCGCTTCCAGCGCATCTGCCGGAACCAGCCAGGAACCACCGATGCACAGCACGCt                 >  1:833236/1‑149 (MQ=255)
                                                                                                                                   aaaaaTGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTACTAATGCGGTCGTAATCGCCCGCTTCCAGCGCATCTGCCGGAACCAGCCAGGAACCACCGATGCACAGCACGCt                 >  1:590155/1‑149 (MQ=255)
                                                                                                                                    aaaaTGCCCGATCGAGGATCGGGCATTTTGACTATTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGAAATTCGGTCGTAATCGCCCGCTTCCAGCGCATCGGCCGGAACCAGCCAGTTACCACCGCTGCACAGCACGCtt                >  2:399777/1‑149 (MQ=255)
                                                                                                                                    aaaaTGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGTAATGCGGTCGTAATCGCCCGCTTCCAGCGCATCTGCCGGAACCAGCCAGGAACCACCGATGCACAGCACGCtt                >  1:1325501/1‑149 (MQ=255)
                                                                                                                                     aaaTGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGTAATGCGGTCGTAATCGCCCGCTTCCAGCGCATCTGCCGGAACCAGCCAGGAACCACCGATGCACAGCACGCttt               <  1:1356193/149‑1 (MQ=255)
                                                                                                                                      aaTGCGCGATCCTAGCTCGGGCATTTTGACTTATACAGCTTACCGCATTCTACAGCTTCACGCGCCAGCTTAGTAATGCGGTCGTAATCGCCCGCTTCCAGCGCATCTGCCGGAACCAGCCAGGAACCACCGATGCACAGCACGCtttt              <  2:846439/149‑1 (MQ=255)
                                                                                                                                      aaTGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGTAATGCGGTCGTAATCGCCCGCTTCCAGCGCATCTGCCGGAACCAGCCAGGAACCACCGATGCACAGCACGCtttt              <  2:164906/149‑1 (MQ=255)
                                                                                                                                        tGCCCGATCCTAGATCGTGCATTTTGACTTTTAGATCTTAGCGCCTTATTCAGCTTCACGTGGCAGCTTATTAATGAGGTAGTAATCGCCCGCTTCCAGCGCATCTGCCGGAACCAGCCAGGAAACACCGATGCACCGCACGCTTTTCa            <  2:1143710/149‑1 (MQ=255)
                                                                                                                                         gCCCGATCGAGGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGTAATGCGGTCGTAATCGCCCGCTTCCAGCGCATCTGCCGGAACCAGCCAGGAACCACCGATGCACAGCACGCTTTTCAg           >  1:591538/1‑149 (MQ=255)
                                                                                                                                          cccGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGTAATGCGGTCGTAATCGCCCGCTTCCAGCGCATCTGCCGGAAccagcca                                         >  1:1337539/1‑118 (MQ=255)
                                                                                                                                          cccGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGTAATGCGGTCGTAATCGCCCGCTTCCAGCGCATCTGCCGGAAccagcca                                         <  2:1337539/118‑1 (MQ=255)
                                                                                                                                                ccTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGTAATGCGGTCGTAATCGCCCGCTTCCAGCGCATCTGCCGGAACCAGCCAGGAACCACCGATGCACAGCACGCTTTTCAGCGCCAgg    <  1:401860/149‑1 (MQ=255)
                                                                                                                                                  tCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGTAATGCGGTCGTAATCGCCCGCTTCCAGCGCATCTGCCGGAACCCGCCAGGAACCACCGATGCACAGCACGCTTTTGAGCGCCAGGtt  >  1:1457692/1‑148 (MQ=255)
                                                                                                                                                  tCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGTAATGCGGTCGTAATCGCCCGCTTCCAGCGCATCTGCCGGAACCAGCCAGGAACCACCGATGCACAGCACGCTTTTCAGCGCCAGGTa  >  1:1303373/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                                    
GAAATTGATGGCAGCCGTCGTCTGGGGGTGGCACTGGCTACTGCAGAGAGTGTTGTTGACGCCATTGAACGCGCGAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGTAATGCGGTCGTAATCGCCCGCTTCCAGCGCATCTGCCGGAACCAGCCAGGAACCACCGATGCACAGCACGCTTTTCAGCGCCAGGTA  >  NC_000913/1931949‑1932243

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