Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,368,101 T→C A21A (GCT→GCC arnA → fused UDP‑L‑Ara4N formyltransferase/UDP‑GlcA C‑4'‑decarboxylase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,368,1010TC99.0% 721.6 / ‑4.6 194A21A (GCT→GCCarnAfused UDP‑L‑Ara4N formyltransferase/UDP‑GlcA C‑4'‑decarboxylase
Reads supporting (aligned to +/- strand):  ref base T (0/2);  new base C (112/80);  total (112/82)
Fisher's exact test for biased strand distribution p-value = 1.77e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 204 total aligned reads displayed.

AAGTTATCCGTCCATCCAGCAAGGAAAATGAATAATGAAAACCGTCGTTTTTGCCTACCACGATATGGGATGCCTCGGTATTGAAGCCCTGCTGGCTGCCGGTTACGAAATTAGCGCCATTTTTACCCATACTGATAATCCCGGTGAAAAAGCCTTTTATGGTTCGGTGGCTCGTCTGGCGGCGGAAAGA  >  NC_000913/2368005‑2368194
                                                                                                |                                                                                             
aaGTTATCCGTCCATCCAGCAAGGAAAATGAATAATGAAAACCGTCGTTTTTGCCTACCACGATATGGGATGCCTCGGTATTGAAGCCCTGCTGgccgcc                                                                                            >  1:1983961/1‑100 (MQ=255)
  gTTATCCGTCCATCCAGCAAGGAAAATGAATAATGAAAACCGTCGTTTTTGCCTACCACGATATGGGATGCCTCGGTATTGAAGCCCTGCTGGCCGCCgg                                                                                          <  1:4304661/100‑1 (MQ=255)
       ccGTCCATCCAGCAAGGAAAATGAATAATGAAAACCGTCGTTTTTGCCTACCACGATATGGGATGCCTCGGTATTGAAGCCCTGCTGGCCGCCGGTTACg                                                                                     >  1:4037613/1‑100 (MQ=255)
        cGTCCATCCAGCAAGGAAAATGAATAATGAAAACCGTCGTTTTTGCCTACCACGATATGGGATGCCTCGGAATTGAAGCCCTGCTGGCCGCCGGTTACGa                                                                                    >  1:3879830/1‑100 (MQ=255)
          tccatccaGCAAGGAAAATGAATAATGAAAACCGTCGTTTTTGCCTACCACGATATGGGATGCCTCGGTATTGAAGCCCTGCTGGCCGCCGGTTACGaaa                                                                                  <  1:767890/100‑1 (MQ=255)
            catccaGCAAGGAAAATGAATAATGAAAACCGTCGTTTTTGCCTACCACGATATGGGATGCCTCGGTATTGAAGCCCTGCTGGCCGCCGGTTACGAAAtt                                                                                >  1:558002/1‑100 (MQ=255)
            catccaGCAAGGAAAATGAATAATGAAAACCGTCGTTTTTGCCTACCACGATATGGGATGCCTCGGTATTGAAGCCCTGCTGGCCGCCGGTTACGAAAtt                                                                                >  1:710229/1‑100 (MQ=255)
            catccaGCAAGGAAAATGAATAATGAAAACCGTCGTTTTTGCCTACCACGATATGAGATGCCTCGGTATTGAAGCCCTGCTGGccgccg                                                                                           >  2:3483813/1‑89 (MQ=255)
             atccaGCAAGGAAAATGAATAATGAAAACCGTCGTTTTTGCCTACCACGATATGGGATGCCTCGGTATTGAAGCCCTGCTGGCCGCCGGTTACGAAATTa                                                                               >  1:3313287/1‑100 (MQ=255)
              tctaGCAAGGAAAATGAATAATGAAAACCGTCGTTTTTGCCTACCACGATATGGGATGCCTCGGTATTGAAGCCCTGCTGGCCGCCGGTTACGAAATTAg                                                                              <  1:4203112/97‑1 (MQ=255)
                caGCAAGGAAAATGAATAATGAAAACCGTCGTTTTTGCCTACCACGATATGGGATGCCTCGGTATTGAAGCCCTGCTGGCCGCCGGTTa                                                                                       >  2:2277452/1‑89 (MQ=255)
                caGCAAGGAAAATGAATAATGAAAACCGTCGTTTTTGCCTACCACGATATGGGATGCCTCGGTATTGAAGCCCTGCTGGCCGCCGGTTACGAAATTAgcg                                                                            >  1:3329012/1‑100 (MQ=255)
                 aGCAAGGAAAATGAATAATGAAAACCGTCGTTTTTGCCTACCACGATATGTGATGCCTCGGTATTGAAGCCCTGCTGGCCGCCGGTTACGAAATTAgcgc                                                                           >  1:4556259/1‑100 (MQ=255)
                   cAAGGAAAATGAATAATGAAAACCGTCGTTTTTGCCTACCACGATATGGGATGCCTCGGTATTGAAGCCCTGCTGGCCGCCGGTTACGAAATTAGCGCCa                                                                         >  1:2518958/1‑100 (MQ=255)
                    aaGGAAAATGAATAATGAAAACCGTCGTTTTTGCCTACCACGATATGGGATGCCTCGGTATTGAAGCCCTGCTGGCCGCCGGTTACGaa                                                                                   >  2:307800/1‑89 (MQ=255)
                     aGGAAAATGAATAATGAAAACCGTCGTTTTTGCCTGCCACGATATGGGATGTCTCGGTATTGAAGCCCTGCTGGCCGCCGGTTACGaaa                                                                                  <  2:2158399/89‑1 (MQ=255)
                     aGGAAAATGAATAATGAAAACCGTCGTTTTTGCCTACCACGATATGGGATGCCTCGGTATTGAAGCCCTGCTGGCCGCCGGTTACGAAATTAGCGCCAtt                                                                       <  1:2992757/100‑1 (MQ=255)
                        aaaaTGAATAATGAAAACCGTCGTTTTTTCCTATCTCGATATGGGATGCCTCGGTATTGAAGCCCTGCTGGCCGCCGGTTACGAAATTa                                                                               >  2:3547366/1‑89 (MQ=255)
                        aaaaTGAATAATGAAAACCGTCGTTTTTGCCTACCACGATATGGGATGCCTCGGTATTGAAGCCCTGCTGGCCGCCGGTTACGAAATTa                                                                               <  2:690213/89‑1 (MQ=255)
                        aaaaTGAATAATGAAAACCGTCGTTTTTGCCTACCACGATATGGGATGCCTCGGTATTGAAGCCCTGCTGGCCGCCGGTTACGAAATTAGCGCCAttttt                                                                    <  1:1121929/100‑1 (MQ=255)
                            tGAATAATGAAAACCGTCGTTTTTGCCTACCACGATATGGGATGCCTCGGTATTGAAGCCCTGCTGGCCGCCGGTTACGAAATTAgcgc                                                                           <  2:4652204/89‑1 (MQ=255)
                            tGAATAATGAAAACCGTCGTTTTTGCCTACCACGATATGGGATGCCTCGGTATTGAAGCCCTGCTGGCCGCCGGTTACGAAATTAGCGCCATTTTTAccc                                                                >  1:4521284/1‑100 (MQ=255)
                             gAATAATGAAAACCGTCGTTTTTGCCTACCACGATATGGGATGCCTCGGTATTGAAGCCCTGCTGGCCGCCGGTTACGAAATTAGCGcc                                                                          >  2:4167613/1‑89 (MQ=255)
                             gAATAATGAAAACCGTCGTTTTTGCCTACCACGATATGGGATGCCTCGGTATTGAAGCCCTGCTGGCCGCCGGTTACGAAATTAGCGCCATTTTTACCCa                                                               >  1:5045129/1‑100 (MQ=255)
                               ataatGAAAACCGTCGTTTTTGCCTACCACGATATGGGATGCCTCGGTATTGAAGCCCTGCTGGCCGCCGGTTACGATATTAGCGCCATTTTTACCCATa                                                             >  1:5166216/1‑100 (MQ=255)
                               ataatGAAAACCGTCGTTTTTGCCTACCACGATATGGGATGCCTCGGTATTGAAGCCCTGCTGGCCGCCGGTTACGAAATTAGCGCCATTTTTACCCATa                                                             >  1:1057270/1‑100 (MQ=255)
                                taatGAAAACCGTCGTTTTTGCCTACCACGATATGGGATGCCTCGGTATTGAAGCCCTGCTGGCCGCCGGTTACGAAATTAGCGCCAtt                                                                       >  2:1211063/1‑89 (MQ=255)
                                  atGAAAACCGTCGTTTTTGCCTATCACGATATGGGATGCCTCGGTATTGAAGTCCTGCTGGCCGCCTGTTACGAAATTATCGCCCtttt                                                                     >  2:4369709/1‑89 (MQ=255)
                                  atGAAAACCGTCGTTTTTGCCTACCACGATATGGGATGCCTCGGTATTGAAGCCCTGCTGGCCGCCGGTTACGAAATTAGCGCCAtttt                                                                     <  2:3173549/89‑1 (MQ=255)
                                     aaaaCCGTCGTTTTTGCCTACCACGATATGGGATGCCTCGGTATTGAAGCCCTGCTGGCCGCCGGTTACGAAATTAGCGCCATTTTTAc                                                                  <  2:4036140/89‑1 (MQ=255)
                                      aaaCCGTCGTTTTTGCCTATCACGATATGGGATGCCTCGGTATTGAAGCCCTGCTGGCCGCCGGTTACGAAATTAGCGCCATTTTTAcc                                                                 >  2:5508132/1‑89 (MQ=255)
                                      aaaCCGTCGTTTTTGCCTACCACGATATGGGATGCCTCGGCATTGAAGCCCTGCTGGCCGCCGGTTACGAAATTAGCGCCATTTTTAcc                                                                 <  2:2787081/89‑1 (MQ=255)
                                       aaCCGTCGTTTTTGCCTACCACGATATGGGATGCCTCGGTATTGAAGCCCTGCTGGCCGCCGGTTACGAAATTAGCGCCATTTTTAccc                                                                <  2:1001605/89‑1 (MQ=255)
                                           gtcgtTTTTGCCTACCACGATATGGGATGCCTCGGTATTGAAGCCCTGCTGGCCGCCGGTTACGAAATTAGCGCCATTTTTACCCATAc                                                            <  2:535884/89‑1 (MQ=255)
                                            tcgtTTTTGCCTACCACGATGTGGGATGCCTCGGTATTGAAGCCCTGCTGGCCGCCGGTTACGAAATTAGCGCCATTTTTACCCATACt                                                           <  2:5441214/89‑1 (MQ=255)
                                            tcgtTTTTGCCTACCACGATATGGGATGCCTCGGTATTGAAGCCCTGCTGGCCGCCGGTTACGAAATTAGCGCCATTTTTACCCATACt                                                           <  2:2431247/89‑1 (MQ=255)
                                              gtTTTTGCCTACCACGATATGGGATGCCTCGGTATTGAAGCCCTGCTGGCCGCCGGTTACGAAATTAGCGCCATTTTTACCCATACTGATAATCCCGGtn                                              <  1:3837136/100‑2 (MQ=255)
                                              gtTTTTGCCTACCACGATATGGGATGCCTCGGTATTGAAGCCCTGCTGGCCGCCGGTTACGAAATTAGCGCCATTTTTACCCATACTGATAATCCCGGTg                                              >  1:3695989/1‑100 (MQ=255)
                                               tttttGCCTACCACGATATGGGATGCCTCGGTATTGAAGCCCTGCTGGCCGCCGGTTACGAAATTAGCGCCATTTTTACCCATACTGAt                                                        <  2:3313287/89‑1 (MQ=255)
                                               tttttGCCTACCACGATATGGGATGCCTCGGTATTGAAGCCCTGCTGGCCGCCGGTTACGAAATTAGCGCCATTTTTACCCATACTGATAATCCCGGTGa                                             <  1:3209815/100‑1 (MQ=255)
                                               tttttGCCTACCACGATATGGGATGCCTCGGTATTGAAGCCCTGCTGGCCGCCGGTTACGAAATTAGCGCCATTTTTACCCATACTGATAATCCCGGTGa                                             >  1:4102260/1‑100 (MQ=255)
                                                 tttGCCTACCACGATATGGGATGCCTCGGTATTGAAGCCCTGCTGGCCGCCGGTTACGAAATTAGCGCCATTTTTACCCATACTGATaa                                                      <  2:2518958/89‑1 (MQ=255)
                                                     ccTACCACGATATGGGATGCCTCGGTATTGAAGCCCTGCTGGCCGCCGGTTACGAAATTAGCGCCATTTTTACCCa                                                               >  2:2477059/1‑76 (MQ=255)
                                                     ccTACCACGATATGGGATGCCTCGGTATTGAAGCCCTGCTGGCCGCCGGTTACGAAATTAGCGCCATTTTTACCCa                                                               <  1:2477059/76‑1 (MQ=255)
                                                     ccTACCACGATATGGGATGCCTCGGTATTGAAGCCCTGCTGGCCGCCGGTTACGAAATTAGCGCCATTTTTACCCATACTGATAATccc                                                  <  2:5131659/89‑1 (MQ=255)
                                                     ccTACCACGATATGGGATGCCTCGGTATTGAAGCCCTGCTGGCCGCCGGTTACGAAATTAGCGCCATTTTTACCCATACTGATAATCCCGGTGAAAAAGc                                       >  1:2197649/1‑100 (MQ=255)
                                                      cTACCACGATATGGGATGCCTCGGTATTGAAGCCCTGCTGGCCGCCGGTTACGAAATTAGCGCCATTTTTACCCATACTGATAATCCCg                                                 >  2:2330480/1‑89 (MQ=255)
                                                      cTACCACGATATGGGATGCCTCGGTATTGAAGCCCTGCTGGCCGCCGGTTACGAAATTAGCGCCATTTTTACCCATACTGATAATCCCGGTGAAAAAGcc                                      >  1:660049/1‑100 (MQ=255)
                                                       tACCACGATATGGGATGCCTCGGTATTGAAGCCCTGCTGGCCGCCGGTTACGAAATTAGCGCCATTTTTACCCATACTGATAATCCCgg                                                <  2:3193699/89‑1 (MQ=255)
                                                        aCCACGATATGGGATGCCTCGGTATTGAAGCCCTGCTGGCTGCCGGTTACGAAATTAGCGCCATTTTTACCCATACTGATAATCCCGGt                                               <  2:4780601/89‑1 (MQ=255)
                                                        aCCACGATATGGGATGCCTCGGTATTGAAGCCCTGCTGGCCGCCGGTTACGAAATTAGCGCCATTTTTACCCATACTGATAATCCCGGt                                               <  2:829372/89‑1 (MQ=255)
                                                         ccACGATATGGGATGCCTCGGTATTGAAGCCCTGCTGGCCGCCGGTTACGAAATTAGCGCCATTTTTACCCATACTGATAATCCCGGTGAAAAAGCCttt                                   >  1:4354035/1‑100 (MQ=255)
                                                         ccACGATATGGGATGCCTCGGTATTGAAGCCCTGCTGGCCGCCGGTTACGAAATTAGCGCCATTTTTACCCATACTGATAATCCCGGTGAAAAAGCCttt                                   >  1:439369/1‑100 (MQ=255)
                                                            cGATATGGGATGCCTCGGTATTGAAGCCCTGCTGGCCGCCGGTTACGAAATTAGCGCCATTTTTACCCATACTGATAATCCCGGTGaaa                                           >  2:222415/1‑89 (MQ=255)
                                                             gATATGGGATGCCTCGGTATTGAAGCCCTGCTGGCCGCCGGTTACGAAATTAGCGCCATTTTTACCCATACTGATAATCCCGGTGAAAAAGCCTTTTATg                               >  1:555649/1‑100 (MQ=255)
                                                              atatGGGATGCCTCGGTATTGAAGCCCTGCTGGCCGCCGGTTACGAAATTAGCGCCATTTTTACCCATACTGATAATCCCGGTGaaaaa                                         >  2:1532846/1‑89 (MQ=255)
                                                               tatGGGATGCCTCGGTATTGAAGCCCTGCTGGCCGCCGGTTACGAAATTAGCGCCATTTTTACCCATACTGATAATCCCGGTGAAAAAg                                        <  2:1340596/89‑1 (MQ=255)
                                                                atGGGATGCCTCGGTATTGAAGCCCTGCTGGCCGCCGGTTACGAAATTAGCGCCATTTTTACCCATACTGATAATCCCGGTGAAAAAGCCTTTTATGGtt                            >  1:1558596/1‑100 (MQ=255)
                                                                 tGGGATGCCTCGGTATTGAAGCCCTGCTGGCCGCCGGTTACGAAATTAGCGCCATTTTTACCCATACTGATAATCCCGGTGAAAAAGcc                                      <  2:3265071/89‑1 (MQ=255)
                                                                 tGGGATGCCTCGGTATTGAAGCCCTGCTGGCCGCCGGTTACGAAATTAGCGCCATTTTTACCCATACTGATAATACCGATGAAAAAGcc                                      <  2:2009311/89‑1 (MQ=255)
                                                                  gggATGCCTCGGTATTGAAGCCCTGCTGGCCGCCGGTTACGAAATTAGCGCCATTTTTACCCATACTGATAACCCCGGTGAAAAAGCCTTTTATGGTTCg                          <  1:5386512/100‑1 (MQ=255)
                                                                   ggATGCCTCGGTATTGAAGCCCTGCTGGCCGCCGGTTACGAAATTAGCGCCATTTTTACCCATACTGATAATCCCGGTGAAAAAGCCtt                                    >  2:519909/1‑89 (MQ=255)
                                                                   ggATGCCTCGGTATTGAAGCCCTGCTGGCCGCCGGTTACGAAATTAGCGCCATTTTTACCCATACTGATAATCCCGGTGAAAAAGCCtt                                    >  2:2975478/1‑89 (MQ=255)
                                                                   ggATGCCTCGGTATTGAAGCCCTGCTGGCCGCCGGTTACGAAATTAGCGCCATTTTTACCCATACTGATAATCCCGGTGAAAAAGCCtt                                    >  2:636160/1‑89 (MQ=255)
                                                                   ggATGCCTCGGTATTGAAGCCCTGCTGGCCGCCGGTTACGAAATTAGCGCCATTTTTACCCATACTGATAATCCCGGTGAAAAAGCCTTTTATGGTTCgg                         <  1:341639/100‑1 (MQ=255)
                                                                   ggATGCCTCGGTATTGAAGCCCTGCTGGCCGCCGGTTACGAAATTAGCGCCATTTTTACCCATACTGATAATCCCGGTGAAAAAGCCTTTTATGGTTCgg                         <  1:1888186/100‑1 (MQ=255)
                                                                    gATGCCTCGGTATTGAAGCCCTGCTGGCCGCCGGTTACGAAATTAGCGCCATTTTTACCCATACTGATAATCCCGGTGAAAAAGCCTTTTATGGTTCGGt                        <  1:2014281/100‑1 (MQ=255)
                                                                     aTGCCTCGGTATTGAAGCCCTGCTGGCCGCCGGTTACGAAATTAGCGCCATTTTTACCCATACTGATAATCCCGGTGAAAAAGCCTTTTATGGTTCGGTg                       >  1:426926/1‑100 (MQ=255)
                                                                      tGCCTCGGTATTGAAGCCCTGCTGGCTGCCGGTTACGAAATTAGCGCCATTTTTACCCATACTGATAATCCCGGTGAAAAAGCCTTTTATGGTTCGGTgg                      <  1:346568/100‑1 (MQ=255)
                                                                       gCCTCGGTATTGAAGCCCTGCTGGCCGCCGGTTACGAAATTAGCGCCATTTTTACCCATACTGATAATCCCGGTGAAAAAGCCTTTTATGGTTCGGTGGc                     <  1:5424305/100‑1 (MQ=255)
                                                                        ccTCGGTATTGAAGCCCTGCTGGCCGCCGGTTACGAAATTAGCGCCATTTTTACCCATACTGATAATCCCGGTGAAAAAGCCTTTTATGGTTCGGTGGCt                    >  1:1592524/1‑100 (MQ=255)
                                                                        ccTCGGTATTGAAGCCCTGCTGGCCGCCGGTTACGAAATTAGCGCCATTTTTACCCATACTGATAATCCCGGTGAAAAAGCCTTTTATGGTTCGGTGGCt                    >  1:4161524/1‑100 (MQ=255)
                                                                        ccTCGGTATTGAAGCCCTGCTGGCCGCCGGTTACGAAATTAGCGCCATTTTTACCCATACTGATAATCCCGGTGAAAAAGCCTTTTATGGTTCGGTGGCt                    >  1:2283743/1‑100 (MQ=255)
                                                                        ccTCGGTATTGAAGCCCTGCTGGCCGCCGGTTACGAAATTAGCGCCATTTTTACCCATACTGATAATCCCGGTGAAAAAGCCTTTTATGGTTCGGTGGCt                    >  1:1281831/1‑100 (MQ=255)
                                                                         cTCGGTATTGAAGCCCTGCTGGCCGCCGGTTACGAAATTAGCGCCATTTTTACCCATACTGATAATCCCGGTGAAAAAGCCTTTTATgg                              >  2:3115363/1‑89 (MQ=255)
                                                                            ggTATTGAAGCCCTGCTGGCCGCCGGTTACGAAATTAGCGCCATTTTTACCCATACTGATAATCCCGGTGAAAAAGCCTTTTATGGTTc                           <  2:1896224/89‑1 (MQ=255)
                                                                               aTTGAAGCCCTGCTGGCCGCCGGTTACGAAATTAGCGCCATTTTTACCCATACTGATAATCCCGGTGAAAAAGCCTTTTATGGTTCGGTGGCTCGTCTgg             <  1:2469716/100‑1 (MQ=255)
                                                                                ttGAAGCCCTGCTGGCCGCCGGTTACGAAATTAGCGCCATTTTTACCCATACTGATAATCCCGGTGAAAAAGCCTTTTATGGTTCGGTg                       >  2:53331/1‑89 (MQ=255)
                                                                                ttGAAGCCCTGCTGGCCGCCGGTTACGAAATTAGCGCCATTTTTACCCATACTGATAATCCCGGTGAAAAAGCCTTTTATGGTTCGGTg                       >  2:1136926/1‑89 (MQ=255)
                                                                                ttGAAGCCCTGCTGGCCGCCGGTTACGAAATTAGCGCCATTTTTACCCATACTGATAATCCCGGTGAAAAAGCCTTTTATGGTTCGGTg                       <  2:2667991/89‑1 (MQ=255)
                                                                                ttGAAGCCCTGCTGGCCGCCGGTTACGAAATTAGCGCCATTTTTACCCATACTGATAATCCCGGTGAAAAAGCCTTTTATGGTTCGGTg                       >  2:4061712/1‑89 (MQ=255)
                                                                                 tGAAGCCCTGCTGGCCGCCGGTTACGAAATTAGCGCCATTTTTACCCATACTGATAATCCCGGTGAAAAAGCCTTTTATGGTTCGGTGGCTCGTCTggcg           >  1:5106418/1‑100 (MQ=255)
                                                                                  gAAGCCCTGCTGGCCGCCGGTTACGAAATTAGCGCCATTTTTACCCATACTGATAATCCCGGTGAAAAAGCCTTTTATGGTTCGGTGGc                     >  2:703422/1‑89 (MQ=255)
                                                                                  gAAGCCCTGCTGGCCGCCGGTTACGAAATTAGCGCCATTTTTACCCATACTGATAATCCCGGTGAAAAAGCCTTTTATGGTTCGGTGGc                     >  2:3820766/1‑89 (MQ=255)
                                                                                  gAAGCCCTGCTGGCCGCCGGTTACGAAATTAGCGCCATTTTTACCCATACTGATAATCCCGGTGAAAAAGCCTTTTATGGTTCGGTGGCTCGTCTggcgg          >  1:2076858/1‑100 (MQ=255)
                                                                                  gAAGCCCTGCTGGCCGCCGGTTACGAAATTAGCGCCATTTTTACCCATACTGATAATCCCGGTGAAAAAGCCTTTTATGGTTCGGTGGCTCGTCTggcgg          >  1:4465546/1‑100 (MQ=255)
                                                                                    aGCCCTGCTGGCCGCCGGTTACGAAATTAGCGCCATTTTTACCCATACTGATAATCCCGGTGAAAAAGCCTTTTATGGTTCGGTGGCTc                   >  2:5271536/1‑89 (MQ=255)
                                                                                    aGCCCTGCTGGCCGCCGGTTACGAAATTAGCGCCATTTTTACCCATACTGATAATCCCGGTGAAAAAGCCTTTTATGGTTCGGTGGCTCGTCTGgcggcg        <  1:1808440/100‑1 (MQ=255)
                                                                                     gCCCTGCTGGCCGCCGGTTACGAAATTAGCGCCATTTTTACCCATACTGATAATCCCGGTGAAAAAGCCTTTTATGGTTCGGTGGCTCg                  >  2:1913793/1‑89 (MQ=255)
                                                                                     gCCCTGCTGGCCGCCGGTTACGAAATTAGCGCCATTTTTACCCATACTGATAATCCCGGTGAAAAAGCCTTTTATGGTTCGGTGGCTCGTCTGGcggcgg       >  1:2174559/1‑100 (MQ=255)
                                                                                      cccTGCTGGCCGCCGGTTACGAAATTAGCGCCATTTTTACCCATACTGATAATCCCGGTGAAAAAGCCTTTTATGGTTCGGTGGCTCGTCTGGCGGCGGa      >  1:5107514/1‑100 (MQ=255)
                                                                                       ccTGCTGGCCGCCGGTTACGAAATTAGCGCCATTTTTACCCATACTGATAATCCCGGTGAAAAAGCCTTTTATGGTTCGGTGGCTCGTCTGGCGGCGGaa     >  1:2615330/1‑100 (MQ=255)
                                                                                        ctgctgGCCGCCGGTTACGAAATTAGCGCCATTTTTACCCATACTGATAATCCCGGTGAAAAAGCCTTTTATGGTTCGGTGGCTCGTCt               <  2:4459766/89‑1 (MQ=255)
                                                                                        ctgctgGCCGCCGGTTACGAAATTAGCGCCATTTTTACCCATACTGATAATCCCGGTGAAAAAGCCTTTTATGGTTCGGTGGCTCGTCt               <  2:1504251/89‑1 (MQ=255)
                                                                                        ctgctgGCCGCCGGTTACGAAATTAGCGCCATTTTTACCCATACTGATAATCCCGGTGAAAAAGCCTTTTATGGTTCGGTGGCTCGTCt               <  2:1769932/89‑1 (MQ=255)
                                                                                        ctgctgGCCGCCGGTTACGAAATTAGCGCCATTTTTACCCATACTGATAATCCCGGTGAAAAAGCCTTTTATGGTTCGGTGGCTCGTCTGGCGGCGGaaa    >  1:5211416/1‑100 (MQ=255)
                                                                                        ctgctgGCCGCCGGTTACGAAATTAGCGCCATTTTTACCCATACTGATAATCCCGGTGAAAAAGCCTTTTATGGTTCGGTGGCTCGTCTGGCGGCGGaaa    >  1:2342057/1‑100 (MQ=255)
                                                                                           ctggccgCCGGTTACGAAATTAGCGCCATTTTTACCCATACTGATAATCCCGGTGAAAAAGCCTTTTATGGTTCGGTGGCTCGTCTggc            >  2:4870691/1‑89 (MQ=255)
                                                                                            tggccgCCGGTTACGAAATTAGCGCCATTTTTACCCATACTGATAATCCCGGTGAAAAAGCCTTTTATGGTTCGGTGGCTCGTCTggcg           <  2:1749079/89‑1 (MQ=255)
                                                                                             ggccgCCGGTTACGAAATTAGCGCCATTTTTACCCATACTGATAATCCCGGTGAAAAAGCCTTTTATGGTTCGGTGGCTCGTCTggcgg          >  2:429383/1‑89 (MQ=255)
                                                                                                |                                                                                             
AAGTTATCCGTCCATCCAGCAAGGAAAATGAATAATGAAAACCGTCGTTTTTGCCTACCACGATATGGGATGCCTCGGTATTGAAGCCCTGCTGGCTGCCGGTTACGAAATTAGCGCCATTTTTACCCATACTGATAATCCCGGTGAAAAAGCCTTTTATGGTTCGGTGGCTCGTCTGGCGGCGGAAAGA  >  NC_000913/2368005‑2368194

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