Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,614,832 T→C 8.7% I389T (ATT→ACT) ‡ nikA → Ni(2(+)) ABC transporter periplasmic binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,614,8320TC 8.7% 228.4 / 10.1 80I389T (ATT→ACT) ‡nikANi(2(+)) ABC transporter periplasmic binding protein
Reads supporting (aligned to +/- strand):  ref base T (31/42);  new base C (3/4);  total (34/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 = 1.00e+00

TTGCTGGAAAAAGCCGGTTGGACGCTGCCTGCGGGCAAAGACATCCGCGAGAAAAATGGTCAGCCGCTGCGCATTGAACTTTCGTTCATCGGCACCGATGCGTTAAGCAAATCGATGGCGGAAATCATTCAGGCTGATATGCGCCAGATTGGCGCAGATGTCTCGCTGATTGGCGAAGAAGAGAGCAGTATCTATGCTCGTCAGCGCGACGGTCGTTTTGGCATGATTTTCCACCGCACCTGGGGCGCGCCATATGATCCACACGCCTTCCTCAGTTCAATGCGCGTACCGTCACA  >  NC_000913/3614684‑3614979
                                                                                                                                                    |                                                                                                                                                   
ttGCTGGAAAAAGCCGGTTGGACGCTGCCTGCGGGCAAAGACATCCGCGAGAAAAATGGTCAGCCGCTGCGCATTGAACTTTCGTTCATCGGCACCGATGCGTTAAGCAAATCGATGGCGGAAATCATTCAGGCTGATATGCGCCAGAtt                                                                                                                                                    <  2:60223/150‑1 (MQ=255)
ttGCTGGAAAAAGCCGGTTGGACGCTGCCTGCGGGCAAAGACATCCGCGAGAAAAATGGTCAGCCGCTGCGCATTGAACTTTCGTTCATCGGCACCGATGCGTTAAGCAAATCGATGGCGGAAATCATTCAGGCTGATATGCGCCAGAtt                                                                                                                                                    <  2:601501/150‑1 (MQ=255)
  gCTGGAAAAAGCCGGTTGGACGCTGCCTGCGGGCAAAGACATCCGCGAGAAAAATGGTCAGCCGCTGCGCATTGAACTTTAGTTCATCGGCACCGATGCGTTAAGCAAATCGATGGCGGAAATCATTCAGGCTGATATGCGCCAGATTgg                                                                                                                                                  >  2:783655/1‑150 (MQ=255)
     ggAAAAAGCCGGTTGGACGCTGCCTGCGGGCAAAGACATCCGCGAGAAAAATGGTCAGCCGCTGCGCATTGAACTTTCGTTCATCGGCACCGATGCGTTAAGCAAATCGATGGCGGAAATCATTCAGGCTGATATGCGCCAGATTGgcgc                                                                                                                                               <  2:225470/150‑1 (MQ=255)
      gAAAAAGCCGGTTGGACGCTGCCTGCGGGCAAAGACATCCGCGAGAAAAATGGTCAGCCGCTGCGCATTGAACTTTCGTTCATCGGCACCGATGCGTTAAGCAAATCGATGGCGGAAATCATTCAGGCTGATATGCGCCAGATTGGCGCa                                                                                                                                              >  1:426459/1‑150 (MQ=255)
        aaaaGCCGGTTGGACGCTGCCTGCGGGCAAAGACATCCGCGAGAAAAATGGTCAGCCGCTGCGCATTGAACTTTCGTTCATCGGCACCGATGCGTTAAGCAAATCGATGGCGGAAATCATTCAGGCTGATATGCGCCAATCTGGCGCata                                                                                                                                            <  1:1320743/150‑3 (MQ=255)
         aaaGCCGGTTGGACGCTGCCTGCGGGCAAAGACATCCGCGAGAAAAATGGTCAGCCGCTGCGCATTGAACTTTCTTTCATCGGCACCGATGCGTTAAGCAAATCGATGGCGGAAATCATTCAGGCTGATATGCGCCAATCTGGCGCatat                                                                                                                                           <  1:1019195/150‑1 (MQ=255)
         aaaGCCGGTTGGACGCTGCCTGCGGGCAAAGACATCCGCGAGAAAAATGGTCAGCCGCTGCGCATTGAACTTTCGTTCATCGGCACCGATGCGTTAAGCAAATCGATGGCGGAAATCATTCAGGCTGATATGCGCCAGATTGGCGCAGAt                                                                                                                                           <  2:703781/150‑1 (MQ=255)
          aaGCCGGTTGGACGCTGCCTGCGGGCAAAGACATCCGCGAGAAAAATGGTCAGCCGCTGCGCATTGAACTTTCGTTCATCGGCACCGATGCGTTAAGCAAATCGATGGCGGAAATCATTCAGGCTGATATGCGCCAATCTGGCGCatatc                                                                                                                                          <  1:392113/150‑2 (MQ=255)
            gCCGGTTGGACGCTGCCTGCGGGCAAAGACATCCGCGAGAAAAATGGTCAGCCGCTGCGCATTGAACTTTCGTTCATCGGCACCGATGCGTTAAGCAAATCGATGGCGGAAATCATTCAGGCTGATATGCGCCAGATTGGCGCAGATGtc                                                                                                                                        <  1:858460/150‑1 (MQ=255)
            gCCGGTTGGACGCTGCCTGCGGGCAAAGACATCCGCGAGAAAAATGGTCAGCCGCTGCGCATTGAACTTTCGTTCATCGGCACCGATGCGTTAAGCAAATCGATGGCGGAAATCATTCAGGCTGATATGCGCCAATCTGGCGCatatcag                                                                                                                                        <  2:1045835/150‑4 (MQ=255)
              cGGTTGGACGCTGCCTGCGGGCAAAGACATCCGCGAGAAAAATGGTCAGCCGCTGCGCATTGAACTTTCGTTCATCGGCACCGATGCGTTAAGCAAATCGATGGCGGAAATCATTCAGGCTGATATGCGCCAGATTGGCGCAGATGtctc                                                                                                                                      >  1:284012/1‑150 (MQ=255)
               ggTTGGACGCTGCCTGCGGGCAAAGACATCCGCGAGAAAAATGGTCAGCCGCTGCGAATTGAACTTTCGTTCATCGGAACCGATGCGTTAAGCAAATCGATGGCGGAAATCATTCAGGCTGATATGCGCCAGATTGGAGCAGATGTATCg                                                                                                                                     >  1:1282949/1‑150 (MQ=255)
                gTTGGACGCTGCCTGCGGGCAAAGACATCCGCGAGAAAAATGGTCAGCCGCTGCGCATTGAACTTTCGTTCATCGGCACCGATGCGTTAAGCAAATCGATGGCGGAAATCATTCAGGCTGATATGCGCCAGATTGGCGCAGATGTCTCGc                                                                                                                                    <  1:1439055/150‑1 (MQ=255)
                  tGGACGCTGCCTGCGGGCAAAGACATCCGCGAGAAAAATGGTCAGCCGCTGCGCATTGAACTTTCGTTCATCGGCACCGATGCGTTAAGCAAATCGATGGCGGAAATCATTCAGGCTGATATGCGCCAGATTGGCGCAGATGTCTCGCTg                                                                                                                                  >  1:9811/1‑150 (MQ=255)
                    gACGCTGCCTGCGGGCAAAGACATCCGCGAGAAAAATGGTCAGCCGCTGCGCATTGAACTTTCGTTCATCGGCACCGATGCGTTAAGCAAATCGATGGCGGAAATCATTCAGGCTGATATGCGCCAGATTGGCGCAGATGTCTCGCTGAt                                                                                                                                <  2:1246226/150‑1 (MQ=255)
                      cGCTGCCTGCGGGCAAAGACATCCGCGAGAAAAATGGTCAGCCGCTGCGCATTGAACTTTCGTTCATCGGCACCGATGCGTTAAGCAAATCGATGGCGGAAATCATTCAGGCTGATATGCGCCAGATTGGCGCAGATGTCTCGCTGATTg                                                                                                                              <  1:1314239/150‑1 (MQ=255)
                        gtgcctgcGGGCAAAGACATCCGCGAGAAAAATGTTCAGCCGCTGCGCATTGAACTTTCGTTCATCGGCACCGATGCGTTAAGCAAATCGATGGCGGAAATCATTCAGGCTGATATGCGCCAGATTGGCGCAGATGTCTCGCTGATTGGc                                                                                                                            <  2:1614236/149‑1 (MQ=255)
                        ctgcctgcGGGCAAAGACATCCGCGAGAAAAATGGTCAGCCGCTGCGCATTGAACTTTCGTTCATCGGCACCGATGCGTTAAGCAAATCGATGGCGGAAATCATTCAGGCTGATATGCGCCAGATTGGCGCAGATGTCTCGCTGATTGGc                                                                                                                            <  2:426459/150‑1 (MQ=255)
                         tgcctgcGGGCAAAGACATCCGCGAGAAAAATGGTCAGCCGCTGCGCATTGAACTTTCGTTCATCGGCACCGATGCGTTAAGCAAATCGATGGCGGAAATCATTCAGGCTGATATGCGCCAGATTGGCGCAGATGTCTCGCTGATTGGCg                                                                                                                           <  2:1302964/150‑1 (MQ=255)
                          gcctgcGGGCAAAGACATCCGCGAGAAAAATGGTCAGCCGCTGCGCATTGAACTTTCGTTCATCGGCACCGATGCGTTAAGCAAATCGATGGCGGAAATCATTCAGGCTGATATGCGCCAGATTGGCGCAGATGTCTCGCTGATTGGCga                                                                                                                          >  1:258599/1‑150 (MQ=255)
                          gcctgcGGGCAAAGACATCCGCGAGAAAAATGGTCAGCCGCTGCGCATTGAACTTTCGTTCATCGGCACCGATGCGTTAAGCAAATCGATGGCGGAAATCATTCAGGCTGATATGCGCCAGATTGGCGCAGATGTCTCGCTGATTGGCga                                                                                                                          <  2:470358/150‑1 (MQ=255)
                           cctgcGGGCAAAGACATCCGCGAGAAAAATGGTCAGCCGCTGCGCATTGAACTTTCGTTCATCGGCACCGATGCGTTAAGCAAATCGATGGCGGAAATCATTCAGGCTGATATGCGCCAGATTGGCGCAGATGTCTCGCTGATTGGCgaa                                                                                                                         >  1:927861/1‑150 (MQ=255)
                                  gCAAAGACATCCGCGAGAAAAATGGTCAGCCGCTGCGCATTGAACTTTCGTTCATCGGCACCGATGCGTTAAGCAAATCGATGGCGGAAATCATTCAGGCTGATATGCGCCAGATTGGCGCAGATGTCTCGCTGATTGGCGAAGAagaga                                                                                                                  <  2:222283/150‑1 (MQ=255)
                                   cAAAGACATCCGCGAGAAAAATGGTCAGCCGCTGCGCATTGAACTTTCGTTCATCGGCACCGATGCGTTAAGCAAATCGATGGCGGAAATCATTCAGGCTGATATGCGCCAGATTGGCGCAGATGTCTCGCTGATTGGCGAAGAagagag                                                                                                                 >  1:124574/1‑150 (MQ=255)
                                     aaGACATCCGCGAGAAAAGTGGTCAGCCGCTGCGCATTGAACTTTCGTTCATCGGCACCGATGCGTTAAGCAAATCGATGGCGGAAATCATTCAGGCTGATATGCGCCAGATTGGCGCAGATGTCTCGCTGATTGGCGAAGAAGAGAGCa                                                                                                               >  2:333154/1‑150 (MQ=255)
                                      aGACATCCGCGAGAAAAATGGTCAGCCGCTGCGCATTGAACTTTCGTTCATCGGCACCGATGCGTTAAGCAAATCGATGGCGGAAATCATTCAGGCTGATATGCGCCAGATTGGCGCAGATGTCTCGCTGATTGGCGAAGAAGAGAGCAg                                                                                                              <  1:777566/150‑1 (MQ=255)
                                        aCATCCGCGAGAAAAATGGTCAGCCGCTGCGCATTGAACTTTCGTTCATCGGCACCGATGCGTTAAGCAAATCGATGGCGGAAATCATTCAGGCTGATATGCGCCAGATTGGCGCAGATGTCTCGCTGATTGGCGAAGAAGAGAGCAGTa                                                                                                            <  1:1148945/150‑1 (MQ=255)
                                        aCATCCGCGAGAAAAATGGTCAGCCGCTGCGCATTGAACTTTCGTTCATCGGCACCGATGCGTTAAGCAAATCGATGGCGGAAATCATTCAGGCTGATATGCGCCAGATTGGCGCAGATGTCTCGCTGATTGGCGAAGAAGAGAGCAGTa                                                                                                            <  1:1143229/150‑1 (MQ=255)
                                          aTCCGCGAGAAAAATGGTCAGCCGCTGCGCATTGAACTTTCGTTCATCGGCACCGATGCGTTAAGCAAATCGATGGCGGAAATCATTCAGGCTGATATGCGCCAGATTGGCGCAGATGTCTCGCTGATTGGCGAAGAAGAGAGCAGTATc                                                                                                          >  2:479667/1‑150 (MQ=255)
                                             cgcgAGAAAAATGGTCAGCCGCTGCGCATTGAACTTTCGTTCATCGGCACCGATGCGTTAAGCAAATCGATGGCGGAAATCATTCAGGCTGATATGCGCCAGATTGGCGCAGATGTCTCGCTGATTGGCGAAGAAGAGAGCAGTATCTAt                                                                                                       <  1:1018918/150‑1 (MQ=255)
                                                  gaAAAATGGTCAGCCGCTGCGCATTGAACTTTCGTTCATCGGCACCGATGCGTTAAGCAAATCGATGGCGGAAATCATTCAGGCTGATATGCGCCAGATTGGCGCAGATGTCTCGCTGATTGGCGAAGAAGAGAGCAGTATCTATGCTCg                                                                                                  >  1:723192/1‑150 (MQ=255)
                                                        tGGTCAGCCGCTGCGCATTGAACTTTCGTTCATCGGCACCGATGCGTTAAGCAAATCGATGGCGGAAATCATTCAGGCTGATATGCGCCAGATTGGCGCAGATGTCTCGCTGATTGGCGAAGAAGAGAGCAGTATCTATGCTCGTCAgcg                                                                                            <  1:1587878/150‑1 (MQ=255)
                                                        tGGTCAGCCGCTGCGCATTGAACTTTCGTTCATCGGCACCGATGCGTTAAGCAAATCGATGGCGGAAATCATTCAGGCTGATATGCGCCAGATTGGCGCAGATGTCTCGCTGATTGGCGAAGAAGAGAGCAGTATCTATGCTCGTCAgcg                                                                                            >  2:54946/1‑150 (MQ=255)
                                                         ggTCAGCCGCTGCGCATTGAACTTTCGTTCATCGGCACCGATGCGTTAAGCAAATCGATGGCGGAAATCATTCAGGCTGATATGCGCCAGATTGGCGCAGATGTCTCGCTGATTGGCGAAGAAGAGAACAGTATCTATGCTCGTCAgcgc                                                                                           <  2:977943/150‑1 (MQ=255)
                                                         ggTCAGCCGCTGCGCATTGAACTTTCGTTCATCGGCACCGATGCGTTAAGCAAATCGATGGCGGAAATCATTCAGGCGGATATGCGCCAGATTGGCGCAGATGTCTCGCTGATTGGCGAAGAAGAGAACAGTATCTATGCTCGTCAgcgc                                                                                           <  2:971918/150‑1 (MQ=255)
                                                            cAGCCGCTGCGCATTTAACTTTCGTTCATCGGCACCGATGCGTTAAGCAAATCGATGGCGGAAATCATTCAGGCTGATATGCGCCAGATTGGCGCAGATGTCTCGCTGATTGGCGAAGAAGAGAGCAGTATCTATGCTCGTCAGCGCGAc                                                                                        <  2:91590/150‑1 (MQ=255)
                                                                      gcATTGAACTTTCGTTCATCGGCACCGATGCGTTAAGCAAATCGATGGCGGAAATCATTCAGGCTGATATGCGCCAGATTGGCGCAGATGTCTCGCTGATTGGCGAAGAAGAGAGCAGTATCTATGCTCGTCAGCGCGACGGTCGTTTTg                                                                              <  1:577521/150‑1 (MQ=255)
                                                                       cATTGAACTTTCGTTCATCGGCACCGATGCGTTAAGCAAATCGATGGCGGAAATCATTCAGGCTGATATGCGCCAGATTGGCGCAGATGTCTCGCTGATTGGCGAAGAAGAGAGCAGTATCTATGCTCGTCAGCGCGACGGTCg                                                                                   >  1:1489887/1‑144 (MQ=255)
                                                                       cATTGAACTTTCGTTCATCGGCACCGATGCGTTAAGCAAATCGATGGCGGAAATCATTCAGGCTGATATGCGCCAGATTGGCGCAGATGTCTCGCTGATTGGCGAAGAAGAGAGCAGTATCTATGCTCGTCAGCGCGACGGTCg                                                                                   <  2:1489887/144‑1 (MQ=255)
                                                                        aTTGAACTTTCGTTCATCGGCACCTATGCGTTAAGCAAATCGATGGCGGAAATCATTCAGGCTGATATGCGCCAGATTGGCGCAGATGTCTCGCTGATTGGCGAAGAAGAGAGCAGTATCTATGCTCGTCAGCGCGACGGTCGTTTTGGc                                                                            <  2:762609/150‑1 (MQ=255)
                                                                        aTTGAACTTTCGTTCATCGGCACCGATGCGTTAAGCAAATCGATGGCGGAAATCATTCAGGCTGATATGCGCCAGATTGGCGCAGATGTCTCGCTGATTGGCGAAGAAGAGAGCAGTATCTATGCTCGTCAGCGCGACGGTCGTTTTGGc                                                                            <  1:176173/150‑1 (MQ=255)
                                                                           gAACTTTCGTTCATCGGCACCGATGCGTTAAGCAAATCGATGGCGGAAATCATTCAGGCTGATATGCGCCAGATTGGCGCAGATGTCTCGCTGATTGGCGAAGAAGAGAGCAGTATCTATGCTCGTCAGCGCGACGGTCGTTTTGGCATg                                                                         >  2:1504686/1‑150 (MQ=255)
                                                                                   gTTCATCGGCACCGATGCGTTAAGCAAATCGATGGCGGAAATCATTCAGGCTGATATGCGCCAGATTGGCGCAGATGTCTCGCTGATTGGCGAAGAAGAGAGCAGTATCTATGCTCGTCAGCGCGACGGTCGTTTTGGCATGATTTTCCa                                                                 >  2:688916/1‑150 (MQ=255)
                                                                                       aTCGGCACCGATGCGTTAAGCAAATCGATGGCGGAAATCATTCAGGCTGATATGCGCCAGATTGGCGCAGATGTCTCGCTGATTGGCGAAGAAGAGAGCAGTATCTATGCTCGTCAGCGCGACGGTCGTTTTGGCATGATTTTCCACCGc                                                             <  1:1592822/150‑1 (MQ=255)
                                                                                        tCGGCACCGATGCGTTAAGCAAATCGATGGCGGAAATCATTCAGGCTGATATGCGCCAGATTGGCGCAGATGTCTCGCTGATTGGCGAAGAAGAGAGCAGTATCTATGCTCGTCAGCGCGACGGTCGTTTTGGCATGATTTTCCACCGCa                                                            <  2:160982/150‑1 (MQ=255)
                                                                                        tCGGCACCGATGCGTTAAGCAAATCGATGGCGGAAATCATTCAGGCTGATATGCGCCAGATTGGCGCAGATGTCTCGCTGATTGGCGAAGAAGAGAGCAGTATCTATGCTCGTCAGCGCGACGGTCGTTTTGGCATGATTTTCCACCGCa                                                            <  1:1325513/150‑1 (MQ=255)
                                                                                         cGGCACCGATGCGTTAAGCAAATCGATGGCGGAAATCATTCAGGCTGATATGCGCCAGATTGGCGCAGATGTCTCGCTGATTGGCGAAGAAGAGAGCAGTATCTATGCTCGTCAGCGCGACGGTCGTTTTGGCATGATTTTCCACCGCAc                                                           >  2:1606359/1‑150 (MQ=255)
                                                                                          ggCACCGATGCGTTAAGCAAATCGATGGCGGAAATCATTCAGGCTGATATGCGCCAGATTGGCGCAGATGTCTCGCTGATTGGCGAAGAAGAGAGCAGTATCTATGCTCGTCAGCGCGACGGTCGTTTTGGCATGATTTTCCACCGCAcc                                                          >  1:170174/1‑150 (MQ=255)
                                                                                           gCACCGATGCGTTAAGCAAATCGATGGCGGAAATCATTCAGGCTGATATGCGCCAGATTGGCGCAGATGTCTCGCTGATTGGCGAAGAAGAGAGCAGTATCTATGCTCGTCAGCGCGACGGTCGTTTTGGCATGATTTTCCAc                                                                >  2:844686/1‑143 (MQ=255)
                                                                                           gCACCGATGCGTTAAGCAAATCGATGGCGGAAATCATTCAGGCTGATATGCGCCAGATTGGCGCAGATGTCTCGCTGATTGGCGAAGAAGAGAGCAGTATCTATGCTCGTCAGCGCGACGGTCGTTTTGGCATGATTTTCCAc                                                                <  1:844686/143‑1 (MQ=255)
                                                                                            cACCGATGCGTTAAGCAAATCGATGGCGGAAATCATTCAGGCTGATATGCGCCAGATTGGCGCAGCTGTCTCGCTGATTGGCGAAGAAGAGAGCAGTATCTATGCTCGTCAGCGCGACGGTCGTTTTGGCATGATTTTCCACCGCACCTg                                                        >  2:1177073/1‑150 (MQ=255)
                                                                                             aCCGATGCGTTAAGCAAATCGATGGCGGAAATCATTCAGGCTGATATGCGCCAGATTGGCGCAGATGTCTCGCTGATTGGCGAAGAAGAGAGCAGTATCTATGCTCGTCAGCGCGACGGTCGTTTTGGCATGATTTTCCACCGCACCTgg                                                       <  2:9811/150‑1 (MQ=255)
                                                                                               cGATGCGTTAAGCAAATCGATGGCGGAAATCATTCAGGCTGATATGCGCCAGATTGGCGCAGATGTCTCGCTGATTGGCGAAGAAGAGAGCAGTATCTATGCTCGTCAGCGCGACGGTCGTTTTGGCATGATTTTCCACCGCACCTgggg                                                     <  1:1504686/150‑1 (MQ=255)
                                                                                                   gCGTTAAGCAAATCGATGGCGGAAATCATTCAGGCTGATATGCGCCAGATTGGCGCAGATGTCTCGCTGATTGGCGAAGAAGAGAGCAGTATCTATGCTCGTCAGCGCGACGGTCGTTTTGGCATGATTTTCCACCGCACCTGGGgcgcg                                                 <  1:996800/150‑1 (MQ=255)
                                                                                                   gCGTTAAGCAAATCGATGGCGGAAATCATTCAGGCTGATATGCGCCAGATTGGCGCAGATGTCTCGCTGATTGGCGAAGAAGAGAGCAGTATCTATGCTCGTCAGCGCGACGGTCGTTTTGGCATGATTTTCCACCGCACCTGGGgcgcg                                                 <  1:1191219/150‑1 (MQ=255)
                                                                                                   gCGTTAAGCAAATCGATGGCGGAAATCATTCAGGCTGATATGCGCCAGATTGGCGCAGATGTCTCGCTGATTGGCGAAGAAGAGAGCAGTATCTATGCTCGTCAGCGCGACGGTCGTTTTGGCATGATTTTCCACCGCACCTGGGgcgcg                                                 <  2:723192/150‑1 (MQ=255)
                                                                                                     gTTAAGCAAATCGATGGCGGAAATCATTCAGGCTGATATGCGCCAGATTGGCGCAGATGTCTCGCTGATTGGCGAAGAAGAGAGCAGTATCTATGCTCGTCAGCGCGACGGTCGTTTTGGCATGATTTTCCACCGCACCTGGGGCGCGcc                                               <  2:284012/150‑1 (MQ=255)
                                                                                                     gTTAAGCAAATCGATGGCGGAAATCATTCAGGCTGATATGCGCCAGATTGGCGCAGATGTCTCGCTGATTGGCGAAGAAGAGAGCAGTATCTATGCTCGTCAGCGCGACGGTCGTTTTGGCATGATTTTCCACCGCACCCGGGGCGCGcc                                               <  1:479667/150‑1 (MQ=255)
                                                                                                      ttAAGCAAATCGATGGCGGAAATCATTCAGGCTGATATGCGCCAGATTGGCGCAGATGTCTCGCTGATTGGCGAAGAAGAGAGCAGTATCTATGCTCGTCAGCGCGACGGTCGTTTTGGCATGATTTTCCACCGCACCTGGGGCGCGCCa                                              <  1:974309/150‑1 (MQ=255)
                                                                                                       tAAGCAAATCGATGGCGGAAATCATTCAGGCTGATATGCGCCAGATTGGCGCAGATGTCTCGCTTATTGGCGAAGAAGAGAGCAGTATCTATGCTCGTCAGCGCGACGGTCGTTTTGGCATGATTTTCCACCGCACCTGGGGCGCGCCat                                             <  2:1202650/150‑1 (MQ=255)
                                                                                                         aGCAAATCGATGGCGGAAATCATTCAGGCTGATATGCGCCAGATTGGCGCAGATGTCTCGCTGATTGGCGAAGAAGAGAGCAGTATCTATGCTCGTCAGCGCGACGGTCGTTTTGGCATGATTTTCCACCGCACCTGGGGCGCGCCatat                                           <  2:1275898/150‑1 (MQ=255)
                                                                                                         aGCAAATCGATGGCGGAAATCATTCAGGCTGATATGCGCCAGATTGGCGCAGATGTCTCGCTGATTGGCGAAGAAGAGAGCAGTATCTATGCTCGTCAGCGCGACGGTCGTTTTGGCATGATTTTCCACCGCACCTGGGGCGCGCCatat                                           >  2:331539/1‑150 (MQ=255)
                                                                                                                     gCGGAAATCATTCAGGCTGATATGCGCCAGATTGGCGCAGATGTCTCGCTGATTGGCGAAGAAGAGAGCAGTATCTATGCTCGTCAGCGCGACGGTCGTTTTGGCATGATTTTCCACCGCACCTGGGGCGCGCCATATGATCCACACGcc                               >  1:1439171/1‑150 (MQ=255)
                                                                                                                           aTCATTCAGGCTGATATGCGCCAGATTGGCGCAGATGTCTCGCTGATTGGCGAAGAAGAGAGCAGTATCTATGCTCGTCAGCGCGACGGTCGTTTTGGCATGATTTTCCACCGCACCTGGGGCGCGCCATATGATCCACACGCCTTCCTc                         <  1:333154/150‑1 (MQ=255)
                                                                                                                             cATTCAGGCTGATATGCGCCAGATTGGCGCAGATGTCTCGCTGATTGGCGAAGAAGAGAGCAGTATCTATGCTCGTCAGCGCGACGGTCGTTTTGGCATGATTTTCCACCGCACCTGGGGCGCGCCATATGATCCACACGCCTTCCTCAg                       >  2:1593705/1‑150 (MQ=255)
                                                                                                                                 cAGGCTGATATGCGCCAGATTGGCGCAGATGTCTCGCTGATTGGCGAAGAAGAGAGCAGTATCTATGCTCGTCAGCGCGACGGTCGTTTTGGCATGATTTTCCACCGCACCTGGGGCGCGCCATATGATCCACACGCCTTCCTCAGTTCa                   >  2:365560/1‑150 (MQ=255)
                                                                                                                                      tGATATGCGCCAGATTGGCGCAGATGTCTCGCTGATTGGCGAAGAAGAGAGCAGTATCTATGCTCGTCAGCGCGACGGTCGTTTTGGCATGATTTTCCACCGCACCTGGGGCGCGCCATATGATCCACACGCCTTCCTCAGTTCAATgcg              >  1:793173/1‑150 (MQ=255)
                                                                                                                                       catctGCGCCAATCTGTCGCAGATGTCTCGCTGATTGGCGAAGAAGAGAGCAGTATCTATGCTCGTCAGCGCGACGGTCGTTTTGGCATGATTTTCCACCGCACCTGGGGCGCGCCATATGATCCACACGCCTTCCTCAGTTCAATgcgc             >  2:262879/5‑150 (MQ=255)
                                                                                                                                        atctGCGCCAATCTGGCGCAGATGTCTCGCTGATTGGCGAAGAAGAGAGCAGTATCTATGCTCGTCAGCGCGACGGTCGTTTTGGCATGATTTTCCACCGCACCTGGGGCGCGCCATATGATCCACACGCCTTCCTCAGTTCAATGcgcg            >  1:1338240/4‑150 (MQ=255)
                                                                                                                                          ctGCGCCAATCTGGCGCAGATGTCTCGCTGATTGGCGAAGAAGAGAGCAGTATCTATGCTCGTCAGCGCGACGGTCGTTTTGGCATGATTTTCCACCGCACCTGGGGCGCGCCATATGATCCACACGCCTTCCTCAGTTCAATGCGCGTa          >  1:1533846/2‑150 (MQ=255)
                                                                                                                                            gcgcCAGATTGGCGCAGATGTCTCGCTGATTGGCGAAGAAGAGAGCAGTATCTATGCTCGTCAGCGCGACGGTCGTTTTGGCATGATTTTCCACCGCCCCTGGGGCGCGCCATATGATCCACACGCCTTCCTGAGTTCAATGCGCGTAcc        >  2:331407/1‑150 (MQ=255)
                                                                                                                                            gcgcCAGATTGGCGCAGATGTCTCGCTGATTGGCGAAGAAGAGAGCAGTATCTATGCTCGTCAGCGCGACGGTCGTTTTGGCATGATTTTCCACCGCACCTGGGGCGCGCCATATGATCCACACGCCTTCCTCAGTTCAATGCGCGTAcc        >  1:1130196/1‑150 (MQ=255)
                                                                                                                                            gcgcCAGATTGGCGCAGATGTCTCGCTGATTGGCGAAGAAGAGAGCAGTATCTATGCTCGTCAGCGCGACGGTCGTTTTGGCATGATTTTCCACCGCACCTGGGGCGCGCCATATGATCCACACGCCTTCCTCAGTTCAATGCGCGTAcc        >  1:256875/1‑150 (MQ=255)
                                                                                                                                             cgcCAGATTGGCGCAGATGTCTCGCTGATTGGCGAAGAAGAGAGCAGTATCTATGCTCGTCAGCGCGACGGTCGTTTTGGCATGATTTTCCACCGCACCTGGGGCGCGCCATATGATCCACACGCCTTCCTCAGTTCAATGCGCGTACCg       >  2:569547/1‑150 (MQ=255)
                                                                                                                                             cgcCAGATTGGCGCAGATGTCTCGCTGATTGGCGAAGAAGAGAGCAGTATCTATGCTCGTCAGCGCGACGGTCGTTTTGGCATGATTTTCCACCGCACCTGGGGCGCGCCATATGATCCACACGCCTTCCTCAGTTCAATGCGCGTACCg       <  2:479712/150‑1 (MQ=255)
                                                                                                                                             cgcCAGATTGGCGCAGATGTCTCGCTGATTCGCGAAGAAGAGAGCAGTATCTATGCTCGTCAGCGCGACGGACGTTTAGGCATGATTTTCCACCGCCCCTCGCGCGCGCCATATGCTCCACACGCCTTCCACACTACACTCCGCGTAccc       >  2:901805/1‑149 (MQ=255)
                                                                                                                                              gcCAGATTGGCGCAGATGTCTCGCTGATTGGCGAAGAAGAGAGCAGTATCTATGCTCGTCAGCGCGACGGTCGTTTTGGCATGATTTTCCACCGCACCTGGGGCGCGCCATATGATCCACACGCCTTCCTCAGTTCAATGCGCGTACCGt      >  2:1451059/1‑150 (MQ=255)
                                                                                                                                              gcCAGATTGGCGCAGATGTCTCGCTGATTGGCGAAGAAGAGAGCAGTATCTATGCTCGTCAGCGCGACGGTCGTTTTGGCATGATTTTCCACCGCACCTGGGGCGCGCCATATGATCCACACGCCTTCCTCAGTTCAATGCGCGTACCGt      >  1:1123614/1‑150 (MQ=255)
                                                                                                                                                 aGATTGGCGCAGATGTCTCGCTGATTGGCGAAGAAGAGAGCAGTATCTATGCTCGTCAGCGCGACGGTCGTTTTGGCATGATTTTCCACCGCACCTGGGGCGCGCCATATGATCCACACGCCTTCCTCAGTTCAATGCGCGTACCGTcac   <  1:604691/150‑1 (MQ=255)
                                                                                                                                                 aGATTGGCGCAGATGTCTCGCTGATTGGCGAAGAAGAGAGCAGTATCTATGCTCGTCAGCGCGACGGTCGTTTTGGCATGATTTTCCACCGCACCTGGGGCGCGCCATATGATCCACACGCCTTCCTCAGTTCAATGCGCGTACCGTcac   <  1:599108/150‑1 (MQ=255)
                                                                                                                                                  gATTGGCGCAGATGTCTCGCTGATTGGCGAAGAAGAGAGCAGTATCTATGCTCGTCAGCGCGACGGTCGTTTTGGCATGATTTTCCACCGCACCTGGGGCGCGCCATATGATCCACACGCCTTCCTCAGTTCAATGCGCGTACCGTcaca  <  2:1296033/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
TTGCTGGAAAAAGCCGGTTGGACGCTGCCTGCGGGCAAAGACATCCGCGAGAAAAATGGTCAGCCGCTGCGCATTGAACTTTCGTTCATCGGCACCGATGCGTTAAGCAAATCGATGGCGGAAATCATTCAGGCTGATATGCGCCAGATTGGCGCAGATGTCTCGCTGATTGGCGAAGAAGAGAGCAGTATCTATGCTCGTCAGCGCGACGGTCGTTTTGGCATGATTTTCCACCGCACCTGGGGCGCGCCATATGATCCACACGCCTTCCTCAGTTCAATGCGCGTACCGTCACA  >  NC_000913/3614684‑3614979

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: