Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,991,539 T→C N45S (AAT→AGT) 
I26T (ATT→ACT) 
cyaY ←
b3808 →
hypothetical protein
hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,991,5390TC99.7% 1183.2 / NA 368N45S (AAT→AGT) 
I26T (ATT→ACT) 
cyaY
b3808
hypothetical protein
hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (0/1);  new base C (175/192);  total (175/193)
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

Only 100 of 387 total aligned reads displayed.

CCTTTCAGGTCAAAATGGTAGCCGCCCTGTTTGGTTGCCAGCCATACCTGGTGCAGCGGCTCCTGGCGGTTGATAATGATTTTGCTGCCATTCTCAAAGGTAATGGTCAGTACGCCGCCGTTGATTTCGCAGTCGATATCGCTGTCGCCATCCCAGTCGTCCAGGCGTTCTTCAATGGTCAG  >  NC_000913/3991449‑3991630
                                                                                          |                                                                                           
ccTTTCAGGTCAAAATGGTAGCCGCCCTGTTTGGTTGCCAGCCATACCTGGTGCAGCGGCTCCTGGCGGTTGATAATGATTTTGCTGCCACtc                                                                                           <  2:3780487/93‑1 (MQ=255)
  tttCAGGTCAAAATGGTAGCCGCCCTGTTTGGTTGCCAGCCATACCTGGTGCAGCGGCTCCTGGCGGTTGATAATGATTTTGCTGCCCCtctc                                                                                         >  2:1515399/1‑93 (MQ=255)
  tttCAGGTCAAAATGGTAGCCGCCCTGTTTGGTTGCCAGCCATACCTGGTGCAGCGGCTCCTGGCGGTTGATAATGATTTTGCTGCCACtctc                                                                                         >  1:5447959/1‑93 (MQ=255)
  tttCAGGTCAAAATGGTAGCCGCCCTGTTTGGTTGCCAGCCATACCTGGTGCAGCGGCTCCTGGCGGTTGATAATGATTTTGCTGCCACtctc                                                                                         >  1:5287817/1‑93 (MQ=255)
  tttCAGGTCAAAATGGTAGCCGCCCTGTTTGGTTGCCAGCCATACCTGGTGCAGCGGCTCCTGGCGGTTGATAATGATTTTGCTGCCACtctc                                                                                         >  2:6255969/1‑93 (MQ=255)
  tttCAGGTCAAAATGGTAGCCGCCCTGTTTGGTTGCCAGCCATACCTGGCGCAGCGGCTCCTGGCGGTTGATAATGATTTTGCTGCCACtctc                                                                                         >  2:74623/1‑93 (MQ=255)
   ttCAGGTCAAAATGGTAGCCGCCCTGTTTGGTTGCCAGCCATACCTGGTGCAGCGGCTCCTGGCGGTTGATAATGATTTTGCTGCCACTCTCa                                                                                        >  2:7430417/1‑93 (MQ=255)
         tCAAAATGGTAGCCGCCCTGTTTGGTTGCCAGCCATACCTGGTGCAGCGGCTCCTGGCGGTTGATAATGATTTTGCTGCCACTCTCAAAGGTa                                                                                  >  1:7635228/1‑93 (MQ=255)
           aaaaTGGTAGCCGCCCTGTTTGGTTGCCAGCCATACCTGGTGCAGCGGCTCCTGGCGGTTGATAATGATTTTGCTGCCACTCTCAAAGGTAAt                                                                                >  2:1777909/1‑93 (MQ=255)
           aaaaTGGTAGCCGCCCTGTTTGGTTGCCAGCCATACCTGGTGCAGCGGCTCCTGGCGGTTGATAATGATTTTGCTGCCACTCTCAAAGGTAAt                                                                                <  2:6078134/93‑1 (MQ=255)
           aaaaTGATAGCCGCCCTGTTGGGTTGCCAGCCATACCTGGTGCAGCGGCTCCTGGCGGTTGATAATGATTTTGCTGCCACTCTCAAAGGTAAt                                                                                <  2:1294978/93‑1 (MQ=255)
             aaTGGTAGCCGCCCTGTTTGGTTGCCAGCCATACCTGGTGCAGCGGCTCCTGGCGGTTGATAATGATTTTGCTGCCACTCTCAAAGGTAATgg                                                                              >  2:2580544/1‑93 (MQ=255)
             aaTGGTAGCCGCCCTGTTTGGTTGCCAGCCATACCTGGTGCAGCGGCTCCTGGCGGTTGATAATGATTTTGCTGCCACTCTCAAAGGTAATgg                                                                              <  1:7745358/93‑1 (MQ=255)
              aTGGTAGCCGCCCTGTTTGGTTGCCAGCCATACCTGGTGCAGCGGCTCCTGGCGGTTGATAATGATTTTGCTGCCACTCTCAAAGGTAATGGt                                                                             <  2:6413385/93‑1 (MQ=255)
              aTGGTAGCCGCCCTGTTTGGTTGCCAGCCATACCTGGTGCAGCGGCTCCTGGCGGTTGATAATGATTTTGCTGCCACTCTCAAAGGTAATGGt                                                                             <  1:4865257/93‑1 (MQ=255)
              aTGGTAGCCGCCCTGTTTGGTTGCCAGCCATACCTGGTGCAGCGGCTCCTGGCGGTTGATAATGATTTTGCTGCCACTCTCAAAGGTAATGGt                                                                             >  1:5857123/1‑93 (MQ=255)
              aTGGTAGCCGCCCTGTTTGGTTGCCAGCCATACCTGGTGCAGCGGCTCCTGGCGGTTGATAATGATTTTGCTGCCACTCTCAAAGGTAATGGt                                                                             >  1:1558259/1‑93 (MQ=255)
               tGGTAGCCGCCTTGTTTGGTTGCCAGCCATGCCTGGTGCAGCGGCTCCTGGCGGTTGATAATGATTTTGCTGCCACTCTCAAAGGTAATGGTc                                                                            <  2:3035491/93‑1 (MQ=255)
               tGGTAGCCGCCCTGTTTGGTTTCCAGGCATACCTGGTGCAGCTGCTCCTGGCGGTTGATTATGACTTTGCTGCCACTCTCTAAGTCTATGGTc                                                                            >  2:619244/1‑93 (MQ=255)
               tGGTAGCCGCCCTGTTTGGTTGCCAGCCATACCTGGTGCAGCGGCTCCTGGCGGTTGATAATGATTTTGCTGCCACTCTCAAAGGTAATGGTc                                                                            <  2:3968956/93‑1 (MQ=255)
               aggTAGCCGCCCTGTTTGGTTGCCAGCCATACCTGGTGCAGCGGCTCCTGGCGGTTGATAATGATTTTGCTGCCACTCTCAAAGGTAATGGTc                                                                            <  1:749719/92‑1 (MQ=255)
                ggTAGCCGCCCTGTTTGGTTGCCAGCCATACCTGGTGCAGCGGCTCCTGGCGGTTGATAATGATTTTGCTGCCACTCTCAAAGGTAATGGTCa                                                                           >  2:3356212/1‑93 (MQ=255)
                 gtcgccgccCTGTTTGGTTGCCAGCCATACCTGGTGCAGCGGCTCCTGGCGGTTGATAATGATTTTGCTGCCACTCTCAAAGGTAATGGTCAg                                                                          <  1:2613574/90‑1 (MQ=255)
                 gTAGCCGCCCTGTTTGGTTGCCAGCCATACCTGGCGCAGCGGCTCCTGGCGGTTGATAATGATTTTGCTGCCACTCTCAAAGGTAATGGTCAg                                                                          >  1:579998/1‑93 (MQ=255)
                   aGCCGCCCTGTTTGGTTGCCAGCCATACCTGGCGCAGCGGCTCCTGGCGGTTGATAATGATTTTGCTGCCACTCTCAAAGGTAATGGTCAGTa                                                                        <  1:7270995/93‑1 (MQ=255)
                    gccgccCTGTTGGGTTGCCAGCCATACCTGGTGCAGCGGCTCCTGGCGGTTGATAATGATTTTGCTGCCACTCTCAAAGGTAATGGTCAGTAc                                                                       <  1:1054560/93‑1 (MQ=255)
                      cgccCTGTTTGGTTGCCAGCCATACCTGGTGCAGCGGCTCCTGGCGGTTGATAATGATTTTGCTGCCACTCTCAAAGGTAATGGTCAGTAcgc                                                                     >  2:6678860/1‑93 (MQ=255)
                      cgccCTGTTTGGTTGCCAGCCATACCTGGCGCAGCGGCTCCTGGCGGTTGATAATGATTTTGCTGCCACTCTCAAAGGTAATGGTCAGTAcgc                                                                     <  1:5545768/93‑1 (MQ=255)
                       gccCTGTTTGGTTGCCAGCCATACCTGGTGCAGCGGCTCCTGGCGGTTGATAATGATTTTGCTGCCACTCTCAAAGGTAATGGTCAGTAcgcc                                                                    >  2:2546312/1‑93 (MQ=255)
                        cccTGTTTGGTTGCCAGCCATACCTGGTGCAGCGGCTCCTGGCGGTTGATAATGATTTTGCTGCCACTCTCAAAGGTAATGGTCAGTAcgccg                                                                   >  1:2665445/1‑93 (MQ=255)
                         cttGTTTGGTGGCCAGCCATACCTGGTGCAGCGGCTCCTGGCGGTTGATAATGATTTTGCTGCCACTCTCAAAGGTAATGGTCAGTAcgccgc                                                                  <  2:780180/91‑1 (MQ=255)
                          cTGTTTGGTTGCCAGCCATACCTGGTGCAGCGGCTCCTGGCGGTTGATAATGATTTTGCTGCCACTCTCAAAGGTAATGGTCAGTACgccgcc                                                                 >  2:6500580/1‑93 (MQ=255)
                           tGTTTGGTTGCCAGCCATACCTGGTGCAGCGGCTCCTGGCGGTTGATAATGATTTTGCTGCCACTCTCAAAGGTAATGGTCAGTACGccgccg                                                                >  1:1967113/1‑93 (MQ=255)
                           tGTTTGGTTGCCAGCCATACCTGGTGCAGCGGCTCCTGGCGGGTGATAATGATTTTGCTGGCACTCTCAAAGGTAATGGTCAGTACGccgccg                                                                >  2:603147/1‑93 (MQ=255)
                            gTTTGGTTGCCAGCCATACCTGGTGCAGCGGCTCCTGGCGGTTGATAATGATTTTGCTGCCACTCTCAAAGGTAATGGTCAGTACGCCGCCGt                                                               >  1:7351614/1‑93 (MQ=255)
                             tttGGTTGCCAGCCATACCTGGTGCAGCGGCTCCTGGCGGTTGATAATGATTTTGCTGCCACTCTCAAAGGTAATGGTCAGTACGCCGCCGtt                                                              >  2:720555/1‑93 (MQ=255)
                               tGGTTGCCAGCCATACCTGGTGCAGCGGCTCCTGGCGGTTGATAATGATTTTGCTGCCACTCTCAAAGGTAATGGTCAGTACGCCGCCGTTGa                                                            <  2:836532/93‑1 (MQ=255)
                               tGGTTGCCAGCCATACCTGGTGCAGCGGCTCCTGGCGGTTGATAATGATTTTGCTGCCACTCTCAAAGGTAATGGTCAGTACGCCGCCGTTGa                                                            <  1:24554/93‑1 (MQ=255)
                                 gTTGCCAGCCATACCTGGTGCAGCGGCTCCTGGCGGTTGATAATGATTTTGCTGCCACTCTCAAAGGTAATGGTCAGTACGCCGCCGTTGAtt                                                          >  2:3878497/1‑93 (MQ=255)
                                 gTTGCCAGCCATACCTGGCGCAGCGGCTCCTGGCGGTTGATAATGATTTTGCTGCCACTCTCAAAGGTAATGGTCAGTACGCCGCCGTTGAtt                                                          <  1:3982993/93‑1 (MQ=255)
                                  ttGCCAGCCATACCTGGTGCAGCGGCTCCTGGCGGTTGATAATGATTTTGCTGCCACTCTCAAAGGTAATGGTCAGTACGCCGCCGTTGAttt                                                         <  2:1389183/93‑1 (MQ=255)
                                   tGCCAGCCATACCTGGTGCAGCGGCTCCTGGCGGTTGATAATGATTTTGCTGCCACTCTCAAAGGTAATGGTCAGTACGCCGCCGTTGATTTc                                                        <  2:5648286/93‑1 (MQ=255)
                                   tGCCAGCCATACCTGGTGCAGCGGCTCCTGGCGGTTGATAATGATTTTGCTGCCACTCTCAAAGGTAATGGTCAGTACGCCGCCGTTGATTTc                                                        <  1:5694895/93‑1 (MQ=255)
                                     ccagccaTACCTGGTGCAGCGGCTCCTGGCGGTTGATAATGATTTTGCTGCCACTCTCAAAGGTAATGGTCAGTACGCCGCCGTTGATTTCGc                                                      >  2:269353/1‑93 (MQ=255)
                                      cagccaTACCTGGTGCAGCGGCTCCTGGCGGTTGATAATTATTTTGCTGCCACTCTCAAAGGTAATGGTCAGTACGCCGCCGTTGATTTCGCa                                                     >  1:4794465/1‑93 (MQ=255)
                                      cagccaTACCTGGTGCAGCGGCTCCTGGCGGTTGATAATGATTTTGCTGCCACTCTCAAAGGTAATGGTCAGTACGCCGCCGTTGATTTCGCa                                                     <  2:4970061/93‑1 (MQ=255)
                                       agccaTACCTGGTGCAGCGGCTCCTGGCGGTTGATAATGATTTTGCTGCCACTCTCAAAGGTAATGGTCAGTACGCCGCCGTTGATTTCGCAg                                                    <  1:6962333/93‑1 (MQ=255)
                                         ccaTACCTGGTGCAGCGGCTCCTGGCGGTTGATAATGATTTTGCTGCCACTCTCAAAGGTAATGGTCAGTACGCCGCCGTTGATTTCGCAGTc                                                  <  2:412664/93‑1 (MQ=255)
                                          caTACCTGGTGCAGCGGCTCCTGGCGGTTGATAATGATTTTGCTGCCACTCTCAAAGGTAATGGTCAGTACGCCGCCGTTGCTTTCGCAGTCg                                                 >  1:6666897/1‑93 (MQ=255)
                                           aTACCTGGTGCAGCGGCTCCTGGCGGTTGATAATGATTTTGCTGCCACTCTCAAAGGTAATGGTCAGTACGCCGCCGTTGATTTCGCAGTCGa                                                <  1:62935/93‑1 (MQ=255)
                                            tACCTGGTGCAGCGGCTCCTGGCGGTTGATAATGATTTTGCTGCCACTCTCAAAGGTAATGGTCAGTACGCCGCCGTTGATTTCGCAGTCGat                                               >  1:2170115/1‑93 (MQ=255)
                                            tACCTGGTGCAGCGGCTCCTGGCGGTTGATAATGATTTTGCTGCCACTCTCAAAGGTAATGGTCAGTACGCCGCCGTTGATTTCGCAGTCGat                                               >  2:7007837/1‑93 (MQ=255)
                                             aCCTGGCGCAGCGGCTCCTGGCGGTTGATAATGATTTTGCTGCCACTCTCAAAGGTAATGGTCAGTACGCCGCCGTTGATTTCGCAGTCGata                                              >  2:437500/1‑93 (MQ=255)
                                               cTGGTGCAGCGGCTCCTGGCGGTTGATAATGATTTTGCTGCCACTCTCAAAGGTAATGGTCAGTACGCCGCCGTTGATTTCGCAGTCGATATc                                            >  2:7579639/1‑93 (MQ=255)
                                                tGGTGCAGCGGCTCCTGGCGGTTGATAATGATTTTGCTGCCACTCTCAAAGGTAATGGTCAGTACGCCGCCGTTGATTTCGCAGTCGATATCg                                           <  2:2363919/93‑1 (MQ=255)
                                                 ggTGCAGCGGCTCCTGGCGGTTGATAATGATTTTGCTGCCACTCTCAAAGGTAATGGTCAGTACGCCGCCGTTGATTTCGCAGTCGATATCGc                                          <  2:545772/93‑1 (MQ=255)
                                                 ggTGCAGCGGCTCCTGGCGGTTGATAATGATTTTGCTGCCACTCTCAAAGGTAATGGTCAGTACGCCGCCGTTGATTTCGCAGTCGATATCGc                                          >  1:5063248/1‑93 (MQ=255)
                                                    gCAGCGGCTCCTGGCGGTTGATAATGATTTTGCTGCCACTCTCAAAGGTAATGGTCAGTACGCCGCCGTTGATTTCGCAGTCGATATCGCTGt                                       <  1:6399452/93‑1 (MQ=255)
                                                      aGCGGCTCCTGGCGGTTGATAATGATTTTGCTGCCACTCTCAAAGGTAATGGTCAGTACGCCGCCGTTGATTTCGCAGTCGATATCGCTGTCg                                     <  1:3199762/93‑1 (MQ=255)
                                                       gCGGCTCCTGGCGGTTGATAATGATTTTGCTGCCACTCTCAAAGGTAATGGTCAGTACGCCGCCGTTGATTTCGCAGTCGATATCGCTGTCGc                                    >  1:4525289/1‑93 (MQ=255)
                                                        cGGCTCCTGGCGGTTGATAATGATTTTGCTGCCACTCTCAAAGGTAATGGTCAGTACGCCGCCGTTGATTTCGCAGTCGATATCGCTGTCGcc                                   <  1:444679/93‑1 (MQ=255)
                                                         ggCTCCTGGCGGTTGATAATGATTTTGCTGCCACTCTCAAAGGTAATGGTCAGTACGCCGCCGTTGATTTCGCAGTCGATATCGCTGTCGCCa                                  <  1:792096/93‑1 (MQ=255)
                                                         ggCTCCTGGCGGTTGATAATGATTTTGCTGCCACTCTCAAAGGTAATGGTCAGTACGCCGCCGTTGATTTCGCAGTCGATATCGCTGTCGCCa                                  <  2:3688533/93‑1 (MQ=255)
                                                         ggCTCCTGGCGGTTAATAATGATTTTGCTGCCACTCTCAAAGGTAATGGTCAGTACGCCGCCGTTGATTTCGCAGTCGATATCGCTGTCGCCa                                  <  1:1827530/93‑1 (MQ=255)
                                                          gCTCCTGGCGGTTGATAATGATTTTGCTGCCACTCTCAAAGGTAATGGTCAGTACGCCGCCGTTGATTTCGCAGTCGATATCGCTGTCGCCAt                                 >  1:4643681/1‑93 (MQ=255)
                                                          gCTCCTGGCGGTTGATAATGATTTTGCTGCCACTCTCAAAGGTAATGGTCAGTACGCCGCCGTTGATTTCGCAGTCGATATCGCTGTCGCCAt                                 >  1:2313714/1‑93 (MQ=255)
                                                            tCCTGGCGGTTGATAATGATTTTGCTGCCACTCTCAAAGGTAATGGTCAGTACGCCGCCGTTGATTTCGCAGTCGATATCGCTGTCGCCATcc                               >  1:2937285/1‑93 (MQ=255)
                                                             ccTGGCGGTTGATAATGATTTTGCTGCCACTCTCAAAGGTAATGGTCAGTACGCCGCCGTTGATTTCGCAGTCGATATCGCTGTCGCCATccc                              <  2:4740954/93‑1 (MQ=255)
                                                              cTGGCGGTTGATAATGATTTTGCTGCCACTCTCAAAGGTAATGGTCAGTACGCCGCCGTTGATTTCGCAGTCGATATCGCTGTCGCCATCCCa                             >  2:6789927/1‑93 (MQ=255)
                                                              cTGGCGGTTGATAATGATTTTGCTGCCACTCTCAAAGGTAATGGTCAGTACGCCGCCGTTGATTTCGCAGTCGATATCGCTGTCGCCATCCCa                             >  1:3585586/1‑93 (MQ=255)
                                                               tGGCGGTTGATAATGATTTTGCTGCCACTCTCAAAGGTAATGGTCAGTACGCCGCCGTTGATTTCGCAGTCGATATCGCTGTCGCCATCCCAg                            <  1:7115540/93‑1 (MQ=255)
                                                                  cGGTTGATAATGATTTTGCTGCCACTCTCAAAGGTAATGGTCAGTACGCCGCCGTTGATTTCGCAGTCGATATCGCTGTCGCCATCCCAgtca                         <  2:5135644/93‑2 (MQ=255)
                                                                    gTTGATAATGATTTTGCTGCCACTCTCAAAGGTAATGGTCAGTACGCCGCCGTTGATTTCGCAGTCGATATCGCTGTCGCCATCCCAgtcatc                       >  2:1969236/1‑93 (MQ=255)
                                                                    gTTGATAATGATTTTGCTGCCACTCTCAAAGGTAATGGTCAGTACGCCGCCGTTGATTTCGCAGTCGATATCGCTGTCGCCATCCCAgtcatc                       >  1:3712675/1‑93 (MQ=255)
                                                                    gTTGATAATGATTTTGCTGCCACTCTCAAAGGTAATGGTCAGTACGCCGCCGTTGATTTCGCAGTCGATATCGCTGTCGCCATCCCAgtcatc                       >  2:4188052/1‑93 (MQ=255)
                                                                    gTTGATAATGATTTTGCTGCCACTCTCAAAGGTAATGGTCAGTACGCCGCCGTTGATTTCGCAGTCGATATCGCTGTCGCCATCCCAgtcatc                       >  2:5773753/1‑93 (MQ=255)
                                                                    gTTGATAATGATTTTGCTGCCACTCTCAAAGGTAATGGTCAGTACGCCGCCGTTGATTTCGCAGTCGATATCGCTGTCGCCATCCCAgtcatc                       >  2:5372913/1‑93 (MQ=255)
                                                                       gATAATGATTTTGCTGCCACTCTCAAAGGTAATGGTCAGTACGCCGCCGTTGATTTCGCAGTCGATATCGCTGTCGCCATCCCAGTCATCCAg                    <  1:971082/93‑1 (MQ=255)
                                                                       gATAATGATTTTGCTGCCACTCTCAAAGGTAATGGTCAGTACGCCGCCGTTGATTTCGCAGTCGATATCGCTGTCGCCATCCCAGTCATCCAg                    >  1:2521524/1‑93 (MQ=255)
                                                                       gATAATGATTTTGCTGCCACTCTCAAAGGTAATGGTCAGTACGCCGCCGTTGATTTCGCAGTCGATATCGCTGTCGCCATCCCAGTCATCCAg                    >  2:2916088/1‑93 (MQ=255)
                                                                        aTAATGATTTTGCTGCCACTCTCAAAGGTAATGGTCAGTACGCCGCCGTTGATTTCGCAGTCGATATCGCTGTCGCCATCCCAGTCATCCAgg                   <  2:2823850/93‑1 (MQ=255)
                                                                        aTAATGATTTTGCTGCCACTCTCAAAGGTAATGGTCAGTACGCCGCCGTTGATTTCGCAGTCGATATCGCTGTCGCCATCCCAGTCATCCAgg                   <  1:4046735/93‑1 (MQ=255)
                                                                        aTAATGATTTTGCTGCCACTCTCAAAGGTAATGGTCAGTACGCCGCCGTTGATTTCGCAGTCGATATCGCTGTCGCCATCCCAGTCATCCAgg                   >  1:2393523/1‑93 (MQ=255)
                                                                          aaTGATTTTGCTGCCACTCTCAAAGGTAATGGTCAGTACGCCGCCGTTGATTTCGCAGTCGATATCGCTGTCGCCATCCCAGTCATCCAGGCg                 >  2:4436826/1‑93 (MQ=255)
                                                                          aaTGATTTTGCTGCCACTCTCAAAGGTAATGGTCAGTACGCCGCCGTTGATTTCGCAGTCGATATCGCTGTCGCCATCCCAGTCATCCAGGCg                 >  2:1625299/1‑93 (MQ=255)
                                                                          aaTGATTTTGCTGCCACTCTCAAAGGTAATGGTAAGTACGCCGCCGTTGATTTCGCAGTCGATATCGCTGTCGCCATCCCAGTCATCCAGGCg                 >  1:4131584/1‑93 (MQ=255)
                                                                           aTGATTTTGCTGCCACTCTCAAAGGTAATGGTCAGTACGCCGCCGTTGATTTCGCAGTCGATATCGCTGTCGCCATCCCAGTCATCCAGGCGt                <  1:6749378/93‑1 (MQ=255)
                                                                           aTGATTTTGCTGCCACTCTCAAAGGTAATGGTCAGTACGCCGCCGTTGATTTCGCAGTCGATATCGCTGTCGCCATCCCAGTCATCCAGGCGt                >  2:3561919/1‑93 (MQ=255)
                                                                           aTGATTTTGCTGCCACTCTCAAAGGTAATGGTCAGTACGCCGCCGTTGATTTCGCAGTCGATATCGCTGTCGCCATCCCAGTCATCCAGGCGt                >  2:3196609/1‑93 (MQ=255)
                                                                            tGATTTTGCTGCCACTCTCAAAGGTAATGGTCAGTACGCCGCCGTTGATTTCGCAGTCGATATCGCTGTCGCCATCCCAGTCATCCAGGCGtt               >  1:4315896/1‑93 (MQ=255)
                                                                            tGATTTTGCTGCCACTCTCAAAGGTAATGGTCAGTACGCCGCCGTTGATTTCGCAGTCGATATCGCTGTCGCCATCCCAGTCATCCAGGCGtt               >  1:7830595/1‑93 (MQ=255)
                                                                              aTTTTGCTGCCACTCTCAAAGGTAATGGTCAGTACGCCGCCGTTGATTTCGCAGTCGATATCGCTGTCGCCATCCCAGTCATCCAGGCGttct             <  1:3275131/93‑1 (MQ=255)
                                                                              aTTTTGCTGCCACTCTCAAAGGTAATGGTCAGTACGCCGCCGTTGATTTCGCAGTCGATATCGCTGTCGCCATCCCAGTCATCCAGGCGttct             >  1:3440804/1‑93 (MQ=255)
                                                                               ttttGCTGCCACTCTCAAAGGTAATGGTCAGTACGCCGCCGTTGATTTCGCAGTCGATATCGCTGTCGCCATCCCAGTCATCCAGGCGttctt            <  2:2171843/93‑1 (MQ=255)
                                                                               ttttGCTGCCACTCTCAAAGGTAATGGTCAGTACGCCGCCGTTGATTTCGCAGTCGATATCGCTGTCGCCATCCCAGTCATCCAGGCGttctt            >  1:6177133/1‑93 (MQ=255)
                                                                                 ttGCTGCCACTCTCAAAGGTAATGGTCAGTACGCCGCCGTTGATTTCGCAGTCGATATCGCTGTCGCCATCCCAGTCATCCAGGCGTTCTTCa          <  2:6565018/93‑1 (MQ=255)
                                                                                    ctgcCACTCTCAAAGGTAATGGTCAGTACGCCGCCGTTGATTTCGCAGTCGATATCGCTGTCGCCATCCCAGTCATCCAGGCGTTCTTCAATg       <  2:490327/93‑1 (MQ=255)
                                                                                      tccACTCTCAAAGGTAATGGTCAGTACGCCGCCGTTGATTTCGCAGTCGATATCGCTGTCGCCATCCCAGTCATCCAGGCGTTCTTCAATGGt     <  2:1209844/92‑1 (MQ=255)
                                                                                       ccACTCTCAAAGGGAATGGTCAGCACGCCGCCGTTGATTTCGCAGTCGATATCCCTGTCGCCATCCCAGCTATCCGGGCGATCTTCAATGGta    >  1:1326343/1‑92 (MQ=255)
                                                                                         aTTCTCAAAGGTAATGCTCAGTAAGCATCCGTAGCTGTCGCATTCGATATCGCTGCCGCCACCCCAGTCTTCCAGGCGTTCTTCAATGGTCAg  <  2:7677986/93‑1 (MQ=255)
                                                                                          |                                                                                           
CCTTTCAGGTCAAAATGGTAGCCGCCCTGTTTGGTTGCCAGCCATACCTGGTGCAGCGGCTCCTGGCGGTTGATAATGATTTTGCTGCCATTCTCAAAGGTAATGGTCAGTACGCCGCCGTTGATTTCGCAGTCGATATCGCTGTCGCCATCCCAGTCGTCCAGGCGTTCTTCAATGGTCAG  >  NC_000913/3991449‑3991630

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