Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,300,041 A→G 5.7% D330D (GAT→GAC mdtP ← outer membrane factor of efflux pump

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,300,0410AG 5.7% 224.9 / 3.7 71D330D (GAT→GACmdtPouter membrane factor of efflux pump
Reads supporting (aligned to +/- strand):  ref base A (31/35);  new base G (2/2);  total (34/37)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.52e-01

GGCGCGCGTGCCTTCGAGATTGGCATTCAACCGTCCACCGTCAAACAGCGGCAATTTCAGACCCGGGATGAAGTTGAACTGGCGACTGGTTTTTTTGAATAAGGTATGCAGATGGATGGAGTCCAGACCGAAAAACGCTTTGATATCAAAGCTCGGATAGAACAACGCCCGCGCGGAATCCACCTGATCTAATGACGCCTGAACATACCAGCGCATGGCTTGCAGATCCGGGCGTCTGGCGAGCAACTCATAAGAGAGTGTCGCCGGAATGCCGGTCTGGACTTGCGGTAAT  >  NC_000913/4299897‑4300188
                                                                                                                                                |                                                                                                                                                   
ggcgCGCGTGCCTTCGAGATTGGCATTCAACCGTCCACCGTCAAACAGCGGCAATTTCAGACCCGGGATGAAGTTGAACTGGCGACTGGTTTTTTTGAATAAGGTATGCAGATGGATGGAGTCCAGACCGAAAAACGCTTTGATATCaa                                                                                                                                                 >  2:447049/1‑149 (MQ=255)
   gcgcgTGCCTTCGAGATTGGCATTCAACCGTCCACCGTCAAACAGCGGCAATTTCAGACCCGGGATGAAGTTGAACTGGCGACTGGTTTTTTTGAATAAGGTATGCAGATGGATGGAGTCCAGACCGAAAAACGCTTTGATATCAAAGc                                                                                                                                              >  4:157066/1‑149 (MQ=255)
       gTGCCTTCGAGATTGGCATTCAACCGTCCACCGTCAAACAGCGGCAATTTCAGACCCGGGATGAAGTTGAACTGGCGACTGGTTTTTTTGAATAAGGTATGCAGATGGATGGAGTCCAGACCGAAAAACGCTTTGATATCAAAGCTCgg                                                                                                                                          <  3:705466/149‑1 (MQ=255)
          ccTTCGAGATTGGCATTCAACCGTCCACCGTCAAACAGCGGCAATTTCAGACCCGGGATGAAGTTGAACTGGCGACTGGTTTTTTTGAATAAGGTATGCAGATGGATGGAGTCCAGACCGAAAAACGCTTTGATATCAAAGCTCGGATa                                                                                                                                       >  4:420056/1‑149 (MQ=255)
            ttCGAGATTGGCATTCAACCGTCCACCGTCAAACAGCGGCAATTTCAGACCCGGGATGAAGTTGAACTGGCGACTGGTTTTTTTGAATAAGGTATGCAGATGGATGGAGTCCAGACCGAAAAACGCTTTGATATCAAAGCTCg                                                                                                                                           >  2:645231/1‑143 (MQ=255)
            ttCGAGATTGGCATTCAACCGTCCACCGTCAAACAGCGGCAATTTCAGACCCGGGATGAAGTTGAACTGGCGACTGGTTTTTTTGAATAAGGTATGCAGATGGATGGAGTCCAGACCGAAAAACGCTTTGATATCAAAGCTCg                                                                                                                                           <  1:645231/143‑1 (MQ=255)
               gagaTTGGCATTCAACCGTCCACCGTCAAACAGCGGCAATTTCAGACCCGGGATGAAGTTGAACTGGCGACTGGTTTTTTTGAATAAGGTATGCAGATGGATGGAGTCCAGACCGAAAAACGCTTTGATATCAAAGCTCGGATAGaaca                                                                                                                                  >  1:786286/1‑149 (MQ=255)
                agaTTGGCATTCAACCGTCCACCGTCAAACAGCGGCAATTTCAGACCCGGGATGAAGTTGAACTGGCGACTGGTTTTTTTGAATAAGGTATGCAGATGGATGGAGTCCAGACCGAAAAACGCTTTGATATCAAAGCTCGGATAGaacaa                                                                                                                                 >  4:137615/1‑149 (MQ=255)
                agaTTGGCATTCAACCGTCCACCGTCAAACAGCGGCAATTTCAGACCCGGGATGAAGTTGAACTGGCGACTGGTTTTTTTGAATAAGGTATGCAGATGGATGGAGTCCAGACCGAAAAACGCTTTGATATCAAAGCTCGGATAGaacaa                                                                                                                                 >  4:731294/1‑149 (MQ=255)
                agaTTGGCATTCAACCGTCCACCGTCAAACAGCGGCAATTTCAGACCCGGGATGAAGTTGAACTGGCGACTGGTTTTTTTGAATAAGGTATGCAGATGGATGGAGTCCAGACCGAAAAACGCTTTGATATCAAAGCTCGGATAGaacaa                                                                                                                                 <  3:97009/149‑1 (MQ=255)
                 gaTTGGCATTCAACCGTCCACCGTCAAACAGCGGCAATTTCAGACCCGGGATGAAGTTGAACTGGCGACTGGTTTTTTTGAATAAGGTATGCAGATGGATGGAGTCCAGACCGAAAAACGCTTTGATATCAAAGCTCGGATAGaacaac                                                                                                                                <  2:723412/149‑1 (MQ=255)
                     ggCATTCAACCGTCCACCGTCAAACAGCGGCAATTTCAGACCCGGGATGAAGTTGAACTGGCGACTGGTTTTTTTGAATAAGGTATGCAGATGGATGGAGTCCAGACCGAAAAACGCTTTGATATCAAAGCTCGGATAGAACAACGccc                                                                                                                            >  3:106348/1‑149 (MQ=255)
                        aTTCAACCGTCCACCGTCAAACAGCGGCAATTTCAGACCCGGGATGAAGTTGAACTGGCGACTGGTTTTTTTGAATAAGGTATGCAGATGGATGGAGTCCAGACCGAAAAACGCTTTGATATCAAAGCTCGGATAGAACAACGCCcgcg                                                                                                                         <  3:82171/149‑1 (MQ=255)
                        aTTCAACCGTCCACCGTCAAACAGCGGCAATTTCAGACCCGGGAAGAAGTTGAACTGGCGACTGGTTTTCTTGAATAAGGTATGCAGATGGATGGAGTCCAGACCGAAAAACGCTTTGATATCAAAGCTCGGATAGAACAACGCCcgcg                                                                                                                         >  2:154574/1‑149 (MQ=255)
                         ttCAACCGTCCACCGTCAAACAGCGGCAATTTCAGACCCGGGATGAAGTTGAACTGGCGACTGGTTTTTTTGAATAAGGTATGCAGATGGATGGAGTCCAGACCGAAAAACGCTTTGATATCAAAGCTCGGATAGAACAACGCCcgcgc                                                                                                                        >  4:540229/1‑149 (MQ=255)
                         ttCAACCGTCCACCGTCAAACAGCGGCAATTTCAGACCCGGGATGAAGTTGAACTGGCGACTGGTTTTTTTGAATAAGGTATGCAGATGGATGGAGTCCAGACCGAAAAACGCTTTGATATCAAAGCTCGGATAGAACAACGCCcgcgc                                                                                                                        >  4:765921/1‑149 (MQ=255)
                           cAACCGTCCACCGTCAAACAGCGGCAATTTCAGACCCGGGATGAAGTTGAACTGGCGACTGGTTTTTTTGAATAAGGTATGCAGATGGATGGAGTCCAGACCGAAAAACGCTTTGATATCAAAGCTCGGATAGAACAACGCCCGCGCgg                                                                                                                      <  4:106348/149‑1 (MQ=255)
                            aaCCGTCCACCGTCAAACAGCGGCAATTTCAGACCCGGGATGAAGTTGAACTGGCGACTGGTTTTTTTGAATAAGGTATGCAGATGGATGGAGTCCAGACCGAAAAACGCTTTGATATCAAAGCTCGGATAGAACAACGCCCGCGCGGa                                                                                                                     <  4:480823/149‑1 (MQ=255)
                                 tCCACCGTCAAACAGCGGCAATTTCAGACCCGGGATGAAGTTGAACTGGCGACTGGTTTTTTTGAATAAGGTATGCAGATGGATGGAGTCCAGACCGAAAAACGCTTTGATATCAAAGCTCGGATAGAACAACGCCCGCGCGGAATCCa                                                                                                                <  3:420056/149‑1 (MQ=255)
                                  ccACCGTCAAACAGCGGCAATTTCAGACCCGGGATGAAGTTGAACTGGCGACTGGTTTTTTTGAATAAGGTATGCAGATGGATGGAGTCCAGACCGAAAAACGCTTTGATATCAAAGCTCGGATAGAACAACGCCCGCGCGGAATCCAc                                                                                                               >  3:479642/1‑149 (MQ=255)
                                    aCCGTCAAACAGCGGCAATTTCAGACCCGGGATGAAGTTGAACTGGCGACTGGTTTTTTTGAATAAGGTATGCAGATGGATGGAGTCCAGACCGAAAAACGCTTTGATATCAAAGCTCGGATAGAACAACGCCCGCGCGGAATCCACCt                                                                                                             >  4:668539/1‑149 (MQ=255)
                                            aCAGCGGCAATTTCAGACCCGGGATGAAGTTGAACTGGCGACTGGTTTTTTTGAATAAGGTATGCAGATGGATGGAGTCCAGACCGAAAAACGCTTTGATATCAAAGCTCGGATAGAACAACGCCCGCGCGGAATCCACCTGATCTAAt                                                                                                     <  4:695861/149‑1 (MQ=255)
                                               gCGGCAATTTGAGACCCGGGATGAAGTTGAACTGGCGACTGGTTTTTTTGAATAAGGTATGCAGATGGATGGAGTCCAGACCGAAAAACGCTTTGATATCAAAGCTCGGATAGAACAACGCCCGCGCGGAATCCACCTGAt                                                                                                          <  3:432332/141‑1 (MQ=255)
                                               gCGGCAATTTGAGACCCGGGATGAAGTTGAACTGGCGACTGGTTTTTTTGAATAAGGTATGCAGATGGATGGAGTCCAGACCGAAAAACGCTTTGATATCAAAGCTCGGATAGAACAACGCCCGCGCGGAATCCACCTGAt                                                                                                          >  4:432332/1‑141 (MQ=255)
                                                   cAATTTCAGACCCGGGATGAAGTTGAACTGGCGACTGGTTTTTTTGAATAAGGTATGCAGATGGATGGAGTCCAGACCGAAAAACGCTTTGATATCAAAGCTCGGATAGAACAACGCCCGCGCGGAATCCACCTGATCTAATGAc                                                                                                  >  2:666093/1‑145 (MQ=255)
                                                   cAATTTCAGACCCGGGATGAAGTTGAACTGGCGACTGGTTTTTTTGAATAAGGTATGCAGATGGATGGAGTCCAGACCGAAAAACGCTTTGATATCAAAGCTCGGATAGAACAACGCCCGCGCGGAATCCACCTGATCTAATGAc                                                                                                  <  1:666093/145‑1 (MQ=255)
                                                         cAGACCCGGGATGAAGTTGAACTGGCGACTGGTTTTTTTGAATAAGGTATGCAGATGGATGGAGTCCAGACCGAAAAACGCTTTGATATCAAAGCTCGGATAGAACAACGCCCGCGCGGAATCCACCTGATCTAATGACGCCTGAACAt                                                                                        <  1:769527/149‑1 (MQ=255)
                                                           gACCCGGGATGAAGTTGAACTGGCGACTGGTTTTTTTGAATAAGGTATGCAGATGGATGGAGTCCAGACCGAAAAACGCTTTGATATCAAAGCTCGGATAGAACAACGCCCGCGCGGAATCCACCTGATCTAATGACGCCTGAACATAc                                                                                      <  2:626295/149‑1 (MQ=255)
                                                           gACCCGGGATGAAGTTGAACTGGCGACTGGTTTTTTTGAATAAGGTATGCAGATGGATGGAGTCCAGACCGAAAAACGCTTTGATATCAAAGCTCGGATAGAACAACGCCCGCGCGGAATCCACCTGATCTAATGACGCCTGAACATAc                                                                                      <  3:137615/149‑1 (MQ=255)
                                                             cccGGGATGAAGTTGAACTGGCGACTGGTTTTTTTGAATAAGGTATGCAGATGGATGGAGTCCAGACCGAAAAACGCTTTGATATCAAAGCTCGGATAGAACAACGCCCGCGCGGAATCCACCTGATCTAATGACGCCTGAACATACCa                                                                                    <  2:233128/149‑1 (MQ=255)
                                                              ccGGGATGAAGTTGAACTGGCGACTGGTTTTTTTGAATAAGGTATGCAGATGGATGGAGTCCAGACCGAAAAACGCTTTGATGTCAAAGCTCGGATAGAACAACGCCCGCGCGGAATCCACCTGATCTAATGACGCCTGAACATACCAg                                                                                   >  3:357709/1‑149 (MQ=255)
                                                              ccGGGATGAAGTTGAACTGGCGACTGGTTTTTTTGAATAAGGTATGCAGATGGATGGAGTCCAGACCGAAAAACGCTTTGATGTCAAAGCTCGGATAGAACAACGCCCGCGCGGAATCCACCTGATCTAATGACGCCTGAACATACCAg                                                                                   >  3:118918/1‑149 (MQ=255)
                                                              ccGGGATGAAGTTGAACTGGCGACTGGTTTTTTTGAATAAGGTATGCAGATGGATGGAGTCCAGACCGAAAAACGCTTTGATATCAAAGCTCGGATAGAACAACGCCCGCGCGGAATCCACCTGATCTAATGACGCCTGAACATACCAg                                                                                   >  3:377844/1‑149 (MQ=255)
                                                              ccGGGATGAAGTTGAACTGGCGACTGGTTTTTTTGAATAAGGTATGCAGATGGATGGAGTCCAGACCGAAAAACGCTTTGATATCAAAGCTCGGATAGAACAACGCCCGCGCGGAATCCACCTGATCTAATGACGCCTGAACATACCAg                                                                                   >  1:14970/1‑149 (MQ=255)
                                                              ccGGGATGAAGTTGAACTGGCGACTGGTTTTTTTGAATAAGGTATGCAGATGGATGGAGTCCAGACCGAAAAACGCTTTGATATCAAAGCTCGGATAGAACAACGCCCGCGCGGAATCCACCTGATCTAATGACGCCTGAACATACCAg                                                                                   >  1:31686/1‑149 (MQ=255)
                                                                 ggATGAAGTTGAACTGGCGACTGGTTTTTTTGAATAAGGTATGCAGATGGATGGAGTCCAGACCGAAAAACGCTTTGATATCAAAGCTCGGATAGAACAACGCCCGCGCGGAATCCACCTGATCTAATGACGCCTGAACATACCAgcgc                                                                                <  2:786286/149‑1 (MQ=255)
                                                                    tGAAGTTGAACTGGCGACTGGTTTTTTTGAATAAGGTATGCAGATGGATGGAGTCCAGACCGAAAAACGCTTTGATATCAAAGCTCGGATAGAACAACGCCCGCGCGGAATCCACCTGATCTAATGACGCCTGAACATACCAgc                                                                                  <  1:738513/144‑1 (MQ=255)
                                                                    tGAAGTTGAACTGGCGACTGGTTTTTTTGAATAAGGTATGCAGATGGATGGAGTCCAGACCGAAAAACGCTTTGATATCAAAGCTCGGATAGAACAACGCCCGCGCGGAATCCACCTGATCTAATGACGCCTGAACATACCAgc                                                                                  >  2:738513/1‑144 (MQ=255)
                                                                       aGTTGAACTGGCGACTGGTTTTTTTGAATAAGGTATGCAGATGGATGGCGTCCAGACCGAAAAACGCTTTGATCTCAAAGCTCGGATAGAGCAACACCCGCGCGGAATCCACCTGATCTAATGACGCCTGAACATACCAGCGCATGGCt                                                                          >  2:249707/1‑149 (MQ=255)
                                                                        gTTGAACTGGCGACTGGTTTTTTTGAATAAGGTATGCAGATGGATGGAGTCCAGACCGAAAAACGCTTTGATATCAAAGCTCGGATAGAACAACGCCCGCGCGGAATCCACCTGATCTAATGACGCCTGAACATACCAGCGCATGGcat                                                                         >  2:443313/1‑147 (MQ=255)
                                                                         ttGAACTGGCGACTGGTTTTTTTGAATAAGGTATGCAGATGGATGGAGTCCAGACCGAAAAACGCTTTGATATCAAAGCTCGGATAGAACAACGCCCGCGCGGAATCCACCTGATCTAATGACGCCTGAACATACCAGCGCATGGCTTg                                                                        >  4:136101/1‑149 (MQ=255)
                                                                             aCTGGCGACTGGTTTTTTTGAATAAGGTATGCAGATGGATGGAGTCCAGACCGAAAAACGCTTTGATATCAAAGCTCGGATAGAACAACGCCCGCGCGGAATCCACCTGATCTAATGACGCCTGAACATACCAGCGCATGGCTTGCAGa                                                                    <  3:668539/149‑1 (MQ=255)
                                                                                  cGACTGGTTTTTTTGAATAAGGTATGCAGATGGATGGAGTCCAGACCGAAAAACGCTTTGATATCAAAGCTCGGATAGAACAACGCCCGCGCGGAATCCACCTGATCTAATGACGCCTGAACATACCAGCGCATGGCTTGCAGATCCgg                                                               <  1:447049/149‑1 (MQ=255)
                                                                                   gACTGGTTTTTTTGAATAAGGTATGCAGATGGATGGAGTCCAGACCGAAAAACGCTTTGATATCAAAGCTCGGATAGAACAACGCCCGCGCGGAATCCACCTGATCTAATGACGCCTGAACATACCAGCGCATGGCTTGCAGATCCggg                                                              >  2:409448/1‑149 (MQ=255)
                                                                                             tttGAATAAGGTATGCAGATGGATGGAGTCCAGACCGAAAAACGCTTTGATATCAAAGCTCGGATAGAACAACGCCCGCGCGGAATCCACCTGATCTAATGACGCCTGAACATACCAGCGCATGGCTTGCAGATCCGGGCGTCTGGCGa                                                    <  1:664222/149‑1 (MQ=255)
                                                                                             tttGAATAAGGTATGCAGATGGATGGAGTCCAGACAGAAAAACGCTTTGATATCAAAGCTCGGATAGAACAACGCCCGCGCGGAATCCACCTGATCTAATGACGCCTGAACATACCAGCGCATGGCTTGCAGATCCGGGCGTCTGGCGa                                                    <  2:767539/149‑1 (MQ=255)
                                                                                                 aaTAAGGTATGCAGATGGATGGAGTCCAGACCGAAAAACGCTTTGATATCAAAGCTCGGATAGAACAACGCCCGCGCGGAATCCACCTGATCTAATGACGCCTGAACATACCAGCGCATGGCTTGCAGATCCGGGCGTCTGGCGAGCaa                                                <  4:161459/149‑1 (MQ=255)
                                                                                                 aaTAAGGTATGCAGATGGATGGAGTCCAGACCGAAAAACGCTTTGATATCAAAGCTCGGATAGAACAACGCCCGCGCGGAATCCACCTGATCTAATGACGCCTGAACATACCAGCGCATGGCTTGCAGATCCGGGCGTCTGGCGAGCaa                                                <  1:314548/149‑1 (MQ=255)
                                                                                                    aaGGTATGCAGATGGATGGAGTCCAGACCGAAAAACGCTTTGATATCAAAGCTCGGATAGAACAACGCCCGCGCGGAATCCACCTGATCTAATGACGCCTGAACATACCAGCGCATGGCTTGCAGATCCGGGCGTCTGGCGAGCAACTc                                             <  2:72450/149‑1 (MQ=255)
                                                                                                        tATGCAGATGGATGGAGTCCAGACCGAAAAACGCTTTGATATCAAAGCTCGGATAGAACAACGCCCGCGCGGAATCCACCTGATCTAATGACGCCTGAACATACCAGCGCATGGCTTGCAGATCCGGGCGTCTGGCGAGCAACTCATaa                                         <  2:31686/149‑1 (MQ=255)
                                                                                                        tATGCAGATGGATGGAGTCCAGACCGAAAAACGCTTTGATATCAAAGCTCGGATAGAACAACGCCCGCGCGGAATCCACCTGATCTAATGACGCCTGAACATACCAGCGCATGGCTTGCAGATCCGGGCGTCTGGCGAGCAACTCATaa                                         <  2:14970/149‑1 (MQ=255)
                                                                                                        tATGCAGATGGATGGAGTCCAGACCGAAAAACGCTTTGATATCAAAGCTCGGATAGAACAACGCCCGCGCGGAATCCACCTGATCTAATGACGCCTGAACATACCAGCGCATGGCTTGCAGATCCGGGCGTCTGGCGAGCAACTCATaa                                         >  1:313405/1‑149 (MQ=255)
                                                                                                         aTGCAGATGGATGGAGTCCAGACCGAAAAACGCTTTGATATCAAAGCTCGGATAGAACAACGCCCGCGCGGAATCCACCTGATCTAATGACGCCTGAACATACCAGCGCATGGCTTGCAGATCCGGGCGTCTGGCGAGCAACTCATaa                                         >  3:553225/1‑148 (MQ=255)
                                                                                                         aTGCAGATGGATGGAGTCCAGACCGAAAAACGCTTTGATATCAAAGCTCGGATAGAACAACGCCCGCGCGGAATCCACCTGATCTAATGACGCCTGAACATACCAGCGCATGGCTTGCAGATCCGGGCGTCTGGCGAGCAACTCATaa                                         <  4:553225/148‑1 (MQ=255)
                                                                                                             aGATGGATGGAGTCCAGACCGAAAAACGCTTTGATATCAAAGCTCGGATAGAACAACGCCCGCGCGGAATCCACCTGATCTAATGACGCCTGAACATACCAGCGCATg                                                                             >  2:497212/1‑108 (MQ=255)
                                                                                                             aGATGGATGGAGTCCAGACCGAAAAACGCTTTGATATCAAAGCTCGGATAGAACAACGCCCGCGCGGAATCCACCTGATCTAATGACGCCTGAACATACCAGCGCATg                                                                             <  1:497212/108‑1 (MQ=255)
                                                                                                             aGATGGATGGAGTCCAGACCGAAAAACGCTTTGATATCAAAGCTCGGATAGAACAACGCCCGCGCGGAATCCACCTGATCTAATGACGCCTGAACATACCAGCGCATGGCTTGCAGATCCgg                                                               >  2:746992/1‑122 (MQ=255)
                                                                                                             aGATGGATGGAGTCCAGACCGAAAAACGCTTTGATATCAAAGCTCGGATAGAACAACGCCCGCGCGGAATCCACCTGATCTAATGACGCCTGAACATACCAGCGCATGGCTTGCAGATCCgg                                                               <  1:746992/122‑1 (MQ=255)
                                                                                                              gatggatgGAGTCCAGACCGAAAAACGCTTTGATATCAAAGCTCGGATAGAACAACGCCCGCGCGGAATCCACCTGATCTAATGACGCCTGAACATACCAGCGCATGGCTTGCAGATCCGGGCGTCt                                                         >  2:415440/1‑127 (MQ=255)
                                                                                                              gatggatgGAGTCCAGACCGAAAAACGCTTTGATATCAAAGCTCGGATAGAACAACGCCCGCGCGGAATCCACCTGATCTAATGACGCCTGAACATACCAGCGCATGGCTTGCAGATCCGGGCGTCt                                                         <  1:415440/127‑1 (MQ=255)
                                                                                                                tggatggaGTCCAGACCGAAAAACGCTTTGATATCAAAGCTCGGATAGAACAACGCCCGCGCGGAATCCACCTGATCTAATGACGCCTGAACATACCAGCGCATGGCTTGCAGAGCCGGGCGTCTGGCGAGCAACTCATAAGAGAgtgt                                 <  1:621673/149‑1 (MQ=255)
                                                                                                                   atggaGTCCAGACCGAAAAACGCTTTGATGTCAAAGCTCGGATAGAACAACGCCCGCGCGGAATCCACCTGATCTAATGACGCCTGAACATACCAGCGCATGGCTTGCAGATCCGGGCGTCTGGCGAGCAACTCATAAGAGAGTGTCGc                              <  4:118918/149‑1 (MQ=255)
                                                                                                                   atggaGTCCAGACCGAAAAACGCTTTGATGTCAAAGCTCGGATAGAACAACGCCCGCGCGGAATCCACCTGATCTAATGACGCCTGAACATACCAGCGCATGGCTTGCAGATCCGGGCGTCTGGCGAGCAACTCATAAGAGAGTGTCGc                              <  4:357709/149‑1 (MQ=255)
                                                                                                                             gACCGAAAAACGCTTTGATATCAAAGCTCGGATAGAACAACGGCCGCGCGGAATCCACCTGATCTAATGACGCCTGAACATACCAGCGCATGGCTTGCAGATCCGGGCGTCTGGCGAGCAACTCATAAGAGAGTGTCGCCGGAATGCCg                    <  1:154574/149‑1 (MQ=255)
                                                                                                                             gACCGAAAAACGCTTTGATATCAAAGCTCGGATAGAACAACGCCCGCGCGGAATCCACCTGATCTAATGACGCCTGAACATACCAGCGCATGGCTTGCAGATCCGGGCGTCTGGCGAGCAACTCATAAGAGAGTGTCGCCGGAATGCCg                    >  2:17502/1‑149 (MQ=255)
                                                                                                                             gACCGAAAAACGCTTTGATATCAAAGCTCGGATAGAACAACGCCCGCGCGGAATCCACCTGATCTAATGACGCCTGAACATACCAGCGCATGGCTTGCAGATCCGGGCGTCTGGCGAGCAACTCATAAGAGAGTGTCGCCGGAATGCCg                    >  1:378469/1‑149 (MQ=255)
                                                                                                                                 gAAAAACGCTTTGATATCAAAGCTCGGATAGAACAACGCCCGCGCGGAATCCACCTGATCTAATGACGCCTGAACATACCAGCGCATGGCTTGCAGATCCGGGCGTCTGGCGAGCAACTCATAAGAGAGTGTCGCCGGAATGCCGGTCt                <  3:765921/149‑1 (MQ=255)
                                                                                                                                  aaaaaCGCTTTGATATCAAAGCTCGGATAGAACAACGCCCGCGCGGAATCCACCTGATCTAATGACGCCTGAACATACCAGCGCATGGCTTGCAGATCCGGGCGTCTGGCGAGCAACTCATAAGAGAGTGTCGCCGGAATGCCGGTCTg               <  1:316593/149‑1 (MQ=255)
                                                                                                                                        gCTTTGATATCAAAGCTCGGATAGAACAACGCCCGCGCGGAATCCACCTGATCTAATGACGCCTGAACATACCAGCGCATGGCTTGCAGATCCGGGCGTCTGGCGAGCAACTCATAAGAGAGTGTCGCCGGAATGCCGGTCTGGACTTg         <  3:731294/149‑1 (MQ=255)
                                                                                                                                          tttGATATCAAAGCTCGGATAGAACAACGCCCGCGCGGAATCCACCTGATCTAATGACGCCTGAACATACCAGCGCATGGCTTGCAGATCCGGGCGTCTGGCGAGCAACTCATAAGAGAGTGTCGCCGGAATGCCGGTCTGGACTTGCg       >  4:289897/1‑149 (MQ=255)
                                                                                                                                          tttGATATCAAAGCTCGGATAGAACAACGCCCGCGCGGAATCCACCTGATCTAATGACGCCTGAACATACCAGCGCATGGCTTGCAGATCCGGGCGTCTGGCGAGCAACTCATAAGAGAGTGTCGCCGGAATGCCGGTCTGGACTTGCg       >  2:683033/1‑149 (MQ=255)
                                                                                                                                               tatCAAAGCTCGGATAGAACAACGCCCGCGCGGAATCCACCTGATCTAATGACGCCTGAACATACCAGCGCATGGCTTGCAGATCCGGGCGTCTGGCGAGCAACTCATAAGAGAGTGTCGCCGGAATGCCGGTCTGGACTTGCGGTAAt  <  4:479642/149‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                                   
GGCGCGCGTGCCTTCGAGATTGGCATTCAACCGTCCACCGTCAAACAGCGGCAATTTCAGACCCGGGATGAAGTTGAACTGGCGACTGGTTTTTTTGAATAAGGTATGCAGATGGATGGAGTCCAGACCGAAAAACGCTTTGATATCAAAGCTCGGATAGAACAACGCCCGCGCGGAATCCACCTGATCTAATGACGCCTGAACATACCAGCGCATGGCTTGCAGATCCGGGCGTCTGGCGAGCAACTCATAAGAGAGTGTCGCCGGAATGCCGGTCTGGACTTGCGGTAAT  >  NC_000913/4299897‑4300188

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