Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,147,987 G→C 6.7% E123Q (GAG→CAG)  plsX → putative phosphate acyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,147,9870GC 6.7% 397.9 / 10.7 116E123Q (GAG→CAG) plsXputative phosphate acyltransferase
Reads supporting (aligned to +/- strand):  ref base G (43/65);  new base C (4/4);  total (47/69)
Fisher's exact test for biased strand distribution p-value = 7.13e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 119 total aligned reads displayed.

CCGGCCTTCGCAAGCTATCCGCGCCAGTCGTGGGAGTTCAATGCGCGTGGCCCTGGAGCTGGTGAAAGAAGGTCGAGCGCAAGCCTGTGTCAGTGCCGGTAATACCGGGGCGCTGATGGGGCTGGCAAAATTATTACTCAAGCCCCTGGAGGGGATTGAGCGTCCGGCGCTGGTGACGGTATTACCACATCAGCAAAAGGGCAAAACGGTGGTCCTTGACTTAGGGGCCAACGTCGATTGTGACAGCACAATGCTGGTGCAATTTGCCATTATGGGCTCAGTTCTGGCTGAAGAGG  >  NC_000913/1147839‑1148134
                                                                                                                                                    |                                                                                                                                                   
ccGGCCTTCGCAAGCTATCCGCGCCAGTCGTGGGAGTTCAATGCGCGTGGCCCTGGAGCTGGTGAAAGAAGGTCGAGCGCAAGCCTGTGTCAGTGCCGGTAATACCGGGGCGCTGATGGGGCTGGCAAAATTATTACTCAAGCCCCTgg                                                                                                                                                     >  4:442766/1‑149 (MQ=255)
   gCCTTCGCAAGCTATCCGCGCCAGTCGTGGGAGTTCAATGCGCGTGGCCCTGGAGCTGGTGAAAGAAGGTCGAGCGCAAGCCTGTGTCAGTGCCGGTAATACCGGGGCGCTGATGGGGCTGGCAAAATTATTACTCAAGCCCCTGGAgg                                                                                                                                                  <  4:588910/149‑1 (MQ=255)
   gCCTTCGCAAGCTATCCGCGCCAGTCGTGGGAGTTCAATGCGCGTGGCCCTGGAGCTGGTGAAAGAAGGTCGAGCGCAAGCCTGTGTCAGTGCCGGTAATACCGGGGCGCTGATGGGGCTGGCAAAATTATTACTCAAGCCCCTGGAgg                                                                                                                                                  <  2:426762/149‑1 (MQ=255)
    ccTTCGCAAGCTATCCGCGCCAGTCGTGGGAGTTCAATGCGCGTGGCCCTGGAGCTGGTGAAAGAAGGTCGAGCGCAAGCCTGTGTCAGTGCCGGTAATACCGGGGCGCTGATGGGGCTGGCAAAATTATTACTCAAGCCCCTGGAggg                                                                                                                                                 <  3:40383/149‑1 (MQ=255)
    ccTTCGCAAGCTATCCGCGCCAGTCGTGGGAGTTCAATGCGCGTGGCCCTGGAGCTGGTGAAAGAAGGTCGAGCGCAAGCCTGTGTCAGTGCCGGTAATACCGGGGCGCTGATGGGGCTGGCAAAATTATTACTCAAGCCCCTGGAggg                                                                                                                                                 <  1:792784/149‑1 (MQ=255)
     cTTCGCAAGCTATCCGCGCCAGTCGTGGGAGTTCAATGCGCGTGGCCCTGGAGCTGGTGAAAGAAGGTCGAGCGCAAGCCTGTGTCAGTGCCGGTAATACCGGGGCGCTGATGGGGCTGGCAAAATTATTACTCAAGCCCCTGGAgggg                                                                                                                                                >  2:567065/1‑149 (MQ=255)
     cTTCGCAAGCTATCCGCGCCAGTCGTGGGAGTTCAATGCGCGTGGCCCTGGAGCTGGTGAAAGAAGGTCGAGCGCAAGCCTGTGTCAGTGCCGGTAATACCGGGGCGCTGATGGGGCTGGCAAAATTATTACTCAAGCCCCTGGAgggg                                                                                                                                                <  2:41296/149‑1 (MQ=255)
      ttCGCAAGCTATCCGCGCCAGTCGTGGGAGTTCAATGCGCGTGGCCCTGGAGCTGGTGAAAGAAGGTCGAGCGCAAGCCTGTGTCAGTGCCGGTAATACCGGGGCGCTGATGGGGCTGGCAAAATTATTACTCAAGCCCCTGGAGGGGa                                                                                                                                               <  4:141384/149‑1 (MQ=255)
      ttCGCAAGCTATCCGCGCCAGTCGTGGGAGTTCAATGCGCGTGGCCCTGGAGCTGGTGAAAGAAGGTCGAGCGCAAGCCTGTGTCAGTGCCGGTAATACCGGGGCGCTGATGGGGCTGGCAAAATTATTACTCAAGCCCCTGGAGGGGa                                                                                                                                               >  1:987896/1‑149 (MQ=255)
       tCGCAAGCTATCCGCGCCAGTCGTGGGAGTTCAATGCGCGTGGCCCTGGAGCTGGTGAAAGAAGGTCGAGCGCAAGCCTGTGTCAGTGCCGGTAATACCGGGGCGCTGATGGGGCTGGCAAAATTATTACTCAAGCCCCTGGAGGGGAt                                                                                                                                              <  4:561066/149‑1 (MQ=255)
          cAAGCTATCCGCGCCAGTCGTGGGAGTTCAATGCGCGTGGCCCTGGAGCTGGTGAAAGAAGGTCGAGCGCAAGCCTGTGTCAGTGCCGGTAATACCGGGGCGCTGATGGGGCTGGCAAAATTATTACTCAATCCCCTCCAGGGGCTTGa                                                                                                                                           <  4:133789/149‑1 (MQ=255)
          cAAGCTATCCGCGCCAGTCGTGGGAGTTCAATGCGCGTGGCCCTGGAGCTGGTGAAAGAAGGTCGAGCGCAAGCCTGTGTCAGTGCCGGTAATACCGGGGCGCTGATGGGGCTGGCAAAATTATTACTCAAGCCCCTGGAGGGGATTGa                                                                                                                                           >  4:307374/1‑149 (MQ=255)
          cAAGCTATCCGCGCCAGTCGTGGGAGTTCAATGCGCGTGGCCCTGGAGCTGGTGAAAGAAGGTCGAGCGCAAGCCTGTGTCAGTGCCGGTAATACCGGGGCGCTGATGGGGCTGGCAAAATTATTACTCAAGCCCCTGGAGGGGATTGa                                                                                                                                           <  1:911036/149‑1 (MQ=255)
          cAAGCTATCCGCGCCAGTCGTGGGAGTTCAATGCGCGTGGCCCTGGAGCTGGTGAAAGAAGGTCGAGCGCAAGCCTGTGTCAGTGCCGGTAATACCGGGGCGCTGATGGGGCTGGCAAAATTATTACTCAAGCCCCTGGAGGGGATTGa                                                                                                                                           <  4:25146/149‑1 (MQ=255)
               tATCCGCGCCAGTCGTGGGAGTTCAATGCGCGTGGCCCTGGAGCTGGTGAAAGAAGGTCGAGCGCAAGCCTGTGTCAGTGCCGGTAATACCGGGGCGCTGATGGGGCTGGCAAAATTATTACTCAAGCCCCTGGAGGGGATTGAGCGTc                                                                                                                                      <  3:587626/149‑1 (MQ=255)
                aTCCGCGCCAGTCGTGGGAGTTCAATGCGCGTGGCCCTGGAGCTGGTGAAAGAAGGTCGAGCGCAAGCCTGTGTCAGTGCCGGTAATACCGGGGCGCTGATGGGGCTGGCAAAATTATTACTCAATCCCCTCCAGGGGCTTGAgtaata                                                                                                                                     <  2:121013/149‑6 (MQ=255)
                aTCCGCGCCAGTCGTGGGAGTTCAATGCGCGTGGCCCTGGAGCTGGTGAAAGAAGGTCGAGCGCAAGCCTGTGTCAGTGCCGGTAATACCGGGGCGCTGATGGGGCTGGCAAAATTATTACTCAAGCCCCTGGAGGGGATTGAGCGTcc                                                                                                                                     <  2:938159/149‑1 (MQ=255)
                 tCCGCGCCAGTCGTGGGAGTTCAATGCGCGTGGCCCTGGAGCTGGTGAAAGAAGGTCGAGCGCAAGCCTGTGTCAGTGCCGGTAATACCGGGGCGCTGATGGGGCTGGCAAAATTATTACTCAATCCCCTCCAGGGGCTTGAgtaa                                                                                                                                       >  4:871796/1‑143 (MQ=255)
                 tCCGCGCCAGTCGTGGGAGTTCAATGCGCGTGGCCCTGGAGCTGGTGAAAGAAGGTCGAGCGCAAGCCTGTGTCAGTGCCGGTAATACCGGGGCGCTGATGGGGCTGGCAAAATTATTACTCAAGCCCCTGGAGGGGATTGAGCGTCCg                                                                                                                                    >  4:363538/1‑149 (MQ=255)
                      gcCAGTCGTGGGAGTTCAATGCGTGTGGCCCTGGAGCTGGTGAAAGAAGGTCGAGCGCAAGCCTGTGTCAGTGCCGGTAATACCGGGGCGCTGATGGGGCTGGCAAAATTATTACTCAAGCCCCTGGAGGGGATTGAGCGTCCGGCGCt                                                                                                                               <  3:218107/149‑1 (MQ=255)
                      gcCAGTCGTGGGAGTTCAATGCGCGTGGCCCTGGAGCTGGTGAAAGAAGGTCGAGCGCAAGCCTGTGTCAGTGCCGGTAATACCGGGGCGCTGATGGGGCTGGCAAAATTATTACTCAAGCCCCTGGAGGGGATTGAGCGTCCGGCGCt                                                                                                                               <  2:757508/149‑1 (MQ=255)
                       ccAGTCGTGGGAGTTCAATGCGCGTGGCCCTGGAGCTGGTGAAAGAAGGTCGAGCGCAAGCCTGTGTCAGTGCCGGTAATACCGGGGCGCTGATGGGGCTGGCAAAATTATTACTCAAGCCCCTGGAGGGGATTGAGCGTCCGGCGCTg                                                                                                                              <  2:311974/149‑1 (MQ=255)
                           tCGTGGGAGTTCAATGCGCGTGGCCCTGGAGCTGGTGAAAGAAGGTCGAGCGCAAGCCTGTGTCAGTGCCGGTAATACCGGGGCGCTGATGGGGCTGGCAAAATTATTACTCAAGCCCCTGGAGGGGATTGAGCGTCCGGCGCTGGTGa                                                                                                                          <  4:235823/149‑1 (MQ=255)
                             gTGGGAGTTCAATGCGCGTGGCCCTGGAGCTGGTGAAAGAAGGTCGAGCGCAAGCCTGTGTCAGTGCCGGTAATACCGGGGCGCTGATGGGGCTGGCAAAATTATTACTCAAGCCCCTGGAGGGGATTGAGCGTCCGGCGCTGGTGACg                                                                                                                        >  1:26221/1‑149 (MQ=255)
                              tGGGAGTTCAATGCGCGTGGCCCTGGAGCTGGTGAAAGAAGGTCGAGCGCAAGCCTGTGTCAGTGCCGGTAATACCGGGGCGCTGATGGGGCTGGCAAAATTATTACTCAAGCCCCTGGAGGGGATTGAGCGTCCGGCGCTGGTGACgg                                                                                                                       <  3:176825/149‑1 (MQ=255)
                              tGGGAGTTCAATGCGCGTGGCCCTGGAGCTGGTGAAAGAAGGTCGAGCGCAAGCCTGTGTCAGTGCCGGTAATACCGGGGCGCTGATGGGGCTGGCAAAATTATTACTCAAGCCCCTGGAGGGGATTGAGCGTCCGGCGCTGGTGACgg                                                                                                                       >  1:721021/1‑149 (MQ=255)
                              tGGGAGTTCAATGCGCGTGGCCCTGGAGCTGGTGAAAGAAGGTCGAGCGCAAGCCTGTGTCAGTGCCGGTAATACCGGGGCGCTGATGGGGCTGGCAAAATTATTACTCAAGCCCCTGGAGGGGATTGAGCGTCCGGCGCTGGTGACgg                                                                                                                       <  1:957491/149‑1 (MQ=255)
                              tGGGAGTTCAATGCGCGTGGCCCTGGAGCTGGTGAAAGAAGGTCGAGCGCAAGCCTGTGTCAGTGCCGGTAATACCGGGGCGCTGATGGGGCTGGCAAAATTATTACTCAAGCCCCTGGAGGGGATTGAGCGTCCGGCGCTGGTGACgg                                                                                                                       <  1:668079/149‑1 (MQ=255)
                               gggAGTTCAATGCGCGTGGCCCTGGAGCTGGTGAAAGAAGGTCGAGCGCAAGCCTGTGTCAGTGCCGGTAATACCGGGGCGCTGATGGGGCTGGCAAAATTATTACTCAAGCCCCTGGAGGGGATTGAGCGTCCGGCGCTGGTGACGGt                                                                                                                      >  4:276300/1‑149 (MQ=255)
                                ggAGTCCAATGCGCGTGGCCCTGGAGCTGGTGAAAGAAGGTCGAGCGCAAGCCTGTGTCAGTGCCGGTAATACCGGGGCGCTGATGGGGCTGGCAAAATTATTACTCAAGCCCCTGGAGGGGATTGAGCGTCCGGCGCTGGTGACGGTa                                                                                                                     <  2:897036/149‑1 (MQ=255)
                                 gAGTTCAATGCGCGTGGCCCTGGAGCTGGTGAAAGAAGGTCGAGCGCAAGCCTGTGTCAGTGCCGGTAATACCGGGGCGCTGATGGGGCTGGCAAAATTATTACTCAAGCCCCTGGAGGGGATTGAGCGTCCGGCGCTGGTGACGGTAt                                                                                                                    >  3:204140/1‑149 (MQ=255)
                                   gTTCAATGCGCGTGGCCCTGGAGCTGGTGAAAGAAGGTCGAGCGCAAGCCTGTGTCAGTGCCGGTAATACCGGGGCGCTGATGGGGCTGGCAAAATTATTACTCAAGCCCCTGGAGGGGATTGAGCGTCCGGCGCTGGTGACGGTATTa                                                                                                                  <  1:988770/149‑1 (MQ=255)
                                    ttCAATGCGCGTGGCCCTGGAGCTGGTGAAAGAAGGTCGAGCGCAAGCCTGTGTCAGTGCCGGTAATACCGGGGCGCTGATGGGGCTGGCAAAATTATTACTCAAGCCCCTGGAGGGGATTGAGCGTCCGGCGCTGGTGACGGTATTAc                                                                                                                 <  2:117878/149‑1 (MQ=255)
                                       aaTGCGCGTGGCCCTGGAGCTGGTGAAAGAAGGTCGAGCGCAAGCCTGTGTCAGTGCCGGTAATACCGGGGCGCTGATGGGGCTGGCAAAATTATTACTCAATCCCCTCCAGGGGCTTGa                                                                                                                                           <  3:890535/120‑1 (MQ=255)
                                       aaTGCGCGTGGCCCTGGAGCTGGTGAAAGAAGGTCGAGCGCAAGCCTGTGTCAGTGCCGGTAATACCGGGGCGCTGATGGGGCTGGCAAAATTATTACTCAATCCCCTCCAGGGGCTTGa                                                                                                                                           >  4:890535/1‑120 (MQ=255)
                                       aaTGCGCGTGGCCCTGGAGCTGGTGAAAGAAGGTCGAGCGCAAGCCTGTGTCAGTGCCGGTAATACCGGGGCGCTGATGGGGCTGGCAAAATTATTACTCAAGCCCCTGGAGGGGATTGAGCGTCCGGCGCTGGTGACGGTATTACcac                                                                                                              >  3:818389/1‑149 (MQ=255)
                                           cgcgTGGCCCTGGAGCTGGTGAAAGAAGGTCGAGCGCAAGCCTGTGTCAGTGCCGGTAATACCGGGGCGCTGATGGGGCTGGCAAAATTATTACTCAAGCCCCTGGAGGGGATTGAGCGTCCGGCGCTGGTGACGGTATTACCACATCa                                                                                                          <  1:695509/149‑1 (MQ=255)
                                                   ccTGGAGCTGGTGAAAGAAGGTCGAGCGCAAGCCTGTGTCAGTGCCGGTAATACCGGGGCGCTGATGGGGCTGGCAAAATTATTACTCAAGCCCCTGGAGGGGATTGAGCGTCCGGCGCTGGTGACGGTATTACCACATCAGCAAAAgg                                                                                                  <  3:319183/149‑1 (MQ=255)
                                                    cTGGAGCTGGTGAAAGAAGGTCGAGCGCAAGCCTGTGTCAGTGCCGGTAATACCGGGGCGCTGATGGGGCTGGCAAAATTATTACTCAAGCCCCTGGAGGGGATTGAGCGTCCGGCGCTGGTGACGGTATTAc                                                                                                                 >  4:481419/1‑133 (MQ=255)
                                                    cTGGAGCTGGTGAAAGAAGGTCGAGCGCAAGCCTGTGTCAGTGCCGGTAATACCGGGGCGCTGATGGGGCTGGCAAAATTATTACTCAAGCCCCTGGAGGGGATTGAGCGTCCGGCGCTGGTGACGGTATTAc                                                                                                                 <  3:481419/133‑1 (MQ=255)
                                                    cTGGAGCTGGTGAAAGAAGGTCGAGCGCAAGCCTGTGTCAGTGCCGGTAATACCGGGGCGCTGATGGGGCTGGCAAAATTATTACTCAAGCCCCTGGAGGGGATTGAGCGTCCGGCGCTGGTGACGGTATTACCACATCAGCAAAAggg                                                                                                 >  4:734434/1‑149 (MQ=255)
                                                    cTGGAGCTGGTGAAAGAAGGTCGAGCGCAAGCCTGTGTCAGTGCCGGTAATACCGGGGCGCTGATGGGGCTGGCAAAATTATTACTCAAGCCCCTGGAGGGGATTGAGCGTCCGGCGCTGGTGACGGTATTACCACATCAGCAAAAggg                                                                                                 >  1:458976/1‑149 (MQ=255)
                                                     tGGAGCTGGTGAAAGAAGGTCGAGCGCAAGCCTGTGTCAGTGCCGGTAATACCGGGGCGCTGATGGGGCTGGCAAAATTATTACTCAAGCCCCTGGAGGGGATTGAGCGTCCGGCGCTGGTGACGGTATTAc                                                                                                                 <  3:824028/132‑1 (MQ=255)
                                                     tGGAGCTGGTGAAAGAAGGTCGAGCGCAAGCCTGTGTCAGTGCCGGTAATACCGGGGCGCTGATGGGGCTGGCAAAATTATTACTCAAGCCCCTGGAGGGGATTGAGCGTCCGGCGCTGGTGACGGTATTAc                                                                                                                 >  4:824028/1‑132 (MQ=255)
                                                       gAGCTGGTGAAAGAAGGTCGAGCGCAAGCCTGTGTCAGTGCCGGTAATACCGGGGCGCTGATGGGGCTGGCAAAATTATTACTCAAGCCCCTGGAGGGGATTGAGCGTCCGGCGCTGGTGACGGTATTACCACATCAGCAAAAGGGCaa                                                                                              <  2:142640/149‑1 (MQ=255)
                                                             gTGAAAGAAGGTCGAGCGCAAGCCTGTGTCAGTGCCGGTAATACCGGGGCGCTGATGGGGCTGGCAAAATTATTACTCAAGCCCCTGGAGGGGATTGAGCGTCCGGCGCTGGTGACGGTATTACCACATCAGCAAAAGGTCAAAACggt                                                                                        <  2:887548/149‑1 (MQ=255)
                                                              tGAAAGAAGGTCGAGCGCAAGGCTGTGTCAGTGCCGGGAATGCCGGGGCGCTGATGGGGCTGGCAAAATTATTACTGACGCCCCTGGAGGGGATTGAGCGTCCGGCGCTGGTGACGGTATTACCACCTCAGGAAAAGGGCAAAACggtg                                                                                       <  4:660225/149‑1 (MQ=255)
                                                              tGAAAGAAGGTCGAGCGCAAGCCTGTGTCAGTGCCGGTAATACCGGGGCGCTGATGGGGCTGGCAAAATTATTACTCAAGCCCCTGGAGGGGATTGAGCGTCCGGCGCTGGTGACGGTATTACCACATCAGCAAAAGGGCAAAACggtg                                                                                       <  2:987896/149‑1 (MQ=255)
                                                                 aagaagGTCGAGCGCAAGCCTGTGTCAGTGCCGGTAATACCGGGGCGCTGATGGGGCTGGCAAAATTATTACTCAAGCCCCTGGAGGGGATTGAGCGTCCGGCGCTGGTGACgg                                                                                                                       >  1:122865/1‑114 (MQ=255)
                                                                 aagaagGTCGAGCGCAAGCCTGTGTCAGTGCCGGTAATACCGGGGCGCTGATGGGGCTGGCAAAATTATTACTCAAGCCCCTGGAGGGGATTGAGCGTCCGGCGCTGGTGACgg                                                                                                                       <  2:122865/114‑1 (MQ=255)
                                                                 aagaagGTCGAGCGCAAGCCTGTGTCAGTGCCGGTAATACCGGGGCGCTGATGGGGCTGGCAAAATTATTACTCAAGCCCCTGGAGGGGATTGAGCGTCCGGCGCTGGTGACGGTATTACCACATCAGCAAAAGGGCAAAACGGTGGTc                                                                                    <  2:492086/149‑1 (MQ=255)
                                                                   gaagGTCGAGCGCAAGCCTGTGTCAGTGCCGGTAATACCGGGGCGCTGATGGGGCTGGCAAAATTATTACTCAAGCCCCTGGAGGGGATTGAGCGTCCGGCGCTGGTGACGGTATTACCACATCAGCAAAAGGGCAAAACGGTGGTCCt                                                                                  <  1:776251/149‑1 (MQ=255)
                                                                    aagGTCGAGCGCAAGCCTGTGTCAGTGCCGGTAATACCGGGGCGCTGATGGGGCTGGCAAAATTATTACTCAAGCCCCTGGAGGGGATTGAGCGTCCGGCGCTGGTGACGGTATTACCACATCAGCAAAAGGGCAAAACGGTGGTCCtt                                                                                 <  1:304890/149‑1 (MQ=255)
                                                                      ggTCGAGCGCAAGCCTGTGTCAGTGCCGGTAATACCGGGGCGCTGATGGGGCTGGCAAAATTATTACTCAAGCCCCTGGAGGGGATTGAGCGTCCGGCGCTGGTGACGGTATTACCACATCAGCAAAAGGGCAAAACGGTGGTCCTTGa                                                                               <  2:906719/149‑1 (MQ=255)
                                                                          gAGCGCAAGCCTGTGTCAGTGCCGGTAATACCGGGGCGCTGATGGGGCTGGCAAAATTATTACTCAAGCCCCTGGAGGGGATTGAGCGTCCGGCGCTGGTGACGGTATTACCACATCAGCAAAAGGGCAAAACGGTGGTCCTTGACTTa                                                                           >  3:1058357/1‑149 (MQ=255)
                                                                            gcgcAAGCCTGTGTCAGTGCCGGTAATACCGGGGCGCTGATGGGGCTGGCAAAATTATTACTCAAGCCCCTGGAGGGGATTGAGCGTCCGGCGCTGGTGACGGTATTACCACATCAGCAAAAGGGCAAAACGGTGGTCCTTGACTTAgg                                                                         <  3:773104/149‑1 (MQ=255)
                                                                             cgcAAGCCTGTGTCAGTGCCGGTAATACCGGGGCGCTGATGGGGCTGGCAAAATTATTACTCAAGCCCCTGGAGGGGATTGAGCGTCCGGCGCTGGTGACGGTATTACCACATCAGCAAAAGGGCAAAACGGTGGTCCTTGACTTAggg                                                                        <  3:825002/149‑1 (MQ=255)
                                                                                aaGCCTGTGTCAGTGCCGGTAATACCGGGGCGCTGATGGGGCTGGCAAAATTATTACTCAAGCCCCTGGAGGGGATTGAGCGTCCGGCGCTGGTGACGGTATTACCACATCAGCAAAAGGGCAAAACGGTGGTCCTTGACTTAGGGGcc                                                                     >  4:38360/1‑149 (MQ=255)
                                                                                  gCCTGTGTCAGTGCCGGTAATACCGGGGCGCTGATGGGGCTGGCAAAATTATTACTCAAGCCCCTGGAGGGGATTGAGCGTCCGGCGCTGGTGACGGTATTACCACATCAGCaaaa                                                                                                    <  1:578360/116‑1 (MQ=255)
                                                                                  gCCTGTGTCAGTGCCGGTAATACCGGGGCGCTGATGGGGCTGGCAAAATTATTACTCAAGCCCCTGGAGGGGATTGAGCGTCCGGCGCTGGTGACGGTATTACCACATCAGCaaaa                                                                                                    >  2:578360/1‑116 (MQ=255)
                                                                                     tgtgTCAGTGCCGGTAATACCGGGGCGCTGATGGGGCTGGCAAAATTATTACTCAAGCCCCTGGAGGGGATTGAGCGTCCGGCGCTGGTGACGGTATTACCACATCAGCAAAAGGGCAAAACGGTGGTCCTTGACTTAGGGGCCAACGt                                                                >  2:792974/1‑149 (MQ=255)
                                                                                      gtgtCAGTGCCGGTAATACCGGGGCGCTGATGGGGCTGGCAAAATTATTACTCAAGCCCCTGGAGGGGATTGAGCGTCCGGCGCTGGTGACGGTATTa                                                                                                                  <  2:812134/98‑1 (MQ=255)
                                                                                          cAGTGCCGGTAATACCGGGGCGCTGATGGGGCTGGCAAAATTATTACTCAAGCCCCTGGAGGGGATTGAGCGTCCGGCGCTGGTGACGGTATTACCACATCAGCAAAAGGGCAAAACGGTGGTCCTTGACTTAGGGGCCAACGTCGAtt                                                           >  1:413414/1‑149 (MQ=255)
                                                                                            gTGCCGGTAATACCGGGGCGCTGATGGGGCTGGCAAAATTATTACTCAAGCCCCTGGAGGGGATTGAGCGTCCGGCGCTGGTGACGGTATTACCACATCAGCAAAAGGGCAAAACGGTGGTCCTTGACTTAGGGGCCAACGTCGATtgt                                                         <  2:1084818/149‑1 (MQ=255)
                                                                                            gTGCCGGTAATACCGGGGCGCTGATGGGGCTGGCAAAATTATTACTCAAGCCCCTGGAGGGGATTGAGCGTCCGGCGCTGGTGACGGTATTACCACATCAGCAAAAGGGCAAAACGGTGGTCCTTGACTTAGGGGCCAACGTCGATtgt                                                         <  3:996150/149‑1 (MQ=255)
                                                                                               ccGGTAATACCGGGGCGCTGATGGGGCTGGCAAAATTATTACTCAAGCCCCTGGAGGGGATTGAGCGTCCGGCGCTGGTGACGGTATTACCACATCAGCAAAAGGGCAAAACGGTGGTCCTTGACTTAGGGGCCAACGTCGATTGTGAc                                                      <  3:525980/149‑1 (MQ=255)
                                                                                                    aaTACCGGGGCGCTGATGGGGCTGGCAAAATTATTACTCAAGCCCCTGGAGGGGATTGAGCGTCCGGCGCTGGTGACGGTATTAcc                                                                                                                >  1:870056/1‑86 (MQ=255)
                                                                                                    aaTACCGGGGCGCTGATGGGGCTGGCAAAATTATTACTCAAGCCCCTGGAGGGGATTGAGCGTCCGGCGCTGGTGACGGTATTAcc                                                                                                                <  2:870056/86‑1 (MQ=255)
                                                                                                             gcgcTGATGGGGCTGGCAAAATTATTACTCAAGCCCCTGGAGGGGATTGAGCGTCCGGCGCTGGTGACGGTATTACCACATCAGCAAAAGGGCAAAACGGTGGTCCTTGACTTAGGGGCCAACGTCGATTGTGACAGCACAATGCTGGt                                        <  4:839868/149‑1 (MQ=255)
                                                                                                             gcgcTGATGGGGCTGGCAAAATTATTACTCAAGCCCCTGGAGGGGATTGAGCGTCCGGCGCTGGTGACGGTATTACCACATCAGCAAAAGGGCAAAACGGTGGTCCTTGACTTAGGGGCCAACGTCGATTGTGACAGCACAATGCTGGt                                        <  3:136345/149‑1 (MQ=255)
                                                                                                                 tGATGGGGCTGTCAAAATTATTACTCAAGCCCCTGGAGGGGATTGAGCGTCCGGCGCTGGTGACGGTATTACCACATCAGCAAAAGGGCAAAACGGTGGTCCTTGACTTAGGGGCCAACGTCGATTGTGACAGCACAATGCTGGTGCaa                                    <  3:701455/149‑1 (MQ=255)
                                                                                                                  gATGGGGCTGGCAAAATTATTACTCAAGCCCCTGGAGGGGATTGAGCGTCCGGCGCTGGTGACGGTATTACCACATCAGCAAAAGGGCAAAACGGTGGTCCTTGACTTAGGGGCCAACGTCGATTGTGACAGCACAATGCTGGTGCAAt                                   >  2:844228/1‑149 (MQ=255)
                                                                                                                      gggCTGGCAAAATTATTACTCAAGCCCCTGGAGGGGATTGAGCGTCCGGCGCTGGTGACGGTATTACCACATCAGCAAAAGGGCAAAACGGTGGTCCTTGACTTAGGGGCCAACGTCGATTGTGACAGCACAATGCTGGTGCAATTTGc                               <  3:363538/149‑1 (MQ=255)
                                                                                                                           ggCAAAATTATTACTCAAGCCCCTGGAGGGGATTGAGCGTCCGGCGCTGGTGACGGTATTACCACATCAGCAAAAGGGCAAAACGGTGGTCCTTGACTTAGGGGCCAACGTCGATTGTGACAGCACAATGCTGGTGCAATTTGCCATTa                          >  3:170108/1‑149 (MQ=255)
                                                                                                                              aaaaTTATTACTCAAGCCCCTGGAGGGGATTGAGCGTCCGGCGCTGGTGACGGTATTACCACATCAGCAAAAGGGCAAAACGGTGGTCCTTGACTTAGGGGCCAACGTCGATTGTGACAGCACAATGCTGGTGCAATTTGCCATTATgg                       >  4:802361/1‑149 (MQ=255)
                                                                                                                              aaaaTTATTACTCAAGCCCCTGGAGGGGATTGAGCGTCCGGCGCTGGTGACGGTATTACCACATCAGCAAAAGGGCAAAACGGTGGTCCTTGACTTAGGGGCCAACGTCGATTGTGACAGCACAATGCTGGTGCAATTTGCCATTATgg                       >  3:426005/1‑149 (MQ=255)
                                                                                                                              aaaaTTATTACTCAAGCCCCTGGAGGGGATTGAGCGTCCGGCGCTGGTGACGGTATTACCACATCAGCAAAAGGGCAAAACGGTGGTCCTTGACTTAGGGGCCAACGTCGATTGTGACAGCACAATGCTGGTGCAATTTGCCATTATgg                       <  2:308304/149‑1 (MQ=255)
                                                                                                                              aaaaTTATTACTCAAGCCCCTGGAGGGGATTGAGCGTCCGGCGCTGGTGACGGTATTACCACATCAGCAAAAGGGCAAAACGGTGGTCCTTGACTTAGGGGCCAACGTCGATTGTGACAGCACAATGCTGGTGCAATTTGCCATTATgg                       >  2:698916/1‑149 (MQ=255)
                                                                                                                              aaaaTTATTACTCAAGCCCCTGGAGGGGATTGAGCGTCCGGCGCTGGTGACGGTATTACCACATCAGCAAAAGGGCAAAACGGTGGTCCTTGACTTAGGGGCCAACGTCGATTGTGACAGCACAATGCTGGTGCAATTTGCCATTATgg                       >  3:39166/1‑149 (MQ=255)
                                                                                                                              aaaaTTATTACTCAAGCCCCTGGAGGGGATTGAGCGTCCGGCGCTGGTGACGGTATTACCACATCAGCAAAAGGGCAAAACGGTGGTCCTTGACTTAGGGGCCAACGTCGATTGTGACAGCACAATGCTGGTGCAATTTGCCATTATgg                       >  2:979869/1‑149 (MQ=255)
                                                                                                                                  ttattaCTCAAGCCCCTGGAGGGGATTGAGCGTCCGGCGCTGGTGACGGTATTACCACATCAGCAAAAGGGCAAAACGGTGGTCCTTGACTTAGGGGCCAACGTCGATTGTGACAGCACAATGCTGGTGCAATTTGCCATTa                          <  4:358662/142‑1 (MQ=255)
                                                                                                                                       gcTCAATCCCCTCCAGGGGCTTGAGCGTCCGGCGCTGGTGACGGTACTACCACATCAGCAAAAGGGCAAAACGGTGGTCCTTGACTTAGGGGCCAACGTCGATTGTGACAGCACAATGCTGGTGCAATTTGCCATTATGGGCTCAGTTc              >  2:937105/2‑149 (MQ=255)
                                                                                                                                         tCAAGCCCCTGGAGGGGATTGAGCGTCCGGCGCTGGTGACGGTATTACCACATCAGCAAAAGGGCAAAACGGTGGTCCTTGACTTAGGGGCCAACGTCGATTGTGACAGCACAATGCTGGTGCAATTTGCCATTATGGGCTCAGTTCTg            >  1:1020387/1‑149 (MQ=255)
                                                                                                                                          cAAGCCCCTGGAGGGGATTGAGCGTCCGGCGCTGGTGACGGTATTACCACATCAGCAAAAGGGCAAAACGGTGGTCCTTGACTTAGGGGCCAACGTCGATTGTGACAGCACAATGCTGGTGCAATTTGCCATTATGGGCTCAGTTCTgg           <  1:844123/149‑1 (MQ=255)
                                                                                                                                          cAAGCCCCTGGAGGGGATTGAGCGTCCGGCGCTGGTGACGGTATTACCACATCAGCAAAAGGGCAAAACGGTGGTCCTTGACTTAGGGGCCAACGTCGATTGTGACAGCACAATGCTGGTGCAATTTGCCATTATGGGCTCAGTTCTgg           >  3:874638/1‑149 (MQ=255)
                                                                                                                                          cAAGCCCCTGGAGGGGATTGAGCGTCCGGCGCTGGTGACGGTATTACCACATCAGCAAAAGGGCAAAACGGTGGTCCTTGACTTAGGGGCCAACGTCGATTGTGACAGCACAATGCTGGTGCAATTTGCCATTATGGGCTCAGTTCTgg           <  3:802361/149‑1 (MQ=255)
                                                                                                                                           aaGCCCCTGGAGGGGATTGAGCGTCCGGCGCTGGTGACGGTATTACCACATCAGCAAAAGGGCAAAACGGTGGTCCTTGACTTAGGGGCCAACGTCGATTGTGACAGCACAATGCTGGTGCAATTTGCCATTATGGGCTCAGTTCTGGc          >  3:348790/1‑149 (MQ=255)
                                                                                                                                           aaGCCCCTGGAGGGGATTGAGCGTCCGGCGCTGGTGACGGTATTACCACATCAGCAAAAGGGCAAAACGGTGGTCCTTGACTTAGGGGCCAACGTCGATTGTGACAGCACAATGCTGGTGCAATTTGCCATTATGGGCTCAGTTCTGGc          <  4:204140/149‑1 (MQ=255)
                                                                                                                                             tccccTCCAGGGGCTTGAGCGTCCGGCGCTGGTGACGGTATTACCACATCAGCAAAAGGGCAAAACGGTGGTCCTTGACTTAGGGGCCAACGTCGATTGTGACAGCACAATGCTGGTGCAATTTGCCATTATGGGCTCAGTTCTGGCTg        >  2:811211/2‑149 (MQ=255)
                                                                                                                                               cccTGGAGGGGATTGAGCGTCCGGCGCTGGTGACGGTATTACCACATCAGCAAAAGGGCAAAACGGTGGTCCTTGACTTAGGGGCCAACGTCGATTGTGATAGCACAATGCTGGTGCAATTTGCCATTATGGGCTCAGTTCTGGCTGaa      <  1:567065/149‑1 (MQ=255)
                                                                                                                                                ccTGGAGGGGATTGAGCGTCCGTCGCTGGTGACGGTATTACCACATCAGCAAAAGGGCAAAACGGTGGCCCTTGACTTAGGGGCCAACGTCGATTGTGACAGCACAATGCTGGTGCAATTTGCCATTATGGGCTCAGTTCTGGCTGAag     >  2:472872/1‑149 (MQ=255)
                                                                                                                                                 cTGGAGGGGATTGAGCGTCCGGCGCTGGTGACGGTATTACCACATCAGCAAAAGGGCAAAACGGTGGTCCTTGACTTAGGGGCCAACGTCGATTGTGACAGCACAATGCTGGt                                        <  4:400467/113‑1 (MQ=255)
                                                                                                                                                 cTGGAGGGGATTGAGCGTCCGGCGCTGGTGACGGTATTACCACATCAGCAAAAGGGCAAAACGGTGGTCCTTGACTTAGGGGCCAACGTCGATTGTGACAGCACAATGCTGGt                                        >  3:400467/1‑113 (MQ=255)
                                                                                                                                                 cTGGAGGGGATTGAGCGTCCGGCGCTGGTGACGGTATTACCACATCAGCAAAAGGGCAAAACGGTGGTCCTTGACTTAGGGGCCAACGTCGATTGTGACAGCACAATGCTGGt                                        <  4:315418/113‑1 (MQ=255)
                                                                                                                                                 cTGGAGGGGATTGAGCGTCCGGCGCTGGTGACGGTATTACCACATCAGCAAAAGGGCAAAACGGTGGTCCTTGACTTAGGGGCCAACGTCGATTGTGACAGCACAATGCTGGt                                        >  3:315418/1‑113 (MQ=255)
                                                                                                                                                 cTGGAGGGGATTGAGCGTCCGGCGCTGGTGACGGTATTACCACATCAGCAAAAGGGCAAAACGGTGGTCCTTGACTTAGGGGCCAACGTCGATTGTGACAGCACAATGCTGGTGCAATTTGCCATTATGGGCTCAGTTCTGGCTGAaga    >  1:1003109/1‑149 (MQ=255)
                                                                                                                                                 cTGGAGGGGATTGAGCGTCCGGCGCTGGTGACGGTATTACCACATCAGCAAAAGGGCAAAACGGTGGTCCTTGACTTAGGGGCCAACGTCGATTGTGACAGCACAATGCTGGTGCAATTTGCCATTATGGGCTCAGTTCTGGCTGAaga    >  3:532965/1‑149 (MQ=255)
                                                                                                                                                 cTGGAGGGGATTGAGCGTCCGGCGCTGGTGACGGTATTACCACATCAGCAAAAGGGCAAAACGGTGGTCCTTGACTTAGGGGCCAACGTCGATTGTGACAGCACAATGCTGGTGCAATTTGCCATTATGGGCTCAGTTCTGGCTGAaga    <  1:571025/149‑1 (MQ=255)
                                                                                                                                                  tGGAGGGGATTGAGCGTCCGGCGCTGGTGACGGTATTACCACATCAGCAAAAGGGCAAAACGGTGGTCCTTGACTTAGGGGCCAACGTCGATTGTGACAGCACAATGCTGGTGCAATTTGCCATTATGGGCTCAGTTCTGGCTGAagag   <  3:38360/149‑1 (MQ=255)
                                                                                                                                                   ggAGGGGATTGAGCGTCCGGCGCTGGTGACGGTATTACCACATCAGCAAAAGGGCAAAACGGTGGTCCTTGACTTAGGGGCCAACGTCGATTGTGACAGCACAATGCTGGTGCAATTTGCCATTATGGGCTCAGTTCTGGCTGAAGAgg  >  1:257056/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
CCGGCCTTCGCAAGCTATCCGCGCCAGTCGTGGGAGTTCAATGCGCGTGGCCCTGGAGCTGGTGAAAGAAGGTCGAGCGCAAGCCTGTGTCAGTGCCGGTAATACCGGGGCGCTGATGGGGCTGGCAAAATTATTACTCAAGCCCCTGGAGGGGATTGAGCGTCCGGCGCTGGTGACGGTATTACCACATCAGCAAAAGGGCAAAACGGTGGTCCTTGACTTAGGGGCCAACGTCGATTGTGACAGCACAATGCTGGTGCAATTTGCCATTATGGGCTCAGTTCTGGCTGAAGAGG  >  NC_000913/1147839‑1148134

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