Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,147,993 A→C 9.3% I125L (ATT→CTT)  plsX → putative phosphate acyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,147,9930AC 9.3% 296.4 / 18.5 95I125L (ATT→CTT) plsXputative phosphate acyltransferase
Reads supporting (aligned to +/- strand):  ref base A (46/40);  new base C (4/5);  total (50/45)
Fisher's exact test for biased strand distribution p-value = 7.31e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.57e-01

TTCGCAAGCTATCCGCGCCAGTCGTGGGAGTTCAATGCGCGTGGCCCTGGAGCTGGTGAAAGAAGGTCGAGCGCAAGCCTGTGTCAGTGCCGGTAATACCGGGGCGCTGATGGGGCTGGCAAAATTATTACTCAAGCCCCTGGAGGGGATTGAGCGTCCGGCGCTGGTGACGGTATTACCACATCAGCAAAAGGGCAAAACGGTGGTCCTTGACTTAGGGGCCAACGTCGATTGTGACAGCACAATGCTGGTGCAATTTGCCATTATGGGCTCAGTTCTGGCTGAAGAGGTGGTGGA  >  NC_000913/1147845‑1148141
                                                                                                                                                    |                                                                                                                                                    
ttCGCAAGCTATCCGCGCCAGTGGTGGGAGTTCAATGCGCGTGGCCCTGGAGCTGGTGAAAGAAGGTCGAGCGCAAGCCTGTGTCAGTGCCGGTAATACCGGGGCGCTGATGGGGCTGGCAAAATTATTACTCAAGCCCCTGGAGGGGa                                                                                                                                                      <  3:435845/149‑1 (MQ=255)
ttCGCAAGCTATCCGCGCCAGTCGTGGGAGTTCAATGCGAGTGGCCCTGGAGCTGGTGAAAGAAGGTCGAGCGCAAGCCTGTGTCAGTGCCGGTAATACCGGGGCGCTGATGGGGCTGGCAAAATTATTACTCAAGCCCCTGGAGGGGa                                                                                                                                                      >  1:604613/1‑149 (MQ=255)
  cGCAAGCTATCCGCGCCAGTCGTGGGAGTTCAATGCGCGTGGCCCTGGAGCTGGTGAAAGAAGGTCGAGCGCAAGCCTGTGTCAGTGCCGGTAATACCGGGGCGCTGATGGGGCTGGCAAAATTATTACTCAAGCCCCTGGAGGGGAtt                                                                                                                                                    <  1:153968/149‑1 (MQ=255)
  cGCAAGCTATCCGCGCCAGTCGTGGGAGTTCAATGCGCGTGGCCCTGGAGCTGGTGAAAGAAGGTCGAGCGCAAGCCTGTGTCAGTGCCGGTAATACCGGGGCGCTGATGGGGCTGGCAAAATTATTACTCAAGCCCCTGGAGGGGAtt                                                                                                                                                    <  2:428280/149‑1 (MQ=255)
  cGCAAGCTATCCGCGCCAGTCGTGGGAGTTCAATGCGCGTGGCCCTGGAGCTGGTGAAAGAAGGTCGAGCGCAAGCCTGTGTCAGTGCCGGTAATACCGGGGCGCTGATGGGGCTGGCAAAATTATTACTCAAGCCCCTGGAGGGGAtt                                                                                                                                                    <  1:640844/149‑1 (MQ=255)
    cAAGCTATCCGCGCCAGTCGTGGGAGTTCAATGCGCGTGGCCCTGGAGCTGGTGAAAGAAGGTCGAGCGCAAGCCTGTGTCAGTGCCGGTAATACCGGGGCGCTGATGGGGCTGGCAAAATTATTACTCAATCCCCTCCAGGGGCTTGa                                                                                                                                                  <  2:582871/149‑1 (MQ=255)
    cAAGCTATCCGCGCCAGTCGTGGGAGTTCAATGCGCGTGGCCCTGGAGCTGGTGAAAGAAGGTCGAGCGCAAGCCTGTGTCAGTGCCGGTAATACCGGGGCGCTGATGGGGCTGGCAAAATTATTACTCAAGCCCCTGGAGGGGATTGa                                                                                                                                                  <  2:758786/149‑1 (MQ=255)
     aaGCTATCCGCGCCAGTCGTGGGAGTTCAATGCGCGTGGCCCTGGAGCTGGTGAAAGAAGGTCGAGCGCAAGCCTGTGTCAGTGCCGGTAATACCGGGGCGCTGATGGGGCTGGCAAAATTATTACTCAATCCCCTCCAGGGGCTTGAg                                                                                                                                                 <  2:414992/149‑1 (MQ=255)
     aaGCTATCCGCGCCAGTCGTGGGAGTTCAATGCGCGTGGCCCTGGAGCTGGTGAAAGAAGGTCGAGCGCAAGCCTGTGTCAGTGCCGGTAATACCGGGGCGCTGATGGGGCTGGCAAAATTATTACTCAAGCCCCTGGAGGGGATTGAg                                                                                                                                                 >  2:497560/1‑149 (MQ=255)
          aTCCGCGCCAGTCGTGGGAGTTCAATGCGCGTGGCCCTGGAGCTGGTGAAAGAAGGTCGAGCGCAAGCCTGTGTCAGTGCCGGTAATACCGGGGCGCTGATGGGGCTGGCAAAATTATTACTCAAGCCCCTGGAGGGGATTGAGCGTcc                                                                                                                                            <  3:120132/149‑1 (MQ=255)
            ccGCGCCAGTCGTGGGAGTTCAATGCGCGTGGCCCTGGAGCTGGTGAAAGAAGGTCGAGCGCAAGCCTGTGTCAGTGCCGGTAATACCGGGGCGCTGATGGGGCTGGCAAAATTATTACTCAAGCCCCTGGAGGGGATTGAGCGTCCgg                                                                                                                                          >  3:30144/1‑149 (MQ=255)
                  cAGTCGTGGGAGTTCAATGCGCGTGGCCCTGGAGCTGGTGAAAGAAGGTCGAGCGCAAGCCTGTGTCAGTGCCGGTAATACCGGGGCGCTGATGGGGCTGGCAAAATTATTACTCAAGCCCCTGGAGGGGATTGAGCGTCCGGCGCTgg                                                                                                                                    >  3:269756/1‑149 (MQ=255)
                   aGTCGTGGGAGTTCAATGCGCGTGGCCCTGGAGCTGGTGAAAGAAGGTCGAGCGCAAGCCTGTGTCAGTGCCGGTAATACCGGGGCGCTGATGGGGCTGGCAAAATTATTACTCAAGCCCCTGGAGGGGATTGAGCGTCCGGCGCTGGt                                                                                                                                   <  2:752595/149‑1 (MQ=255)
                     tCGTGGGAGTTCAATGCGCGTGGCCCTGGAGCTGGTGAAAGAAGGTCGAGCGCAAGCCTGTGTCAGTGCCGGTAATACCGGGGCGCTGATGGGGCTGGCAAAATTATTACTCAAGCCCCTGGAGGGGATTGAGCGTCCGGCGCTGGTGa                                                                                                                                 >  4:596645/1‑149 (MQ=255)
                      cGTGGGAGTTCAATGCGCGTGGCCCTGGAGCTGGTGAAAGAAGGTCGAGCGCAAGCCTGTGTCAGTGCCGGTAATACCGGGGCGCTGATGGGGCTGGCAAAATTATTACTCAAGCCCCTGGAGGGGATTGAGCGACCGGCGCTGGTGAc                                                                                                                                <  3:57231/149‑1 (MQ=255)
                        tGGGAGTTCAATGCGCGTGGCCCTGGAGCTGGTGAAAGAAGGTCGAGCGCAAGCCTGTGTCAGTGCCGGTAATACCGGGGCGCTGATGGGGCTGGCAAAATTATTACTCAAGCCCCTGGAGGGGATTGAGCGTCCGGCGCTGGTGACgg                                                                                                                              >  2:315352/1‑149 (MQ=255)
                          ggtgTTCAATGCGCGTGGCCCTGGAGCTGGTGAAAGAAGGTCGAGCGCAAGCCTGTGTCAGTGCCGGTAATACCGGGGCGCTGATGGGGCTGGCAAAATTATTACTCAAGCCCCTGGAGGGGATTGAGCGTCCGGCGCTGGTGACGGTa                                                                                                                            >  1:446246/4‑149 (MQ=255)
                            aGTTCAATGCGCGTGGCCCTGGAGCTGGTGAAAGAAGGTCGAGCGCAAGCCTGTGTCAGTGCCGGTAATACCGGGGCGCTGATGGGGCTGGCAAAATTATTACTCAAGCCCCTGGAGGGGATTGAGCGTCCGGCGCTGGTGACGGTAtt                                                                                                                          >  4:62342/1‑149 (MQ=255)
                             gTTCAATGCGCGTGGCCCTGGAGCTGGTGAAAGAAGGTCGAGCGCAAGCCTGTGTCAGTGCCGGTAATACCGGGGCGCTGATGGGGCTGGCAAAATTATTACTCAAGCCCCTGGAGGGGATTGAGCGTCCGGCGCTGGTGACGGTATTa                                                                                                                         <  3:172519/149‑1 (MQ=255)
                              ttCAATGCGCGTGGCCCTGGAGCTGGTGAAAGAAGGTCGAGCGCAAGCCTGTGTCAGTGCCGGTAATACCGGGGCGCTGATGGGGCTGGCAAAATTATTACTCAAGCCCCTGGAGGGGATTGAGCGTCCGGCGCTGGTGACGGTATTAc                                                                                                                        >  3:312669/1‑149 (MQ=255)
                              ttCAATGCGCGTGGCCCTGGAGCTGGTGAAAGAAGGTCGAGCGCAAGCCTGTGTCAGTGCCGGTAATACCGGGGCGCTGATGGGGCTGGCAAAATTATTACTCAAGCCCCTGGAGGGGATTGAGCGTCCGGCGCTGGTGACGGTATTAc                                                                                                                        <  4:312669/149‑1 (MQ=255)
                              ttCAATGCGCGTGGCCCTGGAGCTGGTGAAAGAAGGTCGAGCGCAAGCCTGTGTCAGTGCCGGTAATACCGGGGCGCTGATGGGGCTGGCAAAATTATTACTCAAGCCCCTGGAGGGGATTGAGCGTCCGGCGCTGGTGACGGTATTAc                                                                                                                        <  4:424850/149‑1 (MQ=255)
                                cAATGCGCGTGGCCCTGGAGCTGGTGAAAGAAGGTCGAGCGCAAGCCTGTGTCAGTGCCGGTAATACCGGGGCGCTGATGGGGCTGGCAAAATTATTACTCAAGCCCCTGGAGGGGATTGAGCGTCCGGCGCTGGTGACGGTATTACca                                                                                                                      >  3:474696/1‑149 (MQ=255)
                                 aaTGCGCGTGGCCCTGGAGCTGGTGAAAGAAGGTCGAGCGCAAGCCTGTGTCAGTGCCGGTAATACCGGGGCGCTGATGGGGCTGGCAAAATTATTACTCAAGCCCCTGGAGGGGATTGAGCGTCCGGCGCTGGTGACGGTATTACcac                                                                                                                     >  3:672900/1‑149 (MQ=255)
                                 aaTGCGCGTGGCCCCGGAGGGGGTGGAAGAAGGTCGAGCGCAAGCCTGTGTGAGTGCCGGTAATACCGGGGCGCTGATGGGGCTGGCAAAATTATTACTCAAGCCCCTGGCGGGGCTTGATCGTGCGGCGCTGGTGACGGTATTACcac                                                                                                                     >  4:387078/1‑149 (MQ=255)
                                             ccTGGAGCTGGTGAAAGAAGGTCGAGCGCAAGCCTGTGTCAGTGCCGGTAATACCGGGGCGCTGATGGGGCTGGCAAAATTATTACTCAATCCCCTCCAGGGGCTTGa                                                                                                                                                  >  2:215093/1‑108 (MQ=255)
                                             ccTGGAGCTGGTGAAAGAAGGTCGAGCGCAAGCCTGTGTCAGTGCCGGTAATACCGGGGCGCTGATGGGGCTGGCAAAATTATTACTCAATCCCCTCCAGGGGCTTGa                                                                                                                                                  <  1:215093/108‑1 (MQ=255)
                                             ccTGGAGCTGGTGAAAGAAGGTCGAGCGCAAGCCTGTGTCAGTGCCGGTAATACCGGGGCGCTGATGGGGCTGGCAAAATTATTACTCAAGCCCCTGGAGGGGATTGAGCGTCCGGCGCTGGTGACGGTATTACCACCTCAGCAAAAgg                                                                                                         >  1:665335/1‑149 (MQ=255)
                                             ccTGGAGCTGGTGAAAGAAGGTCGAGCGCAAGCCTGTGTCAGTGCCGGTAATACCGGGGCGCTGATGGGGCTGGCAAAATTATTACTCAAGCCCCTGGAGGGGATTGAGCGTCCGGCGCTGGTGACGGTATTACCACATCAGCAAAAgg                                                                                                         >  1:23983/1‑149 (MQ=255)
                                             ccTGGAGCTGGTGAAAGAAGGTCGAGCGCAAGCCTGTGTCAGTGCCGGTAATACCGGGGCGCTGATGGGGCTGGCAAAATTATTACTCAAGCCCCTGGAGGGGATTGAGCGTCCGGCGCTGGTGACGGTATTACCACATCAGCAAAAgg                                                                                                         >  1:378655/1‑149 (MQ=255)
                                              cTGGAGCTGGTGAAAGAAGGTCGAGCGCAAGCCTGTGTCAGTGCCGGTAATACCGGGGCGCTGATGGGGCTGGCAAAATTATTACTCAAGCCCCTGGAGGGGATTGAGCGTCCGGCGCTGGTGACGGTATTACCACATCAGCAAAAggg                                                                                                        >  3:226493/1‑149 (MQ=255)
                                              cTGGAGCTGGTGAAAGAAGGTCGAGCGCAAGCCTGTGTCAGTGCCGGTAATACCGGGGCGCTGATGGGGCTGGCAAAATTATTACTCAAGCCCCTGGAGGGGATTGAGCGTCCGGCGCTGGTGACGGTATTACCACATCAGCAAAAggg                                                                                                        >  3:299955/1‑149 (MQ=255)
                                                   gCTGGTGAAAGAAGGTCGAGCGCAAGCCTGTGTCAGTGCCGGTAATACCGGGGCGCTGATGGGGCTGGCAAAATTATTACTCAAGCCCCTGGAGGGGATTGAGCGTCCGGCGCTGGTGACGGTATTACCACATCAGCAAAAGGGCaaaa                                                                                                   <  4:303547/149‑1 (MQ=255)
                                                   gCTGGTGAAAGAAGGTCGAGCGCAAGCCTGTGTCAGTGCCGGTAATACCGGGGCGCTGATGGGGCTGGCAAAATTATTACTCAAGCCCCTGGAGGGGATTGAGCGTCCGGCGCTGGTGACGGTATTACCACATCAGCAAAAGGGCaaaa                                                                                                   >  4:314925/1‑149 (MQ=255)
                                                       gTGAAAGAAGGTCGAGCGCAAGCCTGTGTCAGTGCCGGTAATACCGGGGCGCTGATGGGGCTGGCAAAATTATTACTCAAGCCCCTGGAGGGGATTGAGCGTCCGGCGCTGGTGACGGTATTACCACATCAGCAAAAGGGCAAAACggt                                                                                               >  4:356662/1‑149 (MQ=255)
                                                        tGAAAGAAGGTGGAGCGCAAGCCTGTGTCAGTGCCGGTAATACCGGGGCGCTGATGGGGCTGGCAAAATTATTACTCAAGCCCCTGGAGGGGATTGAGCGTCCGGCGCTGGTGACGGTATTACCACATCAGCAAAAGGGCAAAACggtg                                                                                              >  1:573825/1‑149 (MQ=255)
                                                        tGAAAGAAGGTCGAGCGCAAGCCTGTGTCAGTGCCGGTAATACCGGGGCGCTGATGGGGCTGGCAAAATTATTACTCAAGCCCCTGGAGGGGATTGAGCGTCCGGCGCTGGTGACGGTATTACCACATCAGCAAAAGGGCAAAACggtg                                                                                              <  2:23983/149‑1 (MQ=255)
                                                        tGAAAGAAGGTCGAGCGCAAGCCTGTGTCAGTGCCGGTAATACCGGGGCGCTGATGGGGCTGGCAAAATTATTACTCAAGCCCCTGGAGGGGATTGAGCGTCCGGCGCTGGTGACGGTATTACCACATCAGCAAAAGGGCAAAACggtg                                                                                              <  4:269756/149‑1 (MQ=255)
                                                        tGAAAGAAGGTCGAGCGCAAGCCTGTGTCAGTGCCGGTAATACCGGGGCGCTGATGGGGCTGGCAAAATTATTACTCAAGCCCCTGGAGGGGATTGAGCGTCCGGCGCTGGTGACGGTATTACCACATCAGCAAAAGGGCAAAACggtg                                                                                              <  3:8923/149‑1 (MQ=255)
                                                        tGAAAGAAGGTCGAGCGCAAGCCTGTGTCAGTGCCGGTAATACCGGGGCGCTGATGGGGCTGGCAAAATTATTACTCAAGCCCCTGGAGGGGATTGAGCGTCCGGCGCTGGTGACGGTATTACCACATCAGCAAAAGGGCAAAACggtg                                                                                              <  3:323483/149‑1 (MQ=255)
                                                         gAAAGAAGGTCGAGCGCAAGCCTGTGTCAGTGCCGGTAATACCGGGGCGCTGATGGGGCTGGCAAAATTATTACTCAAGCCCCTGGAGGGGATTGAGCGTCCGGCGCTGGTGACGGTATTACCACATCAGCAAAAGGGCAAAACggtgg                                                                                             >  3:610073/1‑149 (MQ=255)
                                                         gAAAGAAGGTCGAGCGCAAGCCTGTGTCAGTGCCGGTAATACCGGGGCGCTGATGGGGCTGGCAAAATTATTACTCAAGCCCCTGGAGGGGATTGAGCGTCCGGCGCTGGTGACGGTATTACCACATAAGCAAAAGGGCAAAACggtgg                                                                                             <  1:230264/149‑1 (MQ=255)
                                                            agaagGTCGAGCGCAAGCCTGTGTCAGTGCCGGTAATACCGGGGCGCTGATGGGGCTGGCAAAATTATTACTCAAGCCCCTGGAGGGGATTGAGCGTCCGGCGCTGGTGACGGTATTACCACATCAGCAAAAGGGCAAAACGGTGGTcc                                                                                          >  3:592391/1‑149 (MQ=255)
                                                             gaagGTCGAGCGCAAGCCTGTGTCAGTGCCGGTAATACCGGGGCGCTGATGGGGCTGGCAAAATTATTACTCAAGCCCCTGGAGGGGATTGAGCGTCCGGCGCTGGTGACGGTATTACCACATCAGCAAAAGGGCAAAACGGTGGTCCt                                                                                         >  3:687252/1‑149 (MQ=255)
                                                             gaagGTCGAGCGCAAGCCTGTGTCAGTGCCGGTAATACCGGGGCGCTGATGGGGCTGGCAAAATTATTACTCAAGCCCCTGGAGGGGATTGAGCGTCCGGCGCTGGTGACGGTATTACCACATCAGCAAAAGGGCAAAACGGTGGTCCt                                                                                         >  4:301298/1‑149 (MQ=255)
                                                              aagGTCGAGCGCAAGCCTGTGTCAGTGCCGGTAATACCGGGGCGCTGATGGGGCTGGCAAAATTATTACTCAAGCCCCTGGAGGGGATTGAGCGTCCGGCGCTGGTGACGGTATTACCACATCAGCAAAAGGGCAAAACGGTGGTCCtt                                                                                        <  2:604613/149‑1 (MQ=255)
                                                              aagGTCGAGCGCAAGCCTGTGTCAGTGCCGGTAATACCGGGGCGCTGATGGGGCTGGCAAAATTATTACTCAAGCCCCTGGAGGGGATTGAGCGTCCGGCGCTGGTGACGGTATTACCACATCAGCAAAAGGGCAAAACGGTGGTCCtt                                                                                        >  1:609254/1‑149 (MQ=255)
                                                              aagGTCGAGCGCAAGCCTGTGTCAGTGCCGGTAATACCGGGGCGCTGATGGGGCTGGCAAAATTATTACTCAAGCCCCTGGAGGGGATTGAGCGTCCGGCGCTGGTGACGGTATTACCACATCAGCAAAAGGGCAAAACGGTGGTCCtt                                                                                        >  4:289757/1‑149 (MQ=255)
                                                                ggTCGAGCGCAAGCCTGTGTCAGTGCCGGTAATACCGGGGCGCTGATGGGGCTGGCAAAGTTATTACTCAAGCCCCTGGAGGGGATTGAGCGTCCGCCGCTGGTGACGGTATTACCACATCAGCAAAAGGGCAAAACGGTGGTCCTTGa                                                                                      <  1:89867/149‑1 (MQ=255)
                                                                ggTCGAGCGCAAGCCTGTGTCAGTGCCGGTAATACCGGGGCGCTGATGGGGCTGGCAAAATTATTACTCAAGCCCCTGGAGGGGATTGAGCGTCCGGCGCTGGTGACGGTATTACCACATCAGCAAAAGGGCAAAACGGTGGTCCTTGa                                                                                      <  2:573825/149‑1 (MQ=255)
                                                                 gTCGAGCGCAAGCCTGTGTCAGTGCCGGTAATACCGGGGCGCTGATGGGGCTGGCAAAATTATTACTCAATCCCCTCCAGGGGCTTGa                                                                                                                                                  <  4:339315/88‑1 (MQ=38)
                                                                 gTCGAGCGCAAGCCTGTGTCAGTGCCGGTAATACCGGGGCGCTGATGGGGCTGGCAAAATTATTACTCAATCCCCTCCAGGGGCTTGa                                                                                                                                                  >  3:339315/1‑88 (MQ=38)
                                                                  tCGAGCGCAAGCCTGTGTCAGTGCCGGTAATACCGGGGCGCTGATGGGGCTGGCAAAATTATTACTCAAGCCCCTGGAGGGGATTGAGCGTCCGGCGCTGGTGACGGTATTACCACATCAGCAAAAGGGCAAAACGGTGGTCCTTGACt                                                                                    <  2:406924/149‑1 (MQ=255)
                                                                      gcgcAAGCCTGTGTCAGTGCCGGTAATACCGGGGCGCTGATGGGGCTGGCAAAATTATTACTCAAGCCCCTGGAGGGGATTGAGCGTCCGGCGCTGGTGACGGTATTACCACATCAGCAAAAGGGCAAAACGGTGGTCCTTGACTTAgg                                                                                <  1:558720/149‑1 (MQ=255)
                                                                       cgcAAGCCTGTGTCAGTGCCGGTAATACCGGGGCGCTGATGGGGCTGGCAAAATTATTACTCAAGCCCCTGGAGGGGATTGAGCGTCCGGCGCTGGTGACGGTATTACCACATCAGCAAAAGGGCAAAACGGTGGTCCTTGACTTAggg                                                                               <  1:451163/149‑1 (MQ=255)
                                                                       cgcAAGCCTGTGTCAGTGCCGGTAATACCGGGGCGCTGATGGGGCTGGCAAAATTATTACTCAAGCCCCTGGAGGGGATTGAGCGTCCGGCGCTGGTGACGGTATTACCACATCAGCAAAAGGGCAAAACGGTGGTCCTTGACTTAggg                                                                               <  4:105711/149‑1 (MQ=255)
                                                                              cTGTGTCAGTGCCGGTAATACCGGGGCGCTGATGGGGCTGGCAAAATTATTACTCAAGCCCCTGGAGGGGATTGAGCGTCCGGCGCTGGTGACGGTATTACCACATCAGCAAAAGGGCAAAACGGTGGTCCTTGACTTAggg                                                                               <  4:15324/142‑1 (MQ=255)
                                                                              cTGTGTCAGTGCCGGTAATACCGGGGCGCTGATGGGGCTGGCAAAATTATTACTCAAGCCCCTGGAGGGGATTGAGCGTCCGGCGCTGGTGACGGTATTACCACATCAGCAAAAGGGCAAAACGGTGGTCCTTGACTTAggg                                                                               >  3:15324/1‑142 (MQ=255)
                                                                                   tCAGTGCCGGTAATACCGGGGCGCTGATGGGGCTGGCAAAATTATTACTCAAGCCCCTGGAGGGGATTGAGCGTCCGGCGCTGGTGACGGTATTACCACATCAGCAAAAGGGCAAAACGGTGGTCCTTGACTTAGGGGCCAACGTCGAt                                                                   >  3:704915/1‑149 (MQ=255)
                                                                                      gTGCCGGTAATACCGGGGCGCTGATGGGGCTGGCAAAATTATTACTCAAGCCCCTGGAGGGGATTGAGCGTCCGGCGCTGGTGACGGTATTACCACATCAGCAAAAGGGCAAAACGGTGGTCCTTGACTTAGGGGCCAACGTCGATtgt                                                                >  3:184905/1‑149 (MQ=255)
                                                                                      gTGCCGGTAATACCGGGGCGCTGATGGGGCTGGCAAAATTATTACTCAAGCCCCTGGAGGGGATTGAGCGTCCGGCGCTGGTGACGGTATTACCACATCAGCAAAAGGGCAAAACGGTGGTCCTTGACTTAGGGGCCAACGTCGATtgt                                                                >  1:315096/1‑149 (MQ=255)
                                                                                      gTGCCGGTAATACCGGGGCGCTGATGGGGCTGGCAAAATTATTACTCAAGCCCCTGGAGGGGATTGAGCGTCCGGCGCTGGTGACGGTATTACCACATCAGCAAAAGGGCAAAACGGTGGTCCTTGACTTAGGGGCCAACGTCGATtgt                                                                >  1:284517/1‑149 (MQ=255)
                                                                                      gTGCCGGTAATACCGGGGCGCTGATGGGGCTGGCAAAATTATTACTCAAGCCCCTGGAGGGGATTGAGCGTCCGGCGCTGGTGACGGTATTACCACATCAGCAAAAGGGCAAAACGGTGGTCCTTGACTTAGGGGCCAACGTCGATtgt                                                                <  2:530329/149‑1 (MQ=255)
                                                                                      gTGCCGGTAATACCGGGGCGCTGATGGGGCTGGCAAAATTATTACTCAAGCCCCTGGAGGGGATTGAGCGTCCGGCGCTGGTGACGGTATTACCACATCAGCAAAAGGGCAAAACGGTGGTCCTTGACTTAGGGGCCAACGTCGATtgt                                                                <  1:771244/149‑1 (MQ=255)
                                                                                        gCCGGTAATACCGGGGCGCTGATGGGGCTGGCAAAATTATTACTCAAGCCCCTGGAGGGGATTGAGCGTCCGGCGCTGGTGACGGTATTACCACATCAGCAAAAGGGCAAAACGGTGGTCCTTGACTTAGGGGCCAACGTCGATTGTGa                                                              <  3:356662/149‑1 (MQ=255)
                                                                                        gCCGGTAATACCGGGGCGCTGATGGGGCTGGCAAAATTATTACTCAAGCCCCTGGAGGGGATTGAGCGTCCGGCGCTGGTGACGGTATTACCACATCAGCAAAAGGACAAAACGGTGGTCCTTGACTTAGGGGCCAACGTCGATTGTGa                                                              <  1:600003/149‑1 (MQ=255)
                                                                                            gTAATACCGGGGCGCTGATGGGGCTGGCAAAATTATTACTCAAGCCCCTGGAGGGGATTGAGCGTCCGGCGCTGGTGACGGTATTACCACATCAGCAAAAGGGCAAAACGGTGGTCCTTGACTTAGGGGCCAACGTCGATTGTGACAGc                                                          >  2:541082/1‑149 (MQ=255)
                                                                                            gTAATACCGGGGCGCTGATGGGGCTGGCAAAATTATTACTCAAGCCCCTGGAGGGGATTGAGCGTCCGGCGCTGGTGACGGTATTACCACATCAGCAAAAGGGCAAAACGGTGGTCCTTGACTTAGGGGCCAACGTCGATTGTGACAGc                                                          >  3:58792/1‑149 (MQ=255)
                                                                                              aaTGCCGGGGCGCTGATGGGGCTGGCAAAATTATTACTCAAGCCCCTGGAGGGGATTGAGCATCCGGCGCTGGTGACGGTATTACCACATCAGCAAAAGGGCAAAACGGTGGTCCTTGACTTAGGGGCCAACGTCGATTGTGACAGcac                                                        >  4:481080/1‑149 (MQ=255)
                                                                                               aTACCGGGGCGCTGATGGGGCTGGCAAAATTATTACTCAAGCCCCTGGAGGGGATTGAGCGTCCGGCGCTGGTGACGGTATTACCACATCAGCAAAAGGGCAAAACGGTGGTCCTTGACTTAGGGGCCAACGTCGATTGTGACAGcaca                                                       >  2:118423/1‑149 (MQ=255)
                                                                                                     gggCGCTGATGGGGCTGGCAAAATTATTACTCAAGCCCCTGGAGGGGATTGAGCGTCCGGCGCTGGTGACGGTATTACCACATCAGCAAAAGGGCAAAACGGTGGTCCTTGACTTAGGGGCCAACGTCGATTGTGACAGCACAATGCTg                                                 <  3:631861/149‑1 (MQ=255)
                                                                                                                gggCTGGCAAAATTATTACTCAAGCCCCTGGAGGGGATTGAGCGTCCGGCGCTGGTGACGGTATTACCACATCAGCAAAAGGGCAAAACGGTGGTCCTTGACTTAGGGGCCAACGTCGATTGTGACAGCACAATGCTGGTGCAATTTGc                                      <  3:537997/149‑1 (MQ=255)
                                                                                                                  gCTGGCAAAATTATTACTCAAGCCCCTGGAGGGGATTGAGCGTCCGGCGCTGGTGACGGTATTACCACATCAGCAAAAGGGCAAAACGGTGGTCCTTGACTTAGGGGCCAACGTCGATTGTGACAGCACAATGCTGGTGCAATTTGCCa                                    <  3:559182/149‑1 (MQ=255)
                                                                                                                  gCTGGCAAAATTATTACTCAAGCCCCTGGAGGGGATTGAGCGTCCGGCGCTGGTGACGGTATTACCACATCAGCAAAAGGGCAAAACGGTGGTCCTTGACTTAGGGGCCAACGTCGATTGTGACAGCACAATGCTGGTGCAATTTGCCa                                    >  2:2254/1‑149 (MQ=255)
                                                                                                                     ggCAAAATTATTACTCAAGCCCCTGGAGGGGATTGAGCGTCCGGCGCTGGTGACGGTATTACCACATCAGCAAAAGGGCAAAACGGTGGTCCTTGACTTAGGGGCCAACGTCGATTGTGACAGCACAATGCTGGTGCAATTTGCCATTa                                 >  1:520410/1‑149 (MQ=255)
                                                                                                                       cAAAATTATTACTCAAGCCCCTGGAGGGGATTGAGCGTCCGGCGCTGGTGACGGTATTACCACATCAGCAAAAGGGCAAAACGGTGGTCCTTGACTTAGGGGCCAACGTCGATTGTGACAGCACAATGCTGGTGCAATTTGCCATTATg                               <  1:315352/149‑1 (MQ=255)
                                                                                                                        aaaaTTATTACTCAAGCCCCTGGAGGGGATTGAGCGGCCGGCGCTGGTGACGGTATTACCACATCAGCAAAAGGGCAAAACGGTGGGCCTTGACTTAGGGGCCAACGTCGATTGTGACAGCACAATGCTGGTGCAATTTGCCATTATgg                              <  1:731295/149‑1 (MQ=255)
                                                                                                                         aaaTTATTACTCAAGCCCCTGGAGGGGATTGAGCGTCCGGCGCTGGTGACGGTATTACCACATCAGCAAAAGGGCAAAACGGTGGTCCTTGACTTAGGGGCCAACGTCGATTGTGACAGCACAATGCTGGTGCAATTTGCCATTATggg                             >  2:406849/1‑149 (MQ=255)
                                                                                                                            ttattaCTCAAGCCCCTGGAGGGGATTGAGCGTCCGGCGCTGGTGACGGTATTACCACATCAGCAAAAGGGCAAAACGGTGGTcc                                                                                          <  2:681318/85‑1 (MQ=255)
                                                                                                                            ttattaCTCAAGCCCCTGGAGGGGATTGAGCGTCCGGCGCTGGTGACGGTATTACCACATCAGCAAAAGGGCAAAACGGTGGTcc                                                                                          >  1:681318/1‑85 (MQ=255)
                                                                                                                             tattaCTCAAGCCCCTGGAGGGGATTGAGCGTCCGGCGCTGGTGACGGTATTACCACATCAGCAAAAGGGCAAAACGGTGGTCCTTGACTTAGGGGCCAACGTCGATTGTGACAGCACAATGCTGGTGCAATTTGCCATTATGGGCTCa                         >  2:621353/1‑149 (MQ=255)
                                                                                                                                  cTCAATCCCCTCCAGGGGCTTGAGCGTCCGGCGCTGGTGACGGTATTACCACATCAGCAAAAGGGCAAAACGGTGGTCCTTGACTTAGGGGCCAACGTCGATTGTGACAGCACAATGCTGGTGCAATTTGCCATTATGGGCTCAGTTCt                    >  3:135395/1‑149 (MQ=255)
                                                                                                                                  cTCAAGCCCCTGGAGGGGATTGAGCGTCCGGCGCTGGTGACGGTATTACCACATCAGCAAAAGGGCAAAACGGTGGTCCTTGACTTAGGGGCCAACGTCGATTGTGACAGCACAATGCTGGTGCAATTTGCCATTATGGGCTCAGTTCt                    <  4:687252/149‑1 (MQ=255)
                                                                                                                                  cTCAAGCCCCTGGAGGGGATTGAGCGTCCGGCGCTGGTGACGGTATTACCACATCAGCAAAAGGGCAAAACGGTGGTCCTTGACTTAGGGGCCAACGTCGATTGTGACAGCACAATGCTGGTGCAATTTGCCATTATGGGCTCAGTTCt                    <  1:118423/149‑1 (MQ=255)
                                                                                                                                  cTCAAGCCCCTGGAGGGGATTGAGCGTCCGGCGCTGGTGACGGTATTACCACATCAGCAAAAGGGCAAAACGGTGGTCCTTGACTTAGGGGCCAACGTCGATTGTGACAGCACAATGCTGGTGCAATTTGCCATTATGGGCTCAGTTCt                    >  2:165908/1‑149 (MQ=255)
                                                                                                                                   tCAAGCCCCTGGAGGGGATTGAGCGTCCGGCGCTGGTGACGGTATTACCACATCAGCAAAAGGGCAAAACGGTGGTCCTTGACTTAGGGGCCAACGTCGATTGTGACAGCACAATGCTGGTGCAATTTGCCATTATGGGCTCAGTTCTg                   >  2:154941/1‑149 (MQ=255)
                                                                                                                                   tCAAGCCCCTGGAGGGGATTGAGCGTCCGGCGCTGGTGACGGTATTACCACATCAGCAAAAGGGCAAAACGGTGGTCCTTGACTTAGGGGCCAACGTCGATTGTGACAGCACAATGCTGGTGCAATTTGCCATTATGGGCTCAGTTCTg                   >  3:274090/1‑149 (MQ=255)
                                                                                                                                   tCAAGCCCCTGGAGGGGATTGAGCGTCCGGCGCTGGTGACGGTATTACCACATCAGCAAAAGGGCAAAACGGTGGTCCTTGACTTAGGGGCCAACGTCGATTGTGACAGCACAATGCTGGTGCAATTTGCCATTATGGGCTCAGTTCTg                   >  2:771556/1‑149 (MQ=255)
                                                                                                                                     aaGCCCCTGGAGGGGATTGAGCGTCCGGCGCTGGTGACGGTATTACCACATCAGCAAAAGGGCAAAACGGTGGTCCTTGACTTAGGGGCCAACGTCGATTGTGACAGCACAATGCTGGTGCAATTTGCCATTATGGGCTCAGTTCTGGc                 >  3:19484/1‑149 (MQ=255)
                                                                                                                                     aaGCCCCTGGAGGGGATTGAGCGTCCGGCGCTGGTGACGGTATTACCACATCAGCAAAAGGGCAAAACGGTGGTCCTTGACTTAGGGGCCAACGTCGATTGTGACAGCACAATGCTGGTGCAATTTGCCATTATGGGCTCAGTTCTGGc                 >  4:386262/1‑149 (MQ=255)
                                                                                                                                      aGCCCCTGGAGGGGATTGAGCGTCCGGCGCTGGTGACGGTATTACCACATCAGCAAAAGGGCAAAACGGTGGTCCTTGACTTAGGGGCCAACGTCGATTGTGACAGCACAATg                                                    >  1:370183/1‑113 (MQ=255)
                                                                                                                                      aGCCCCTGGAGGGGATTGAGCGTCCGGCGCTGGTGACGGTATTACCACATCAGCAAAAGGGCAAAACGGTGGTCCTTGACTTAGGGGCCAACGTCGATTGTGACAGCACAATg                                                    <  2:370183/113‑1 (MQ=255)
                                                                                                                                         cccTGGAGGGGATTGAGCGTCCGGCGCTGGTGACGGTATTACCACATCAGCAAAAGGGCAAAACGGTGGTCCTTGACTTAGGGGCCAACGTCGATTGTGACAGCACAATGCTGGTGCAATTTGCCATTATGGGCTCAGTTCTGGCTGaa             >  4:187261/1‑149 (MQ=255)
                                                                                                                                              gAGGGGATTGAGCGTCCGGCGCTGGTGACGGTATTACCACATCAGCAAAAGGGCAAAACGGTGGTCCTTGACTTAGGGGCCAACGTCGATTGTGACAGCACAATGCTGGTGCAATTTGCCATTATGGGCTCAGTTCTGGCTGAAGAggt        <  2:315096/149‑1 (MQ=255)
                                                                                                                                              gAGGGGATTGAGCGTCCGGCGCTGGTGACGGTATTACCACATCAGCAAAAGGGCAAAACGGTGGTCCTTGACTTAGGGGCCAACGTCGATTGTGACAGCACAATGCTGGTGCAATTTGCCATTATGGGCTCAGTTCTGGCTGAAGAggt        <  3:726524/149‑1 (MQ=255)
                                                                                                                                               aGGGGCTTGAGCGTCCGGCGCTGGTGACGGTATTACCACATCAGCAAAAGGGCAAAACGGTGGTCCTTGACTTAGGGGCCAACGTCGATTGTGACAGCACAATGCTGGTGCAATTTGCCATTATGGGCTCAGTTCTGGCTGAAGAggtg       <  4:135395/149‑1 (MQ=255)
                                                                                                                                                 gggATTGAGCGTCCGGCGCTGGTGACGGTATTACCACATCAGCAAAAGGGCAAAACGGTGGTCCTTGACTTAGGGGCCAACGTCGATTGTGACAGCACAATGCTGGTGCAATTTGCCATTATGGGCTCAGTTCTGGCTGAAGAggtggt     <  2:750995/149‑1 (MQ=255)
                                                                                                                                                    aTTGAGCGTCCGGCGCTGGTGACGGTATTACCACATCAGCAAAAGGGCAAAACGGTGGTCCTTGACTTAGGGGCCAACGTCGATTGTGACAGCACAATGCTGGTGCAATTTGCCATTATGGGCTCAGTTCTGGCTGAAGAGGTGGTGGa  >  2:473729/1‑149 (MQ=255)
                                                                                                                                                    aTTGAGCGTCCGGCGCTGGTGACGGTATTACCACATCAGCAAAAGGGCAAAACGGTGGTCCTTGACTTAGGGGCCAACGTCGATTGTGACAGCACAATGCTGGTGCAATTTGCCATTATGGGCTCAGTTCTGGCTGAAGAGGTGGTGGa  <  2:706626/149‑1 (MQ=255)
                                                                                                                                                    aTTGAGCGTCCGGCGCTGGTGACGGTATTACCACATCAGCAAAAGGGCAAAACGGTGGTCCTTGACTTAGGGGCCAACGTCGATTGTGACAGCACAATGCTGGTGCAATTTGCCATTATGGGCTCAGTTCTGGCTGAAGAGGTGGTGGa  >  2:267151/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
TTCGCAAGCTATCCGCGCCAGTCGTGGGAGTTCAATGCGCGTGGCCCTGGAGCTGGTGAAAGAAGGTCGAGCGCAAGCCTGTGTCAGTGCCGGTAATACCGGGGCGCTGATGGGGCTGGCAAAATTATTACTCAAGCCCCTGGAGGGGATTGAGCGTCCGGCGCTGGTGACGGTATTACCACATCAGCAAAAGGGCAAAACGGTGGTCCTTGACTTAGGGGCCAACGTCGATTGTGACAGCACAATGCTGGTGCAATTTGCCATTATGGGCTCAGTTCTGGCTGAAGAGGTGGTGGA  >  NC_000913/1147845‑1148141

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: