Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,298,174 A→C L29R (CTG→CGG)  napC ← quinol dehydrogenase, electron source for NapAB

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,298,1740AC99.4% 588.3 / NA 164L29R (CTG→CGG) napCquinol dehydrogenase, electron source for NapAB
Reads supporting (aligned to +/- strand):  ref base A (1/0);  new base C (79/84);  total (80/84)
Fisher's exact test for biased strand distribution p-value = 4.88e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 169 total aligned reads displayed.

TGGCAGCTAATGCAGAACTCTTCGGTATTGGCTTTTTCCATCCCGGTGTTAAAGCCACCCCAGAAGACGATGCCGCCAACAAAACCGATCAACAGCAGCGTCCCCAGCGCCAGACGGCTGGGGGTACGCCACCATTTCCACAGGCGCTTAATCAGACCAGGCTTACGGTCAGAATTTCCCATAATAACCTCT  >  NC_000913/2298078‑2298269
                                                                                                |                                                                                               
tGGCAGCTAATGCAGAACTCTTCGGTATTGGCTTTTTCCATCCCGGTGTTAAAGCCACCCCAGAAGACGATGCCGCCAACAAAACCGATCAACAGCcgcg                                                                                              <  1:37056/100‑1 (MQ=255)
tGGCAGCTAATGCAGAACTCTTCGGTATTGGCTTTTTCCATCCCGGTGTTAAAGCCACCCCAGAAGACGATGCCGCCAACAAAACCGATCAACAGCcgcg                                                                                              >  1:2624409/1‑100 (MQ=255)
 ggCAGCTAATGCAGAACTCTTCGGTATTGGCTTTTTCCATCCCGGTGTTCAAGCCACCCCAGAAGACGATGCCGCCACCAAAACCGATCAACAGCCGCGt                                                                                             >  2:870656/1‑100 (MQ=255)
 ggCAGCTAATGCAGAACTCTTCGGTATTGGCTTTTTCCATCCCGGTGTTAAAGCCACCCCAGAAGACGACGCCGCCAACAAAACCGATCAACAGCCGCGt                                                                                             >  1:2453697/1‑100 (MQ=255)
     gCTAATGCAGAACTCTTCGGTATTGGCTTTTTCCATCCCGGTGTTAAAGCCACCCCAGCAGCCGATGCCGCCAACAAAACCGATCAACAGCCGCGTcccc                                                                                         >  1:5149000/1‑100 (MQ=255)
      cTAATGCAGAACTCTTCGGTATTGGCTTTTTCCATCCCGGTGTTAAAGCCACCCCAGAAGCCGATGCCGCCAACAAAACCGATCAGCAGCCGCGTCCCCa                                                                                        <  2:2589885/100‑1 (MQ=255)
           gCAGAACTCTTCGGTATTGGCTTTTTTTATCTCTGTGTTAAAGCCACCCCAGAAGACGATGCCGCCAACAAAACCGATCAACAGCCGCGTCCCCAGCGcc                                                                                   >  2:1478997/1‑100 (MQ=255)
           gCAGAACTCTTCGGTATTGGCTTTTTTCTTTTTTTTGTTAAAGCCACCCCAGAAGACGATGCCGCCAACAAAACCGATCAACAGCCGCGTCCCCAGCGcc                                                                                   >  2:2779714/1‑100 (MQ=255)
           gCAGAACTCTTCGGTATTGGCTTTTTCCATCCCGGTGTTAAAGCCACCCCAGAAGACGATGCCGCCAACAAAACCGATCAACAGCCGCGTCCCCAGCGcc                                                                                   >  2:2377649/1‑100 (MQ=255)
           gCAGAACTCTTCGGTATTGGCTTTTTCCATCCCGGTGTTAAAGCCACCCCAGAAGACGATGCCGCCAACAAAACCGATCAACAGCCGCGTCCCCAGCGcc                                                                                   >  2:4823264/1‑100 (MQ=255)
            cAGAACTCTTCGGTATTGGCTTTTTCCATCCCGGTGTTAAAGCCACCCCAGAAGACGATGCCGCCAACAAAACCGATCAACAGCCGCGTCCCCAGCGCCa                                                                                  >  1:4112331/1‑100 (MQ=255)
             aGAACTCTTCGGTATTGGCTTTTTCCATCCCGGTGTTAAAGCCACCCCAGAAGACGATGCCGCCAACAAAACCGATCAACAGCCGCGTCCCCAGCGCCAg                                                                                 <  1:1616295/100‑1 (MQ=255)
             aGAACTCTTCGGTATTGGCTTTTTCCATCCCGGTGTTAAAGCCACCCCAGAAGACGATGCCGCCAACAAAACCGATCAACAGCCGCGTCCCCAGCGCCAg                                                                                 <  1:4489773/100‑1 (MQ=255)
             aGAACTCTTCGGTATTGGCTTTTTCCATCCCGGTGTTAAAGCCACCCCAGAAGACGATGCCGCCAACAAAACCGATCAACAGCCGCGTCCCCAGCGCCAg                                                                                 <  1:3772518/100‑1 (MQ=255)
             aGAACTCTTCGGTATTGGCTTTTTCCATCCCGGTGTTAAAGCCACCCCAGAAGACGATGCCGCCAACAAAACCGATCAACAGCCGCGTCCCCAGCGCCAg                                                                                 <  1:3503798/100‑1 (MQ=255)
               aaCTCTTCGGTATTGGCTTTTTCCATCCCGGTGTTAAAGCCACCCCAGAAGACGATGCCGCCAACAAAACCGATCAACAGCCGCGTCCCCAGCGCCAGAc                                                                               <  1:1789394/100‑1 (MQ=255)
               aaCTCTTCGGTATTGGCTTTTTCCATCCCGGTGTTAAAGCCACCCCAGAAGACGATGCCGCCAACAAAACCGATCAACAGCCGCGTCCCCAGCGCCAGAc                                                                               <  1:3346579/100‑1 (MQ=255)
                   ctTCGGTATTGGCTTTTTCCATCCCGGTGTTAAAGCCACCCCAGAAGACGATGCCGCCCACAAAACCGATCAACAGCCGCGTCCCCAGCGCCAGACGGCt                                                                           >  2:2158567/1‑100 (MQ=255)
                     tCGGTATTGGCTTTTTCCATCCCGGTGTTAAAGCCACCCCAGAAGACGATGCCGCCAACAAAACCGATCAACAGCCGCGTCCCCAGCGCCAGACGGCTgg                                                                         <  2:1050216/100‑1 (MQ=255)
                     tCGGTATTGGCTTTTTCCATCCCGGTGTTAAAGCCACCCCAGAAGACGATGCCGCCAACAAAACCGATCAACAGCCGCGTCCCCAGCGCCAGACGGCTgg                                                                         >  1:2652747/1‑100 (MQ=255)
                      cGGTATTGGCTTTTTCCATCCCGGTGTTAAAGCCACCCCAGAAGACGATGCCGCCAACAAAACCGATCAACAGCCGCGTCCCCAGCGCCAGACGGCTggg                                                                        <  1:3023218/100‑1 (MQ=255)
                      cGGTATTGGCTTTTTCCATCCCGGTGTTAAAGCCACCCCAGAAGACGATGCCGCCAACAAAACCGATCAACAGCCGCGTCCCCAGCGCCAGACGGCTggg                                                                        <  2:1517717/100‑1 (MQ=255)
                      cGGTATTGGCTTTTTCCATCCCGGTGTTAAAGCCACCCCAGAAGACGATGCCGCCAACAAAACCGATCAACAGCCGCGTCCCCAGCGCCAGACGGCTggg                                                                        <  2:1272589/100‑1 (MQ=255)
                       ggTATTGGCTTTTTCCATCCCGGTGTTAAAGCCACCCCAGAAGACGATGCCGCCAACAAAACCGATCAACAGCCGCGTCCCCAGCGCCAGACGGCTgggg                                                                       >  1:712067/1‑100 (MQ=255)
                       ggTATTGGCTTTTTCCATCCCGGTGTTAAAGCCACCCCAGAAGACGATGCCGCCAACAAAACCGATCAACAGCCGCGTCCCCAGCGCCAGACGGCTgggg                                                                       >  1:4787620/1‑100 (MQ=255)
                       ggTATTGGCTTTTTCCATCCCGGTGTTAAAGCCACCCCAGAAGACGATGCCGCCAACAAAACCGATCAACAGCCGCGTCCCCAGCGCCAGACGGCTgggg                                                                       >  1:4676875/1‑100 (MQ=255)
                       agTATTGGCTTTTTCCATCCCGGTGTTAAAGCCACCCCAGAAGACGATGCCGCCAACAAAACCGATCAACAGCCGCGTCCCCAGCGCCAGACGGCTgggg                                                                       <  1:4678175/99‑1 (MQ=255)
                        gTATTGGCTTTTTCCATCCCGGTGTTAAAGCCACCCCAGAAGACGATGCCGCCAACAAAACCGATCAACAGCCGCGTCCCCAGCGCCAGACGGCTggggg                                                                      >  2:4827839/1‑100 (MQ=255)
                          aTTGGCTTTTTCCATCCCGGTGTTAAAGCCACCCCAGAAGACGATGCCGCCCACAAAACAGATCAGCAGCAGCGTCCCCAGCGCCGGACGGATGGGGATa                                                                    >  1:5563899/1‑100 (MQ=255)
                          aTTGGCTTTTTCCATCCCGGTGTTAAAGCCACCCCAGAAGACGATGCCGCCAACAAAACCGATCAACAGCCGCGTCCCCAGCGCCAGACGGCTGGGGGTa                                                                    >  1:3612199/1‑100 (MQ=255)
                          aTTGGCTTTTTCCATCCCGGTGTTAAAGCCACCCCAGAAGACGATGCCGCCAACAAAACCGATCAACAGCCGCGTCCCCAGCGCCAGACGGCTGGGGGTa                                                                    <  2:2214183/100‑1 (MQ=255)
                          aTTGGCTTTTTCCATCCCGGTGTTAAAGCCACCCCAGAAGACGATGCCGCCAACAAAACCGATCAACAGCCGCGTCCCCAGCGCCAGACGGCTGGGGGTa                                                                    >  2:1179516/1‑100 (MQ=255)
                           ttGGCTTTTTCCATCCCGGTGTTAAAGCCACCCCAGAAGACGATGCCGCCAACAAAACCGATCAACAGCCGCGTCCCCAGCGCCAGACGGCTGGGGGTAc                                                                   <  2:4542755/100‑1 (MQ=255)
                           ttGGCTTTTTCCATCCCGGTGTTAAAGCCACCCCAGAAGACGATGCCGCCAACAAAACCGATCAACAGCCGCGTCCCCAGCGCCAGACGGCTGGGGGTAc                                                                   <  2:4428672/100‑1 (MQ=255)
                           ttGGCTTTTTCCATCCCGGTGTTAAAGCCACCCCAGAAGACGATGCCGCCAACAAAACCGATCAACAGCCGCGTCCCCAGCGCCAGACGGCTGGGGGTAc                                                                   <  1:2981275/100‑1 (MQ=255)
                           ttGGCTTTTTCCATCCCGGTGTTAAAGCCACCCCAGAAGACGATGCCGCCAACAAAACCGATCAACAGCCGCGTCCCCAGCGCCAGACGGCTGGGGGTAc                                                                   <  1:1352387/100‑1 (MQ=255)
                            tGGCTTTATCCATCCGGGAGTTAAGGCCACCCCAGAAGACGATGCCGCCAACAAAACCGATCAACAGCCGCGTCGCAAGCGGCAGACGGCTGGGGGTACg                                                                  <  2:1466533/100‑1 (MQ=255)
                              gCTTTTTCCATCCCGGTGTTAAAGCCACCCCAGAAGACGATGCCGCCAACAAAACCGATCAACAGCCGCGTCCCCAGCGCCAGACGGCTGGGGGTACGcc                                                                <  1:4493916/100‑1 (MQ=255)
                              gCTTTTTCCATCCCGGTGTTAAAGCCACCCCAGAAGACGATGCCGCCAACAAAACCGATCAACAGCCGCGTCCCCAGCGCCAGACGGCTGGGGGTACGcc                                                                <  1:4056640/100‑1 (MQ=255)
                              gCTTTTTCCATCCCGGTGTTAAAGCCACCCCAGAAGACGATGCCGCCAACAAAACCGATCAACAGCCGCGTCCCCAGCGCCAGACGGCTGGGGGTACGcc                                                                >  2:4706169/1‑100 (MQ=255)
                              gCTCTTTCCATCCCGGTGTTAGAGCCACCCCAGAAGACGATGCCGCCAACAAAACCGATCAACAGCCGAGTCCTCAGCGCCAGACGGCTGGGGGTACGcc                                                                <  2:1970024/100‑1 (MQ=255)
                                   ttCCATCCCGGTGTTAAAGCCACCCCAGAAGACGATGCCGCCAACAAAACCGATCAACAGCCGCGTCCCCAGCGCCAGACGGCTGGGGGTACGCCACCAt                                                           >  2:5004603/1‑100 (MQ=255)
                                   ttCCATCCCGGTGTTAAAGCCACCCCAGAAGACGATGCCGCCAACAAAACCGATCAACAGCCGCGTCCCCAGCGCCAGACGGCTGGGGGTACGCCACCAt                                                           >  2:2279545/1‑100 (MQ=255)
                                   ttCCATCCCGGTGGTAAAGCCACCCCAGAAGACGATGCCGCCAACAAAACCGATCAACAGCCGCGTCCCCAGCGCCAGACGGCTGGGGGTACGCCACCAt                                                           >  2:2831779/1‑100 (MQ=255)
                                    tCCATCCCGGTGTTAAAGCCACCCCAGAAGACGATGCCGCCAACAAAACCGATCAACAGCCGCGTCCCCAGCGCCAGACGGCTGGGGGTACGCCACCAtt                                                          >  1:4519767/1‑100 (MQ=255)
                                     ccATCCCGGTGTTAAAGCCACCCCAGAAGACGATGCCGCCAACAAAACCGATCAACAGCCGCGTCCCCAGCGCCAGACGGCTGGGGGTACGCCACCAttt                                                         <  1:3491110/100‑1 (MQ=255)
                                     ccATCCCGGTGTTAAAGCCACCCCAGAAGACGATGCCGCCAACAAAACCGATCAACAGCCGCGTCCCCAGCGCCAGACGGCTGGGGGTACGCCACCAttt                                                         <  1:3935243/100‑1 (MQ=255)
                                      cATCCCGGTGTTAAAGCCACCCCAGAAGACGATGCCGCCAACAAAACCGATCAACAGCCGCGTCCCCAGCGCCAGACGGCTGGGGGTACGCCACCATTTc                                                        >  2:4532547/1‑100 (MQ=255)
                                      cATCCCGGTGTTAAAGCCACCCCAGAAGACGATGCCGCCAACAAAACCGATCAACAGCCGCGTCCCCAGCGCCAGACGGCTGGGGGTACGCCACCATTTc                                                        >  2:4455382/1‑100 (MQ=255)
                                      cATCCCGGTGTTAAAGCCACCCCAGAAGACGATGCCGCCAACAAAACCGATCAACAGCCGCGTCCCCAGCGCCAGACGGCTGGGGGTACGCCACCATTTc                                                        <  1:3911684/100‑1 (MQ=255)
                                       aTCCCGGTGTTAAAGCCACCCCAGAAGACGATGCCGCCAACAAAACCGATCAACAGCCGCGTCCCCAGCGCCAGACGGCTGGGGGTACGCCACCATTTcc                                                       <  1:1247633/100‑1 (MQ=255)
                                       aTCCCGGTGTTAAAGCCACCCCAGAAGACGATGCCGCCAACAAAACCGATCAACAGCCGCGTCCCCAGCGCCAGACGGCTGGGGGTACGCCACCATTTcc                                                       <  2:2883002/100‑1 (MQ=255)
                                         cccGGTGTTAAAGCCACCCCAGAAGACGATGCCGCCAACAAAACCGATCAACAGCCGCGTCCCCAGCGCCAGACGGCTGGGGGTACGCCACCATTTCcac                                                     >  1:4864735/1‑100 (MQ=255)
                                         cccGGTGTTAAAGCCACCCCAGAAGACGATGCCGCCAACAAAACCGATCAACAGCCGCGTCCCCAGCGCCAGACGGCTGGGGGTACGCCACCATTTCcac                                                     >  2:3206355/1‑100 (MQ=255)
                                         cccGGTGTTAAAGCCACCCCAGAAGACGATGCCGCCAACAAAACCGATCAACAGCCGCGTCCCCAGCGCCAGACGGCTGGGGGTACGCCACCATTTCcac                                                     >  2:1766413/1‑100 (MQ=255)
                                          ccGGTGTTAAAGCCACCCCAGAAGACGATGCCGCCAACAAAACCGATCAACAGCCGCGTCCCCAGCGCCAGACGGCTGGGGGTACGCCACCATTTCcaca                                                    >  1:4745455/1‑100 (MQ=255)
                                            tctgtTAAAGCCACCCCAGAAGACGATGCCGCCAACAAAACCGATCAACAGCCGCGTCCCCAGCGCCAGACGGCTGGGGGTACGCCACCATTTCCACAgg                                                  <  1:2206282/98‑1 (MQ=255)
                                            ggTGTTAAAGCCACCCCAGAACACGATGCCGCCAACAAAACCGATCAACAGCCGCGTCCCCAGCGCCAGACGGCTGGGGGTACGCCACCATTTCCACAgg                                                  <  1:5190130/100‑1 (MQ=255)
                                             gtgtTAAAGCCACCCCAGAAGACGATGCCGCCAACAAAACCGATCAACAGCCGCGTCCCCAGCGCCAGACGGCTGGGGGTACGCCACCATTCCCACAGgc                                                 >  2:3198985/1‑100 (MQ=255)
                                              tgtTAAAGCCACCCCAGAAGACGATGCCGCCAACAAAACCGATCAACAGCCGCGTCCCCAGCGCCAGACGGCTGGGGGTACGCCACCATTTCCACAGgcg                                                >  1:1186392/1‑100 (MQ=255)
                                               gtTAAAGCCACCCCAGAAGACGATGCCGCCAACAAAACCGATCAACAGCCGCGTCCCCAGCGCCAGACGGCTGGGGGTACGCCACCATTTCCACAggaga                                               >  2:2183455/1‑97 (MQ=255)
                                                ttAAAGCCACCCCAGAAGACGATGCCGCCAACAAAACCGATCAACAGCCGCGTCCCCAGCGCCAGACGGCTGGGGGTACGCCACCATTTCCACAGGCGCt                                              <  1:1514036/100‑1 (MQ=255)
                                                 tAAAGCCACCCCAGAAGACGATGCCGCCAACAAAACCGATCAACAGCCGCGTCCCCAGCGCCAGACGGCTGGGGGTACGCCACCATTTCCACAGGCGCtt                                             >  2:1793520/1‑100 (MQ=255)
                                                 tAAAGCCACCCCAGAAGACGATGCCGCCAACAAAACCGATCAACAGCCGCGTCCCCAGCGCCAGACGGCTGGGGGTACGCCACCATTTCCACAGGCGCtt                                             <  2:3753958/100‑1 (MQ=255)
                                                    aGCCACCCCAGAAGACGATGCCGCCAACAAAACCGATCAACAGCCGCGTCCCCAGCGCCAGACGGCTGGGGGTACGCCACCATTTCCACAGGCGCTTAAt                                          >  1:5402956/1‑100 (MQ=255)
                                                     gCCACCCCAGAAGACGATGCCGCCAACAAAACCGATCAACAGCCGCGTCCCCAGCGCCAGACGGCTGGGGGTACGCCACCATTTCCACAGGCGCTTAATc                                         >  2:4012367/1‑100 (MQ=255)
                                                      ccACCCCAGAAGACGATGCCGCCAACAAAACCGATCAACAGCCGCGTCCCCAGCGCCAGACGGCTGGGGGTACGCCAGCATTTCCACAGGCGCTTAATCa                                        <  1:4659710/100‑1 (MQ=255)
                                                      ccACCCCAGAAGACGATGCCGCCAACAAAACCGATCAACAGCCGCGTCCCCAGCGCCAGACGGCTGGGGGTACGCCACCATTTCCACAGGCGCTTAATCa                                        <  1:5074815/100‑1 (MQ=255)
                                                      ccACCCCAGAAGACGATGCCGCCAACAAAACCGATCAACAGCCGCGTCCCCAGCGCCAGACGGCTGGGGGTACGCCACCATTTCCACAGGCGCTTAATCa                                        <  1:2303032/100‑1 (MQ=255)
                                                       cACCCCAGAAGACGATGCCGCCAACAAAACCGATCAACAGCCGCGTCCCCAGCGCCAGACGGCTGGGGGTACGCCACCATTTCCACAGGCGCTTAATCAg                                       <  2:4003092/100‑1 (MQ=255)
                                                       cACCCCAGAAGACGATGCCGCCAACAAAACCGATCAACAGCCGCGTCCCCAGCGCCAGACGGCTGGGGGTACGCCACCATTTCCACAGGCGCTTAATCAg                                       >  1:817252/1‑100 (MQ=255)
                                                         ccccAGAAGACGATGCCGCCAACAAAACCGATCAACAGCCGCGTCCCCAGCGCCAGACGGCTGGGGGTACGCCACCATTTCCACAGGCGCTTAATCAGAc                                     >  2:2172099/1‑100 (MQ=255)
                                                         ccccAGAAGACGATGCCGCCAACAAAACCGATCAACAGCCGCGTCCCCAGCGCCAGACGGCTGGGGGTACGCCACCATTTCCACAGGCGCTTAATCAGAc                                     >  1:3700445/1‑100 (MQ=255)
                                                          cccAGAGGACGATGCCGCCGACAAAACCGATCAACAGCCGCGTCCCCAGCGCCAGACTGCTGGGGGTACGCCACCATTTCCACAGGCGCTTAATCAGAcc                                    <  1:1119192/100‑1 (MQ=255)
                                                               tgtaCGCTGCCGCAAACAACACCGATCAACACCCTCGTTCCCAGCGCAAGACGGCTGTCGGTACGCCACCATTTCCACAGGCGCTTAAACAGACCAGGCt                               <  2:4492274/97‑1 (MQ=255)
                                                                        ccgccAACAAAACCGATCAACAGCCGCGTCCCCAGCGCCAGACGGCTGGGGGTACGCCACCATTTCCACAGGCGCTTAATCAGACCAGGCTTACGGTCAg                      >  2:2455068/1‑100 (MQ=255)
                                                                        ccgccAACAAAACCGATCAACAGCCGCGTCCCCAGCGCCAGACGGCTGGGGGTACGCCACCATTTCCACAGGCGCTTAATCAGACCAGGCTTACGGTCAg                      <  2:1667687/100‑1 (MQ=255)
                                                                          gccAACAAAACCGATCAACAGCCGCGTCCTCAGCGCCAGACGGCTGGGGGTACGCCACCATTTCCACAGGCGCTTAATCAGACCAGGCTTACGGTCAGaa                    <  2:1121844/100‑1 (MQ=255)
                                                                               caaAACCGATCAACAGCCGCGTCCCCAGCGCCAGACGGCTGGGGGTACGCCACCATTTCCACAGGCGCTTAATCAGACCAGGCTTACGGTCAGAATTTcc               >  2:3051869/1‑100 (MQ=255)
                                                                               caaAACCGATCAACAGCCGCGTCCCCAGCGCCAGACGGCTGGGGGTACGCCACCATTTCCACAGGCGCTTAATCAGACCAGGCTTACGGTCAGAATTTcc               >  2:1892494/1‑100 (MQ=255)
                                                                               caaAACCGATCAACAGCCGCGTCCCCAGCGCCAGACGGCTGGGGGTACGCCACCAATTCCACAGGCGCTTAATCAGACCAGGCTTACGGTCAGAATTTcc               <  1:2901723/100‑1 (MQ=255)
                                                                                aaaaCCGATCAACAGCCGCGTCCCCAGCGCCAGACGGCTGGGGGTACGCCACCATTTCCACAGGCGCTTAATCAGACCAGGCTTACGGTCAGAATTTccc              <  2:1007233/100‑1 (MQ=255)
                                                                                  aaCCGATCAACAGCCGCGTCCCCAGCGCCAGACGGCTTGGGGTACGCCACCATTTCCACAGGCGCTTAATCAGACCAGGCTTACGGTCAGAATTGCCCat            <  2:2233252/100‑1 (MQ=255)
                                                                                  aaCCGATCAACAGCCGCGTCCCCAGCGCCAGACGGCTGGGGGTACGCCACCATTTCCACAGGCGCTTAATCAGACCAGGCTTACGGTCAGAATTTCCCat            >  1:3191255/1‑100 (MQ=255)
                                                                                  aaCCGATCAACAGCCGCGTCCCCAGCGCCAGACGGCTGGGGGTACGCCACCATTTCCACAGGCGCTTAATCAGACCAGGCTTACGGTCAGAATTTCCCat            >  2:1862665/1‑100 (MQ=255)
                                                                                    ccGATCAACAGCCGCGTCCCCAGCGCCAGACGGCTGGGGGTACGCCACCATTTCCACAGGCGCTTAATCAGACCAGGCTTACGGTCAGAATTTCCCataa          >  2:3911804/1‑100 (MQ=255)
                                                                                     cGATCAACAGCCGCGTCCCCAGCGCCAGACGGCTGGGGGTACGCCACCATTTCCACAGGCGCTTAATCAGACCAGGCTTACGGTCAGAATTTCCCataat         <  2:4160449/100‑1 (MQ=255)
                                                                                     cGATCAACAGCCGCGTCCCCAGCGCCAGACGGCTGGGGGTACGCCACCATTTCCACAGGCGCTTAATCAGACCAGGCTTACGGTCAGAATTTCCCataat         <  1:1039101/100‑1 (MQ=255)
                                                                                     cGATCAACAGCCGCGTCCCCAGCGCCAGACGGCTGGGGGTACGCCACCATTTCCACAGGCGCTTAATCAGACCAGGCTTACGGTCAGAATTTCCCataat         <  2:332989/100‑1 (MQ=255)
                                                                                      gATCAACAGCCGCGTCCCCAGCGCCAGACGGCTGGGGGTACGCCACCATTTCCACAGGCGCTTAATCAGACCAGGCTTACGGTCAGAATTTCCCataata        >  1:4374635/1‑100 (MQ=255)
                                                                                        tCAACAGCCGCGTCCCCAGCGCCAGACGGCTGGGGGTACGCCACCATTTCCACAGGCGCTTAATCAGACCAGGCTTACGGTCAGAATTTCCCATAATAAc      >  1:5554460/1‑100 (MQ=255)
                                                                                         tctcagccgCGTCCCCAGCGCCAGACGGCTGGGGGTACGCCACCATTTCCACAGGCGCTTAATCAGACCAGGCTTACGGTCAGAATTTCCCATAATAAcc     <  1:4256420/97‑1 (MQ=255)
                                                                                         cAACAGCCGCGTCCCCAGCGCCAGACGGCTGGGGGTACGCCACCATTTCCACAGGCGCTTAATCAGACCAGGCTTACGGTCAGAATTTCCCATAATAAcc     <  2:3763689/100‑1 (MQ=255)
                                                                                         cAACAGCCGCGTCCCCAGCGCCAGACGGCTGGGGGTACGCCACCATTTCCACAGGCGCTTAATCAGACCAGGCTTACGGTCAGAATTTCCCATAATAAcc     <  2:4774132/100‑1 (MQ=255)
                                                                                          aaCAGCCGCGTCCCCAGCGCCAGACGGCTGGGGGTACGCCACCATTTCCACAGGTGCTTAATCAGACCAAGCTTACGGTCAGAATTTCCCATAATAACct    >  2:5175927/1‑100 (MQ=255)
                                                                                          aaCAGCCGCGTCCCCAGCGCCAGACGGCTGGGGGTACGCCACCATTTCCACAGGCGCTTAATCAGACCAGGCTTACGGTCAGAATTTCCCATAATAACct    >  1:5135706/1‑100 (MQ=255)
                                                                                          aaCAGCCGCGTCCCCAGCGCCAGACGGCTGGGGGTACGCCACCATTTCCACAGGCGCTTAATCAGACCAGGCTTACGGTCAGAATTTCCCATAATAACct    >  2:660588/1‑100 (MQ=255)
                                                                                           aCAGCCGCGTCCCCAGCGCCAGACGGCTGGGGGTACGCCACCATTTCCACAGGCGCTTAATCAGACCAGGCTTACGGTCAGAATTTCCCATAATAACctc   >  2:2735780/1‑100 (MQ=255)
                                                                                           aCAGCCGCGTCCCCAGCGCCAGACGGCTGGGGGTACGCCACCATTTCCACAGGCGCTTAATCAGACCAGGCTTACGGTCAGAATTTCCCATAATAACctc   >  1:1243943/1‑100 (MQ=255)
                                                                                            cagccgCGTCCCCAGCGCCAGACGGCTGGGGGTACGCCACCATTTCCACAGGCGCTTAATCAGACCAGGCTTACGGTCAGAATTTCCCATAATAAccaga  >  2:4217539/1‑97 (MQ=255)
                                                                                                |                                                                                               
TGGCAGCTAATGCAGAACTCTTCGGTATTGGCTTTTTCCATCCCGGTGTTAAAGCCACCCCAGAAGACGATGCCGCCAACAAAACCGATCAACAGCAGCGTCCCCAGCGCCAGACGGCTGGGGGTACGCCACCATTTCCACAGGCGCTTAATCAGACCAGGCTTACGGTCAGAATTTCCCATAATAACCTCT  >  NC_000913/2298078‑2298269

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