Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,147,993 A→C 8.9% 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 8.9% 248.9 / 17.7 79I125L (ATT→CTT) plsXputative phosphate acyltransferase
Reads supporting (aligned to +/- strand):  ref base A (38/34);  new base C (3/4);  total (41/38)
Fisher's exact test for biased strand distribution p-value = 7.05e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.82e-01

TTCGCAAGCTATCCGCGCCAGTCGTGGGAGTTCAATGCGCGTGGCCCTGGAGCTGGTGAAAGAAGGTCGAGCGCAAGCCTGTGTCAGTGCCGGTAATACCGGGGCGCTGATGGGGCTGGCAAAATTATTACTCAAGCCCCTGGAGGGGATTGAGCGTCCGGCGCTGGTGACGGTATTACCACATCAGCAAAAGGGCAAAACGGTGGTCCTTGACTTAGGGGCCAACGTCGATTGTGACAGCACAATGCTGGTGCAATTTGCCATTATGGGCTCAGTTCTGGCTGAAGAGGTGGTGG  >  NC_000913/1147845‑1148140
                                                                                                                                                    |                                                                                                                                                   
ttCGCAAGCTATCCGCGCCAGTCGTGGGAGTTCAATGCGCGTGGCCCTGGAGCTGGTGAAAGAAGGTCGAGCGCAAGCCTGTGTCAGTGCCGGTAATACCGGGGCGCTGATGGGGCTGGCAAAATTATTACTCAAGCCCCTGGAGGGGa                                                                                                                                                     >  2:544870/1‑149 (MQ=255)
ttCGCAAGCTATCCGCGCCAGTCGTGGGAGTTCAATGCGCGTGGCCCTGGAGCTGGTGAAAGAAGGTCGAGCGCAAGCCTGTGTCAGTGCCGGTAATACCGGGGCGCTGATGGGGCTGGCAAAATTATTACTCAAGCCCCTGGAGGGGa                                                                                                                                                     >  2:23896/1‑149 (MQ=255)
  cGCAAGCTATCCGCGCCAGTCGTGGGAGTTCAATGCGCGTGGCCCTGGAGCTGGTGAAAGAAGGTCGAGCGCAAGCCTGTGTCAGTGCCGGTAATACCGGGGCGCTGATGGGGCTGGCAAAATTATTACTCAATCCCCTCCAGGGGCtt                                                                                                                                                   <  3:392184/149‑1 (MQ=255)
    cAAGCTATCCGCGCCAGTCGTGGGAGTTCAATGCGCGTGGCCCTGGAGCTGGTGAAAGAAGGTCGAGCGCAAGCCTGTGTCAGTGCCGGTAATACCGGGGCGCTGATGGGGCTGGCAAAATTATTACTCAAGCCCCTGGAGGGGATTGa                                                                                                                                                 <  4:787383/149‑1 (MQ=255)
     aaGCTATCCGCGCCAGTCGTGGGAGTTCAATGCGCGTGGCCCTGGAGCTGGTGAAAGAAGGTCGAGCGCAAGCCTGTGTCAGTGCCGGTAATACCGGGGCGCTGATGGGGCTGGCAAAATTATTACTCAAGCCCCTGGAGGGGATTGAg                                                                                                                                                <  1:544870/149‑1 (MQ=255)
        cTATCCGCGCCAGTCGTGGGAGTTCAATGCGCGTGGCCCTGGAGCTGGTGAAAGAAGGTCGAGCGCAAGCCTGTGTCAGTGCCGGTAATACCGGGGCGCTGATGGGGCTGGCAAAATTATTACTCAATCCCCTCCAGGGGCTTGAgtaa                                                                                                                                             <  3:635241/149‑4 (MQ=255)
         tATCCGCGCCAGTCGTGGGAGTTCAATGCGCGTGGCCCTGGAGCTGGTGAAAGAAGGTCGAGCGCAAGCCTGTGTCAGTGCCGGTAATACCGGGGCGCTGATGGGGCTGGCAAAATTATTACTCAAGCCCCTGGAGGGGATTGAGCGTc                                                                                                                                            >  1:336244/1‑149 (MQ=255)
           tCCGCGCCAGTCGTGGGAGTTCAATGCGCGTGGCCCTGGAGCTGGTGAAAGAAGGTCGAGCGCAAGCCTGTGTCAGTGCCGGTAATACCGGGGCGCTGATGGGGCTGGCAAAATTATTACTCAAGCCCCTGGAGGGGATTGAGCGTCCg                                                                                                                                          >  2:260952/1‑149 (MQ=255)
                         gggAGTTCAATGCGCGTGGCCCTGGAGCTGGTGAAAGAAGGTCGAGCGCAAGCCTGTGTCAGTGCCGGTAATACCGGGGCGCTGATGGGGCTGGCAAAATTATTACTCAAGCCCCTGGAGGGGATTGAGCGTCCGGCGCTGGTGACGGt                                                                                                                            >  1:215842/1‑149 (MQ=255)
                          ggAGTTCAATGCGCGTGGCCCTGGAGCTGGTGAAAGAAGGTCGAGCGCAAGCCTGTGTCAGTGCCGGTAATACCGGGGCGCTGATGGGGCTGGCAAAATTATTTCTCAATCCCCTCCAGGGGCTTGAgtaataa                                                                                                                                          <  1:285232/134‑7 (MQ=255)
                          ggAGTTCAATGCGCGTGGCCCTGGAGCTGGTGAAAGAAGGTCGAGCGCAAGCCTGTGTCAGTGCCGGTAATACCGGGGCGCTGATGGGGCTGGCAAAATTATTACTCAAGCCCCTGGAGGGGATTGAGCGTCCGGCGCTGGTGACGGTa                                                                                                                           <  1:616573/149‑1 (MQ=255)
                          ggAGCTCAATGCGCGTCGCCCTGGAGCTGGTGAAAGAAGGTCGAGCGCAAGCCTGTGTCAGTGCCGGTAATACCGGGGCGCTGATGGGGCTGGCAAAATTATTACTCAATCCCCTCCAGGGGCTTGAgtaataa                                                                                                                                          >  2:285232/1‑128 (MQ=255)
                              ttCAATGCGCGTGGCCCTGGAGCTGGTGAAAGAAGGTCGAGCGCAAGCCTGTGTCAGTGCCGGTAATACCGGGGCGCTGATGGGGCTGGCAAAATTATTACTCAAGCCCCTGGAGGGGATTGAGCGTCCGGCGCTGGTGACGGTATTAc                                                                                                                       <  2:588010/149‑1 (MQ=255)
                                cAATGCGCGTGGCCCTGGAGCTGGTGAAAGAAGGTCGAGCGCAAGCCTGTGTCAGTGCCGGTAATACCGGGGCGCTGATGGGGCTGGCAAAATTATTACTCAAGCCCCTGGAGGGGATTGAGCGTCCGGCGCTGGTGACGGTATTACca                                                                                                                     >  2:321782/1‑149 (MQ=255)
                                            cccTGGAGCTGGTGAAAGAAGGTCGAGCGCAAGCCTGTGTCAGTGCCGGTAATACCGGGGCGCTGATGGGGCTGGCAAAATTATTACTCAAGCCCCTGGAGGGGATTGAGCGTCCGGCGCTGGTGACGGTATTACCACATCAGCAAAAg                                                                                                         <  3:786404/149‑1 (MQ=255)
                                            cccTGGAGCTGGTGAAAGAAGGTCGAGCGCAAGCCTGTGTCAGTGCCGGTAATACCGGGGCGCTGATGGGGCTGGCAAAATTATTACTCAAGCCCCTGGAGGGGATTGAGCGTCCGGCGCTGGTGACGGTATTACCACATCAGCAAAAg                                                                                                         <  1:23896/149‑1 (MQ=255)
                                             ccTGGAGCTGGTGAAAGAAGGTCGAGCGCAAGCCTGTGTCAGTGCCGGTAATACCGGGGCGCTGATGGGGCTGGCAAAATTATTACTCAAGCCCCTGGAGGGGATTGAGCGTCCGGCGCTGGTGACGGTATTACCACATCAGCAAAAgg                                                                                                        >  3:256141/1‑149 (MQ=255)
                                              cTGGAGCTGGTGAAAGAAGGTCGAGCGCAAGCCTGTGTCAGTGCCGGTAATACCGGGGCGCTGATGGGGCTGGCAAAATTATTACTCACGCCCCTGGAGGGGATTGAGCGTCCGGCGCTGGTGACGGTAtt                                                                                                                         >  3:554210/1‑131 (MQ=255)
                                              cTGGAGCTGGTGAAAGAAGGTCGAGCGCAAGCCTGTGTCAGTGCCGGTAATACCGGGGCGCTGATGGGGCTGGCAAAATTATTACTCAAGCCCCTGGAGGGGATTGAGCGTCCGGCGCTGGTGACGGTATTACCACATCAGCAAAAggg                                                                                                       >  1:720531/1‑149 (MQ=255)
                                              cTGGAGCTGGTGAAAGAAGGTCGAGCGCAAGCCTGTGTCAGTGCCGGTAATACCGGGGCGCTGAAGGGGCTGGCAAAATTATTACTCAAGCCCCTGGAGGGGATTGAGCGTCCGGCGCTGGTGACGGTAtt                                                                                                                         <  4:554210/131‑1 (MQ=255)
                                               tGGAGCTGGTGAAAGAAGGTCGAGCGCAAGCCTGTGTCAGTGCCGGTAATACCGGGGCGCTGATGGGGCTGGCAAAATTATTACTCAAGCCCCTGGAGGGGATTGAGCGTCCGGCGCTGGTGACGGTATTACCACATCAGCAAAAGGGc                                                                                                      >  1:528558/1‑149 (MQ=255)
                                               tGGAGCTGGTGAAAGAAGGTCGAGCGCAAGCCTGTGTCAGTGCCGGTAATACCGGGGCGCTGATGGGGCTGGCAAAATTATTACTCAAGCCCCTGGAGGGGATTGAGCGTCCGGCGCTGGTGACGGTATTACCACATCAGCAAAAGGGc                                                                                                      <  2:336244/149‑1 (MQ=255)
                                                   gCTGGTGAAAGAAGGTCGAGCGCAAGCCTGTGTCAGTGCCGGTAATACCGGGGCGCTGATGGGGCTGGCAAAATTATTACTCAAGCCCCTGGAGGGGATTGAGCGTCCGGCGCTGGTGACGGTATTACCACATCAGCAAAAGGGCaaaa                                                                                                  >  3:172852/1‑149 (MQ=255)
                                                     tGGTGAAAGAAGGTCGAGCGCAAGCCTGTGTCAGTGCCGGTAATACCGGGGCGCTGATGGGGCTGGCAAAATTATTACTCAAGCCCCTGGAGGGGATTGAGCGTCCGGCGCTGGTGACGGTATTACCACATCAGCAAAAGGGCAAAACg                                                                                                <  4:328706/149‑1 (MQ=255)
                                                      ggTGAAAGAAGGTCGAGCGCAAGCCTGTGTCAGTGCCGGTAATACCGGGGCGCTGATGGGGCTGGCAAAATTATTACTCAAGCCCCTGGAGGGGATTGAGCGTCCGGCGCTGGTGACGGTATTACCACATCAGCAAAAGGGCAAAACgg                                                                                               <  4:172852/149‑1 (MQ=255)
                                                       gTGAAAGAAGGTCGAGCGCAAGCCTGTGTTAGTGCCGGTAATACCGGGGCGCTGATGGGGCTGGCAAAATTATTACTCAAGCCCCTGGAGGGGATTGAGCGTCCGGCGCTGGTGACGGTATTACCACATCAGCAAAAGGGCAAAACggt                                                                                              >  3:652956/1‑149 (MQ=255)
                                                        tGAAAGAAGGTCGAGCGCAAGCCTGTGTCAGTGCCGGTAATACCGGGGCGCTGATGGGGCTGGCAAAATTATTACTCAAGCCCCTGGAGGGGATTGAGCGTCCGGCGCTGGTGACGGTATTACCACATCAGCAAAAGGGCAAAACggtg                                                                                             >  1:298211/1‑149 (MQ=255)
                                                         gAAAGAAGGTCGAGCGCAAGCCTGTGTCAGTGCCGGTAATACCGGGGCGCTGATGGGGCTGGCAAAATTATTACTCAAGCCCCTGGAGGGGATTGAGCGTCCGGCGCTGGTGACGGTATTACCACATCAGCAAAAGGGCAAAACggtgg                                                                                            <  2:528558/149‑1 (MQ=255)
                                                          aaaGAAGGTCGAGCGCAAGCCTGTGTCAGTGCCGGTAATACCGGGGCGCTGATGGGGCTGGCAAAATTATTACTCAAGCCCCTGGAGGGGATTGAGCGTCCGGCGCTGGTGACGGTATTACCACATCAGCAAAAGGGCAAAACggtggt                                                                                           >  3:185672/1‑149 (MQ=255)
                                                            agaagGTCGAGCGCAAGCCTGTGTCAGTGCCGGTAATACCGGGGCGCTGATGGGGCTGGCAAAATTATTACTCAAGCCCCTGGAGGGGATTGAGCGTCCGGCGCTGGTGACGGTATTACCACATCAGCAAAAGGGCAAAACGGTGGTcc                                                                                         >  1:94987/1‑149 (MQ=255)
                                                             gaagGTCGAGCGCAAGCCTGTGTCAGTGCCGGTAATACCGGGGCGCTGATGGGGCTGGCAAAATTATTACTCAAGCCCCTGGAGGGGATTGAGCGTCCGGCGCTGGTGACGGTATTACCACATCAGCAAAAGGGCAAAACGGTGGTCCt                                                                                        >  3:279436/1‑149 (MQ=255)
                                                                ggTCGAGCGCAAGCCTGTGTCAGTGCCGGTAATACCGGGGCGCTGATGGGGCTGGCAAAATTATTACTCAAGCCCCTGGAGGGGATTGAGCGTCCGGCGCTGGTGACGGTATTACCACATCAGCAAAAGGGCAAAACGGTGGTCCTTGa                                                                                     <  2:376643/149‑1 (MQ=255)
                                                                  tCGAGCGCAAGCCTGTGTCAGTGCCGGTAATACCGGGGCGCTGATGGGGCTGGCAAAATTATTACTCAAGCCCCTGGAGGGGATTGAGCGTCCGGCGCTGGTGACGGTATTACCACATCAGCAAAAGGGCAAAACGGTGGTCCTTGACt                                                                                   <  4:217829/149‑1 (MQ=255)
                                                                   cGAGCGCAAGCCTGTGTCAGTGCCGGTAATACCGGGGCGCTGATGGGGCTGGCAAAATTATTACTCAAGCCCCTGGAGGGGATTGAGCGTCCGGCGCTGGTGACGGTATTACCACATCAGCAAAAGGGCAAAACGGTGGTCCTTGACtt                                                                                  <  4:279436/149‑1 (MQ=255)
                                                                      gcgcAAGCCTGTGTCAGTGCCGGTAATACCGGGGCGCTGATGGGGCTGGCAAAATTATTACTCAAGCCCCTGGAGGGGATTGAGCGTCCGGCGCTGGTGACGGTATTACCACATCAGCAAAAGGGCAAAACGGTGGTCCTTGACTTAgg                                                                               <  1:269700/149‑1 (MQ=255)
                                                                      gcgcAAGCCTGTGTCAGTGCCGGTAATACCGGGGCGCTGATGGGGCTGGCAAAATTATTACTCAAGCCCCTGGAGGGGATTGAGCGTCCGGCGCTGGTGACGGTATTACCACATCAGCAAAAGGGCAAAACGGTGGTCCTTGACTTAgg                                                                               <  1:321782/149‑1 (MQ=255)
                                                                      gcgcAAGCCTGTGTCAGTGCCGGTAATACCGGGGCGCTGATGGGGCTGGCAAAATTATTACTCAAGCCCCTGGAGGGGATTGAGCGTCCGGCGCTGGTGACGGTATTACCACATCAGCAAAAGGGCAAAACGGTGGTCCTTGACTTAgg                                                                               <  2:19556/149‑1 (MQ=255)
                                                                          aaGCCTGTGTCAGTGCCGGTAATACCGGGGCGCTGATGGGGCTGGCAAAATTATTACTCAAGCCCCTGGAGGGGATTGAGCGTCCGGCGCTGGTGACGGTATTACCACATCAGCAAAAGGGCAAAACggtg                                                                                             >  4:300333/1‑131 (MQ=255)
                                                                          aaGCCTGTGTCAGTGCCGGTAATACCGGGGCGCTGATGGGGCTGGCAAAATTATTACTCAAGCCCCTGGAGGGGATTGAGCGTCCGGCGCTGGTGACGGTATTACCACATCAGCAAAAGGGCAAAACggtg                                                                                             <  3:300333/131‑1 (MQ=255)
                                                                              cTGTGTCAGTGCCGGTAATACCGGGGCGCTGATGGGGCTGGCAAAATTATTACTCAAGCCCCTGGAGGGGATTGAGCGTCCGGCGCTGGTGACGGTATTACCACATCAGCAAAAGGGCAAAACGGTGGTCCTTGACTTAGGGGCCAACg                                                                       >  3:483366/1‑149 (MQ=255)
                                                                                   tCAGTGCCGGTAATACCGGGGCGCTGATGGGGCTGGCAAAATTATTACTCAAGCCCCTGGAGGGGATTGAGCGTCCGGCGCTGGTGACGGTATTACCACATCAGCAAAAGGGCAAAACGGTGGTCCTTGACTTAGGGGTCAACGTCGag                                                                  >  3:500342/1‑148 (MQ=255)
                                                                                        gCCGGTAATACCGGGGCGGTGATGGGGCTGGCAAAATTATTACTCAAGCCCCTGGAGGGGATTGAGCGTCCGGCGCTGGTGACGGTATTACCACATCAGCAAAAGGGCAAAACGGTGGTCCTTGACTTAGGGGCCAACGTCGATTGTGa                                                             <  4:166758/149‑1 (MQ=255)
                                                                                        gCCGGTAATACCGGGGCGCTGATGGGGCTGGCAAAATTATTACTCAAGCCCCTGGAGGGGATTGAGCGTCCGGCGCTGGTGACGGTATTACCACATCAGCAAAAGGGCAAAACGGTGGTCCTTGACTTAGGGGCCAACGTCGATTGTGa                                                             <  4:652956/149‑1 (MQ=255)
                                                                                        gCCGGTAATACCGGGGCGCTGATGGGGCTGGCAAAATTATTACTCAAGCCCCTGGAGGGGATTGAGCGTCCGGCGCTGGTGACGGTATTACCACATCAGCAAAAGGGCAAAACGGTGGTCCTTGACTTAGGGGCCAACGTCGATTGTGa                                                             <  2:446254/149‑1 (MQ=255)
                                                                                          cGGTAATACCGGGGCGCTGATGGGGCTGGCAAAATTATTACTCAAGCCCCTGGAGGGGATTGAGCGTCCGGCGCTGGTGACGGTATTACCACATCAGCAAAAGGGCAAAACGGTGGTCCTTGACTTAGGGGCCAACGTCGATTGTGACa                                                           <  2:215842/149‑1 (MQ=255)
                                                                                           ggTAATACCGGGGCGCTGATGGGGCTGGCAAAATTATTACTCAAGCCCCTGGAGGGGATTGAGCGTCCGGCGCTGGTGACGGTATTACCACATCAGCAAAAGGGCAAAACggtg                                                                                             >  4:514289/1‑114 (MQ=255)
                                                                                           ggTAATACCGGGGCGCTGATGGGGCTGGCAAAATTATTACTCAAGCCCCTGGAGGGGATTGAGCGTCCGGCGCTGGTGACGGTATTACCACATCAGCAAAAGGGCAAAACggtg                                                                                             <  3:514289/114‑1 (MQ=255)
                                                                                              aaTACCGGGGCGCTGATGGGGCTGGCAAAATTATTACTCAAGCCCCTGGAGGGGATTGAGCGTCCGGCGCTGGTGACGGTATTACCACATCAGCAAAAGGGCAAAACGGTGGTCCTTGACTTAGGGGCCAACGTCGATTGTGACAGcac                                                       <  3:372441/149‑1 (MQ=255)
                                                                                                 aCCGGGGCGCTGATGGGGCTGGCAAAATTATTACTCAAGCCCCTGGAGGGGATTGAGCGTCCGGCGCTGGTGACGGTATTACCACATCAGCAAAAGGGCAAAACGGTGGTCCTTGACTTAGGGGCCAACGTCGATTGTGACAGCACAAt                                                    >  3:797985/1‑149 (MQ=255)
                                                                                                 aCCGGGGCGCTGATGGGGCTGGCAAAATTATTACTCAAGCCCCTGGAGGGGATTGAGCGTCCGGCGCTGGTGACGGTATTACCACATCAGCAAAAGGGCAAAACGGTGGTCCTTGACTTAGGGGCCAACGTCGATTGTGACAGCACAAt                                                    >  3:378098/1‑149 (MQ=255)
                                                                                                  ccGGGGCGCTGATGGTGCTGGCAAAATTATTACTCAAGCCCCTGGAGGGGATTGAGCGTCCGGCGCTGGTGACGGTATTACCACATCAGCAAAAGGGCAAAACGGTGGTCCTTGACTTAGGGGCCAACGTCGATTGTGACAGCACAATg                                                   >  4:173527/1‑149 (MQ=255)
                                                                                                   cGGGGCGCTGATGGGGCTGGCAAAATTATTACTCAAGCCCCTGGAGGGGATTGAGCGTCCGGCGCTGGTGACGGTATTACCACATCAGCAAAAGGGCAAAACGGTGGTCCTTGACTTAGGGGCCAACGTCGATTGTGACAGCACAATGc                                                  >  4:126218/1‑149 (MQ=255)
                                                                                                    ggggCGCTGATGGGGCTGGCAAAATTATTACTCAAGCCCCTGGAGGGGATTGAGCGTCCGGCGCTGGTGACGGTATTACCACATCAGCAAAAGGGCAAAACGGTGGTCCTTGACTTAGGGGCCAACGTCGa                                                                   >  1:367312/1‑131 (MQ=255)
                                                                                                    ggggCGCTGATGGGGCTGGCAAAATTATTACTCAAGCCCCTGGAGGGGATTGAGCGTCCGGCGCTGGTGACGGTATTACCACATCAGCAAAAGGGCAAAACGGTGGTCCTTGACTTAGGGGCCAACGTCGa                                                                   <  2:367312/131‑1 (MQ=255)
                                                                                                            gATGGGGCTGGCAAAATTATTACTCAAGCCCCTGGAGGGGATTGAGCGTCCGGCGCTGGTGACGGTATTACCACATCAGCAAAAGGGCAAAACGGTGGTCCTTGACTTAGGGGCCAACGTCGATTGTGACAGCACAATGCTGGTGCAAt                                         <  2:207292/149‑1 (MQ=255)
                                                                                                              tGGGGCTGGCAAAATTATTACTCAAGCCCCTGGAGGGGATTGAGCGTCCGGCGCTGGTGACGGTATTACCACATCAGCAAAAGGGCAAAACGGTGGTCCTTGACTTAGGGGCCAACGTCGATTGTGACAGCACAATGCTGGTGCAAttt                                       >  4:581700/1‑149 (MQ=255)
                                                                                                                gggCTGGCAAAATTATTACTCAAGCCCCTGGAGGGGATTGAGCGTCCGGCGCTGGTGACGGTATTACCACATCAGCAAAAGGGCAAAACGGTGGTCCTTGACTTAGGGGCCAACGTCGATTGTGACAGCACAATGCTGGTGCAATTTGc                                     >  1:348647/1‑149 (MQ=255)
                                                                                                                gggCTGGCAAAATTATTACTCAAGCCCCTGGAGGGGATTGAGCGTCCGGCGCTGGTGACGGTATTACCACATCAGCAAAAGGGCAAAACGGTGGTCCTTGACTTAGGGGCCAACGTCGATTGTGACAGCACAATGCTGGTGCAATTTGc                                     <  4:558862/149‑1 (MQ=255)
                                                                                                                  gCTGGCAAAATTATTACTCAAGCCCCTGGAGGGGATTGAGCGTCCGGCGCTGGTGACGGTATTACCACATCAGCAAAAGGGCAAAACGGTGGTCCTTGACTTAGGGGCCAACGTCGATTGTGACAGCACAATGCTGGt                                              <  4:712793/138‑1 (MQ=255)
                                                                                                                  gCTGGCAAAATTATTACTCAAGCCCCTGGAGGGGATTGAGCGTCCGGCGCTGGTGACGGTATTACCACATCAGCAAAAGGGCAAAACGGTGGTCCTTGACTTAGGGGCCAACGTCGATTGTGACAGCACAATGCTGGt                                              >  3:712793/1‑138 (MQ=255)
                                                                                                                      gCAAAATTATTACTCAAGCCCCTGGAGGGGATTGAGCGTCCGGCGCTGGTGACGGTATTACCACATCAGCAAAAGGGCAAAACGGTGGTCCTTGACTTAGGGGCCAACGTCGATTGTGACAGCACAATGCTGGTGCAATTTGCCATTAt                               <  4:256141/149‑1 (MQ=255)
                                                                                                                       cAAAATTATTACTCAAGCCCCTGGAGGGGATTGAGCGTCCGGCGCTGGTGACGGTATTACCACATCAGCAAAAGGGCAAAACggt                                                                                              >  3:514101/1‑85 (MQ=255)
                                                                                                                       cAAAATTATTACTCAAGCCCCTGGAGGGGATTGAGCGTCCGGCGCTGGTGACGGTATTACCACATCAGCAAAAGGGCAAAACggt                                                                                              <  4:514101/85‑1 (MQ=255)
                                                                                                                        aaaaTTATTACTCAAGCCCCTGGAGGGGATTGAGCGTCCGGCGCTGGTGACGGTATTACCACATCAGCAAAAGGGCAAAACGGTGGTCCTTGACTTAGGGGCCAACGTCGATTGTGACAGCACAATGCTGGTGCAATTTGCCATTATgg                             >  2:24537/1‑149 (MQ=255)
                                                                                                                           attattACTCAAGCCCCTGGAGGGGATTGAGCGTCCGGCGCTGGTGATGGTATTACCACATCAGCAAAAGGGCAAAACGGTGGTCCTTGACTTAGGGGCCAACGTCGATTGTGACAGCACAATGCTGGTGCAATTTGCCATTATGGGCt                          >  4:78633/1‑149 (MQ=255)
                                                                                                                           attattACTCAAGCCCCTGGAGGGGATTGAGCGTCCGGCGCTGGTGACGGTATTACCACATCAGCAAAAGGGCAAAACGGTGGTCCTTGACTTAGGGGCCAACGTCGATTGTGACAGCACAATGCTGGTGCAATTTGCCATTATGGGCt                          >  2:782142/1‑149 (MQ=255)
                                                                                                                             tattaCTCAAGCCCCTGGAGGGGATTGAGCGTCCGGCGCTGGTGACGGTATTACCACATCAGCAAAAGGGCAAAACGGTGGTCCTTGACTTAGGGGCCAACGTCGATTGTGACAGCACAATGCTGGTGCAATTTGCCATTATGGGCTCa                        <  2:94987/149‑1 (MQ=255)
                                                                                                                              gacgcTCAATCCCCTCCAGGGGCTTGAGCGTCCGGCGCTGGTGACGGTATTACCACATCAGCAAAAGGGCAAAACGGTGGTCCTTGACTTAGGGGCCAACGTCGATTGTGACAGCACAATGCTGGTGCAATTTGCCATTATGGGCTCAg                       >  1:313512/5‑149 (MQ=255)
                                                                                                                                  cTCAAGCCCCTGGAGGGGATTGAGCGTCCGGCGCTGGTGACGGTATTACCACATCAGCAAAAGGGCAAAACGGTGGTCCTTGACTTAGGGGCCAACGTCGATTGTGACAGCACAATGCTGGTGCAATTTGCCATTATGGGCTCAGTTCt                   >  4:139085/1‑149 (MQ=255)
                                                                                                                                  cTCAAGCCCCTGGAGGGGATTGAGCGTCCGGCGCTGGTGACGGTATTACCACATCAGCAAAAGGGCAAAACGGTGGTCCTTGACTTAGGGGCCAACGTCGATTGTGACAGCACAATGCTGGTGCAATTTGCCATTATGGGCTAAGTTCt                   >  2:311751/1‑149 (MQ=255)
                                                                                                                                   tCAATCCCCTCCAGGGGCTTGAGCGTCCGGCGCTGGTGACGGTATTACCACATCAGCAAAAGGGCAAAACGGTGGTCCTTGACTTAGGGGCCAACGTCGATTGTGACAGCACAATGCTGGTGCAATTTGCCATTATGGGCTCAGTTCTg                  >  4:43751/1‑149 (MQ=255)
                                                                                                                                   tCAAGCCCCTGGAGGGGATTGAGCGTCCGGCGCTGGTGATGGTATTACCACATCAGCAAAAGGGCAAAACGGTGGTCCTTGACTTAGGGGCCAACGTCGATTGTGACAGCACAATGCTGGTGCAATTTGCCATTATGGGCTCAGTTCTg                  <  3:78633/149‑1 (MQ=255)
                                                                                                                                      atccccTCCAGGGGCTTGAGCGTCCGGCGCTGGTGACGGTATTACCACATCAGCAAAAGGGCAAAACGGTGGTCCTTGACTTAGGGGCCAACGTCGATTGTGACAGCACAATGCTGGTGCAATTTGCCATTATGGGCTCAGTTCTGGCt               <  3:43751/147‑1 (MQ=255)
                                                                                                                                         cccTGGAGGGGATTGAGCGTCCGGCGCTGGTGACGGTATTACCACATCAGCAAAAGGGCAAAACGGTGGTCCTTGACTTAGGGGCCAACGTCGAt                                                                  >  2:254325/1‑95 (MQ=255)
                                                                                                                                         cccTGGAGGGGATTGAGCGTCCGGCGCTGGTGACGGTATTACCACATCAGCAAAAGGGCAAAACGGTGGTCCTTGACTTAGGGGCCAACGTCGAt                                                                  <  1:254325/95‑1 (MQ=255)
                                                                                                                                         cccTGGAGGGGATTGAGCGTCCGGCGCTGGTGACGGTATTACCACATCAGCAAAAGGGCAAAACGGTGGTCCTTGACTTAGGGGCCAACGTCGATTGTGACAGCACAATGCTGGTGCAATTTGCCATTATGGGCTCAGTTCTGGCTGaa            >  2:744243/1‑149 (MQ=255)
                                                                                                                                          ccTGGAGGGGATTGAGCGTCCGGCGCTGGTGACGGTATTACCACATCAGCAAAAGGGCAAAACGGTGGTCCTTGACTTAGGGGCCAACGTCGATTGTGACAGCACAATGCTGGTGCAATTTGCCATTATGGGCTCAGTTCTGGCTGAag           >  4:515615/1‑149 (MQ=255)
                                                                                                                                                 gggATTGAGCGTCCGGCGCTGGTGACGGTATTACCACATCAGCAAAAGGGCAAAACGGTGGTCCTTGACTTAGGGGCCAACGTCGATTGTGACAGCACAATGCTGGTGCAATTTGCCATTATGGGCTCAGTTCTGGCTGAAGAggtggt    >  3:247482/1‑149 (MQ=255)
                                                                                                                                                  ggATTGAGCGTCCGGCGCTGGTGACGGTATTACCACATCAGCAAAAGGGCAAAACGGTGGTCCTTGACTTAGGGGCCAACGTCGATTGTGACAGCACAATGCTGGTGCAATTTGCCATTATGGGCTCAGTTCTGGCTGAAGAGgtggtg   <  2:298211/149‑1 (MQ=255)
                                                                                                                                                   gATTGAGCGTCCGGCGCTGGTGACGGTATTACCACATCAGCAAAAGGGCAAAACGGTGGTCCTTGACTTAGGGGCCAACGTCGATTGTGACAGCACAATGCTGGTGCAATTTGCCATTATGGGCTCAGTTCTGGCTGAAGAGGtggtgg  >  2:20948/1‑149 (MQ=255)
                                                                                                                                                   gATTGAGCGTCCGGCGCTGGTGACGGTATTACCACATCAGCAAAAGGGCAAAACGGTGGTCCTTGACTTAGGGGCCAACGTCGATTGTGACAGCACAATGCTGGTGCAATTTGCCATTATGGGCTCAGTTCTGGCTGAAGAGGtggtgg  >  4:802330/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
TTCGCAAGCTATCCGCGCCAGTCGTGGGAGTTCAATGCGCGTGGCCCTGGAGCTGGTGAAAGAAGGTCGAGCGCAAGCCTGTGTCAGTGCCGGTAATACCGGGGCGCTGATGGGGCTGGCAAAATTATTACTCAAGCCCCTGGAGGGGATTGAGCGTCCGGCGCTGGTGACGGTATTACCACATCAGCAAAAGGGCAAAACGGTGGTCCTTGACTTAGGGGCCAACGTCGATTGTGACAGCACAATGCTGGTGCAATTTGCCATTATGGGCTCAGTTCTGGCTGAAGAGGTGGTGG  >  NC_000913/1147845‑1148140

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