Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,614,830 G→A 8.2% Q388Q (CAG→CAA)  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,8300GA 8.2% 183.3 / 5.2 61Q388Q (CAG→CAA) nikANi(2(+)) ABC transporter periplasmic binding protein
Reads supporting (aligned to +/- strand):  ref base G (28/28);  new base A (3/2);  total (31/30)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.46e-01

CTGGAAAAAGCCGGTTGGACGCTGCCTGCGGGCAAAGACATCCGCGAGAAAAATGGTCAGCCGCTGCGCATTGAACTTTCGTTCATCGGCACCGATGCGTTAAGCAAATCGATGGCGGAAATCATTCAGGCTGATATGCGCCAGATTGGCGCAGATGTCTCGCTGATTGGCGAAGAAGAGAGCAGTATCTATGCTCGTCAGCGCGACGGTCGTTTTGGCATGATTTTCCACCGCACCTGGGGCGCGCCATATGATCCACACGCCTTCCTCAGTTCAATGCGCGTACCGTCA  >  NC_000913/3614687‑3614977
                                                                                                                                               |                                                                                                                                                   
cTGGAAAAAGCCGGTTGGACGCTGCCTGCGGGCAAAGACATCCGCGAGAAAAATGGTCAGCCGCTGCGCATTGAACTTTCGTTCATCGGCACCGATGCGTTAAGCAAATCGATGGCGGAAATCATTCAGGCTGATATGCGCCAGATTGgc                                                                                                                                               >  1:1189630/1‑150 (MQ=255)
 tGGAAAAAGCCGGTTGGACGCTGCCTGCGGGCAAAGACATCCGCGAGAAAAATGGTCAGCCGCTGCGCATTGAACTTTCGTTCATCGGCACCGATGCGTTAAGCAAATCGATGGCGGAAATCATTCAGGCTGATATGCGCCAGATTGgcg                                                                                                                                              >  1:235688/1‑150 (MQ=255)
       aaGCCGGTTGGACGCTGCCTGCGGGCAAAGACATCCGCGAGAAAAATGGTCAGCCGCTGCGCATTGAACTTTCGTTCATCGGCACCGATGCGTTAAGCAAATCGATGGCGGAAATCATTCAGGCTGATATGCGCCAGATTGGCGCAGATg                                                                                                                                        <  1:990531/150‑1 (MQ=255)
       aaGCCGGTTGGACGCTGCCTGAGGGGAAAGACATCCGCGAGAAAAATGGTCAGCCGCTGCGCATTGAACTTTCGTTCATCGGCACCGATGCGTTAAGCAAATCGATGGCGGAAATCATTCAGGCTGATATGCGCGAATCTGGCGCatatc                                                                                                                                        <  2:1030177/150‑2 (MQ=255)
         gCCGGTTGGACGCTGCCTGCGGGCAAAGACATCCGCGAGAAAAATGGTCAGCCGCTGCGCATTGAACTTTCGTTCATCGGCACCGATGCGTTAAGCAAATCGATGGCGGAAATCATTCAGGCTGATATGCGCCAATCTGGCGCatatcag                                                                                                                                      <  1:513377/150‑4 (MQ=255)
            ggTTGGACGCTGCCTGCGGGCAAAGACATCCGCGAGAAAAATGGTCAGCCGCTGCGCATTGAACTTTCGTTCATCGGCACCGATGCGTTAAGCAAATCGATGGCGGAAATCATTCAGGCTGATATGCGCCAGATTGGCGCAGATGTCTCg                                                                                                                                   >  2:1395840/1‑150 (MQ=255)
              ttGGACGCTGCCTGCGGGCAAAGACATCCGCGAGAAAAATGGTCAGCCGCTGCGCATTGAACTTTCGTTCATCGGCACCGATGCGTTAAGCAAATCGATGGCGGAAATCATTCAGGCTGATATGCGCCAGATTGGCGCAGATGTCTCGCt                                                                                                                                 >  1:486025/1‑150 (MQ=255)
               tGGACGCTGCCTGCGGGCAAAGACATCCGCGAGAAAAATGGTCAGCCGCTGCGCATTGAACTTTCGTTCATCGGCACCGATGCGTTAAGCAAATCGATGGCGGAAATCATTCAGGCTGATATGCGCCAGATTGGCGCAGATGTCTCGCTg                                                                                                                                <  1:188677/150‑1 (MQ=255)
                ggACGCTGCCTGCGGGCAAAGACATCCGCGAGAAAAATCGTCAGCCGCTGCGCATTGAACTTTCGTTCATAGGCACCGATGCGGTAAGCAAATCGAAGGAGGAAATCATTCAGGCTGATATGCGCCAGATTGGCGCAGATGTCTCGCTGa                                                                                                                               >  2:1349840/1‑150 (MQ=255)
                    gCTGCCTGCGGGCAAAGACATCCGCGAGAAAAATGGTCAGCCGCTGCGCATTGAACTTTCGTTCATCGGCACCGATGCGTTAAGCAAATCGATGGCGGAAATCATTCAGGCTGATATGCGCCAGATTGGCGCAGATGTCTCGCTGATTgg                                                                                                                           <  2:1099117/150‑1 (MQ=255)
                    gCTGCCTGCGGGCAAAGACATCCGCGAGAAAAATGGTCAGCCGCTGCGCATTGAACTTTCGTTCATCGGCACCGATGCGTTAAGCAAATCGATGGCGGAAATCATTCAGGCTGATATGCGCCAGATTGGCGCAGATGTCTCGCTGATTgg                                                                                                                           <  1:533659/150‑1 (MQ=255)
                     ctgcctgcGGGCAAAGACATCCGCGAGAAAAATGGTCAGCCGCTGCGCATTGAACTTTCGTTCATCGGCACCGATGCGTTAAGCAAATCGATGGCGGAAATCATTCAGGCTGATATGCGCCAGATTGGCGCAGATGTCTCGCTGATTGGc                                                                                                                          <  1:333883/150‑1 (MQ=255)
                       gcctgcGGGCAAAGACATCCGCGAGAAAAATGGTCAGCCGCTGCGCATTGAACTTTCGTTCATCGGCACCGATGCGTTAAGCAAATCGATGGCGGAAATCATTCAGGCTGATATGCGCCAGATTGGCGCAGATGTCTCGCTGATTGGCga                                                                                                                        >  1:293618/1‑150 (MQ=255)
                         ctgcGGGCAAAGACATCCGCGAGAAAAATGGTCAGCCGCTGCGCATTGAACTTTCGTTCATCGGCACCGATGCGTTAAGCAAATCGATGGCGGAAATCATTCAGGCTGATATGCGCCAGATTGGCGCAGATGTCTCGCTGATTGGCgaag                                                                                                                      >  2:1040201/1‑150 (MQ=255)
                          tgcGGGCAAAGACATCCGCGAGAAAAATGGTCAGCCGCTGCGCATTGAACTTTCGTTCATCGGCACCGATGCGTTAAGCAAATCGATGGCGGAAATCATTCAGGCTGATATGCGCCAGATTGGCGCAGATGTCTCGCTGATTGGCgaaga                                                                                                                     <  1:1208462/150‑1 (MQ=255)
                               gCAAAGACATCCGCGAGAAAAATGGTCAGCCGCTGCGCATTGAACTTTCGTTCATCGGCACCGATGCGTTAAGCAAATCGATGGCGGAAATCATTCAGGCTGATATGCGCCAGATTGGCGCAGATGTCTCGCTGATTGGCGAAGAagaga                                                                                                                <  1:956672/150‑1 (MQ=255)
                                 aaaGACATCCGCGAGAAAAATGGTCAGCCGCTGCGCATTGAACTTTCGTTCATCGGCACCGATGCGTTAAGCAAATCGATGGCGGAAATCATTCAGGCTGATATGCGCCAGATTGGCGCAGATGTCTCGCTGATTGGCGAAGAAGAGAGc                                                                                                              <  2:440506/150‑1 (MQ=255)
                                 aaaGACATCCGCGAGAAAAATGGTCAGCCGCTGCGCATTGAACTTTCGTTCATCGGCACCGATGCGTTAAGCAAATCGATGGCGGAAATCATTCAGGCTGATATGCGCCAGATTGGCGCAGATGTCTCGCTGATTGGCGAAGAAGAGAGc                                                                                                              <  2:1261667/150‑1 (MQ=255)
                                    gACATCCGCGAGAAAAATGGTCAGCCGCTGCGCATTGAACTTTCGTTCATCGGCACCGATGCGTTAAGCAAATCGATGGCGGAAATCATTCAGGCTGATATGCGCCAGATTGGCGCAGATGTCTCGCTGATTGGCGAAGAAGAGAGCAGt                                                                                                           <  2:713388/150‑1 (MQ=255)
                                    gACATCCGCGAGAAAAATGGTCAGCCGCTGCGCATTGAACTTTCGTTCATCGGCACCGATGCGTTAAGCAAATCGATGGCGGAAATCATTCAGGCTGATATGCGCCAGATTGGCGCAGATGTCTCGCTGATTGGCGAAGAAGAGAGCAGt                                                                                                           >  1:261692/1‑150 (MQ=255)
                                        tCCGCGAGAAAAATGGTCAGCCGCTGCGCATTGAACTTTCGTTCATCGGCACCGATGCGTTAAGCAAATCGATGGCGGAAATCATTCAGGCTGATATGCGCCAGATTGGCGCAGATGTCTCGCTGATTGGCGAAGAAGAGAGCAGTATCt                                                                                                       >  1:1381232/1‑150 (MQ=255)
                                        tCCGCGAGAAAAATGGTCAGCCGCTGCGCATTGAACTTTCGTTCATCGGCACCGATGCGTTAAGCAAATCGATGGCGGAAATCATTCAGGCTGATATGCGCCAGATTGGCGCAGATGTCTCGCTGATTGGCGAAGAAGAGAGCAGTATCt                                                                                                       >  2:152722/1‑150 (MQ=255)
                                         ccGCGAGAAAAATGGTCAGCCGCTGCGCATTGAACTTTCGTTCATCGGCACCGATGCGTTAAGCAAATCGATGGCGGAAATCATTCAGGCTGATATGCGCCAGATTGGCGCAGATGTCTCGCTGATTGGCGAAGAAGAGAGCAGTATCTa                                                                                                      <  1:1395840/150‑1 (MQ=255)
                                          cgcgAGAAAAATGGTCAGCCGCTGCGCATTGAACTTTCGTTCATCGGCACCGATGCGTTAAGCAAATCGATGGCGGAAATCATTCAGGCTGATATGCGCCAGATTGGCGCAGATGTCTCGCTGATTGGCGAAGAAGAGAGCAGTATCTAt                                                                                                     <  1:1112988/150‑1 (MQ=255)
                                           gcgAGAAAAATGGTCAGCCGCTGCGCATTGAACTTTCGTTCATCGGCACCGATGCGTTAAGCAAATCGATGGCGGAAATCATTCAGGCTGATATGCGCCAGATTGGCGCAGATGTCTCGCTGATTGGCGAAGAAGAGAGCAGTATCTATg                                                                                                    >  2:1287908/1‑150 (MQ=255)
                                             gagaAAAATGGTCAGCCGCTGCGCATTGAACTTTCGTTCATCGGCACCGATGCGTTAAGCAAATCGATGGCGGAAATCATTCAGGCTGATATGCGCCAGATTGGCGCAGATGTCTCGCTGATTGGCGAAGAAGAGAGCAGTATCTATGCt                                                                                                  <  1:538499/150‑1 (MQ=255)
                                             gagaAAAATGGTCAGCCGCTGCGCATTGAACTTTCGTTCATCGGCACCGATGCGTTAAGCAAATCGATGGCGGAAATCATTCAGGCTGATATGCGCCAGATTGGCGCAGATGTCTCGCTGATTGGCGAAGAAGAGAGCAGTATCTATGCt                                                                                                  >  2:411874/1‑150 (MQ=255)
                                              agaAAAATGGTCAGCCGCTGCGCATTGAACTTTCGTTCATCGGCACCGATGCGTTAAGCAAATCGATGGCGGAAATCATTCAGGCTGATATGCGCCAGATTGGCGCAGATGTCTCGCTGATTGGCGAAGAAGAGAGCAGTATCTATGCTc                                                                                                 <  1:60145/150‑1 (MQ=255)
                                                 aaaaTGGTCAGCCGCTGCGCATTGAACTTTCGTTCATCGGCACCGATGCGTTAAGCAAATCGATGGCGGAAATCATTCAGGCTGATATGCGCCAGATTGGCGCAGATGTCTCGCTGATTGGCGAAGAAGAGAGCAGTATCTATGCTCGTc                                                                                              <  2:770631/150‑1 (MQ=255)
                                                   aaTGGTCAGCCGCTGCGCATTGAACTTTCGTTCATCGGCACCGATGCGTTAAGCAAATCGATGGCGGAAATCATTCAGGCTGATATGCGCCAGATTGGCGCAGATGTCTCGCTGATTGGCGAAGAAGAGAGCAGTATCTATGCTCGTCAg                                                                                            <  1:744281/150‑1 (MQ=255)
                                                         cagcaacTGCTCATTTAACTTTAGTTCATCGGCAACGATGCGTTAAGCAAATAGATGACGGAAATCATTCAGGCTGATATGCGCCAGATTGGCGCAGATTTCTCGATTATTGGCGAAGAAGAGAGCAGTATCTATGCTCGTCAGCGCGAc                                                                                      >  2:1145080/7‑150 (MQ=255)
                                                           gCCGCTGCGCATTGAACTTTCGTTCATCGGCACCGATGCGTTAAGCAAATCGATGGCGGAAATCATTCAGGCTGATATGCGCCAGATTGGCGCAGATGTCTCGCTGATTGGCGAAGAAGAGAGCAGTATCTATGCTCGTCAGCGCGACgg                                                                                    >  1:682661/1‑150 (MQ=255)
                                                                tGCGCATTGAACTTTCGTTCATCGGCACCGATGCGTTAAGCAAATCGATGGCGGAAATCATTCAGGCTGATATGCGCCAGATTGGCGCAGATGTCTCGCTGATTGGCGAAGAAGAGAGCAGTATCTATGCTCGTCAGCGCGACGGTCGtt                                                                               >  2:522244/1‑150 (MQ=255)
                                                                   gcATTGAACTTTCGTTCATCGGCACCGATGCGTTAAGCAAATCGATGGCGGAAATCATTCAGGCTGATATGCGCCAGATTGGCGCAGATGTCTCGCTGATTGGCGAAGAAGAGAGCAGTATCTATGCTCGTCAGCGCGACGGTCGTTTTg                                                                            >  2:124534/1‑150 (MQ=255)
                                                                           cTTTCGTTCATCGGCACCGATGCGTTAAGCAAATCGATGGCGGAAATCATTCAGGCTGATATGCGCCAGATTGGCGCAGATGTCTCGCTGATTGGCGAAGAAGAGAGCAGTATCTATGCTCGTCAGCGCGACGGTCGTTTTGGCATGAtt                                                                    >  1:1328543/1‑150 (MQ=255)
                                                                                 ttCATCGGCACCGATGCGTTAAGCAAATCGATGGCGGAAATCATTCAGGCTGATATGCGCCAGATTGGCGCAGATGTCTCGCTGATTGGCGAAGAAGAGAGCAGTATCTATGCTCGTCAGCGCGACGGTCGTTTTGGCATGATTTTCCAc                                                              <  2:682661/150‑1 (MQ=255)
                                                                                                   ttAAGCAAATCGATGGCGGAAATCATTCAGGCTGATATGCGCCAGATTGGCGCAGATGTCTCGCTGATTGGCGAAGAAGAGAGCAGTATCTATGCTCGTCAGCGCGACGGTCGTTTTGGCATGATTTTCCACCGCACCTGGGGCGCGCCa                                            <  2:56938/150‑1 (MQ=255)
                                                                                                    tAAGCAAATCGATGGCGGAAATCATTCAGGCTGATATGCGCCAGATTGGCGCAGATGTCTCGCTGATTGGCGAAGAAGAGAGCAGTATCTATGCTCGTCAGCGCGACGGTCGTTTTGGCATGATTTTCCACCGCACCTGGGGCGCGCCat                                           <  1:1287908/150‑1 (MQ=255)
                                                                                                      aGCAAATCGATGGCGGAAATCATTCAGGCTGATATGCGCCAGATTGGCGCAGATGTCTCGCTGATTGGCGAAGAAGAGAGCAGTATCTATGCTCGTCAGCGCGACGGTCGTTTTGGCATGATTTTCCACCGCACCTGGGGCGCGCCatat                                         >  1:44552/1‑150 (MQ=255)
                                                                                                       gCAAATCGATGGCGGAAATCATTCAGGCTGATATGCGCCAGATTGGCGCAGATGTCTCGCTGATTGGCGAAGAAGAGAGCAGTATCTATGCTCGTCAGCGCGACGGTCGTTTTGGCATGATTTTCCACCGCACCTGGGGCGCGCCATATg                                        >  2:50641/1‑150 (MQ=255)
                                                                                                          aaTCGATGGCGGAAATCATTCAGGCTGATATGCGCCAGATTGGCGCAGATGTCTCGCTGATTGGCGAAGAAGAGAGCAGTATCTATGCTCGTCAGCGCGACGGACGTTTTGGCATGATTTTCCACCGAACCTGGGGCGCGCCATATGATc                                     >  1:1138799/1‑150 (MQ=255)
                                                                                                            tCGATGGCGGAAATCATTCAGGCTGATATGCGCCAGATTGGCGCAGATGTCTCGCTGATTGGCGAAGAAGAGAGCAGTATCTATGCTCGTCAGCGCGACGGTCGTTTTGGCATGATTTTCCACCGCACCTGGGGCGCGCCATATGATCca                                   <  1:1074324/150‑1 (MQ=255)
                                                                                                             cGATGGCGGAAATCATTCAGGCTGATATGCGCCAGATTGGCGCAGATGTCTCGCTGATTGGCGAAGAAGAGAGCAGTATCTATGCTCGTCAGCGCGACGGTCGTTTTGGCATGATTTTCCACCGCACCTGGGGCGCGCCATATGATCcac                                  <  1:1040201/150‑1 (MQ=255)
                                                                                                                 ggCGGAAATCATTCAGGCTGATATGCGCCAGATTGGCGCAGATGTCTCGCTGATTGGCGAAGAAGAGAGCAGTATCTATGCTCGTCAGCGCGACGGTCGTTTTGGCATGATTTTCCACCGCACCTGGGGCGCGCCATATGATCCACACGc                              >  2:684057/1‑150 (MQ=255)
                                                                                                                      aaaTCATTCAGGCTGATATGCGCCAGATTGGCGCAGATGTCTCGCTGATTGGCGAAGAAGAGAGCAGTATCTATGCTCGTCAGCGCGACGGTCGTTTTGGCATGATTTTCCACCGCACCTGGGGCGCGCCATATGATCCACACGCCTTcc                         <  2:519310/150‑1 (MQ=255)
                                                                                                                            ttCAGGCTGATATGCGCCAGATTGGCGCAGATGTCTCGCTGATTGGCGAAGAAGAGAGCAGTATCTATGCTCGTCAGCGCGACGGTCGTTTTGGCATGATTTTCCACCGCACCTGGGGCGCGCCATATGATCCACACGCCTTCCTCAGtt                   >  2:396914/1‑150 (MQ=255)
                                                                                                                             tCAGGCTGATATGCGCCAGATTGGCGCAGATGTCTCGCTGATTGGCGAAGAAGAGAGCAGTATCTATGCTCGTCAGCGCGACGGTCGTTTTGGCATGATTTTCCACCGCACCTGGGGCGCGCCATATGATCCACACGCCTTCCTCAGTTc                  <  2:235688/150‑1 (MQ=255)
                                                                                                                             tCAGGCTGATATGCGCCAGATTGGCGCAGATGTCTCGCTGATTGGCGAAGAAGAGAGCAGTATCTATGCTCGTCAGCGCGACGGTCGTTTTGGCATGATTTTCCACCGCACCTGGGGCGCGCCATATGATCCACACGCCTTCCTCAGTTc                  <  2:843124/150‑1 (MQ=255)
                                                                                                                               aGGCTGATATGCGCCAGATTGGCGCAGATGTCTCGCTGATTGGCGAAGAAGAGAGCAGTATCTATGCTCGTCAGCGCGACGGTCGTTTTGGCATGATTTTCCACCGCACCTGGGGCGCGCGATATGATCCACACGCCTTCCTCAGTTCaa                <  2:293618/150‑1 (MQ=255)
                                                                                                                                  cTGATATGCGCCAGATTGGCGCAGATGTCTCGCTGATTGGCGAAGAAGAGAGCAGTATCTATGCTCGTCAGCGCGACGGTCGTTTTGGCATGATTTTCCACCGCACCTGGGGCGCGCCATATGATCCACACGCCTTCCTCAGTTCAATgc             >  1:692126/1‑150 (MQ=255)
                                                                                                                                  cTGATATGCGCCAGATTGGCGCAGATGTCTCGCTGATTGGCGAAGAAGAGAGCAGTATCTATGCTCGTCAGCGCGACGGTCGTTTTGGCATGATTTTCCACCGCACCTGGGGCGCGCCATATGATCCACACGCCTTCCTCAGTTCAATgc             >  2:868575/1‑150 (MQ=255)
                                                                                                                                    catctGCGCCAATCTGGCGCAGATGTCTCGCTGATTGGCGAAGAAGAGAGCAGTATCTATGCTCGTCAGCGCGACGGTCGTTTTGGCATGATTTTCCACCGCACCTGGGGCGCGCCATATGATCCACACGCCTTCCTCAGTTCAATgcgc           >  1:160301/5‑150 (MQ=255)
                                                                                                                                     atctGCGCCAATCTGGCGCAGATGTCTCGCTGATTGGCGAAGAAGAGAGCAGTATCTATGCTCGTCAGCGCGACGGTCGTTTTGGCATGATTTTCCACCGCACCTGGGGCGCGCCATATGATCCACACGCCTTCCTCAGTTCAATGcgcg          >  1:321422/4‑150 (MQ=255)
                                                                                                                                     atctGCGCCAATCTGGCGCAGATGTCTCGCTGATTGGCGAAGAAGAGAGCAGTATCTATGCTCGTCAGCGCGACGGTCGTTTTGGCATGATTTTCCACCGCACCTGGGGCGCGCCATATGATCCACACGCCTTCCTCAGTTCAATGcgcg          >  1:626907/4‑150 (MQ=255)
                                                                                                                                     atatGCGCCAGATTGGCGCAGATGTCTCGCTGATTGGCGAAGAAGAGAGCAGTATCTATGCTCGTCAGCGCGACGGTCGTTTTGGCATGATTTTCCACCGCACCTGGGGCGCGCCATATGATCCACACGCCTTCCTCAGTTCAATGcgcg          <  1:684057/150‑1 (MQ=255)
                                                                                                                                      tatGCGCCAGATTGGCGCAGATGTCTCGCTGATTGGCGAAGAAGAGAGCAGTATCTATGCTCGTCAGCGCGACGGTCGTTTTGGCATGATTTTCCACCGCACCTGGGGCGCGCCATATGATCCACACGCCTTCCTCAGTTCAATGCGCGt         >  1:142501/1‑150 (MQ=255)
                                                                                                                                        tGCGCCAGAATGGCGCAGATCTCTCGATGATTGGAGAAGAAGACAGAAGTATCTATGCTCGTCAGCGCGACGGTCGTTTTGGCATGATTTTCCACCGCACCTGGGGCGCGCCATATGATCCACACGCCTTCCTCAGTTCAATGCGCGTAc       >  2:895897/1‑150 (MQ=255)
                                                                                                                                         gcgcCAGATTGGCGCAGATGTCTCGCTGATTGGCGAAGAAGAGAGCAGTATCTATGCTCGTCAGCGCGACGGTCGTTTTGGCATGATTTTCCACCGCACCTGGGGCGCGCCATATGATCCACACGCCTTCCTCAGTTCAATGCGCGTAcc      >  1:382840/1‑150 (MQ=255)
                                                                                                                                            ccAGATTGGCGCAGATGTCTCGCTGATTGGCGAAGAAGAGAGCAGTATCTATGCTCGTCAGCGCGACGGTCGTTTTGGCATGATTTTCCACCGCACCTGGGGCGCGCCATATGATCCACACGCCTTCCTCAGTTCAATgcg            <  1:835669/141‑1 (MQ=255)
                                                                                                                                            ccAGATTGGCGCAGATGTCTCGCTGATTGGCGAAGAAGAGAGCAGTATCTATGCTCGTCAGCGCGACGGTCGTTTTGGCATGATTTTCCACCGCACCTGGGGCGCGCCATATGATCCACACGCCTTCCTCAGTTCAATgcg            >  2:835669/1‑141 (MQ=255)
                                                                                                                                             cAGATTGGCGCAGATGTCTCGCTGATTGGCGAAGAAGAGAGCAGTATCTATGCTCGTCAGCGCGACGGTCGTTTTGGCATGATTTTCCACCGCACCTGGGGCGCGCCATATGATCCACACGCCTTCCTCAGTTCAATGCGCGTACCGTca  <  1:1312084/150‑1 (MQ=255)
                                                                                                                                               |                                                                                                                                                   
CTGGAAAAAGCCGGTTGGACGCTGCCTGCGGGCAAAGACATCCGCGAGAAAAATGGTCAGCCGCTGCGCATTGAACTTTCGTTCATCGGCACCGATGCGTTAAGCAAATCGATGGCGGAAATCATTCAGGCTGATATGCGCCAGATTGGCGCAGATGTCTCGCTGATTGGCGAAGAAGAGAGCAGTATCTATGCTCGTCAGCGCGACGGTCGTTTTGGCATGATTTTCCACCGCACCTGGGGCGCGCCATATGATCCACACGCCTTCCTCAGTTCAATGCGCGTACCGTCA  >  NC_000913/3614687‑3614977

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