Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,147,986 G→C 7.8% L122L (CTG→CTC)  plsX → putative phosphate acyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,147,9860GC 7.8% 755.7 / 44.2 217L122L (CTG→CTC) plsXputative phosphate acyltransferase
Reads supporting (aligned to +/- strand):  ref base G (88/112);  new base C (7/10);  total (95/122)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.34e-01

Only 100 of 220 total aligned reads displayed.

CCGGCCTTCGCAAGCTATCCGCGCCAGTCGTGGGAGTTCAATGCGCGTGGCCCTGGAGCTGGTGAAAGAAGGTCGAGCGCAAGCCTGTGTCAGTGCCGGTAATACCGGGGCGCTGATGGGGCTGGCAAAATTATTACTCAAGCCCCTGGAGGGGATTGAGCGTCCGGCGCTGGTGAC‑GGTATTACCACATCAGCAAAAGGGCAAAACGGTGGTCCTTGACTTAGGGGCCAACGTCGATTGTGACAGCACAATGCTGGTGCAATTTGCCATTATGGGCTCAGTTCTGGCTGAAGA  >  NC_000913/1147839‑1148132
                                                                                                                                                   |                                                                                                                                                   
ccGGCCTTCGCAAGCTATCCGCGCCAGTCGTGGGAGTTCAATGCGCGTGGCCCTGGAGCTGGTGAAAGAAGGTCGAGCGCAAGCCTGTGTCAGTGCCGGTAATACCGGGGCGCTGATGGGGCTGGCAAAATTATTACTCAAGCCCCTgg                                                                                                                                                    <  1:287425/149‑1 (MQ=255)
   gCCTTCGCAAGCTATCCGCGCCAGTCGTGGGAGTTCAATGCGCGTGGCCCTGGAGCTGGTGAAAGAAGGTCGAGCGCAAGCCTGTGTCAGTGCCGGTAATACCGGGGCGCTGATGGGGCTGGCAAAATTATTACTCAAGCCCCTGGAgg                                                                                                                                                 <  2:142316/149‑1 (MQ=255)
    ccTTCGCAAGCTATCCGCGCCAGTCGTGGGAGTTCAATGCGCGTGGCCCTGGAGCTGGTGAAAGAAGGTCGAGCGCAAGCCTGTGTCAGTGCCGGTAATACCGGGGCGCTGATGGGGCTGGCAAAATTATTACTCAAGCCCCTGGAggg                                                                                                                                                <  1:844695/149‑1 (MQ=255)
        cGCAAGCTATCCGCGCCAGTCGTGGGAGTTCAATGCGCGTGGCCCTGGAGCTGGTGAAAGAAGGTCGAGCGCAAGCCTGTGTCAGTGCCGGTAATACCGGGGCGCTGATGGGGCTGGCAAAATTATTACTCAATCCCCTCCAGGGGCtt                                                                                                                                            <  2:833369/149‑1 (MQ=255)
         gCAAGCTATCCGCGCCAGTCGTGGGAGTTCAATGCGCGTGGCCCTGGAGCTGGTGAAAGAAGGTCGAGCGCAAGCCTGTGTCAGTGCCGGTAATACCGGGGCGCTGATGGGGCTGGCAAAATTATTACTCAAGCCCCTGGAgg                                                                                                                                                 <  1:982692/143‑1 (MQ=255)
         gCAAGCTATCCGCGCCAGTCGTGGGAGTTCAATGCGCGTGGCCCTGGAGCTGGTGAAAGAAGGTCGAGCGCAAGCCTGTGTCAGTGCCGGTAATACCGGGGCGCTGATGGGGCTGGCAAAATTATTACTCAAGCCCCTGGAgg                                                                                                                                                 >  2:982692/1‑143 (MQ=255)
          cAAGCTATCCGCGCCAGTCGTGGGAGTTCAATGCGCGTGGCCCTGGAGCTGGTGAAAGAAGGTCGAGCGCAAGCCTGTGTCAGTGCCGGTAATACCGGGGCGCTGATGGGGCTGGCAAAATTATTACTCAATCCCCTCCAGGGGCTTGa                                                                                                                                          <  1:342328/149‑1 (MQ=255)
          cAAGCTATCCGCGCCAGTCGTGGGAGTTCAATGCGCGTGGCCCTGGAGCTGGTGAAAGAAGGTCGAGCGCAAGCCTGTGTCAGTGCCGGTAATACCGGGGCGCTGATGGGGCTGGCAAAATTATTACTCAATCCCCTCCAGGGGCTTGa                                                                                                                                          <  2:706594/149‑1 (MQ=255)
          cAAGCTATCCGCGCCAGTCGTGGGAGTTCAATGCGCGTGGCCCTGGAGCTGGTGAAAGAAGGTCGAGCGCAAGCCTGTGTCAGTGCCGGTAATACCGGGGCGCTGATGGGGCTGGCAAAATTATTACTCAATCCCCTCCAGGGGCTTGa                                                                                                                                          <  3:794979/149‑1 (MQ=255)
          cAAGCTATCCGCGCCAGTCGTGGGAGTTCAATGCGCGTGGCCCTGGAGCTGGTGAAAGAAGGTCGAGCGCAAGCCTGTGTCAGTGCCGGTAATACCGGGGCGCTGATGGGGCTGGCAAAATTATTACTCAAGCCCCTGGAGGGGATTGa                                                                                                                                          <  2:1383765/149‑1 (MQ=255)
          cAAGCTATCCGCGCCAGTCGTGGGAGTTCAATGCGCGTGGCCCTGGAGCTGGTGAAAGAAGGTCGAGCGCAAGCCTGTGTCAGTGCCGGTAATACCGGGGCGCTGATGGGGCTGGCAAAATTATTACTCAAGCCCCTGGAGGGGATTGa                                                                                                                                          <  4:1270671/149‑1 (MQ=255)
          cAAGCTATCCGCGCCAGTCGTGGGAGTTCAATGCGCGTGGCCCTGGAGCTGGTGAAAGAAGGTCGAGCGCAAGCCTGTGTCAGTGCCGGTAATACCGGGGCGCTGATGGGGCTGGCAAAATTATTACTCAAGCCCCTGGAGGGGATTGa                                                                                                                                          <  4:1315056/149‑1 (MQ=255)
            aGCTATCCGCGCCAGTCGTGGGAGTTCAATGCGCGTGGCCCTGGAGCTGGTGAAAGAAGGTCGAGCGCAAGCCTGTGTCAGTGCCGGTAATACCGGGGCGCTGATGGGGCTGGCAAAATTATTACTCAAGCCCCTGGAGGGGATTGAGc                                                                                                                                        <  1:15713/149‑1 (MQ=255)
              cTATCCGCGCCAGTCGTGGGAGTTCAATGCGCGTGGCCCTGGAGCTGGTGAAAGAAGGTCGAGCGCAAGCCTGTGTCAGTGCCGGTAATACCGGGGCGCTGATGGGGCTGGCAAAATTATTACTCAAGCCCCTGGAg                                                                                                                                                  <  2:252798/137‑1 (MQ=255)
               tATCCGCGCCAGTCGTGGGAGTTCAATGCGCGTGGCCCTGGAGCTGGTGAAAGAAGGTCGAGCGCAAGCCTGTGTCAGTGCCGGTAATACCGGGGCGCTGATGGGGCTGGCAAAATTATTACTCAAGCCCCTGGAGGGGATTGAGCGTc                                                                                                                                     <  3:1360917/149‑1 (MQ=255)
                 tCCGCGCCAGTCGTGGGAGTTCAATGCGCGTGGCCCTGGAGCTGGTGAAAGAAGGTCGAGCGCAAGCCTGTGTCAGTGCCGGTAATACCGGGGCGCTGATGGGGCTGGCAAAATTATTACTCAAGCCCCTGGAGGGGATTGAGCGTCCg                                                                                                                                   >  3:502251/1‑149 (MQ=255)
                    gcgcCAGTCGTGGGAGTTCAATGCGCGTGGCCCTGGAGCTGGTGAAAGAAGGTCGAGCGCAAGCCTGTGTCAGTGCCGGTAATACCGGGGCGCTGATGGGGCTGGCAAAATTATTACTCAAGCCCCTGGAGGGGATTGAGCGTCCGgcg                                                                                                                                >  1:813573/1‑149 (MQ=255)
                           tCGTGGGAGTTCAATGCGCGTGGCCCTGGAGCTGGTGAAAGAAGGTCGAGCGCAAGCCTGTGTCAGTGCCGGTAATACCGGGGCGCTGATGGGGCTGGCAAAATTATTACTCAAGCCCCTGGAGGGGATTGAGCGTCCGGCGCTGGTGa                                                                                                                         <  1:911254/149‑1 (MQ=255)
                           tCGTGGGAGTTCAATGCGCGTGGCCCTGGAGCTGGTGAAAGAAGGTCGAGCGCAAGCCTGTGTCAGTGCCGGTAATACCGGGGCGCTGATGGGGCTGGCAAAATTATTACTCAAGCCCCTGGAGGGGATTGAGCGTCCGGCGCTGGTGa                                                                                                                         <  3:824340/149‑1 (MQ=255)
                             gTGGGAGTTCAATGCGCGTGGCCCTGGAGCTGGTGAAAGAAGGTCGAGCGCAAGCCTGTGTCAGTGCCGGTAATACCGGGGCGCTGATGGGGCTGGCAAAATTATTACTCAAGCCCCTGGAGGGGATTGAGCGTCCGGCGCTGGTGAC‑g                                                                                                                      <  1:1441661/149‑1 (MQ=255)
                               gggAGTTCAATGCGCGTGGCCCTGGAGCTGGTGAAAGAAGGTCGAGCGCAAGCCTGTGTCAGTGCCGGTAATACCGGGGCGCTGATGGGGCTGGCAAAATTATTACTCAAGCCCCTGGAGGGGATTGAGCGTCCGGCGCTGGTGAC‑GGt                                                                                                                    >  2:944076/1‑149 (MQ=255)
                                    ttCAATGCGCGTGGCCCTGGAGCTGGTGAAAGAAGGTCGAGCGCAAGCCTGTGTCAGTGCCGGTAATACCGGGGCGCTGATGGGGCTGGCAAAATTATTACTCAAGCCCCTGGAGGGGATTGAGCGTCCGGCGCTGGTGAC‑GGTa                                                                                                                   <  1:1007271/145‑1 (MQ=255)
                                    ttCAATGCGCGTGGCCCTGGAGCTGGTGAAAGAAGGTCGAGCGCAAGCCTGTGTCAGTGCCGGTAATACCGGGGCGCTGATGGGGCTGGCAAAATTATTACTCAAGCCCCTGGAGGGGATTGAGCGTCCGGCGCTGGTGAC‑GGTATTAc                                                                                                               >  3:580176/1‑149 (MQ=255)
                                       aaTGCGCGTGGCCCTGGAGCTGGTGAAAGAAGGTCGAGCGCAAGCCTGTGTCAGTGCCGGTAATACCGGGGCGCTGATGGGGCTGGCAAAATTATTACTCAAGCCCCTGGAGGGGATTGa                                                                                                                                          <  4:757195/120‑1 (MQ=255)
                                       aaTGCGCGTGGCCCTGGAGCTGGTGAAAGAAGGTCGAGCGCAAGCCTGTGTCAGTGCCGGTAATACCGGGGCGCTGATGGGGCTGGCAAAATTATTACTCAAGCCCCTGGAGGGGATTGa                                                                                                                                          >  3:757195/1‑120 (MQ=255)
                                       aaTGCGCGTGGCCCTGGAGCTGGTGAAAGAAGGTCGAGCGCAAGCCTGTGTCAGTGCCGGTAATACCGGGGCGCTGATGGGGCTGGCAAAATTATTACTCAAGCCCCTGGAGGGGATTGAGCGTCCGGCGCTGGTGAC‑GGTATTACcac                                                                                                            >  1:140143/1‑149 (MQ=255)
                                       aaTGCGCGTGGCCCTGGAGCTGGTGAAAGAAGGTCGAGCGCAAGCCTGTGTCAGTGCCGGTAATACCGGGGCGCTGATGGGGCTGGCAAAATTATTACTCAAGCCCCTGGAGGGGATTGAGCGTCCGGCGCTGGTGAC‑GGTATTACcac                                                                                                            >  3:213481/1‑149 (MQ=255)
                                        aTGCGCGTGGCCCTGGAGCTGGTGAAAGAAGGTCGAGCGCAAGCCTGTGTCAGTGCCGGTAATACCGGGGCGCTGATGGGGCTGGCAAAATTATTACTCAAGCCCCTGGAGGGGATTGAGCGTCCGGCGCTGGTGAC‑GGTATTACcaca                                                                                                           >  4:1230690/1‑149 (MQ=255)
                                                 gCCCTGGGGCTGGTGAAAGAAGGTCGAGCGCAAGCCTGTGTCAGTGCCTGTAATACCGGGGCGCTGATGGGGCTGGCAAAATTATTACTCAAGCCCCTGGAGGGGATTGAGCGTCCGGCGCTGGTGAC‑GGTATTACCACATCAGCaaaa                                                                                                  <  1:1322552/149‑1 (MQ=255)
                                                  cccTGGAGCTGGTGAAAGAAGGTCGAGCGCAAGCCTGTGTCAGTGCCGGTAATACCGGGGCGCTGATGGGGCTGGCAAAATTATTACTCAAGCCCCTGGAGGGGATTGAGCGTCCGGCGCTGGTGAC‑GGTATTACCACATCAGCAAAAg                                                                                                 >  1:1087935/1‑149 (MQ=255)
                                                   ccTGGAGCTGGTGAAAGAAGGTCGAGCGCAAGCCTGTGTCAGTGCCGGTAATACCGGGGCGCTGATGGGGCTGGCAAAATTATTACTCAAGCCCCTGGAGGGGATTGAGCGTCCGGCGCTGGTGAC‑GGTATTACCACATCAGCAAAAgg                                                                                                <  3:613507/149‑1 (MQ=255)
                                                   ccTGGAGCTGGTGAAAGAAGGTCGAGCGCAAGCCTGTGTCAGTGCCGGTAATACCGGGGCGCTGATGGGGCTGGCAAAATTATTACTCAAGCCCCTGGAGGGGATTGAGCGTCCGGCGCTGGTGAC‑GGTATTACCACATCAGCAAAAgg                                                                                                <  1:1507447/149‑1 (MQ=255)
                                                   ccTGGAGCTGGTGAAAGAAGGTCGAGCGCAAGCCTGTGTCAGTGCCGGTAATACCGGGGCGCTGATGGGGCTGGCAAAATTATTACTCAAGCCCCTGGAGGGGATTGAGCGTCCGGCGCTGGTGAC‑GGTATTACCACATCAGCAAAAgg                                                                                                <  4:1384936/149‑1 (MQ=255)
                                                    cTGGAGCTGGTGAAAGAAGGTCGAGCGCAAGCCTGTGTCAGTGCCGGTAATACCGGGGCGGTGATGGGGCTGGCAAAATTATTACTCAAGCCCATGGAGGGGATTGAGCGTCCGGCGCTGGTGAC‑GGTATTACCACATCAGCAAAAggg                                                                                               >  3:1158445/1‑149 (MQ=255)
                                                    cTGGAGCTGGTGAAAGAAGGTCGAGCGCAAGCCTGTGTCAGTGCCGGTAATACCGGGGCGCTGATGGGGCTGGCAAAATTATTACTCAAGCCCCTGGAGGGGATTGAGCGTCCGGCGCTGGTGAC‑GGTATTACCACATCAGCAAAAggg                                                                                               >  3:303617/1‑149 (MQ=255)
                                                         gCTGGTGAAAGAAGGTCGAGCGCAAGCCTGTGTCAGTGCCGGTAATACCGGGGCGCTGATGGGGCTGGCAAAATTATTACTCAAGCCCCTGGAGGGGATTGAGCGTCCGGCGCTGGTGAC‑GGTATTACCACATCAGCAAAAGGGCaaaa                                                                                          <  4:35356/149‑1 (MQ=255)
                                                            ggTGAAAGAAGGTCGAGCGCAAGCCTGTGTCAGTGCCGGTAATACCGGGGCGCTGATGGGGCTGGCAAAATTATTACTCAAGCCCCTGGAGGGGATTGAGCGTCCGGCGCTGGTGAC‑GGTATTACCACATCAGCAAAAGGGCAAAACgg                                                                                       <  2:67118/149‑1 (MQ=255)
                                                             gTGAAAGAAGGTCGAGCGCAAGCCTGTGTCAGTGCCGGTAATACCGGGGCGCTGATGGGGCTGGCAAAATTATTACTCAAGCCCCTGGAGGGGATTGAGCGTCCGGCGCTGGTGAC‑GGTATTACCACATCAGCAAAAGGGCAAAACggt                                                                                      >  2:1087314/1‑149 (MQ=255)
                                                             gTGAAAGAAGGTCGAGCGCAAGCCTGTGTCAGTGCCGGTAATACCGGGGCGCTGATGGGGCTGGCAAAATTATTACTCAAGCCCCTGGAGGGGATTGAGCGTCCGGCGCTGGTGAC‑GGTATTACCACATCAGCAAAAGGGCAAAACggt                                                                                      <  1:1355062/149‑1 (MQ=255)
                                                               gAAAGAAGGTCGAGCGCAAGCCTGTGTCAGTGCCGGTAATACCGGGGCGCTGATGGGGCTGGCAAAATTATTACTCAAGCCCCTGGAGGGGAtt                                                                                                                                            <  1:1268280/94‑1 (MQ=255)
                                                               gAAAGAAGGTCGAGCGCAAGCCTGTGTCAGTGCCGGTAATACCGGGGCGCTGATGGGGCTGGCAAAATTATTACTCAAGCCCCTGGAGGGGATTGAGCGTCCGGCGCTGGTGAC‑GGTATTACCACATCAGCAAAAGGGCAAAACggtgg                                                                                    <  3:495016/149‑1 (MQ=255)
                                                                aaaGAAGGTCGAGCGCAAGCCTGTGTCAGTGCCGGTAATACCGGGGCGCTGATGGGGCTGGCAAAATTATTACTCAAGCCCCTGGAGGGGATTGAGCGTCCGGCGCTGGTGAC‑GGTATTACCACATCAGCAAAAGGGCAAAACggtggt                                                                                   >  3:858002/1‑149 (MQ=255)
                                                                aaaGAAGGTCGAGCGCAAGCCTGTGTCAGTGCCGGTAATACCGGGGCGCTGATGGGGCTGGCAAAATTATTACTCAAGCCCCTGGAGGGGATTGAGCGTCCGGCGCTGGTGAC‑GGTATTACCACATCAGCAAAAGGGCAAAACggtggt                                                                                   <  2:390745/149‑1 (MQ=255)
                                                                 aagaagGTCGAGCGCAAGCCTGAGTCAGTGCCGGTAATACCGGGGCGCTGAGGGGGCTGGCAAAATTATTACTCAAGCCCCTGGAGGGGATTGAGCGTCCGGCGCTGGTGAC‑GGTATTACCACATCAGCAAAAGGGCAAAACGGTGGTc                                                                                  >  4:1682498/1‑149 (MQ=255)
                                                                  agaagGTCGAGCGCAAGCCTGTGTCAGTGCCGGTAATACCGGGGCGCTGATGGGGCTGGCAAAATTATTACTCAAGCCCCTGGAGGGGATTGAGCGTCCGGCGCTGGTGAC‑GGTATTACCACATCAGCAAAAGGGCAAAACGGTGGTcc                                                                                 >  2:1276416/1‑149 (MQ=255)
                                                                  agaagGTCGAGCGCAAGCCTGTGTCAGTGCCGGTAATACCGGGGCGCTGATGGGGCTGGCAAAATTATTACTCAAGCCCCTGGAGGGGATTGAGCGTCCGGCGCTGGTGAC‑GGTATTACCACATCAGCAAAAGGGCAAAACGGTGGTcc                                                                                 >  4:261972/1‑149 (MQ=255)
                                                                   gaagGTCGAGCGCAAGCCTGTGTCAGTGCCGGTAATACCGGGGCGCTGATGGGGCTGGCAAAATTATTACTCAAGCCCCTGGAGGGGATTGAGCGTCCGGCGCTGGTGAC‑GGTATTACCACATCAGCAAAAGGGCAAAACGGTGGTCCt                                                                                >  3:1548330/1‑149 (MQ=255)
                                                                    aagGTCGAGCGCAAGCCTGTGTCAGTGCCGGTAATACCGGGGCGCTGATGGGGCTGGCAAAATTATTACTCAAGCCCCTGGAGGGGATTGAGCGTCCGGCGCTGGTGAC‑GGTATTACCACATCAGCAAAAGGGCAAAACGGTGGTCCtt                                                                               <  4:1490447/149‑1 (MQ=255)
                                                                       gTCGAGCGCAAGCCTGTGTCAGTGCCGGTAATACCGGGGCGCTGATGGGGCTGGCAAAATTATTACTCAAGCCCCTGGAGGGGATTGAGc                                                                                                                                        >  1:1068359/1‑90 (MQ=255)
                                                                       gTCGAGCGCAAGCCTGTGTCAGTGCCGGT‑ATTCCGGGGCGGCGGAGGGGCTGGCAAAAATTATACTTAAGGCCCCGGGGGGGAATGGGCGGTCGGGGGCGGTGACGGGAATTACAAC‑TCCGCAAAAGGGGCAAACGGGTGTCCtttgct                                                                           >  1:1695433/1‑145 (MQ=255)
                                                                          gAGCGCAAGCCTGTGTCAGTGCCGGTAATACCGGGGCGCTGATGGGGCTGGCAAAATTATTACTCAAGCCCCTGGAGGGGATTGAGCGTCCGGCGCTGGTGAC‑GGTATTACCACATCAGCAAAAGGGCAAAACGGTGGTCCTTGACTTa                                                                         <  1:607253/149‑1 (MQ=255)
                                                                           aGCGCAAGCCTGTGTCAGTGCCGGTAATACCGGGGCGCTGATGGGGCTGGCAAAATTATTACTCAAGCCCCTGGAGGGGATTGAGCGTCCGGCGCTGGTGAC‑GGTATTACCACATCAGCAAAAGGGCAAAACGGTGGTCCTTGACTTAg                                                                        <  2:1658989/149‑1 (MQ=255)
                                                                           aGCGCAAGCCTGTGTCAGTGCCGGTAATACCGGGGCGCTGATGGGGCTGGCAAAATTATTACTCAAGCCCCTGGAGGGGATTGAGCGTCCGGCGCTGGTGAC‑GGTATTACCAAATCAGCAAAAGGGCAAAACGGTGGTCCTTGACTTAg                                                                        <  4:1596378/149‑1 (MQ=255)
                                                                            gcgcAAGCCTGTGTCAGTGCCGGTAATACCGGGGCGCTGATGGGCCTGGCAAAATTATTACTCAAGCCCCTGGAGGGGATTGAGCGTCCGGCGCTGGTGAC‑GGTATTACCACATCAGCAAAAGGGCAAAACGGTGGTCCTTGACTTAgg                                                                       <  3:1230690/149‑1 (MQ=255)
                                                                             cgcAAGCCTGTGTCAGTGCCGGTAATACCGGGGCGCTGATGGGGCTGGCAAAATTATTACTCAAGCCCCTGGAGGGGATTGAGCGGCCGGCGCTGGTGAC‑GGGATTACCACATCAGCAAAAGGGCAAAACGGTGGTCCTTAACTTAggg                                                                      <  2:759686/149‑1 (MQ=255)
                                                                              gcAAGCCTGTGTCAGTGCCGGTAATACCGGGGCGCTGATGGGGCTGGCAAAATTATTACTCAAGCCCCTGGAGGGGATTGAGCGTCCGGCGCTGGTGAC‑GGTATTACCACATCAGCAAAAGGGCAAAACGGTGGTCCTTGACTTAgggg                                                                     >  2:578409/1‑149 (MQ=255)
                                                                                aaGCCTGTGTCAGTGCCGGTAATACCGGGGCGCTGATGGGGCTGGCAAAATTATTACTCAAGCCCCTGGAGGGGATTGAGCGTCCGGCGCTGGTGAC‑GGTATTACCACATCAGCAAAAGGGCAAAACGGTGGTCCTTGACTTAGGGGcc                                                                   >  3:978333/1‑149 (MQ=255)
                                                                                  gCCTGTGTCAGTGCCGGTAATACCGGGGCGCTGATGGGGCTGGCAAAATTATTACTCAATCCCCTCCAGGGGCTTGa                                                                                                                                          <  1:579234/77‑1 (MQ=38)
                                                                                  gCCTGTGTCAGTGCCGGTAATACCGGGGCGCTGATGGGGCTGGCAAAATTATTACTCAAGCCCCTGGa                                                                                                                                                   <  1:198668/68‑1 (MQ=255)
                                                                                      gtgtCAGTGCCGGTAATACCGGGGCGCTGATGGGGCTGGCAAAATTATTACTCAAGCCCCTGGAGGGGATTGAGCGTCCGGCGCTGGTGAC‑GGTATTACCACATCAGCAAAAGGGCAAAACGGTGGTCCTTGACTTAGGGGCCAACGTc                                                             <  3:261972/149‑1 (MQ=255)
                                                                                        gtCAGTGCCGGTAATACCGGGGCGCTGATGGGGCTGGCAAAATTATTACTCAAGCCCCTGGAGGGGATTGAGCGTCCGGCGCTGGTGAC‑GGTATTACCACATCAGCAAAAGGGCAAAACGGTGGTCCTTGACTTAGGGGCCAACGTCGa                                                           <  3:741821/149‑1 (MQ=255)
                                                                                         tCAGTGCCGGTAATACCGGGGCGCTGATGGGGCTGGCAAAATTATTACTCAAGCCCCTGGAGGGGATTAAGCGTCCGGCGCTGGTGAC‑GGTATTACCACATCAGCAAAAGGGCAAAACGGTGGTCCTTGACTTAGGGGCCAACGTCGAt                                                          <  4:252137/149‑1 (MQ=255)
                                                                                         tCAGTGCCGGTAATACCGGGGCGCTGATGGGGCTGGCAAAATTATTACTCAAGACCCTGGAGGGGATTGTGCGCCCGGCGCTGGTGAC‑GGTATTACCCCATCAGCAAAAGGGCAAAACGGTGGTCCTTGAGTTATGGGCCAACGTCGAt                                                          >  4:629057/1‑149 (MQ=255)
                                                                                            gTGCCGGTAATACCGGGGCGCTGATGGGGCTGGCAAAATTATTACTCAAGCCCCTGGAGGGGATTGAGCGTCCGGCGCTGGTGAC‑GGTATTACCACATCAGCAAAAGGGCAAAACGGTGGTCCTTGACTTAGGGGCCAACGTCGATtgt                                                       <  4:303617/149‑1 (MQ=255)
                                                                                            gTGCCGGTAATACCGGGGCGCTGATGGGGCTGGCAAAATTATTACTCAAGCCCCTGGAGGGGATTGAGCGTCCGGCGCTGGTGAC‑GGTATTACCACATCAGCAAAAGGGCAAAACGGTGGTCCTTGACTTAGGGGCCAACGTCGATtgt                                                       <  4:1338164/149‑1 (MQ=255)
                                                                                               ccGGTAATACCGGGGCGCTGATGGGGCTGGCAAAATTATTACTCAAGCCCCTGGAGGGGATTGAGCGTCCGGCGCTGGTGAC‑GGTATTACCACATCAGCAAAAGGGCAAAACGGTGGTCCTTGACTTAGGGGCCAACGTCGATTGTGAc                                                    <  4:580176/149‑1 (MQ=255)
                                                                                                cGGTAATACCGGGGCGCTGATGGGGCTGGCAAAATTATTACTCAAGCCCCTGGAGGGGATTGAGCGTCCGGCGCTGGTGAC‑GGTATTACCACATCAGCAAAAGGGCAAAACGGTGGTCCTTGACTTAGGGGCCAACGTCGATTGTGACa                                                   <  1:1390769/149‑1 (MQ=255)
                                                                                                 ggTAATACCGGGGCGCTGATGGGGCTGGCAAAATTATTACTCAAGCCCCTGGAGGGGATTGAGCGTCCGGCGCTGGTGAC‑GGTATTACCACATCAGCAAAAGGGCAAAACggtggt                                                                                   <  1:864091/116‑1 (MQ=255)
                                                                                                 ggTAATACCGGGGCGCTGATGGGGCTGGCAAAATTATTACTCAAGCCCCTGGAGGGGATTGAGCGTCCGGCGCTGGTGAC‑GGTATTACCACATCAGCAAAAGGGCAAAACggtggt                                                                                   >  2:864091/1‑116 (MQ=255)
                                                                                                  gTAATACCGGGGCGCTGATGGGGCTGGCAAAATTATTACTCAAGCCCCTGGAGGGGATTGAGCGTCCGGCGCTGGTGAC‑GGTATTACCACATCAGCAAAAGGGCAAAACGGTGGTCCTTGACTTAGGGGCCAACGTCGATTGTGACAGc                                                 >  2:1478360/1‑149 (MQ=255)
                                                                                                   tAATACCGGGGCGCTGATGGGGCTGGCAAAATTATTACTCAAGCCCCTGGAGGGGATTGAGCGTCCGGCGCTGGTGAC‑GGTATTACCACATCAGCAAAAGGGCAAAACGGTGGTCCTTGACTTAGGGGCCAACGTCGATTGTGACAGca                                                <  1:747142/149‑1 (MQ=255)
                                                                                                      tACCGGGGCGCTGATGGGGCTGGCAAAATTATTACTCAAGCCCCTGGAGGGGATTGAGCGTCCGGCGCTGGTGAC‑GGTATTACCACATCAGCAAAAGGGCAAAACGGTGGTCCTTGACTTAGGGGCCAACGTCGATTGTGACAGCACaa                                             <  4:502251/149‑1 (MQ=255)
                                                                                                       aCCGGGGCGCTGATGGGGCTGGCAAAATTATTACTCAAGCCCCTGGAGGGGATTGAGCGTCCGGCGCTGGTGAC‑GGTATTACCACATCAGCAAAAGGGCAAAACGGTGGTCCTTGACTTAggg                                                                      >  2:334206/1‑123 (MQ=255)
                                                                                                           gggCGCTGATGGGGCTGGCAAAATTATTACTCAAGCCCCTGGAGGGGATTGAGCGTCCGGCGCTGGTGAC‑GGTATTACCACATCAGCAAAAGGGCAAAACGGTGGTCCTTGACTTAGGGGCCAACGTCGATTGTGACAGCACAATGCTg                                        <  4:217788/149‑1 (MQ=255)
                                                                                                                     ggggCTGGCAAAATTATTACTCAAGCCCCTGGAGGGGATTGAGCGTCCGGCGCTGGTGAC‑GGTATTACCACATCAGCAAAAGGGCAAAACGGTGGTCCTTGACTTAGGGGCCAACGTCGATTGTGACAGCACAATGCTGGTGCAATTTg                              <  2:1614114/149‑1 (MQ=255)
                                                                                                                           ggCAAAATTATTACTCAAGCCCCTGGAGGGGATTGAGCGTCCGGCGCTGGTGAC‑GGTATTACCACATCAGCAAAAGGGCAAAACGGTGGTCCTTGACTTAGGGGCCAACGTCGATTGTGACAGCACAATGCTGGTGCAATTTGCCATTa                        >  3:1381058/1‑149 (MQ=255)
                                                                                                                            gCAAAAGTATAACTCAAGCCCCTGGAGGGGATTGAGCGTGCGGCGCTGGTGAC‑GGTATTACCAAATCAGCAAAAGGGCAAAACGGTGGTCTTCGACTTAGGGGCCAACGTCGATTGTGACAGCACAATCCTGGTGCAATTTGCGATTAt                       <  4:455723/149‑1 (MQ=255)
                                                                                                                             cAAAATTATTACTCAAGCCCCTGGAGGGGATTGAGCGTCCGGCGCTGGTGAC‑GGTATTACCACATCAGCAAAAGGGCAAAACGGTGGTCCTTGACTTAGGGGCCAACGTCGATTGTGACAGCACAATGCTGGTGCAATTTGCCATTATg                      <  4:1121359/149‑1 (MQ=255)
                                                                                                                              aaaaTTATTACTCAAGCCCCTGGAGGGGATTGAGCGTCCGGCGCTGGTGAC‑GGTATTACCACATCAGCAAAAGGGCAAAACGGTGGTCCTTGACTTAGGGGCCAACGTCGATTGTGACAGCACAATGCTGGTGCAATTTGc                             <  3:25099/141‑1 (MQ=255)
                                                                                                                               aaaTTATTACTCAAGCCCCTGGAGGGGATTGAGCGTCCGGCGCTGGTGAC‑GGTATTACCACATCAGCAAAAGGGCAAAACGGTGGTCCTTGACTTAGGGGCCAACGTCGATTGTGACAGCACAATGCTGGTGCAATTTGCCATTATggg                    >  1:1153778/1‑149 (MQ=255)
                                                                                                                                 attattACTCAAGCCCCTGGAGGGGATTGAGCGTCCGGCGCTGGTGAC‑GGTATTACCACATCAGCAAAAGGGCAAAACGGTGGTCCTTGACTTAGGGGCCAACGTCGATTGTGACAGCACAATGCTGGTGCAATTTGCCATTATGGGCt                  >  3:1622344/1‑149 (MQ=255)
                                                                                                                                  ttattaCTCAAGCCCCTGGAGGGGATTGAGCGTCCGGCGCTGGTGAC‑GGTATTACCACATCAGCAAAAGGGCAAAACGGTGGTCCTTGACTTAGGGGCCAACGTCGATTGTGACAGCACAATGCTGGTGCAATTTGCCATTAt                       >  2:588963/1‑143 (MQ=255)
                                                                                                                                  ttattaCTCAAGCCCCTGGAGGGGATTGAGCGTCCGGCGCTGGTGAC‑GGTATTACCACATCAGCAAAAGGGCAAAACGGTGGTCCTTGACTTAGGGGCCAACGTCGATTGTGACAGCACAATGCTGGTGCAATTTGCCATTAt                       <  1:588963/143‑1 (MQ=255)
                                                                                                                                    attaCTCAAGCCCCTGGAGGGGATTGAGCGTCCGGCGCTGGTGAC‑GGTATTACcaca                                                                                                           <  4:1026434/57‑1 (MQ=255)
                                                                                                                                     ttaCTCAAGCCCCTGGAGGGGATTGAGCGTCCGGCGCTGGTGAC‑GGTATTACCACATCAGCAAAAGGGCAAAACGGTGGTCCTTGACTTAGGGGCCAACGTCGATTGTGACAGCACAATGCTGGTGCAATTTGCCATTATGGGCTCAGt              >  1:115113/1‑149 (MQ=255)
                                                                                                                                     ttaCTCAAGCCCCTGGAGGGGATTGAGCGTCCGGCGCTGGTGAC‑GGTATTACCACATCAGCAAAAGGGCAAAACGGTGGTCCTTGACTTAGGGGCCAACGTCGATTGTGACAGCACAATGCTGGTGCAATTTGCCATTATGGGCTCAGt              >  3:482783/1‑149 (MQ=255)
                                                                                                                                       gcTCAATCCCCTCCAGGGGCTTGAGCGTCCGGCGCTGGTGAC‑GGTATTACCACATCAGCAAAAGGGCAAAACGGTGGTCCTTGACTTAGGGGCCAACGTCGATTGTGACAGCACAATGCTGGTGCAATTTGCCATTATGGGCTCAGTTc            >  4:183533/2‑149 (MQ=255)
                                                                                                                                       aCTCAAGCCCCTGGAGGGGATTGAGCGTCCGGCGCTGGTGAC‑GGTATTACCACATCAGCAAAAGGGCAAAACGGTGGTCCTTGACTTAGGGGCCAACGTCGATTGTGACAGCACAATGCTGGTGCAATTTGCCATTATGGGCTCAGTTc            <  4:482783/149‑1 (MQ=255)
                                                                                                                                         tCAAGCCCCTGGAGGGGATTGAGCGTCCGGCGCTGGTGAC‑GGTATTACCACATCAGCAAAAGGGCAAAACGGTGGTCCTTGACTTAGGGGCCAACGTCGATTGTGACAGCACAATGCTGGTGCAATTTGCCATTATGGGCTCAGTTCTg          >  2:1697650/1‑149 (MQ=255)
                                                                                                                                         tCAAGCCCCTGGAGGGGATTGAGCGTCCGGCGCTGGTGAC‑GGTATTACCACATCAGCAAAAGGGCAAAACGGTGGTCCTTGACTTAGGGGCCAACGTCGATTGTGACAGCACAATGCTGGTGCAATTTGCCATTATGGGCTCAGTTCTg          >  3:1647554/1‑149 (MQ=255)
                                                                                                                                         tCAAGCCCCTGGAGGGGATTGAGCGTCCGGCGCTGGTGAC‑GGTATTACCACATCAGCAAAAGGGCAAAACGGTGGTCCTTGACTTAGGGGCCAACGTCGATTGTGACAGCACAATGCTGGTGCAATTTGCCATTATGGGCTCAGTTCTg          >  1:412503/1‑149 (MQ=255)
                                                                                                                                         tCAAGCCCCTGGAGGGGATTGAGCGTCCGGCGCTGGTGAC‑GGTATTACCACATCAGCAAAAGGGCAAAACGGTGGTCCTTGACTTAGGGGCCAACGTCGATTGTGACAGCACAATGCTGGTGCAATTTGCCATTATGGGCTCAGTTCTg          >  1:257594/1‑149 (MQ=255)
                                                                                                                                           aaGCCCCTGGAGGGGATTGAGCGTCCGGCGCTGGTGAC‑GGTATTACCACATCAGCAAAAGGGCAAAACGGTGGTCCTTGACTTAGGGGCCAACGTCGATtgtg                                                      <  2:1153726/103‑1 (MQ=255)
                                                                                                                                           aaGCCACTGGAGGGGATTGAGCGTCCAGCGCTGGTGAC‑GGTATTACCACATCAGCAAAAGGGCAAAACGGTGGTCCTTGACTTAGGGGCCAACGTCGATTGTGACAGCACAATGCTGGTGCAATTTGCCATTATGGGCTCAGTTCTGGc        <  1:231340/149‑1 (MQ=255)
                                                                                                                                               cccTGGAGGGGATTGAGCGTCCGGCGCTGGTGAC‑GGTATTACCACATCAGCAAAAgg                                                                                                >  4:1028775/1‑57 (MQ=255)
                                                                                                                                               cccTGGAGGGGATTGAGCGTCCGGCGCTGGTGAC‑GGTATTACCACATCAGCAAAAgg                                                                                                <  3:1028775/57‑1 (MQ=255)
                                                                                                                                               cccTGGAGGGGATTGAGCGTCCGGCGCTGGTGAC‑GGTATTACCACATCAGCAAAAGGGCAAAACGGTGGTCCTTGACTTAGGGGCCAACGTCGATTGTGACAGCACAATGCTGGTGCAATTTGCCATTATGGGCTCAGTTCTGGCTGaa    >  1:1281514/1‑149 (MQ=255)
                                                                                                                                                 cTGGAGGGGATTGAGCGTCCGGCGCTGGTGAC‑GGTATTACCACATCAGCAAAAGGGCAAAACGGTGGTCCTTGACTTAGGGGCCAACGTCGATTGTGACAGCACAATGCTGGt                                      >  3:727874/1‑113 (MQ=255)
                                                                                                                                                 cTGGAGGGGATTGAGCGTCCGGCGCTGGTGAC‑GGTATTACCACATCAGCAAAAGGGCAAAACGGTGGTCCTTGACTTAGGGGCCAACGTCGATTGTGACAGCACAATGCTGGt                                      <  4:727874/113‑1 (MQ=255)
                                                                                                                                                 cTGGAGGGGATTGAGCGTCCGGCGCTGGTGAC‑GGTATTACCACATCAGCAAAAGGGCAAAACGGTGGTCCTTGACTTAGGGGCCAACGTCGATTGTGACAGCACAATGCTGGTGCAATTTGCCATTATGGGCTCAGTTCTGGCTGAaga  <  2:1062880/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
CCGGCCTTCGCAAGCTATCCGCGCCAGTCGTGGGAGTTCAATGCGCGTGGCCCTGGAGCTGGTGAAAGAAGGTCGAGCGCAAGCCTGTGTCAGTGCCGGTAATACCGGGGCGCTGATGGGGCTGGCAAAATTATTACTCAAGCCCCTGGAGGGGATTGAGCGTCCGGCGCTGGTGAC‑GGTATTACCACATCAGCAAAAGGGCAAAACGGTGGTCCTTGACTTAGGGGCCAACGTCGATTGTGACAGCACAATGCTGGTGCAATTTGCCATTATGGGCTCAGTTCTGGCTGAAGA  >  NC_000913/1147839‑1148132

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: ‑