Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,116,250 G→A 5.8% I111I (ATC→ATT)  lpxL ← lauryl‑acyl carrier protein (ACP)‑dependent acyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,116,2500GA 5.8% 227.7 / 3.6 69I111I (ATC→ATT) lpxLlauryl‑acyl carrier protein (ACP)‑dependent acyltransferase
Reads supporting (aligned to +/- strand):  ref base G (32/33);  new base A (2/2);  total (34/35)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.42e-01

CGTTCGGGCGATAAACGCCAATACCCGGTTCCTGCATACCAAACTGCCGCGCACCCAGCTCCAGTGTCAGAAAATGGATGCCAACTAACAGGATGCCGCGTTTTTGCGCCTGCACGTCACGAATGTGTTCCATGCCGATCACTTCCGTCCAGCGGGCGATTCGGCGGTCCGGCCAGAACCACGCCATGCCGGTTTCCATCAGGCCCATGCCAACGGATTCGAAATTCTTCACCACCATTTTACGGCGTTCTTGTTCGCTCATTTCCGGGAAGCACAGTTCCAGG  >  NC_000913/1116114‑1116397
                                                                                                                                        |                                                                                                                                                   
cGTTCGGGCGATAAACGCCAATACCCGGTTCCTGCATACCAAACTGCCGCGCACCCAGCTCCAGTGTCAGAAAATGGATGCCAACTAACAGGATGCCGCGTTTTTGCGCCTGCACGTCACGAATGTGTTCCATGCCGATCACTTCCGTc                                                                                                                                         <  2:747584/149‑1 (MQ=255)
  ttCGGGCGATAAACGCCAATACCCGGTTCCTGCATACCAAACTGCCGCGCACCCAGCTCCAGTGTCAGAAAATGGATGCCAACTAACAGGATGCCGCGTTTTTGCGCCTGCACGTCACGAATGTGTTCCATGCCGATCACTTCCGTCCa                                                                                                                                       >  4:158641/1‑149 (MQ=255)
  ttCGGGCGATAAACGCCAATACCCGGTTCCTGCATACCAAACTGCCGCGCACCCAGCTCCAGTGTCAGAAAATGGATGCCAACTAACAGGATGCCGCGTTTTTGCGCCTGCACGTCACGAATGTGTTCCATGCCAATCACTTCCGTCCa                                                                                                                                       >  2:403129/1‑149 (MQ=255)
  ttCGGGCGATAAACGCCAATACCCGGTTCCTGCATACCAAACTGCCGCGCACCCAGCTCCAGTGTCAGAAAATGGATGCCAACTAACAGGATGCCGCGTTTTTGCGCCTGCACGTCACGAATGTGTTCCATGCCAATCACTTCCGTCCa                                                                                                                                       >  4:756967/1‑149 (MQ=255)
    cGGGCGATAAACGCCAATACCCGGTTCCTGCATACCAAACTGCCGCGCACCCAGCTCCAGTGTCAGAAAATGGATGCCAACTAACAGGATGCCGCGTTTTTGCGCCTGCACGTCACGAATGTGTTCCATGCCGATCACTTCCGTCCAGc                                                                                                                                     >  3:7205/1‑149 (MQ=255)
       gCGATAAACGCCAATACCCGGTTCCTGCATACCAAACTGCCGCGCACCCAGCTCCAGTGTCAGAAAATGGATGCCAACTAACAGGATGCCGCGTTTTTGCGCCTGCACGTCACGAATGTGTTCCATGCCGATCACTTCCGTCCAGCggg                                                                                                                                  <  4:7205/149‑1 (MQ=255)
          aTAAACGCCAATACCCGGTTCCTGCATACCAAACTGCCGCGCACCCAGCTCCAGTGTCAGAAAATGGATGCCAACTAACAGGATGCCGCGTTTTTGCGCCTGCACGTCACGAATGTGTTCCATGCCGATCACTTCCGTCCAGCGGGCGa                                                                                                                               <  3:153415/149‑1 (MQ=255)
              aCGCCAATACCCGGTTCCTGCATACCAAACTGCCGCGCACCCAGCTCCAGTGTCAGAAAATGGATGCCAACTAACAGGATGCCGCGTTTTTGCGCCTGCACGTCACGAATGTGTTCCATGCCGATCACtt                                                                                                                                              <  3:409555/130‑1 (MQ=255)
              aCGCCAATACCCGGTTCCTGCATACCAAACTGCCGCGCACCCAGCTCCAGTGTCAGAAAATGGATGCCAACTAACAGGATGCCGCGTTTTTGCGCCTGCACGTCACGAATGTGTTCCATGCCGATCACtt                                                                                                                                              >  4:409555/1‑130 (MQ=255)
                gCCAATACCCGGTTCCTGCATACCAAACTGCCGCGCACCCAGCTCCAGTGTCAGAAAATGGATGCCAACTAACAGGATGCCGCGTTTTTGCGCCTGCACGTCACGAATGTGTTCCATGCCGATCACTTCCGTCCAGCGGGCGATTcggc                                                                                                                         >  4:92948/1‑149 (MQ=255)
                gCCAATACCCGGTTCCTGCATACCAAACTGCCGCGCACCCAGCTCCAGTGTCAGAAAATGGATGCCAACTAACAGGATGCCGCGTTTTTGCGCCTGCACGTCACGAATGTGTTCCATGCCGATCACTTCCGTCCAGCGGGCGATTcggc                                                                                                                         >  2:159825/1‑149 (MQ=255)
                  cAATACCCGGTTCCTGCATACCAAACTGCCGCGCACCCAGCTCCAGTGTCAGAAAATGGATGCCAACTAACAGGATGCCGCGTTTTTGCGCCTGCACGTCACGAATGTGTTCCATGCCGATCACTTCCGTCCAGCGGGCGATTcggcgg                                                                                                                       <  4:641877/149‑1 (MQ=255)
                    aTACCCGGTTCCTGCATACCAAACTGCCGCGCACCCAGCTCCAGTGTCAGAAAATGGATGCCAACTAACAGGATGCCGCGTTTTTGCGCCTGCACGTCACGAATGTGTTCCATGCCGATCACTTCCGTCCAGCGGGCGATTCGGCGGTc                                                                                                                     >  3:688085/1‑149 (MQ=255)
                     tACCCGGTTCCTGCATACCAAACTGCCGCGCACCCAGCTCCAGTGTCAGAAAATGGATGCCAACTAACAGGATGCCGCGTTTTTGCGCCTGCACGTCACGAATGTGTTCCATGCCGATCACTTCCGTCCAGCGGGCGATTCGGCGGTcc                                                                                                                    >  2:594989/1‑149 (MQ=255)
                     tACCCGGTTCCTGCATACCAAACTGCCGCGCACCCAGCTCCAGTGTCAGAAAATGGATGCCAACTAACAGGATGCCGCGTTTTTGCGCCTGCACGTCACGAATGTGTTCCATGCCGATCACTTCCGTCCAGCGGGCGATTCGGCGGTcc                                                                                                                    <  3:174509/149‑1 (MQ=255)
                        ccGGTTCCTGCATACCAAACTGCCGCGCACCCAGCTCCAGTGTCAGAAAATGGATGCCAACTAACAGGATGCCGCGTTTTTGCGCCTGCACGTCACGAATGTGTTCCATGCCAATCACTTCCGTCCAGCGGGCGATTCGGCGGTCCGGc                                                                                                                 <  3:756967/149‑1 (MQ=255)
                        ccGGTTCCTGCATACCAAACTGCCGCGCACCCAGCTCCAGTGTCAGAAAATGGATGCCAACTAACAGGATGCCGCGTTTTTGCGCCTGCACGTCACGAATGTGTTCCATGCCAATCACTTCCGTCCAGCGGGCGATTCGGCGGTCCGGc                                                                                                                 <  1:403129/149‑1 (MQ=255)
                            ttCCTGCATACCAAACTGCCGCGCACCCAGCTCCAGTGTCAGAAAATGGATGCCAACTAACAGGATGCCGCGTTTTTGCGCCTGCACGTCACGAATGTGTTCCATGCCGATCACTTCCGTCCAGCGGGCGATTCGGCGGTCCGGCCAGa                                                                                                             >  3:584208/1‑149 (MQ=255)
                             tCCTGCATACCAAACTGCCGCGCACCCAGCTCCAGTGTCAGAAAATGGATGCCAACTAACAGGATGCCGCGTTTTTGCGCCTGCACGTCACAAATGTGTTCCATGCCGATCACTTCCGTCCAGCGGGCGATTCGGCGGTCCGGCCAGaa                                                                                                            <  2:664243/149‑1 (MQ=255)
                                 gCATACCAAACTGCCGCGCACCCAGCTCCAGTGTCAGAAAATGGATGCCAACTAACAGGATGCCGCGTTTTTGCGCCTGCACGTCACGAATGTGTTCCATGCCGATCACTTCCGTCCAGCGGGCGATTCGGCGGTCCGGCCAGAACCAc                                                                                                        <  2:54249/149‑1 (MQ=255)
                                   aTACCAATCTGCCGCGCACCCAGCTCCAGTGTCAGAAAATGGATGCCAACTAACAGGATGCCGCGTTTTTGCGCCTGCACGTCACGAATGAGTTCCATGCCGATCACTTCCGTCCAGCGGTCGATTCGGCGGTCCGGCCAGAACCACGc                                                                                                      <  3:688268/149‑1 (MQ=255)
                                       cAAACTGCCGCGCACCCGGCTCCAGTGTCAGAAAATGGATGCCAACTAACAGGATGCCGCGTTTTTGCGCCTGCACGTCACGAATGTGTTCCATGCCGATCACTTCCGTCCAGCGGGCGATTCGGCGGTCCGGCCAGAACCACg                                                                                                       <  2:317212/144‑1 (MQ=255)
                                       cAAACTGCCGCGCACCCAGCTCCAGTGTCAGAAAATGGATGCCAACTAACAGGATGCCGCGTTTTTGCGCCTGCACGTCACGAATGTGTTCCATGCCGATCACTTCCGTCCAGCGGGCGATTCGGCGGTCCGGCCAGAACCACg                                                                                                       <  2:371759/144‑1 (MQ=255)
                                       cAAACTGCCGCGCACCCAGCTCCAGTGTCAGAAAATGGATGCCAACTAACAGGATGCCGCGTTTTTGCGCCTGCACGTCACGAATGTGTTCCATGCCGATCACTTCCGTCCAGCGGGCGATTCGGCGGTCCGGCCAGAACCACg                                                                                                       >  1:317212/1‑144 (MQ=255)
                                       cAAACTGCCGCGCACCCAGCTCCAGTGTCAGAAAATGGATGCCAACTAACAGGATGCCGCGTTTTTGCGCCTGCACGTCACGAATGTGTTCCATGCCGATCACTTCCGTCCAGCGGGCGATTCGGCGGTCCGGCCAGAACCACg                                                                                                       >  1:371759/1‑144 (MQ=255)
                                       cAAACTGCCGCGCACCCAGCTCCAGTGTCAGAAAATGGATGCCAACTAACAGGATGCCGCGTCTTTGCGCCTGCACGTCACGAATGTGTTCCATGCCGATCACTTCCGTCCAGCGGGCGATTCGGCGGTCCGGCCAGAACCACg                                                                                                       >  3:227591/1‑144 (MQ=255)
                                       cAAACTGCCGCGCACCCAGCTCCAGTGTCAGAAAATGGATGCCAACTAACAGGATGCCGCGTCTTTGCGCCTGCACGTCACGAATGTGTTCCATGCCGATCACTTCCGTCCAGCGGGCGATTCGGCGGTCCGGCCAGAACCACg                                                                                                       <  4:227591/144‑1 (MQ=255)
                                             gCCGCGCACCCAGCTCCAGTGTCAGAAAATGGATGCCAACTAACAGGATGCCGCGTTTTTGCGCCTGCACGTCACGAATGTGTTCCATGCCGATCACTTCCGTCCAGCGGGCGATTCGGCGGTCCGGCCAGAACCACGCCATGCCGGtt                                                                                            >  3:391802/1‑149 (MQ=255)
                                                gcgcACCCAGCTCCAGTGTCAGAAAATGGATGCCGACTAACAGGATGCCGCGTTTTTGCGCCTGCACGTCACGAATGTGTTCCATGCCGATCACTTCCGTCCAGCGGGCGa                                                                                                                               >  2:170271/1‑111 (MQ=255)
                                                gcgcACCCAGCTCCAGTGTCAGAAAATGGATGCCGACTAACAGGATGCCGCGTTTTTGCGCCTGCACGTCACGAATGTGTTCCATGCCGATCACTTCCGTCCAGCGGGCGa                                                                                                                               <  1:170271/111‑1 (MQ=255)
                                                  gcACCCAGCTCCAGTGTCAGAAAATGGATGACAACTAACAGGATGCCGCGTTTTTGCGCCTGCACGTCACGAATGTGTTCCATGCCGATCACTTCCGTCCAGCGGGCGATTCGGCGGTCCGGCCAGAACCACGCCATGCCGGTTTCCAt                                                                                       <  1:527610/149‑1 (MQ=255)
                                                    aCCCAGCTCCAGTGTCAGAAAATGGATGCCAACTAACAGGATGCCGCGTTTTTGCGCCTGCACGTCACGAATGTGTTCCATGCCGATCACTTCCGTCCAGCGGGCGATTCGGCGGTCCGGCCAGAACCACGCCATGCCGGTTTCCATCa                                                                                     <  3:92948/149‑1 (MQ=255)
                                                        aGCTCCAGTGTCAGAAAATGGATGCCAACTAACAGGATGCCGCGTTTTTGCGCCTGCACGTCACGAATGTGTTCCATGCCGATCACTTCCGTCCAGCGGGCGATTCGGCGGTCCGGCCAGAACCACGCCATGCCGGTTTCCATCAGGcc                                                                                 >  4:709739/1‑149 (MQ=255)
                                                         gCTCCAGTGTCAGAAAATGGATGCCAACTAACAGGATGCCGCGTTTTTGCGCCTGCACGTCACGAATGTGTTCCATGCCGATCACTTCCGTCCAGCGGGCGATTCGGCGGTCCGGCCAGAACCACGCCATGCCGGTTTCCATCAGGccc                                                                                <  3:158641/149‑1 (MQ=255)
                                                           tCCAGTGTCAGAAAATGGATGCCAACTAACAGGATGCCGCGTTTTTGCGCCTGCACGTCACGAATGTGTTCCATGCCGATCACTTCCGTCCAGCGGGCGATTCGGCGGTCCGGCCAGAACCACGCCATGCCGGTTTCCATCAGGCCCAt                                                                              <  3:699955/149‑1 (MQ=255)
                                                           tCCAGTGTCAGAAAATGGATGCCAACTAACAGGATGCCGCGTTTTTGCGCCTGCACGTCACGAATGTGTTCCATGCCGATCACTTCCGTCCAGCGGGCGATTCGGCGGTCCGGCCAGAACCACGCCATGCCGGTTTCCATCAGGCCCAt                                                                              >  2:588774/1‑149 (MQ=255)
                                                                 gtCAGAAAATGGATGCCAACTAACAGGATGCCGCGTTTTTGCGCCTGCACGTCACGAATGTGTTCCATGCCGATCACTTCCGTCCAGCGGGCGATTCGGCGGTCCGGCCAGAACCACGCCATGCCGGTTTCCATCAGGCCCATGCCAAc                                                                        <  2:812206/149‑1 (MQ=255)
                                                                 gtCAGAAAATGGATGCCAACTAACAGGATGCCGCGTTTTTGCGCCTGCACGTCACGAATGTGTTCCATGCCGATCACTTCCGTCCAGCGGGCGATTCGGCGGTCCGGCCAGAACCACGCCATGCCGGTTTCCATCAGGCCCATGCCAAc                                                                        <  2:441881/149‑1 (MQ=255)
                                                                  tCAGAAAATGGATGCCAACTAACAGGATGCCGCGTTTTTGCGCCTGCACGTCACGAATGTGTTCCATGCCGATCACTTCCGTCCAGCGGGCGATTCGGCGGTCCGGCCAGAACCACGCCATGCCGGTTTCCATCAGGCCCATGCCAACg                                                                       <  2:136708/149‑1 (MQ=255)
                                                                     gAAAATGGATGCCAACTAACAGGATGCCGCGTTTTTGCGCCTGCACGTCACGAATGTGTTCCATGCCGATCACTTCCGTCCAGCGGGCGATTCGGCGGTCCGTCCAGAACCACGCCATGCCGGTTTCCATCAGGCCCATGCCAACGGAt                                                                    >  1:237301/1‑149 (MQ=255)
                                                                     gAAAATGGATGCCAACTAACAGGATGCCGCGTTTTTGCGCCTGCACGTCACGAATGTGTTCCATGCCGATCACTTCCGTCCAGCGGGCGATTCGGCGGTCCGGCCAGAACCACGCCATGCCGGTTTCCATCAGGCCCATGCCAACGGAt                                                                    >  3:242432/1‑149 (MQ=255)
                                                                     gAAAATGGATGCCAACTAACAGGATGCCGCGTTTTTGCGCCTGCACGTCACGAATGTGTTCCATGCCGATCACTTCCGTCCAGCGGGCGATTCGGCGGTCCGGCCAGAACCACGCCATGCCGGTTTCCATCAGGCCCATGCCAACGGAt                                                                    >  2:387769/1‑149 (MQ=255)
                                                                       aaaTGGATGCCAACTAACAGGATGCCGCGTTTTTGCGCCTGCACGTCACGAATGTGTTCCATGCCGATCACTTCCGTCCAGCGGGCGATTCGGCGGTCCGGCCAGAACCACGCCATGCCGGTTTCCATCAGGCCCATGCCAACGGATTc                                                                  <  1:512200/149‑1 (MQ=255)
                                                                                 cAACTAACAGGATGCCGCGTTTTTGCGCCTGCACGTCACGAATGTGTTCCATGCCGATCACTTCCGTCCAGCGGGCGATTCGGCGGTCCGGCCAGAACCACGCCATGCCGGTTTCCATCAGGCCCATGCCAACGGATTCGAAAttcttc                                                        <  4:79070/149‑1 (MQ=255)
                                                                                 cAACTAACAGGATGCCGCGTTTTTGCGCCTGCACGTCACGAATGTGTTCCATGCCGATCACTTCCGTCCAGCGGGCGATTCGGCGGTCCGGCCAGAACCACGCCATGCCGGTTTCCATCAGGCCCATGCCAACGGATTCGAAAttcttc                                                        <  1:387769/149‑1 (MQ=255)
                                                                                  aaCTAACAGGATGCCGCGTTTTTGCGCCTGCACGTCACGAATGTGTTCCATGCCGATCACTTCCGTCCAGCGGGCGATTCGGCGGTCCGGCCAGAACCACGCCATGCCGGTTTCCATCAGGCCCATGCCAACGGATTCGAAATTCTTca                                                       <  3:369493/149‑1 (MQ=255)
                                                                                              gCCGCGTTTTTGCGCCTGCACGTCACGAATGTGTTCCATGCCGATCACTTCCGTCCAGCGGGCGATTCGGCGGTCCGGCCAGAACCACGCCATGCCGGTTTCCATCAGGCCCATGCCAACGGATTCGAAATTCTTCACCACCATTTTAc                                           <  3:731088/149‑1 (MQ=255)
                                                                                               ccGCGTTTTTGCGCCTGCACGTCACGAATGTGTTCCATGCCGATCACTTCCGTCCAGCGGGCGATTCGGCGGTCCGGCCAGAACCACGCCATGCCGGTTTCCATCAGGCCCATGCCAACGGATTCGAAATTCTTCACCACCATTTTACg                                          >  4:224914/1‑149 (MQ=255)
                                                                                                cgcgTTTTTGCGCCTGCACGTCACGAATGTGTTCCATGCCGATCACTTCCGTCCAGCGGGCGATTCGGCGGTCCGGCCAGAACCACGCCATGCCGGTTTCCATCAGGCCCATGCCAACGGATTCTAAATTCTTCACCACCATTTTACg                                          >  4:776782/1‑148 (MQ=255)
                                                                                                cgcgTTTTTGCGCCTGCACGTCACGAATGTGTTCCATGCCGATCACTTCCGTCCAGCGGGCGATTCGGCGGTCCGGCCAGAACCACGCCATGCCGGTTTCCATCAGGCCCATGCCAACGGATTCGAAATTCTTCACCACCATTTTACgg                                         >  1:262403/1‑149 (MQ=255)
                                                                                                cgcgTTTTTGCGCCTGCACGTCACGAATGTGTTCCATGCCGATCACTTCCGTCCAGCGGGCGATTCGGCGGTCCGGCCAGAACCACGCCATGCCGGTTTCCATCAGGCCCATGCCAACGGATTCGAAATTCTTCACCACCATTTTACgg                                         >  4:316870/1‑149 (MQ=255)
                                                                                                cgcgTTTTTGCGCCTGCACGTCACGAATGTGTTCCATGCCGATCACTTCCGTCCAGCGGGCGATTCGGCGGTCCGGCCAGAACCACGCCATGCCGGTTTCCATCAGGCCCATGCCAACGGATTCGAAATTCTTCACCACCATTTTACg                                          <  3:776782/148‑1 (MQ=255)
                                                                                                 gcgTTTTTGCGCCTGCACGTCACGAATGTGTTCCATGCCGATCACTTCCGTCCAGCGGGCGATTCGGCGGTCCGGCCAGAACCACGCCATGCCGGTTTCCATCAGGCCCATGCCAACGGATTCGAAATTCTTCACCACCATTTTACGGc                                        >  1:580663/1‑149 (MQ=255)
                                                                                                 gcgTTTTTGCGCCTGCACGTCACGAATGTGTTCCATGCCGATCACTTCCGTCCAGCGGGCGATTCGGCGGTCCGGCCAGAACCACGCCATGCCGGTTTCCATCAGGCCCATGCCAACGGATTCGAAATTCTTCACCACCATTTTACGGc                                        >  2:805642/1‑149 (MQ=255)
                                                                                                       ttGCGCCTGCACGTCACGAATGTGTTCCATGCCGATCACTTCCGTCCAGCGGGCGATTCGGCGGTCCGGCCAGAACCACGCCATGCCGGTTTCCATCAGGCCCATGCCAACGGATTCGAAATTCTTCACCACCATTTTACGGCGTTCtt                                  <  3:316870/149‑1 (MQ=255)
                                                                                                         gcgcCTGCACGTCACGAATGTGTTCCATGCCGATCACTTCCGTCCAGCGGGCGATTCGGCGGTCCGGCCAGAACCACGCCATGCCGGTTTCCATCAGGCCCATGCCAACGGATTCGAAATTCTTCACCACCATTTTACGGCGTTCTTGt                                <  3:184542/149‑1 (MQ=255)
                                                                                                           gcCTGCACGTCACGAATGTGTTCCATGCCGATCACTTCCGTCCAGCGGGCGATTCGGCGGTCCGGCCAGAACCACGCCATGCCGGTTTCCATCAGGCCCATGCCAACGGATTCGAAATTCTTCACCACCATTTTACGGCGTTCTTGTTc                              >  3:143568/1‑149 (MQ=255)
                                                                                                                   gTCACGAATGTGTTCCATGCCGATCACTTCCGTCCAGCGGGCGATTCGGCGGTCCGGCCAGAACCACGCCATGCCGGTTTCCATCAGGCCCATGCCAACGGATTCGAAATTCTTCACCACCATTTTACGGCGTTCTTGTTCGCTCAttt                      >  4:781492/1‑149 (MQ=255)
                                                                                                                        gAATGTGTTCCATGCCGATCACTTCCGTCCAGCGGGCGATTCGGCGGTCCGGCCAGAACCACGCCATGCCGGTTTCCATCAgg                                                                                   >  2:239465/1‑83 (MQ=255)
                                                                                                                        gAATGTGTTCCATGCCGATCACTTCCGTCCAGCGGGCGATTCGGCGGTCCGGCCAGAACCACGCCATGCCGGTTTCCATCAgg                                                                                   <  1:239465/83‑1 (MQ=255)
                                                                                                                        gAATGTGTTCCATGCCGATCACTTCCGTCCAGCGGGCGATTCGGCGGTCCGGCCAGAACCACGCCATGCCGGTTTCCATCAGGCCCATGCCAACGGATTCGAAATTCTTCACCACCATTTTACGGCGTTCTTGTTCGCTCATTTCCggg                 >  1:449738/1‑149 (MQ=255)
                                                                                                                          aTGTGTTCCATGCCGATCACTTCCGTCCAGCGGGCGATTCGGCGGTCCGGCCAGAACCACGCCATGCCGGTTTCCATCAGGCCCATGCCAACGGATTCGAAATTCTTCACCACCATTTTACGGCGTTCTTGTTCGCTCATTTCCGGGaa               <  3:224914/149‑1 (MQ=255)
                                                                                                                           tgtgTTCCATGCCGATCACTTCCGTCCAGCGGGCGATTCGGCGGTCCGGCCAGAACCACGCCATGCCGGTTTCCATCAGGCCCATGCCAACGGATTCGAAATTCTTCACCACCATTTTACGGCGTTCTTGTTCGCTCATTTCCGGGAAg              >  2:659027/1‑149 (MQ=255)
                                                                                                                            gtgtTCCATGCCGATCACTTCCGTCCAGCGGGCGATTCGGCGGTCCGGCCAGAACCACGCCATGCCGGTTTCCATCAGGCCCATGCCAACGGATTCGAAATTCTTCACCACCATTTTACGGCGTTCTTGTTCGCTCATTTCCGGGAAGc             >  3:513279/1‑149 (MQ=255)
                                                                                                                               ttCCATGCCGATCACTTCCGTCCAGCGGGCGATTCGGCGGTCCGGCCAGAACCACGCCATGCCGGTTTCCATCAGGCCCATGCCAACGGATTCGAAATTCTTCACCACCATTTTACg                                          >  1:103557/1‑117 (MQ=255)
                                                                                                                               ttCCATGCCGATCACTTCCGTCCAGCGGGCGATTCGGCGGTCCGGCCAGAACCACGCCATGCCGGTTTCCATCAGGCCCATGCCAACGGATTCGAAATTCTTCACCACCATTTTACg                                          <  2:103557/117‑1 (MQ=255)
                                                                                                                                 ccATGCCGATCACTTCCGTCCAGCGGGCGATTCGGCGGTCCGGCCAGAACCACGCCATGCCGGTTTCCATCAGGCCCATGCCAACGGATTCGAAATTCTTCACCACCATTTTACGGCGTTCTTGTTCGCTCATTTCCGGGAAGCACAGt        <  1:159825/149‑1 (MQ=255)
                                                                                                                                      ccGATCACTTCCGTCCAGCGGGCGATTCGGCGGTCCGGCCAGAACCACGCCATGCCGGTTTCCATCAGGCCCATGCCAACGGATTCGAAATTCTTCACCACCATTTTACGGCGTTCTTGTTCGCTCATTTCCGGGAAGCACAGTTCCAg   <  4:584208/149‑1 (MQ=255)
                                                                                                                                       cGATCACTTCCGTCCAGCGGGCGATTCGGCGGTCCGGCCAGAACCACGCCATGCCGGTTTCCATCAGGCCCATGCCAACGGATTCGAAATTCTTCACCACCATTTTACGGCGTTCTTGTTCGCTCATTTCCGGGAAGCACAGTTCCAgg  >  2:782806/1‑149 (MQ=255)
                                                                                                                                        gATCACTTCCGTCCAGCGGGCGATTCGGCGGTCCGGCCAGAACCACGCCATGCCGGTTTCCATCAGGCCCATGCCAACGGATTCGAAATTCTTcaccac                                                   <  2:567317/99‑1 (MQ=255)
                                                                                                                                        gATCACTTCCGTCCAGCGGGCGATTCGGCGGTCCGGCCAGAACCACGCCATGCCGGTTTCCATCAGGCCCATGCCAACGGATTCGAAATTCTTcaccac                                                   >  1:567317/1‑99 (MQ=255)
                                                                                                                                        |                                                                                                                                                   
CGTTCGGGCGATAAACGCCAATACCCGGTTCCTGCATACCAAACTGCCGCGCACCCAGCTCCAGTGTCAGAAAATGGATGCCAACTAACAGGATGCCGCGTTTTTGCGCCTGCACGTCACGAATGTGTTCCATGCCGATCACTTCCGTCCAGCGGGCGATTCGGCGGTCCGGCCAGAACCACGCCATGCCGGTTTCCATCAGGCCCATGCCAACGGATTCGAAATTCTTCACCACCATTTTACGGCGTTCTTGTTCGCTCATTTCCGGGAAGCACAGTTCCAGG  >  NC_000913/1116114‑1116397

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