Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 469,577 G→T 5.7% R236L (CGG→CTG)  mdlA → ABC transporter family protein MdlA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913469,5770GT 5.7% 236.7 / 2.9 70R236L (CGG→CTG) mdlAABC transporter family protein MdlA
Reads supporting (aligned to +/- strand):  ref base G (33/33);  new base T (2/2);  total (35/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 = 9.73e-01

CTGGCACAGGCGGCGTTTTCCAGTCTTAATGACCGCACCCAGGAAAGCCTCACCAGTATCCGCATGATCAAAGCCTTTGGTCTGGAAGATCGCCAGTCGGCGTTATTTGCCGCGGATGCCGAAGATACCGGCAAAAAAAACATGCGGGTGGCGCGTATTGATGCTCGTTTCGACCCGACCATCTATATCGCGATTGGTATGGCGAACTTGCTGGCGATTGGCGGCGGTAGTTGGATGGTGGTGCAGGGCAGTTTAACGCTGGGCCAGCTCACCAGTTTTATGATGTATTT  >  NC_000913/469432‑469721
                                                                                                                                                 |                                                                                                                                                
cTGGCACAGGCGGCGTTTTCCAGTCTTAATGACCGCACCCAGGAAAGCCTCACCAGTATCCGCATGATCAAAGCCTTTGGTCTGGAAGATCGCCAGTCGGCGTTATTTGCCGCGGATGCCGAAGATACCGGCAAAAAAAACATGCGGGTg                                                                                                                                              <  2:101763/150‑1 (MQ=255)
 tGGCACAGGCGGCGTTTTCCAGTCTTAATGACCGCACCCAGGAAAGCCTCACCAGTATCCGCATGATCAAAGCCTTTGGTCTGGAAGATCGCCAGTCGGCGTTATTTGCCGCGGATGCCGAAGATACCGGCAAAAAAAACATGCGGGTgg                                                                                                                                             >  2:1741276/1‑150 (MQ=255)
     acaGGCGGCGTTTTCCAGTCTTAATGACCGCACCCAGGAAAGCCTCACCAGTATCCGCATGATCAAAGCCTTTGGTCTGGAAGATCGCCAGTCGGCGTTATTTGCCGCGGATGCCGAAGATACCGGCAAAAAAAACATGCGGGTGGcgcg                                                                                                                                         <  2:664961/150‑1 (MQ=255)
      caGGCGGCGTTTTCCAGTCTTAATGACCGCACCCAGGAAAGCCTCACCAGTATCCGCATGATCAAAGCCTTTGGTCTGGAAGATCGCCAGTCGGCGTTATTTGCCGCGGATGCCGAAGATACCGGCAAAAAAAACATGCGGGTGGCGCGt                                                                                                                                        >  2:1813682/1‑150 (MQ=255)
        ggcggcGTTTTCCAGTCTTAATGACCGCACCCAGGAAAGCCTCACCAGTATCCGCATGATCAAAGCCTTTGGTCTGGAAGATCGCCAGTCGGCGTTATTTGCCGCTGATGCCGAAGATACCGGCAAAAAAAACATTCGGGTGGCGCGTAt                                                                                                                                      <  2:1760374/150‑1 (MQ=255)
         gcggcgTTTTCCAGTCTTAATGACCGCACCCAGGAAAGCCTCACCAGTATCCGCATGATCAAAGCCTTTGGTCTGGAAGATCGCCAGTCGGCGTTATTTGCCGCGGATGCCGAAGATACCGGCAAAAAAAACATGCGGGTGGCGCGTAtt                                                                                                                                     <  2:1667990/150‑1 (MQ=255)
          cggcgTTTTCCAGTCTTAATGACCGCACCCAGGAAAGCCTCACCAGTATCCGCATGATCAAAGCCTTTGGTCTGGAAGATCGCCAGTCGGCGTTATTTGCCGCGGATGCCGAAGATACCGGCAAAAAAAACATGCGGGTGGCGCGTATTg                                                                                                                                    <  1:1703927/150‑1 (MQ=255)
           ggcgTTTTCCAGTCTTAATGACCGCACCCAGGAAAGCCTCACCAGTATCCGCATGATCAAAGCCTTTGGTCTGGAAGATCGCCAGTCGGCGTTATTTGCCGCGGATGCCGAAGATACCGGCAAAAAAAACATGCGGGTGGCGCGTATTGa                                                                                                                                   <  1:227490/150‑1 (MQ=255)
               ttttCCAGTCTTAATGACCGCACCCAGGAAAGCCTCACCAGTATCCGCATGATCAAAGCCTTTGGTCTGGAAGATCGCCAGTCGGCGTTATTTGCCGCGGATGCCGAAGATACCGGCAAAAAAAACATGCGGGTGGCGCGTATTGATGCt                                                                                                                               >  1:958386/1‑150 (MQ=255)
                tttCCAGTCTTAATGACCGCACCCAGGAAAGCCTCACCAGTATCCGCATGATCAAAGCCTTTGGTCTGGAAGATCGCCAGTCGGCGTTATTTGCCGCGGATGCCGAAGATACCGGCAAAAAAAACATGCGGGTGGCGCGTATTGATGCTc                                                                                                                              <  1:1375849/150‑1 (MQ=255)
                   ccAGTCTTAATGACCGCACCCAGGAAAGCCTCACCAGTATCCGCATGATCAAAGCCTTTGGTCTGGAAGATCGCCAGTCGGCGTTATTTGCCGCGGATGCCGAAGATACCGGCAAAAAAAACATGCGGGTGGCGCGTATTGATGCTCGtt                                                                                                                           >  1:1870744/1‑150 (MQ=255)
                     aGTCTTAATGACCGCACCCAGGAAAGCCTCACCAGTATCCGCATGATCAAAGCCTTTGGTCTGGAAGATCGCCAGTCGGCGTTATTTGCCGCGGATGCCGAAGATACCGGCAAAAAAAACATGCGGGTGGCGCGTATTGATGCTCGTTTc                                                                                                                         <  2:1870744/150‑1 (MQ=255)
                     aGTCTTAATGACCGCACCCAGGAAAGCCTCACCAGTATCCGCATGATCAAAGCCTTTGGTCTGGAAGATCGCCAGTCGGCGTTATTTGCCGCGGATGCCGAAGATACCGGCAAAAAAAACATGCGGGTGGCGCGTATTGATGCTCGTTTc                                                                                                                         <  1:511587/150‑1 (MQ=255)
                     aGTCTTAATGACCGCACCCAGGAAAGCCTCACCAGTATCCGCATGATCAAAGCCTTTGGTCTGGAAGATCGCCAGTCGGCGTTATTCGCCGCGGATGCCGAAGATACCGGCAAAAAAAACATGCGGGTGGCGCGTATTGATGCTCGTTTc                                                                                                                         >  2:613608/1‑150 (MQ=255)
                          ttaTGACCGCACCCAGGAAAGCCTCACCAGTATCCGCATGATCAAAGCCTTTGGTCTGGAAGATCGCCAGTCGGCGTTATTTGCCGCGGATGCCGAAGATACCGGCAAAAAAAACATGCGGGTGGCGCGTATTGATGCTCGTTTcgacct                                                                                                                    <  2:920231/148‑2 (MQ=255)
                                  gCACCCAGGAAAGCCTCACCAGTATCCGCATGATCAAAGCCTTTGGTCTGGAAGATCGCCAGTCGGCGTTATTTGCCGCGGATGCCGAAGATACCGGCAAAAAAAACATGCGGGTGGCGCGTATTGATGCTCGTTTCGACCCGACCATCt                                                                                                            <  2:958386/150‑1 (MQ=255)
                                      ccAGGAAAGCCTCACCAGTATCCGCATGATCAAAGCCTTTGGTCTGGAAGATCGCCAGTCGGCGTTATTTGCCGCGGATGCCGAAGATACCGGCAAAAAAAACATGCGGGTGGCGCGTATTGATGCTCGTTTCGACCCGACCATCTatat                                                                                                        <  1:1861185/150‑1 (MQ=255)
                                       cAGGAAAGCCTCACCAGTATCCGCATGATCAAAGCCTTTGGTCTGGAAGATCGCCAGTCGGCGTTATTTGCCGCGGATGCCGAAGATACCGGCAAAAAAACATGCGGGTGGCGCGTATTGATGCTCGTTTCGACCCGACCATCTATATcg                                                                                                      >  2:1064441/1‑150 (MQ=255)
                                         ggAAAGCCTCACCAGTATCCGCATGATCAAAGCCTTTGGTCTGGAAGATCGCCAGTCGGCGTTATTTGCCGCGGATGCCGAAGATACCGGCAAAAAAAACATGCGGGTGGCGCGTATTGATGCTCGTTTCGACCCGACCATCTATATcgc                                                                                                     <  1:833248/150‑1 (MQ=255)
                                              gCCTCACCAGTATCCGCATGATCAAAGCCTTTGGTCTGGAAGATCGCCAGTCGGCGTTATTTGCCGCGGATGCCGAAGATACCGGCAAAAAAAACATGCGGGTGGCGCGTATTGATGCTCGTTTCGACCCGACCATCTATATCGCGATTg                                                                                                >  2:1706016/1‑150 (MQ=255)
                                               ccTCACCAGTATCCGCATGATCAAAGCCTTTGGTCTGGAAGATCGCCAGTCGGCGTTATTTGCCGCGGATGCCGAAGATACCGGCAAAAAAAACATGCGGGTGGCGCGTATTGATGCTCGTTTCGACCCGACCATCTATATCGCGATTgg                                                                                               <  1:978768/150‑1 (MQ=255)
                                                  cACCAGTATCCGCATGATCAAAGCCTTTGGTCTGGAAGATCGCCAGTCGGCGTTATTTGCCGCGGATGCCGAAGATACCGGCAAAAAAAACATGCGGGTGGCGCGTATTGATGCTCGTTTCGACCCGACCATCTATATCGCGATTGGTAt                                                                                            >  2:84392/1‑150 (MQ=255)
                                                   aCCAGTATCCGCATGATCAAAGCCTTTGGTCTGGAAGATCGCCAGTCGGCGTTATTTGCCGCGGATGCCGAAGATACCGGCAAAAAAAACATGCGGGTGGCGCGTATTGATGCTCGTTTCGACCCGACCATCTATATCGCGATTGGTATg                                                                                           <  1:1741276/150‑1 (MQ=255)
                                                      attATCCGCATCATCAAAGCATTTGCTCTGGAAGATCTCCAGTCGGCGTTATTTGCCGCGGATGCCGAAGATACCCGAAAAAAAAACATCCGTGAGCCGATTATTTATGCACCTATCGACCCGACCATCTATATCCCGATTGGTATggag                                                                                        >  2:1823879/3‑148 (MQ=255)
                                                         aTCCGCATGATCAAAGCCTTTGGTCTGGAAGATCGCCAGTCGGCGTTATTTGCCGCGGATGCCGAAGATACCGGCAAAAAAAACATGCGGGTGGCGCGTATTGATGCTCGTTTCGACCCGACCATCTATATCGCGATTGGTATGGCGAAc                                                                                     >  1:661319/1‑150 (MQ=255)
                                                         aTCCGCATGATCAAAGCCTTTGGTCTGGAAGATCGCCAGTCGGCGTTATTTGCCGCGGATGCCGAAGATACCGGCAAAAAAAACATGCGGGTGGCGCGTATTGATGCTCGTTTCGACCCGACCATCTATATCGCGATTGGTATGGCGAAc                                                                                     >  2:206244/1‑150 (MQ=255)
                                                           ccGCATGATCAAAGCCTTTGGTCTGGAAGATCGCCAGTCGGCGTTATTTGCCGCGGATGCCGAAGATACCGGCAAAAAAAACATGCGGGTGGCGCGTATTGATGCTCGTTTCGACCCGACCATCTATATCGCGATTGGTATGGCGAACtt                                                                                   >  2:795906/1‑150 (MQ=255)
                                                            cGCATGATCAAAGCCTTTGGTCTGGAAGATCGCCAGTCGGCGTTATTTGCCGCGGATGCCGAAGATACCGGCAAAAAAAACATGCGGGTGGCGCGTATTGATGCTCGTTTCGACCCGACCATCTATATCGCGATTGGTATGGCGAACTTg                                                                                  >  1:1034351/1‑150 (MQ=255)
                                                            cGCATGATCAAAGCCTTTGGTCTGGAAGATCGCCAGTCGGCGTTATTTGCCGCGGATGCCGAAGATACCGGCAAAAAAAACATGCGGGTGGCGCGTATTGATGCTCGTTTCGACCCGACCATCTATATCGCGATTGGTATGGCGAACTTg                                                                                  >  1:1041259/1‑150 (MQ=255)
                                                              cATGATCAAAGCCTTTGGTCTGGAAGATCGCCAGTCGGCGTTATTCGCCGCGGATGCCGAAGATACCGGCAAAAAAAACATGCGGGTGGCGCGTATTGATGCTCGTTTCGACCCGACCATCTATATCGCGATTGGTATGGCGAACTTGCt                                                                                <  1:613608/150‑1 (MQ=255)
                                                               aTGATCAAAGCCTTTGGTCTGGAAGATCGCCAGTCGGCGTTATTTGCCGCGGATGCCGAAGATACCGGCAAAAAAAACATGCGGGTGGCGCGTATTGATGCTCGTTTCGACCCGACCATCTATATCGCGATTGGTATGGCGAACTTGCTg                                                                               <  2:1915637/150‑1 (MQ=255)
                                                                tGATCAAAGCCTTTGGTCTGGAAGATCGCCAGTCGGCGTTATTTGCCGCGGATGCCGAAGATACCGGCAAAAAAAACATGCGGGTGGCGCGTATTGATGCTCGTTTCGACCCGACCATCTATATCGCGATTGGTATGGCGAACTTGCTgg                                                                              <  1:1128062/150‑1 (MQ=255)
                                                                     aaaGCCTTTGGTCTGGAAGATCGCCAGTCGGCGTTATTTGCCGCGGATGCCGAAGATACCGGCAAAAAAAACATGCGGGTGGCGCGTATTGATGCTCGTTTCGACCCGACCATCTATATCGCGATTGGTATGGCGAACTTGCTGGCGAtt                                                                         >  2:882365/1‑150 (MQ=255)
                                                                     aaaGCCTTTGGTCTGGAAGATCGCCAGTCGGCGTTATTTGCCGCGGATGCCGAAGATACCGGCAAAAAAAACATGCGGGTGGCGCGTATTGATGCTCGTTTCGACCCGACCATCTATATCGCGATTGGTATGGCGAACTTGCTGGCGAtt                                                                         >  1:1306906/1‑150 (MQ=255)
                                                                             tGGTCTGGAAGATCGCCAGTCGGCGTTATTTGCCGCGGATGCCGAAGATACCGGCAAAAAAAACATGCGGGTGGCGCGTATTGATGCTCGTTTCGACCCGACCATCTATATCGCGATTGGTATGGCGAACTTGCTGGCGATTGGcggcgg                                                                 <  1:867766/150‑1 (MQ=255)
                                                                                 cTGGAAGATCGCCAGTCGGCGTTATTTGCCGCGGATGCCGAAGATACCGGCAAAAAAAACATGCGGGTGGCGCGTATTGATGCTCGTTTCGACCCGACCATCTATATCGCGATTGGTATGGCGAACTTGCTGGCGATTGGCGGCGGTAGt                                                             >  2:1855947/1‑150 (MQ=255)
                                                                                   ggAAGATCGCCAGTCGGCGTTATTTGCCGCGGATGCCGAAGATACCGGCAAAAAAAACATGCGGGTGGCGCGTATTGATGCTCGTTTCGACCCGACCATCTATATCGCGATTGGTATGGCGAACTTGCTGGCGATTGGCGGCGGTAGTTg                                                           >  2:1785556/1‑150 (MQ=255)
                                                                                   ggAAGATCGCCAGTCGGCGTTATTTGCCGCGGATGCCGAAGATACCGGCAAAAAAAACATGCGGGTGGCGCGTATTGATGCTCGTTTCGACCCGACCATCTATATCGCGATTGGTATGGCGAACTTGCTGGCGATTGGCGGCGGTAGTTg                                                           >  1:1148231/1‑150 (MQ=255)
                                                                                    gAAGATCGCCAGTCGGCGTTATTTGCCGCGGATGCCGAAGATACCGGCAAAAAAAACATGCGGGTGGCGCGTATTGATGCTCGTTTCGACCCGACCATCTATATCGCGATTGGTATGGCGAACTTGCTGGCGATTGGCGGCGGTAGTTgg                                                          >  1:944222/1‑150 (MQ=255)
                                                                                    gAAGATCGCCAGTCGGCGTTATTTGCCGCGGATGCCGAAGATACCGGCAAAAAAAAAATGCGGGTGGCGCGTATTGATGTTCGTTTCGACCCCACCATCTATATAGCGATTGGTATGGCGAACTTGATGGCGATTGGCGGCGGTAGTTgg                                                          >  2:663924/1‑150 (MQ=255)
                                                                                      aGATCGCCAGTCGGCGTTATTTGCCGCGGATGCCGAAGATACCGGCAAAAAAAACATGCGGGTGGCGCGTATTGATGCTCGTTTCGACCCGACCATCTATATCGCGATTGGTATGGCGAACTTGCTGGCGATTGGCGGCGGTAGTTGGAt                                                        <  2:1195498/150‑1 (MQ=255)
                                                                                          cGCCAGTCGGCGTTATTTGCCGCGGATGCCGAAGATACCGGCAAAAAAAACATGCGGGTGGCGCGTATTGATGCTCGTTTCGACCCGACCATCTATATcgc                                                                                                     <  2:1327395/101‑1 (MQ=255)
                                                                                          cGCCAGTCGGCGTTATTTGCCGCGGATGCCGAAGATACCGGCAAAAAAAACATGCGGGTGGCGCGTATTGATGCTCGTTTCGACCCGACCATCTATATcgc                                                                                                     >  1:1327395/1‑101 (MQ=255)
                                                                                           gCCAGTCGGCGTTATTTGCCGCGGATGCCGAAGATACCGGCAAAAAAAACATGCGGGTGGCGCGTATTGATGCTCGTTTCGACCCGACCATCTATATCGCGATTGGTATGGCGAACTTGCTGGCGATTGGCGGCGGTAGTTGGAtggtgg                                                   <  2:1306906/150‑1 (MQ=255)
                                                                                           gCCAGTCGGCGTTATTTGCCGCGGATGCCGAAGATACCGGCAAAAAAAACATGCGGGTGGCGCGTATTGATGCTCGTTTCGACCCGACCATCTATATCGCGATTGGTATGGCGAACTTGCTGGCGATTGGCGGCGGTAGTTGGAtggtgg                                                   <  1:1823879/150‑1 (MQ=255)
                                                                                                   gCGTTATTTGCCGCGGATGCCGAAGATACCGGCAAAAAAAACATGCGGGTGGCGCGTATTGATGCTCGTTTCGACCCGACCATCTATATCGCGATTGGTATGGCGAACTTGCTGGCGATTGGCGGCGGTAGTTGGATGGTGGTGCAGGGc                                           <  1:1813682/150‑1 (MQ=255)
                                                                                                   gCGTTATTTGCCGCGGATGCCGAAGATACCGGCAAAAAAAACATGCGGGTGGCGCGTATTGATGCTCGTTTCGACCCGACCATCTATATCGCGATTGGTATGGCGAACTTGCTGGCGATTGGCGGCGGTAGTTGGATGGTGGTGCAGGGc                                           <  2:506071/150‑1 (MQ=255)
                                                                                                   gCGTTATTTGCCGCGGATGCCGAAGATACCGGCAAAAAAAACATGCGGGTGGCGCGTATTGATGCTCGTTTCGACCCGACCATCTATATCGCGATTGGTATGGCGAACTTGCTGGCGATTGGCGGCGGTAGTTGGATGGTGGTGCAGGGc                                           <  1:84392/150‑1 (MQ=255)
                                                                                                       tATTTGCCGCGGATGCCGAAGATACCGGCAAAAAAAACATGCGGGTGGCGCGTATTGATGCTCGTTTCGACCCGACCATCTATATCGCGATTGGTATGGCGAACTTGCTGGCGATTGGCGGCGGTAGTTGGATGGTGGTGCAGGGCAGtt                                       >  1:1855914/1‑150 (MQ=255)
                                                                                                          ttGCCGCGGATGCCGAAGATACCGGCAAAAAAAACATGCGGGTGGCGCGTATTGATGCTCGTTTCGACCCGACCATCTATATCGCGATTGGTATGGCGAACTTGCTGGCGATTGGCGGCGGTAGTTGGATGGTGGTGCAGGGCAGTTTaa                                    >  1:1303763/1‑150 (MQ=255)
                                                                                                            gCCGCGGATGCCGAAGATACCGGCAAAAAAAACATGCTGGTGGCGCGTATTGATGCTCGTTTCGACCCGACCATCTATATCGCGATTGGTATGGCGAACTTGCTGGCGATTGGCGGCGGTAGTTGGATGGTGGTGCAGGGCAGTTTAACg                                  >  2:1304265/1‑150 (MQ=255)
                                                                                                            gCCGCGGATGCCGAAGATACCGGCAAAAAAAACATGCTGGTGGCGCGTATTGATGCTCGTTTCGACCCGACCATCTATATCGCGATTGGTATGGCGAACTTGCTGGCGATTGGCGGCGGTAGTTGGATGGTGGTGCAGGGCAGTTTAACg                                  >  2:1310869/1‑150 (MQ=255)
                                                                                                             ccGCGGATGCCGAAGATACCGGCAAAAAAAACATGCTGGTGGCGCGTATTGATGCTCGTTTCGACCCGACCATCTATATCGCGATTGGTATGGCGAACTTGCTGGCGATTGGCGGCGGTAGTTGGATGGTGGTGCAGGGCAGTTTAACGc                                 <  1:1304265/150‑1 (MQ=255)
                                                                                                             ccGCGGATGCCGAAGATACCGGCAAAAAAAACATGCTGGTGGCGCGTATTGATGCTCGTTTCGACCCGACCATCTATATCGCGATTGGTATGGCGAACTTGCTGGCGATTGGCGGCGGTAGTTGGATGGTGGTGCAGGGCAGTTTAACGc                                 <  1:1310869/150‑1 (MQ=255)
                                                                                                                 ggATGCCGAAGATACCGGCAAAAAAAACATGCGGGTGGCGCGTATTGATGCTCGTTTCGACCCGACCATCTATATCGCGATTGGTATGGCGAACTTGCTGGCGATTGGCGGCGGTAGTTGGATGGTGGTGCAGGGCAGTTTAACGCTggg                             <  2:1148231/150‑1 (MQ=255)
                                                                                                                   aTGCCGAAGATACCGGCAAAAAAAACATGCGGGTGGCGCGTATTGATGCTCGTTTCGACCCGACCATCTATATCGCGATTGGTATGGCGAACTTGCTGGCGATTGGCGGCGGTAGTTGGATGGTGGTGCAGGGCAGTTTAACGCTGGGcc                           >  2:1818824/1‑150 (MQ=255)
                                                                                                                    tGCCGAAGATACCGGCAAAAAAAACATGCGGGTGGCGCGTATTGATGCTCGTTTCGACCCGACCATCTATATCGCGATTGGTATGGCGAATTTGCTGGCGATTGGCGGCGGTAGTTGGATGGTGGTGCAGGGCAGTTTAACGCTGGGCCa                          >  2:1848103/1‑150 (MQ=255)
                                                                                                                     gCCGAAGATACCGGCAAAAAAAACATGCGGGTGGCGCGTATTGATGCTCGTTTCGACCCGACCATCTATATCGCGATTGGTATGGCGAACTTGCTGGCGATTGGCGGCGGTAGTTGGATGGTGGTGCAGGGCAGTTTAACGCTGGGCCAg                         <  1:882365/150‑1 (MQ=255)
                                                                                                                      ccGAAGATACCGGCAAAAAAAACATGCGGGTGGCGCGTATTGATGCTCGTTTCGACCCGACCATCTATATCGCGATTGGTATGGCGAACTTGCTGGCGATTGGCGGCGGTAGTTGGATGGTGGTGCAGGGCAGTTTAACGCTGGGCCAGc                        >  2:882138/1‑150 (MQ=255)
                                                                                                                      ccGAAGATACCGGCAAAAAAAACATGCGGGGGGCGCGTATTGATGCTCGTTTCGACCCGACCATCTATATCGCGATTGGTATGGCGAACTTGCTGGCGATTGGCGGCGGTAGTTGGAGGGTGGTGCAGGGCAGTTTAACGCTGGGCCAGc                        <  1:1602776/150‑1 (MQ=255)
                                                                                                                      ccGAACATACCGGCAAAAAAAACATGCGGGTGGAGCGTATTGATGATCGTTTCGACCCGACCAACTATATCGCGATTGGTATGGCGAACTTGCTGGCGATTGGCGGCGGTAGATGGATGGTGGTGCAGGGCAGTTTAACGCTGGGCCAGc                        >  2:1658849/1‑150 (MQ=255)
                                                                                                                       cGAAGATACCGGCAAAAAAAACATGCGGGTGGCGCGTATTGATGCTCGTTTCGACCCGACCATCTATATCGCGATTGGTATGGCGAACTTGCTGGCGATTGGCGGCGGTAGTTGGATGGTGGTGCAGGGCAGTTTAACGCTGGGCCAGCt                       >  2:509730/1‑150 (MQ=255)
                                                                                                                               ccGGCAAAAAAAACATGCGGGTGGCGCGTATTGATGCTCGTTTCGACCCGACCATCTATATCGCGATTGGTATGGCGAACTTGCTGGCGATTGGCGGCGGTAGTTGGATGGTGGTGCAGGGCAGTTTAACGCTGGGCCAGCTCACCAGtt               >  2:985170/1‑150 (MQ=255)
                                                                                                                               ccGGCAAAAAAAACATGCGGGTGGCGCGTATTGATGCTCGTTTCGACCCGACCATCTATATCGCGATTGGTATGGCGAACTTGCTGGCGATTGGCGGCGGTAGTTGGATGGTGGTGCAGGGCAGTTTAACGCTGGGCCAGCTCACCAGtt               >  1:822304/1‑150 (MQ=255)
                                                                                                                                cGGCAAAAAAAACATGCGGGTGGCGCGTATTGATGCTCGTTTCGACCCGACCATCTATATCGCGATTGGTATGGCGAACTTGCTGGCGATTGGCGGCGGTAGTTGGATGGTGGTGCAGGGCAGTTTAACGCTGGGCCAGCTCACCAGttt              >  2:601765/1‑150 (MQ=255)
                                                                                                                                 ggCAAAAAAAACATGCGGGTGGCGCGTATTGATGCTCGTTTCGACCCGACCATCTATATCGCGATTGGTATGGCGAACTTGCTGGCGATTGGCGGCGGTAGTTGGATGGTGGTGCAGGGCAGTTTAACGCTGGGCCAGCTCACCAGtttt             >  2:381740/1‑150 (MQ=255)
                                                                                                                                     aaaaaaaCATGCGGGTGGCGCGTATTGATGCTCGTTTCGACCCGACCATCTATATCGCGATTGGTATGGCGAACTTGCTGGCGATTGGCGGCGGTAGTTGGATGGTGGTGCAGGGCAGTTTAACGCTGGGCCAGCTCACCAGTTTTatga         <  1:158860/150‑1 (MQ=255)
                                                                                                                                     aaaaaaaCATGCGGGGGGCGCGTCTTGATGCTCGTTTCGACCCGACCATCTATATCGCGATTGGTATGGCGAACTTGCTGGCGATTGGCGGCGGTAGTTGGATGGTGGTGCAGGGCAGTTTAACGCTGGGCCAGCTCACCAGTTTTatga         <  1:1938748/150‑1 (MQ=255)
                                                                                                                                            cATGCGGGTGGCGCGTATTGATGCTCGTTTCGACCCGACCATCTATATCGCGATTGGTATGGCGAACTTGCTGGCGATTGGCGGCGGTAGTTGGATGGTGGTGCAGGGCAGTTTAACGCTGGGCCAGCTCACCAGTTTTATGATGTAttt  <  2:944222/150‑1 (MQ=255)
                                                                                                                                            cATGCGGGTGGCGCGTATTGATGCTCGTTTCGACCCGACCATCTATATCGCGATTGGTATGGCGAACTTGCTGGCGATTGGCGGCGGTAGTTGGATGGTGGTGCAGGGCAGTTTAACGCTGGGCCAGCTCACCAGTTTTATGATGTAttt  <  2:1303763/150‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                
CTGGCACAGGCGGCGTTTTCCAGTCTTAATGACCGCACCCAGGAAAGCCTCACCAGTATCCGCATGATCAAAGCCTTTGGTCTGGAAGATCGCCAGTCGGCGTTATTTGCCGCGGATGCCGAAGATACCGGCAAAAAAAACATGCGGGTGGCGCGTATTGATGCTCGTTTCGACCCGACCATCTATATCGCGATTGGTATGGCGAACTTGCTGGCGATTGGCGGCGGTAGTTGGATGGTGGTGCAGGGCAGTTTAACGCTGGGCCAGCTCACCAGTTTTATGATGTATTT  >  NC_000913/469432‑469721

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: