Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,660,594 A→T 5.4% E13D (GAA→GAT ynfF → S‑ and N‑oxide reductase, A subunit, periplasmic

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,660,5940AT 5.4% 244.7 / 2.2 84E13D (GAA→GATynfFS‑ and N‑oxide reductase, A subunit, periplasmic
Reads supporting (aligned to +/- strand):  ref base A (36/43);  new base T (2/3);  total (38/46)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.55e-01

GAAGGGAAACCCGCAGCACACTAATCTGGTCGAGATCGAAAAAATATAACCCACGACAACCATAAATTCTGGCATGACATTTTGTTTGAAAAGCAATAAGTGAGTAATGATGAAAATCCATACCACAGAGGCGCTCATGAAGGCTGAAATCAGCCGTAGAAGTCTGATGAAAACCTCCGCACTTGGCAGTCTTGCGCTGGCAAGCAGTGCTTTCACTCTGCCATTTTCCCAAATGGTCCGGGCGGCAGAGGCTCCGGTAGAAGAGAAAGCGGTCTGGAGTTCCTGCACCGTTAACTG  >  NC_000913/1660447‑1660743
                                                                                                                                                   |                                                                                                                                                     
gAAGGGAAACCCGCAGCACACTAATCTGGTCGAGATCGAAAAAATATAACCCACGACAACCATAAATTCTGGCATGACATTTTGTTTGAAAAGCAATAAGTGAGTAATGATGAAAATCCATACCACAGAGGCGCTCATGAAGGCTGAAAt                                                                                                                                                     >  2:118616/1‑150 (MQ=255)
gAAGGGAAACCCGCAGCACACTAATCTGGTCGAGATCGAAAAAATATAACCCACGACAACCATAAATTCTGGCATGACATTTTGTTTGAAAAGCAATAAGTGAGTAATGATGAAAATCCATACCACAGAGGCGCTCATGAAGGCTGAAAt                                                                                                                                                     <  1:118616/150‑1 (MQ=255)
   ggtaaaCCCGCAGCACACTAATCTGGTCGAGATCGAAAAAATATAACCCACGACAACCATAAATTCTGGCATGACATTTTGTTTGAAAAGCAATAAGTGAGTAATGATGAAAATCCATACCACAGAGGCGCTCATGAAGGCTGAAATCAg                                                                                                                                                  <  1:895409/147‑1 (MQ=255)
    ggAAACCCGCAGCACACTAATCTGGTCGAGATCGAAAAAATATAACCCACGACAACCATAAATTCTGGCATGACATTTTGTTTGAAAAGCAATAAGTGAGTAATGATGAAAATCCATACCACAGAGGCGCTCATGAAGGCTGAAATCAGc                                                                                                                                                 <  2:551907/150‑1 (MQ=255)
      aaaCCCGCAGCACACTAATCTGGTCGAGATCGAAAAAATATAACCCACGACAACCATAAATTCTGGCATGACATTTTGTTTGAAAAGCAATAAGTGAGTAATGATGAAAATCCATACCACAGAGGCGCTCATGAAGGCTGAAATCAGCCg                                                                                                                                               <  2:1302442/150‑1 (MQ=255)
        aCCCGCAGCACACTAATCTGGTCGAGATCGAAAAAATATAACCCACGAAAACCATAAATTCTTGCATGACATTTTGTTTGAAAAGCAATAAGTGAGTAATGATGAAAATCCATACCACAGAGGCGCTCATGAAGGCTGAAATCAGCCGTa                                                                                                                                             <  2:141352/150‑1 (MQ=255)
         cccGCAGCACACTAATCTGGTCGAGATCGAAAAAATATAACCCACGACATCCATAAATTCTGGCATGACATTTTGTTTGAAAAGCAATAAGTGAGTAATGATGAAAATCCATACCACAGAGGCGCTCATGAAGGCTGAAATCAGCCGTAg                                                                                                                                            <  2:1501067/150‑1 (MQ=255)
           cGCAGCTCACTAATCTGGTCGAGATCGAAAAAATATAACCCACGATAACCATAAATTCTGGCATGACATTTTGTTTGAAAAGCAATAAGTGAGTAATGATGAAAATCCATACCACAGAGGCGCTCATGAAGGCTGATTTCAGccttcatg                                                                                                                                          <  1:313303/150‑7 (MQ=255)
              agcaCACTAATCTGGTCGAGATCGAAAAAATATAACCCACGACAACCATAAATTCTGGCATGACATTTTGTTTGAAAAGCAATAAGTGAGTAATGATGAAAATCCATACCACAGAGGCGCTCATGAAGGCTGAAATCAGCCGTAGAAGTc                                                                                                                                       <  2:61874/150‑1 (MQ=255)
               gcaCACTAATCTGGTCGAGATCGAAAAAATATAACCCACGACAACCATAAATTCTGGCATGACATTTTGTTTGAAAAGCAATAAGTGAGTAATGATGAAAATCCATACCACAGAGGCGCTCATGAAGGCTGAAATCAGCCGTAGAAGTCt                                                                                                                                      <  1:1147632/150‑1 (MQ=255)
                cacaCTAATCTGGTCGAGATCGAAAAAATATAACCCACGACAACCATAAATTCTGGCATGACATTTTGTTTGAAAAGCAATAAGTGAGTAATGATGAAAATCCATACCACAGAGGCGCTCATGAAGGCTGAAATCAGCCGTAGAAGTCtg                                                                                                                                     <  2:1372461/150‑1 (MQ=255)
                 acacTAATCTGGTCGAGATCGAAAAAATATAACCCACGACAACCATAAATTCTGGCATGACATTTTGTTTGAAAAGCAATAAGTGAGTAATGATGAAAATCCATACCACAGAGGCGCTCATGAAGGCTGAAATCAGCCGTAGAAGTCtga                                                                                                                                    <  2:1023999/150‑1 (MQ=255)
                 acacTAATCTGGTCGAGATCCAAAAAATATAACCCACGACAACCATAAATTCTGGCATGCCATTTTGTTTGAAAAGCAATAAGTGACTAATAATGAAAATCCATACCCCAGAGGCCCTCATGAAGGCTGAAATCAGCCGAAGAAGACtga                                                                                                                                    >  1:1026112/1‑150 (MQ=255)
                  cacTAATCTGGTCGAGATCGAAAAAATATAACCCACGACAACCATAAATTCTGGCATGACATTTTGTTTGAAAAGCAATAAGTGAGTAATGATGAAAATCCATACCACAGAGGCGCTCATGAAGGCTGAAATCAGCCGTAGAAGTCtcaa                                                                                                                                   >  1:300201/1‑147 (MQ=255)
                   acTAATCTGGTCGAGATCGAAAAAATATAACCCACGACAACCATAAATTCTGGCATGACATTTTGTTTGAAAAGCAATAAGTGAGTAATGATGAAAATCCATACCACAGAGGCGCTCATGAAGGCTGAAATCAGCCGTAGAAGTCtgatg                                                                                                                                  <  1:317777/150‑1 (MQ=255)
                   acTAATCTGGTCGAGATCGAAAAAATATAACCCACGACAACCATAAATTCTGGCATGACATTTTGTTTGAAAAGCAATAAGTGAGTAATGATGAAAATCCATACCACAGAGGCGCTCATGAAGGCTGAAATCAGCCGTAGAAGTCtgatg                                                                                                                                  <  1:1700278/150‑1 (MQ=255)
                    cTAATCTGGTCGAGATCGAAAAAATATAACCCACGACAACCATAAATTCTGGCATGACATTTTGTTTGAAAAGAAATAAGTGAGTAATGATGAAAATCCATACCACAGAGGCGCTCATGAAGAATGAAATCAGACGTAGAAGACagatga                                                                                                                                 >  1:132241/1‑150 (MQ=255)
                      aaTCTGGTCGAGATCGAAAAAATATAACCCACGACAACCATAAATTCTGGCATGACATTTTGTTTGAAAAGCAATAAGTGAGTAATGATGAAAATCCATACCACAGAGGCGCTCATGAAGGCTGAAATCAGCCGTAGAAGTCTGATGaaa                                                                                                                               <  1:738515/150‑1 (MQ=255)
                               gagaTCGAAAAAATATAACCCACGACAACCATAAATTCTGGCATGACATTTTGTTTGAAAAGCAATAAGTGAGTAATGATGAAAATCCATACCACAGAGGCGCTCATGAAGGCTGAAATCAGCCGTAGAAGTCTGATGAAAACCTCCGCa                                                                                                                      >  1:1522884/1‑150 (MQ=255)
                                     taaaaaaTATAACCCACGACAACCATAAATTCTGGCATGACATTTTGTTTGAAAAGCAATAAGTGAGTAATGATGAAAATCCATACCACAGAGGCGCTCATGAAGGCTGAAATCAGCCGTAGAAGTCTGATGAAAACCTCCGCACTTGGc                                                                                                                <  1:1392941/149‑1 (MQ=255)
                                     gAAAAAATATAACCCACGACAACCATAAATTCTGGCATGACATTTTGTTTGAAAAGCAATAAGTGAGTAATGATGAAAATCCATACCACAGAGGCGCTCATGAAGGCTGAAATCAGCCGTAGAAGTCTGATGAAAACCTCCGCACTTGGc                                                                                                                <  1:1401258/150‑1 (MQ=255)
                                          aaTATAACCCACGACAACCATAAATTCTGGCATGACATTTTGTTTGAAAAGCAATAAGTGAGTAATGATGAAAATCCATACCACAGAGGCGCTCATGAAGGCTGAAATCAGCCGTAGAAGTCTGATGAAAACCTCCGCACTTGGCAGTCt                                                                                                           >  2:1419696/1‑150 (MQ=255)
                                           atatAACCCACGACAACCATAAATTCTGGCATGACATTTTGTTTGAAAAGCAATAAGTGAGTAATGATGAAAATCCATACCACAGAGGCGCTCATGAAGGCTGAAATCAGCCGTAGAAGTCTGATGAAAACCTCCGCACTTGGCAGTCtt                                                                                                          <  1:690146/150‑1 (MQ=255)
                                              taACCCACGACAACCATAAATTCTGGCATGACATTTTGTTTGAAAAGCAATAAGTGAGTAATGATGAAAATCCATACCACAGAGGCGCTCATGAAGGCTGAAATCAGCCGTAGAAGTCTGATGAAAACCTCCGCACTTGGCAGTCTTgcg                                                                                                       <  2:1211438/150‑1 (MQ=255)
                                                aCCCACGACAACCATAAATTCTGGCATGACATTTTGTTTGAAAAGCAATAAGTGAGTAATGATGAAAATCCATACCACAGAGGCGCTCATGAAGGCTGAAATCAGCCGTAGAAGTCTGATGAAAACCTCCGCACTTGGCAGTCTTGCGCt                                                                                                     <  1:1545093/150‑1 (MQ=255)
                                                aCCCACGACAACCATAAATTCTGGCATGACATTTTGTTTGAAAAGCAATAAGTGAGTAATGATGAAAATCCATACCACAGAGGCGCTCATGAAGGCTGAAATCAGCCGTAGAAGTCTGATGAAAACCTCCGCACTTGGCAGTCTTGCGCt                                                                                                     >  2:21528/1‑150 (MQ=255)
                                                 cccACGACAACCATAAATTCTGGCATGACATTTTGTTTGAAAAGCAATAAGTGAGTAATGATGAAAATCCATACCACAGAGGCGCTCATGAAGGCTGATTTCAGccttc                                                                                                                                             >  2:1108987/1‑106 (MQ=255)
                                                 cccACGACAACCATAAATTCTGGCATGACATTTTGTTTGAAAAGCAATAAGTGAGTAATGATGAAAATCCATACCACAGAGGCGCTCATGAAGGCTGATTTCAGccttc                                                                                                                                             <  1:1108987/109‑4 (MQ=255)
                                                 cccACGACAACCATAAATTCTGGCATGACATTTTGTTTGAAAAGCAATAAGTGAGTAATGATGAAAATCCATACCACAGAGGCGCTCATGAAGGCTGAAATCAGCCGTAGAAGTCTGATGAAAACCTCCGCACTTGGCAGTCTTGCGCTg                                                                                                    >  2:60323/1‑150 (MQ=255)
                                                 cccACGACAACCATAAATTCTGGCATGACATTTTGTTTGAAAAGCAATAAGTGAGTAATGATGAAAATCCATACCACAGAGGCGCTCATGAAGGCTGAAATCAGCCGTAGAAGTCTGATGAAAACCTCCGCACTTGGCAGTCTTGCGCTg                                                                                                    <  2:568652/150‑1 (MQ=255)
                                                 cccACGACAACCATAAATTCTCGCATGACATTTTGTTTGAAAAGCAATAAGTGAGTAATGATGAAAATCCATACCACAGAGGCGCTCATGAAGGCTGAAATCAGCAGTAGAAGTCTGATGAAAACCTCCGCACTTGGCAGTCTTGCTCTg                                                                                                    >  1:36248/1‑150 (MQ=255)
                                                  ccACGACAACCATAAATTCTGGCATGACATTTTGTTTGAAAAGCAATAAGTGAGTAATGATGAAAATCCATACCACAGAGGCGCTCATGAAGGCTGAAATCAGCCGTAGAAGTCTGATGAAAACCTCCGCACTTGGCAGTCTTGCGCTgg                                                                                                   >  1:14338/1‑150 (MQ=255)
                                                  ccACGACAACCATAAATTCTGGCATGACATTTTGTTTGAAAAGCAATAAGTGAGTAATGATGAAAATCCATACCACAGAGGCGCTCATGAAGGCTGAAATCAGCCGTAGAAGTCTGATGAAAACCTCCGCACTTGGCAGTCTTGCGCTgg                                                                                                   >  1:1256560/1‑150 (MQ=255)
                                                   cACGACAACCATAAATTCTGGCATGACAGTTTGTTTGAAAAGCAATAAGTGAGTAATGATGAAAATCCATACCACAGAGGCGCTCATGAAGGCTGAAATCAGCCGTAGAAGTCTGATGAAAACCTCCGCACTTGGCAGTCATGAGCTGGc                                                                                                  >  2:324728/1‑150 (MQ=255)
                                                           ccATAAATTCTGGCATGACATTTTGTTTGAAAAGCAATAAGTGAGTAATGATGAAAATCCATACCACAGAGGCGCTCATGAAGGCTGAAATCAGCCGTAGAAGTCTGATGAAAACCTCCGCACTTGGCAGTCTTGCGCTGGCAAGCAGTg                                                                                          <  2:1314720/150‑1 (MQ=255)
                                                            cATAAATTCTGGCATGACATTTTGATTGAAAAGCAATAAGTGAGTAATGCTGAAAATCCATACCACAGAGGCGCTCATGAAGGCTGAAATCAGCCGTAGAAGTCTGATGAAAACCTCCGCACTTGGCAGTCTTGCGCTGGCAAGCAGTGc                                                                                         <  1:1580600/150‑1 (MQ=255)
                                                                aaatCTGGCAAGACCTTGTGTTTGAACAGCACTACGTGACTAATGATGACAATCCATACCCCCGCGGCGCGCCGGCCGGCTGATATCAGCCGTAGAAGTCTGATGAAAACCTCCGCACTTGGCAGTCTTGCGCTGGCAAGCAGTGCTTTc                                                                                     <  1:89942/147‑1 (MQ=255)
                                                                aaTTCTGGCATGACATTTTGTTTGAAAAGCAATAAGTGAGTAATGATGAAAATCCATACCACAGAGGCGCTCATGAAGGCTGAAATCAGCCGTAGAAGTCTGATGAAAACCTCCGCACTTGGCAGTCTTGCGCTGGCAAGCAGTGCTTTc                                                                                     >  2:1606686/1‑150 (MQ=255)
                                                                aaTTCTGGCATGACATTTTGTATGAAAAGCAATAAGTGAGTAATGATGAAAATCCATACCACAGAGGCGCTCATGAAGGCTGAAATCAGCCGAAGAAGTCTGATCAAAACCTCCGCACTTGGCAGTCTTGCGCTGGCAAGCAGTGCTTTc                                                                                     <  2:1530084/150‑1 (MQ=255)
                                                                      ggCATGACATTTTGTTTGAAAAGCAATAAGTGAGTAATGATGAAAATCCATACCACAGAGGCGCTCATGAAGGCGGAAATCAGCCGTAGAAGTCTGATGAAAACCTCCGCACTTGGCAGTCTTGCGCTGGCAAGCAGTGCTTTCACTCTg                                                                               <  1:444779/150‑1 (MQ=255)
                                                                       gCATGACATTTTGTTTGAAAAGCAATAAGTGAGTAATGATGAAAATCCATACCACAGAGGCGCTCATGAAGGCTGAAATCAGCCGTAGAAGTCTGATGAAAACCTCCGCACTTGGCAGTCTTGCGCTGGCAAGCAGTGCTTTCACTCTGc                                                                              <  2:129481/150‑1 (MQ=255)
                                                                        gaTGACATTTTGTGTGAAAAGAAATAATTGAGTAATGATGAAAATCCATACCACAGCGGCGAGCAAGACGGCTGAAATCAGCCGTAGAAGTCTTATGAAAACCTACGCACTTTGCAGTATTTCGATGGCAAGCAGTGCTTTCACTCTGcc                                                                             <  2:881085/149‑1 (MQ=255)
                                                                        cATGACATTTTGTTTGAAAAGCAATAAGTGAGTAATGATGAAAATCCATACCACAGAGGCGCTCATGAAGGCTGAAATCAGCCGTAGAAGTCTGATGAAAACCTCCGCACTTGGCAGTCTTGCGCTGGCAAGCAGTGCTTTCACTCTGcc                                                                             >  2:203810/1‑150 (MQ=255)
                                                                         aTGACATTTTGTTTGAAAAGCAATAAGTGAGTAATGATGAAAATCCATACCACAGAGGCGCTCATGAAGGCTGAAATCAGCCGTAGAAGTCTGATGAAAACCTCCGCACTTGGCAGTCTTGCGCTGGCAAGCAGTGCTTTCACTCTGCCa                                                                            >  2:1300179/1‑150 (MQ=255)
                                                                           gACATTTTGTTTGAAAAGCAATAAGTGAGTAATGATGAAAATCCATACCACAGAGGCGCTCATGAAGGCTGAAATCAGCCGTAGAAGTCTGATGAAAACCTCCGCACATGGCAGTCTTGCGCTGGCAAGCAGTGCTTTCACTCTGCCAtt                                                                          <  2:1237471/150‑1 (MQ=255)
                                                                              aTTTTGTTTGAAAAGCAATAAGTGAGTAATGATGAAAATCCATACCACAGAGGCGCTCATGAAGGCTGAAATCAGCCGTAGAAGTCTGATGAAAACCTCCGCACTTGGCAGTCTTGCGCTGGCAAGCAGTGCTTTCACTCTGCCATTTTc                                                                       <  1:1357216/150‑1 (MQ=255)
                                                                              aTTTTGTTTGAAAAGCAATAAGTGAGTAATGATGAAAATCCATACCACAGAGGCGCTCATGAAGGCTGAAATCAGCCGTAGAAGTCTGATGAAAACCTCCGCACTTGGCAGTCTTGCGCTGGCAAGCAGTGCTTTCACTCTGCCATTTTc                                                                       <  1:1350454/150‑1 (MQ=255)
                                                                              aTTTTGTTTGAAAAGCAATAAGTGAGTAATGATGAAAATCCATACCACAGAGGCGCTCATGAAGGCTGAAATCAGCCGTAGAAGTCTGATGAAAACCTCCGCACTTGGCAGTCTTGCGCTGGCAAGCAGTGCTTTCACTCTGCCATTTTc                                                                       <  1:1619887/150‑1 (MQ=255)
                                                                                tttgtttgAAAAGCAATAAGTGAGTAATGATGAAAATCCATACCACAGAGGCGCTCATGAAGGCTGAAATCAGCCGTAGAAGTCTGATGAAAACCTCCGCACTTGGCAGTCTTGCGCTGGCAAGCAGTGCTTTCACTCTGCCATTTTccc                                                                     >  1:669814/1‑150 (MQ=255)
                                                                                tttgtttgAAAAGCAATAAGTGAGTAATGATGAAAATCCATACCACAGAGGCGCTCATGAAGGCTGAAATCAGCCGTAGAAGTCTGATGAAAACCTCCGCACTTGGCAGTCTAGCGCTGCCAAGCAGTGCTTTCACTCTGCCATTATccc                                                                     >  1:551657/1‑150 (MQ=255)
                                                                                  tgtttgAAAAGCAATAAGTGAGTAATGATGAAAATCCATACCACAGAGGCGCTCATGAAGGCTGAAATCAGCCGTAGAAGTCTGATGAAAACCTCCGCACTTGGCAGTCTTGCGCTGGCAAGCAGTGCTTTCACTCTGCCATTTTCCCaa                                                                   <  1:771933/150‑1 (MQ=255)
                                                                                  tgtttgAAAAGCAATAAGTGAGTAATGATGAAAATCCATACCACAGAGGCGCTCATGAAGGATGAAATCAGCCGTAGAAGTCTGATGAAAACCTCCGCACTTGGCAGTCTTGCGCTGGCAAGCAGTGCTTACACACTGCCATGTTCCCaa                                                                   >  1:1314686/1‑150 (MQ=255)
                                                                                      tgAAAAGCAATAAGTGAGTAATGATGAAAATCCATACCACAGAGGCGCTCATGAAGGCTGAAATCAGCCGTAGAAGTCTGATGAACACCTCCGCACTTGGCAGTCTTCCGCTGGCAAGAAGTGCATTCACTCTGCAATTTTCCCAAATgg                                                               >  2:102400/1‑150 (MQ=255)
                                                                                      tgAAAAGCAATAAGTGAGTAATGATGAAAATCCATACCACAGAGGCGCTCATGAAGGCTGAAATCAGCCGTAGAAGTCTGATGAAAACCTCCGCACTTGGCAGTCTTGCGCTGGCAAGCAGTGCTTTCACTCTGCCATTTTCCCAAATgg                                                               <  1:864608/150‑1 (MQ=255)
                                                                                       gAAAAGCAATAAGTGAGTAATGATGAAAATCCATACCACAGAGGCGCTCATGAAGGCTGAAATCAGCCGTAGAAGTCTGATGAAAACCTCCGCACTTGGCAGTCTTGCGCTGGCAAGCAGTGCTTTCACTCTGCCATTTTCCCAAATGGt                                                              >  1:224810/1‑150 (MQ=255)
                                                                                       gAAAAGCAATAAGTGAGTAATGATGAAAATCCATACCACAGAGGCGCTCATGAAGGCTGAAATCAGCCGTAGAAGTCTGATGAAAACCTCCGCACTTGGCAGTCTTGCGCTGGCAAGCAGTGCTTTCACTCTGCCATTTTCCCAAATGGt                                                              >  2:633938/1‑150 (MQ=255)
                                                                                        aaaaGCAATAAGTGAGTAATGATGAAAATCCATACCACAGAGGCGCTCATGAAGGCTGAAATCAGCCGTAGAAGTCTGATGAAAACCTCCGCACTTGGCAGTCTTGCGCTGGCAAGCAGTGCTTTCACTCTGCCATTTTCCCAAATGGTc                                                             >  2:902042/1‑150 (MQ=255)
                                                                                        aaaaGCAATAAGTGAGTAATGATGAAAATCCATACCACAGAGGCGCTCATGAAGGCTGAAATCAGCCGTAGAAGTCTGATGAAAACCTCCGCACTTGGCAGTCTTGCGCTGGCAAGCAGTGCTTTCACTCTGCCATTTTCCCAAATGGTc                                                             >  2:888599/1‑150 (MQ=255)
                                                                                          aaGCAATAAGTGAGTAATGATGAAAATCCATACCACAGAGGCGCTCATGAAGGCTGAAATCAGCCGTAGAAGTCTGATGAAAACCTCCGCACTTGGCAGTCTTGCGCTGGCAAGCAGTGCTTTCACTCTGCCATTTTCCCAAATGGTCCg                                                           >  1:704366/1‑150 (MQ=255)
                                                                                           aGCAATAATTTAGTAATGATGAAAATCCATACCACAGAGGCGCTCATGAAGGCTGAAATCAGCCGTAGAAGTCTGATGAAAACCTCCGCACTTGGCAGTCTTGCGCTGGCAAGCAGTGCTTTCACTCTGCCATTTTCCCAAATGGTCCgg                                                          <  1:437763/150‑1 (MQ=255)
                                                                                              aaTAAGTGAGTAATGATGAAAATCCATACCACAGAGGCGCTCATGAAGGCTGAAATCAGCCGTAGAAGTCTGATGAAAACCTCCGCACTTGGCAGTCTTGCGCTGGCAAGCAGTGCTTTCACTCTGCCATTTTCCCAAATGGTCCGggcg                                                       >  2:178982/1‑150 (MQ=255)
                                                                                                  aGTGAGTAATGATGAAAATCCATACCACAGAGGCGCTCATGAAGGCTGAAATCAGCCGTAGAAGTCTGATGAAAACCTCCGCACTTGGCAGTCTTGCGCTGGCAAGCAGTGCTTTCACTCTGCCATTTTCCCAAATGGTCCGGGCGGCag                                                   >  2:1247802/1‑150 (MQ=255)
                                                                                                   gTGAGTAATGATGAAAATCCATACCACAGAGGCGCTCATGAAGGCTGAAATCAGCCGTAGAAGTCTGATGAAAACCTCCGCACTTGGCAGTCTTGCGCTGGCAAGCAGTGCTTTCACTCTGCCATTTTCCCAAATGGTCCGGGCGGCaga                                                  <  1:564595/150‑1 (MQ=255)
                                                                                                      aGTAATGATGAAAATCCATACCACAGAGGCGCTCATGAAGGCTGAAATCAGCCGTAGCAGTCTGATGAAAACCTCCGCACTTGGCAGTCTTGCGCTGGCAAGCAGTGCTTTCACTCTGCCATTTTCCCAAATGGTCCGGGCGGCAGAGGc                                               <  1:681128/150‑1 (MQ=255)
                                                                                                       gTAATGATGAAAATCCATACCACAGAGGCGCTCATGAAGGCTGAAATCAGCCGTAGAAGTCTGATGAAAACCTCCGCACTTGGCAGTCTTGCGCTGGCAAGCAGTGCTTTCACTCTGCCATTTTCCCAAATGGTCCGGGCGGCAGAGGCt                                              >  2:509813/1‑150 (MQ=255)
                                                                                                         aaTGATGAAAATCCATACCACAGAGGCGCTCATGAAGGCTGAAATCAGCCGTAGAAGTCTGATGAAAACCTCCGCACTTGGCAGTCTTGCGCTGGCAAGCAGTGCTTTCACTCTGCCATTTTCCCAAATGGTCCGGGCGGCAGAGGCTcc                                            <  2:179814/150‑1 (MQ=255)
                                                                                                           ttatgaAAATCCATACCACAGAGGCGCTCATGATGGCTTAAATCAGCCGTAGAAGTCTGATGAAAACCTCCGCACTTGGCAGTCTTGCGCTGGCAAGCAGTGCTTTCACTCTGCCATTTTCACAAATGGTCCGGGCGGCAGAGGCTCCgg                                          <  1:958907/148‑1 (MQ=255)
                                                                                                           tgatgaAAATCCATACCACAGAGGCGCTCATGAAGGCTGAAATCAGCCGTAGAAGTCTGATGAAAACCTCCGCACTTGGCAGTCTTGCGCTGGCAAGCAGTGCTTTCACTCTGCCATTTTCCCAAATGGTCCGGGCGGCAGAGGCTCCgg                                          >  2:884822/1‑150 (MQ=255)
                                                                                                               gaAAATCCATACCACAGCGGCGCTCATGAAGGCTGAAATCAGCCGTAGAAGTCTGATGAAAACCTCCGCACTTGGCAGTCTTGCGCTGGCAAGCAGTGCTTTCACTCTGCCATTTTCCCAAATGGTCCGGGCGGCAGAGGCTCCGGTaga                                      >  2:1305421/1‑150 (MQ=255)
                                                                                                                 aaaTCCATACCACAGAGGCGCTCATGAAGGCTGAAATCAGCCGTAGAAGTCTGATGAAAACCTCCGCACTTGGCAGTCTTGCGCTGGCAAGCAGTGCTTTCACTCTGCCATTTTCCCAAATGGTCCGGGCGGCAGAGGCTCCGGTagaag                                    <  1:1306897/150‑1 (MQ=255)
                                                                                                                 aaaTCCATACCACAGAGGCGCTCATGAAGGCTGAAATCAGCCGTAGAAGTCTGATGAAAACCTCCGCACTTGGCAGTCTTGCGCTGGCAAGCAGTGCTTTCACTCTGCCATTTTCCCAAATGGTCCGGGCGGCAGAGGCTCCGGTagaag                                    <  1:26160/150‑1 (MQ=255)
                                                                                                                   aTCCATACCACAGAGGCGCTCATGAAGGCTGAAATCAGCCGTAGAAGTCTGATGAAAACCTCCGCACTTGGCAGTCTTGCGCTGGCAAGCAGTGCTTTCACTCTGCCATTTTCCCAAATGGTCCGGGCGGCAGAGGCTCCGGTAGAagag                                  >  2:1278653/1‑150 (MQ=255)
                                                                                                                   aTCCATACCACAGAGGCGCTCATGAAGGCTGAAATCAGCCGTAGAAGTCTGATGAAAACCTCCGCACTTGGCAGTCTTGCGCTGGCAAGCAGTGCTTTCACTCTGCCATTTTCCCAAATGGTCCGGGCGGCAGAGGCTCCGGTAGAagag                                  <  1:1621841/150‑1 (MQ=255)
                                                                                                                       aTACCACAGTGGCGCTCCTGACGGCTTAAATCAGCCGTATAAGTCTGATGAAAACCTCCGCACTTGTCAGTCTTGCGCTGGCAAGCAGTGCTATCACTCTGCCATTTTCCCAAATGGTCCGGGCGGCAGAGGCTCCGGTAGAAGAGAAAg                              <  2:236373/150‑1 (MQ=255)
                                                                                                                       aTACCACAGAGGCGCTCATGAAGGCGGAAATCAGCCGTAGAAGTCTGATGAAAACCTCCGCACTGGGCAGTCTTGCGCTGGCAAGCAGTGCTTTCACTCTGCCATTTTCCCAAATGGTCCGGGCGGCAGAGGCTCCGGTAGAAGAGAAAg                              <  2:94611/150‑1 (MQ=255)
                                                                                                                        tACCACAGAGGCGCTCATGAAGGCTGAAATCAGCCGTAGAAGTCTGATGAAAACCTCCGCACTTGGCAGGCTTGCGCTGGCAAGCAGTGCTTTCACTCTGCCATTTTCCCAAATGGTCCGGGCGGCAGAGGCTCCGGTAGAAGAGAAAGc                             <  1:773931/150‑1 (MQ=255)
                                                                                                                         aCCACAGAGGCGCTCATGAAGGCTGAAATCAGCCGTAGAAGTCTGATGAAAACCTCCGCACTTGGCAGTCTTGCGCTGGCAAGCAGTGCTTTCACTCTGCCATTTTCCCAAATGGTCCGGGCGGCAGAGGCTCCGGTAGAAGAGAAAGCg                            >  1:1717990/1‑150 (MQ=255)
                                                                                                                          ccACAGAGGCGCTCATGAAGGCTGAAATCAGCCGTAGAAGTCTGATGAAAACCTCCGCACTTTGCAGTCTTGCGCTGGCAAGCAGGGCTTTCACTCTGCCATTTTCGCAAATGGTCCGGGCGGCAGAGGCTCAGGTAGAAGAGAAAGCgg                           >  2:220774/1‑150 (MQ=255)
                                                                                                                                agGCGCTCATGAAGGCTGAAATCAGCCGTAGAAGTCTGATGAAAACCTCCGCACTTGGCAGTCTTGCGCTGGCAAGCAGTGCTTTCACTCTGCCATTTTCCCAAATGGTCCGGGCGGCAGAGGCTCCGGTAGAAGAGAAAGCGGTCTggg                     >  1:599351/1‑149 (MQ=255)
                                                                                                                                     cTCATGAAGGCTGAAATCAGCCGTAGAAGTCTGATGAAAACCTCCGCACTTGGCAGTCTTGCGCTGGCAAGCAGTGCTTTCACTCTGCCATTTTCCCAAATGGTCCGGGCGGCAGAGGCTCCGGTAGAAGAGAAAGCGGTCTGGAGTTcc                >  2:1625506/1‑150 (MQ=255)
                                                                                                                                        tctacggCTGATTTCAGCCGTAGAAGTCTGATGAAAACCTCCGCACTTGGCAGTCTTGCGCTGGCAAGCAGTGCTTTCACTCTGCCATTTTCCCAAATGGTCCGGGCGGCAGAGGCTCCGGTAGAAGAGAAAGCGGTCTGGAGTTCCTGc             >  2:1533833/6‑150 (MQ=255)
                                                                                                                                        aTGAAGGCGGAAATCAGCCGTAGAAGTCTGATGAAAACCTCCGCACTTGGCAGTCTTGCGCTGGCAAGCAGTGCTTTCACTCTGCCATTTTCCCAAATGGTCCGGGCGGCAGAGGCTCCGGTAGAAGAGAAAGCGGTCTGGAGTTCCTGc             <  1:927241/150‑1 (MQ=255)
                                                                                                                                            atgatGAAATCAGACGTAGAAGTCTGATGAAAACCTCAGCACTTGGCAGTCTTGCGATGGCAATCAGTGATTTAACTCTGACATTTTCCCAAATGGTCAGGGCGGCAGAGGCTCCGGTAGAATATAAAGCGGTCTGGAGTTCATGCACCg         >  2:1021269/5‑150 (MQ=255)
                                                                                                                                               cTGAAATCAGCCGTAGAAGTCTGATGAAAACCTCCGCACTTGGCAGTCTTGCGCTGGCAAGCAGTGCTTTCACTCTGCCATTTTCCCAAATGGTCCGGGCGGCAGAGGCTCCGGTAGAAGAGAAAGCGGTCTGGAGTTCCTGCACCGTTa      >  2:887750/1‑150 (MQ=255)
                                                                                                                                                   aaTCAGCCGTAGAAGTCTGATGAAAACCTCCGCACTTGGCAGTCTTGCGCTGGCAAGCAGTGCTTTCACTCTGCCATTTTCCCAAATGGTCCGGGCGGCAGAGGCTCCGGTAGAAGAGAAAGCGGTCTGGAGTTCCTGCACCGTTAACTg  >  1:1463253/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
GAAGGGAAACCCGCAGCACACTAATCTGGTCGAGATCGAAAAAATATAACCCACGACAACCATAAATTCTGGCATGACATTTTGTTTGAAAAGCAATAAGTGAGTAATGATGAAAATCCATACCACAGAGGCGCTCATGAAGGCTGAAATCAGCCGTAGAAGTCTGATGAAAACCTCCGCACTTGGCAGTCTTGCGCTGGCAAGCAGTGCTTTCACTCTGCCATTTTCCCAAATGGTCCGGGCGGCAGAGGCTCCGGTAGAAGAGAAAGCGGTCTGGAGTTCCTGCACCGTTAACTG  >  NC_000913/1660447‑1660743

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