Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,342,604 G→A A10A (GCC→GCT b4523 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,342,6040GA84.4% 627.7 / 75.8 231A10A (GCC→GCTb4523hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (25/11);  new base A (113/82);  total (138/93)
Fisher's exact test for biased strand distribution p-value = 2.67e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.77e-02

Only 100 of 234 total aligned reads displayed.

CGCAAGAGTTCTGCTTTAAGTATCTCGAGATTATGGATCGCAGAGTGGTAATCGCCTGCCACAAGGATATCTAACACAATATCAATACGTTCGGCAACGCGCTGGGCATCAAATTCAGACATCAAGGACATCCTTATAGCAAAAGGGGGAAAATGAATAATGCAAAAATTGCCGCTAAAAGA  >  NC_000913/1342512‑1342693
                                                                                            |                                                                                         
cgcAAGAGTTCTGCTTTAAGTATCTCGAGATTATGGATCGCAGAGTGGTAATCGCCTGCCACAAGGATATCTAACACAATATCAATACGTTCg                                                                                           >  2:5527911/1‑93 (MQ=255)
   aaGAGTTCTGCTTTAAGTATCTCGAGATTATGGATCGCAGAGTGGTAATCGCCTGCCACAAGGATATCTAACACAATATCAATACGTTCAGCa                                                                                        >  1:1148969/1‑93 (MQ=255)
   aaGAGTTCTGCTTTAAGTATCTCGAGATTATGGATCGCAGAGTGGTAATCGCCTGCCACAAGGATATCTAACACAATATCAATACGTTCAGCa                                                                                        >  1:2997405/1‑93 (MQ=255)
    agagTTCTGCTTTAAGTATCTCGAGATTATGGATCGCAGAGTGGTAATCGCCTGCCACAAGGATATCTAACACAATATCAATACGTTCGGCaa                                                                                       >  1:5112005/1‑93 (MQ=255)
    agagTTCTGCTTTAAGTATCTCGAGATTATGGATCGCAGAGTGGTAATCGCCTGCCACAAGGATATCTAACACAATATCAATACGTTCAGCaa                                                                                       >  1:5928624/1‑93 (MQ=255)
    agagTTCTGCTTTAAGTATCTCGAGATTATGGATCGCAGAGTGGTAATCGCCTGCCACAAGGATATCTAACACAATATCAATACGTTCAGCaa                                                                                       >  2:1254554/1‑93 (MQ=255)
     gagTTCTGCTTTAAGTATCTCGAGATTATGGATCGCAGAGTGGTAATCGCCTGCCACAAGGATATCTAACACAATATCAATACGTTCGGCAAc                                                                                      >  1:6988114/1‑93 (MQ=255)
     gagTTCTGCTTTAAGTATCTCGAGATTATGGATCGCAGAGTGGTAATCGCCTGCCACAAGGATATCTAACACAATATCAATACGTTCGGCAAc                                                                                      <  2:4943779/93‑1 (MQ=255)
      agTTCTGCTTTAAGTATCTCGAGATTATGGATCGCAGAGTGGTAATCGCCTGCCACAAGGATATCTAACACAATATCAATACGTTCAGCAAcg                                                                                     >  2:3846394/1‑93 (MQ=255)
       gTTCTGCTTTAAGTATCTCGAGATTATGGATCGCAGAGTGGTAATCGCCTGCCACAAGGATATCTAACACAATATCAATACGTTCAGCAAcgc                                                                                    >  1:5355333/1‑93 (MQ=255)
       gTTCTGCTTTAAGTATCTCGAGATTATGGATCGCAGAGTGGTAATCGCCTGCCACAAGGATATCTAACACAATATCAATACGTTCAGCAAcgc                                                                                    >  1:2384792/1‑93 (MQ=255)
       gTTCTGCTTTAAGTATCTCGAGATTATGGATCGCAGAGTGGTAATCGCCTGCCACAAGGATATCTAACACAATATAAATACGTTCGGCAAcgc                                                                                    >  1:4266499/1‑93 (MQ=255)
        ttCTGCTTTAAGTATCTCGAGATTATGGATCGCAGAGTGGTAATCGCCTGCCACAAGGATATCTAACACAATATCAATACGTTCAGCAAcgcg                                                                                   >  2:7166337/1‑93 (MQ=255)
         tCTGCTTTAAGTATCTCGAGATTATGGATCGCAGAGTGGTAATCGCCTGCCACAAGGATATCTAACACAATATCAATACGTTCAGCAACgcgc                                                                                  >  1:2677434/1‑93 (MQ=255)
         tCTGCTTTAAGTATCTCGAGATTATGGATCGCAGAGTGGTAATCGCCTGCCACAAGGATATCTAACACAATATCAATACGTTCAGCAACgcgc                                                                                  >  1:895289/1‑93 (MQ=255)
          cTGCTTTAAGTATCTCGAGATTATGGATCGCAGAGTGGTAATCGCCTGCCACAAGGATATCTAACACAATATCAATACGTTCAGCAACGCGCt                                                                                 <  2:5575501/93‑1 (MQ=255)
          cTGCTTTAAGTATCTCGAGATTATGGATCGCAGAGTGGTAATCGCCTGCCACAAGGATATCTAACACAATATCAATACGTTCAGCAACGCGCt                                                                                 <  2:616596/93‑1 (MQ=255)
            gCTTTAAGTATCTCGAGATTATGGATCGCAGAGTGGTAATCGCCTGCCACAAGGATATCTAACACAATATCAATACGTTCAGCAACGCGCTgg                                                                               >  2:3373020/1‑93 (MQ=255)
             cTTTAAGTATCTCGAGATTATGGATCGCAGAGTGGTAATCGCCTGCCACAAGGATATCTAACACAATATCAATACGTTCAGCAACGCGCTggg                                                                              >  1:5522442/1‑93 (MQ=255)
              tttAAGTATCTCGAGATTATGGATCGCAGAGTGGTAATCGCCTGCCACAAGGATATCTAACACAATATCAATACGTTCAGCAACGCGCTGGGc                                                                             >  2:1860907/1‑93 (MQ=255)
              tttAAGTATCTCGAGATTATGGATCGCAGAGTGGTAATCGCCTGCCACAAGGATATCTAACACAATATCAATACGTTCAGCAACGCGCTGGGc                                                                             >  1:7453973/1‑93 (MQ=255)
                 aaGTATCTCGAGATTATGGATCGCAGAGTGGTAATCGCCTGCCACAAGGATATCTAACACAATATCAATACGTTCAGCAACGCGCTGGGCATc                                                                          <  2:7393921/93‑1 (MQ=255)
                  aGTATCTCGAGATTATGGATCGCAGAGTGGTAATCGCCTGCCACAAGGATATCTAACACAATATCAATACGTTCAGCAACGCGCTGGGCATCa                                                                         >  2:4224213/1‑93 (MQ=255)
                  aGTATCTCGAGATTATGGATCGCAGAGTGGTAATCGCCTGCCACAAGGATATCTAACACAATATCAATACGTTCAGCAACGCGCTGGGCATCa                                                                         >  2:3785846/1‑93 (MQ=255)
                  aGTATCTCGAGATTATGGATCGCAGAGTGGTAATCGCCTGCCACAAGGATATCTAACACAATATCAATACGTTCAGCAACGCGCTGGGCATCa                                                                         >  1:62557/1‑93 (MQ=255)
                   gTATCTCGAGATTATGGATCGCAGAGTGGTAATCGCCTGCCACAAGGATATCTAACACAATATCAATACGTTCAGCAACGCGCTGGGCATCaa                                                                        >  2:1661874/1‑93 (MQ=255)
                     aTCTCGAGATTATGGATCGCAGAGTGGTAATCGCCTGCCACAAGGATATCTAACACAATATCAATACGTTCGGCAACGCGCTGGGCATCAAAt                                                                      >  2:4676331/1‑93 (MQ=255)
                     aTCTCGAGATTATGGATCGCAGAGTGGTAATCGCCTGCCACAAGGATATCTAACACAATATCAATACGTTCAGCAACGCGCTGGGCATCAAAt                                                                      >  1:5433008/1‑93 (MQ=255)
                     aTCTCGAGATTATGGATCGCAGAGTGGTAATCGCCTGCCACAAGGATATCTAACACAATATCAATACGTTCAGCAACGCGCTGGGCATCAAAt                                                                      >  2:7182193/1‑93 (MQ=255)
                     aTCTCGAGATTATGGATCGCAGAGTGGTAATCGCCTGCCACAAGGATATCTAACACAATATCAATACGTTCAGCAACGCGCTGGGCATCAAAt                                                                      >  1:84836/1‑93 (MQ=255)
                      tctcGAGATTATGGATCGCAGAGTGGTAATCGCCTGCCACAAGGATATCTAACACAATATCAATACGTTCAGCAACGCGCTGGGCATCAAAtt                                                                     <  2:2344061/93‑1 (MQ=255)
                          gagaTTATGGATCGCAGAGTGGTAATCGCCTGCCACAAGGATATCTAACACAATATCAATACGTTCGGCAACGCGCTGGGCATCAAATTCAGa                                                                 <  2:7072328/93‑1 (MQ=255)
                          gagaTTATGGATCGCAGAGTGGTAATCGCCTGCCACAAGGATATCTAACACAATATCAATACGTTCAGCAACGCGCTGGGCATCAAATTCAGa                                                                 <  2:146365/93‑1 (MQ=255)
                          gagaTTATGGATCGCAGAGTGGTAATCGCCTGCCACAAGGATATCTAACACAATATCAATACGTTCAGCAACGCGCTGGGCATCAAATTCAGa                                                                 >  1:7156477/1‑93 (MQ=255)
                            gaTTATGGATCGCAGAGTGGTAATCGCCTGCCACAAGGATATCTAACACAATATCAATACGTTCAGCAACGCGCTGGGCATCAAATTCAGACa                                                               >  2:1141884/1‑93 (MQ=255)
                              ttATGGATCGCAGAGTGGTAATCGCCTGCCACAAGGATATCTAACACAATATCAATACGTTCAGCAACGCGCTGGGCATCAAATTCAGACATc                                                             >  2:3246775/1‑93 (MQ=255)
                              ttATGGATCGCAGAGTGGTAATCGCCTGCCACAAGGATATCTAACACAATATCAATACGTTCAGCAACGCGCTGGGCATCAAATTCAGACATc                                                             >  2:6254018/1‑93 (MQ=255)
                               tATGGATCGCAGAGTGGTAATCGCCTGCCACAAGGATATCTAACACAATATCAATACGTTCAGCAACGCGCTGGGCATCAAATTCAGACATCa                                                            >  1:7282138/1‑93 (MQ=255)
                               tATGGATCGCAGAGTGGTAATCGCCTGCCACAAGGATATCTAACACAATATCAATACGTTCAGCAACGCGCTGGGCATCAAATTCAGACATCa                                                            >  2:1844207/1‑93 (MQ=255)
                                aTGGATCGCAGAGTGGTAATCGCCTGCCACAAGGATATCTAACACAATATCAATACGTTCAGCAACGCGCTGGGCATCAAATTCAGACATCaa                                                           <  2:5827017/93‑1 (MQ=255)
                                 tGGATCGCAGAGTGGTAATCGCCTGCCACAAGGATATCTAACACAATATCAATACGTTCGGCAACGCGCTGGGCATCAAATTCAGACATCAAg                                                          <  1:455773/93‑1 (MQ=255)
                                 tGGATCGCAGAGTGGTAATCGCCTGCCACAAGGATATCTAACACAATATCAATACGTTCAGCAACGCGCTGGGCATCAAATTCAGACATCAAg                                                          >  2:5988033/1‑93 (MQ=255)
                                 tGGATCGCAGAGTGGTAATCGCCTGCCACAAGGATATCTAACACAATATCAATACGTTCAGCAACGCGCTGGGCATCAAATTCAGACATCAAg                                                          >  1:4470495/1‑93 (MQ=255)
                                   gATCGCAGAGTGGTAATCGCCTGCCACAAGGATATCTAACACAATATCAATACGTTCAGCAACGCGCTGGGCATCAAATTCAGACATCAAGGa                                                        >  1:6369038/1‑93 (MQ=255)
                                      cGCAGAGTGGTAATCGCCTGCCACAAGGATATCTAACACAATATCAATACGTTCAGCAACGCGCTGGGCATCAAATTCAGACATCAAGGACAt                                                     <  1:5721405/93‑1 (MQ=255)
                                       gCAGAGTGGTAATCGCCTGCCACAAGGATATCTAACACAATATCAATACGTTCAGCAACGCGCTGGGCATCAAATTCAGACATCAAGGACATc                                                    >  1:2460051/1‑93 (MQ=255)
                                         agagTGGTAATCGCCTGCCACAAGGATATCTAACACAATATCAATACGTTCGGCAACGCGCTGGGCATCAAATTCAGACATCAAGGACATCCt                                                  >  2:1762774/1‑93 (MQ=255)
                                          gagTGGTAATCGCCTGCCACAAGGATATCTAACACAATATCAATACGTTCAGCAACGCGCTGGGCATCAAATTCAGACATCAAGGACATCCtt                                                 <  1:4392450/93‑1 (MQ=255)
                                           agTGGTAATCGCCTGCCACAAGGATATCTAACACAATATCAATACGTTCAGCAACGCGCTGGGCATCAAATTCAGACATCAAGGACATCCTta                                                >  1:4102270/1‑93 (MQ=255)
                                            gTGGTAATCGCCTGCCACAAGGATATCTAACACAATATCAATACGTTCAGCAACGCGCTGGGCATCAAATTCAGACATCAAGGACATCCTtat                                               <  2:2478922/93‑1 (MQ=255)
                                            gTGGTAATCGCCTGCCACAAGGATATCTAACACAATATCAATACGTTCAGCAACGCGCTGGGCATCAAATTCAGACATCAAGGACATCCTtat                                               >  2:5792359/1‑93 (MQ=255)
                                             tGGTAATCGCCTGCCACAAGGATATCTAACACAATATCAATACGTTCAGCAACGCGCTGGGCATCAAATTCAGACATCAAGGACATCCTtata                                              >  2:4316352/1‑93 (MQ=255)
                                             tGGTAATCGCCTGCCACAAGGATATCTAACACAATATCAATACGTTCAGCAACGCGCTGGGCATCAAATTCAGACATCAAGGACATCCTtata                                              <  2:7424411/93‑1 (MQ=255)
                                              ggTAATCGCCTGCCACAAGGATATCTAACACAATATCAATACGTTCAGCAACGCGCTGGGCATCAAATTCAGACATCAAGGACATCCTTATAg                                             <  1:1386418/93‑1 (MQ=255)
                                                 aaTCGCCTGCCACAAGGATATCTAACACAATATCAATACGTTCAGCAACGCGCTGGGCATCAAATTCAGACATCAAGGACATCCTTATAGCaa                                          >  2:7341770/1‑93 (MQ=255)
                                                  aTCGCCTGCCACAAGGATATCTAACACAATATCAATACGTTCGGCAACGCGCTGGGCATCAAATTCAGACATCAAGGACATCCTTATAGCaaa                                         >  1:6851528/1‑93 (MQ=255)
                                                   tCGCCTGCCACAAGGATATCTAACACAATATCAATACGTTCAGCAACGCGCTGGGCATCAAATTCAGACATCAAGGACATCCTTATAGCaaaa                                        <  2:5487050/93‑1 (MQ=255)
                                                   tCGCCTGCCACAAGGATATCTAACACAATATCAATACGTTCAGCAACGCGCTGGGCATCAAATTCAGACATCAAGGACATCCTTATAGCaaaa                                        <  1:552026/93‑1 (MQ=255)
                                                   tCGCCTGCCACAAGGATATCTAACACAATATCAATACGTTCAGCAACGCGCTGGGCATCAAATTCAGACATCAAGGACATCCTTATAGCaaaa                                        >  2:2101010/1‑93 (MQ=255)
                                                    cGCCTGCCACAAGGATATCTAACACAATATCAATACGTTCAGCAACGCGCTGGGCATCAAATTCAGACATCAAGGACATCCTTATAGCAAAAg                                       <  1:3950489/93‑1 (MQ=255)
                                                       cTGCCACAAGGATATCTAACACAATATCAATACGTTCAGCAACGCGCTGGGCATCAAATTCAGACATCAAGGACATCCTTATAGCAAAAgggg                                    <  2:4532554/93‑1 (MQ=255)
                                                         gCCACAAGGATATCTAACACAATATCAATACGTTCAGCAACGCGCTGGGCATCAAATTCAGACATCAAGGACATCCTTATAGCAAAAGGGGGa                                  <  2:6712290/93‑1 (MQ=255)
                                                            acaAGGATATCTAACACAATATCAATACGTTCAGCAACGCGCTGGGCATCAAATTCAGACATCAAGGACATCCTTATAGCAAAAGGGGGaaaa                               >  2:4791693/1‑93 (MQ=255)
                                                             caAGGATATCTAACACAATATCAATACGTTCAGCAACGCGCTGGGCATCAAATTCAGACATCAAGGACATCCTTATAGCAAAAGGGGGAAAAt                              >  2:6320042/1‑93 (MQ=255)
                                                             caAGGATATCTAACACAATATCAATACGTTCAGCAACGCGCTGGGCATCAAATTCAGACATCAAGGACATCCTTATAGCAAAAGGGGGAAAAt                              >  1:5163661/1‑93 (MQ=255)
                                                             caAGGATATCTAACACAATATCAATACGTTCAGCAACGCGCTGGGCATCAAATTCAGACATCAAGGACATCCTTATAGCAAAAGGGGGAAAAt                              <  1:205113/93‑1 (MQ=255)
                                                             caAGGATATCTAACACAATATCAATACGTTCAGCAACGCGCTGGGCATCAAATTCAGACATCAAGGACATCCTTATAGCAAAAGGGGGAAAAt                              <  1:491575/93‑1 (MQ=255)
                                                               aGGATATCTAACACAATATCAATACGTTCGGCAACGCGCTGGGCATCAAATTCAGACATCAAGGACATCCTTATAGCAAAAGGGGGAAAATGa                            <  1:4676331/93‑1 (MQ=255)
                                                               aGGATATCTAACACAATATCAATACGTTCAGCAACGCGCTGGGCATCAAATTCAGACATCAAGGACATCCTTATAGCAAAAGGGGGAAAATGa                            >  1:2681943/1‑93 (MQ=255)
                                                                ggATATCTAACACAATATCAATACGTTCAGCAACGCGCTGGGCATCAAATTCAGACATCAAGGACATCCTTATAGCAAAAGGGGGAAAATGaa                           <  2:2737788/93‑1 (MQ=255)
                                                                ggATATCTAACACAATATCAATACGTTCAGCAACGCGCTGGGCATCAAATTCAGACATCAAGGACATCCTTATAGCAAAAGGGGGAAAATGaa                           <  1:3284759/93‑1 (MQ=255)
                                                                  atatCTAACACAATATCAATACGTTCAGCAACGCGCTGGGCATCAAATTCAGACATCAAGGACATCCTTATAGCAAAAGGGGGAAAATGaata                         <  2:7505495/93‑1 (MQ=255)
                                                                  atatCTAACACAATATCAATACGTTCAGCAACGCGCTGGGCATCAAATTCAGACATCAAGGACATCCTTATAGCAAAAGGGGGAAAATGaata                         >  1:6460301/1‑93 (MQ=255)
                                                                   tatCTAACACAATATCAATACGTTCAGCAACGCGCTGGGCATCAAATTCAGACATCAAGGACATCCTTATAGCAAAAGGGGGAAAATGaataa                        >  2:7478541/1‑93 (MQ=255)
                                                                    atCTAACACAATATCAATACGTTCGGCAACGCGCTGGGCATCAAATTCAGACATCAAGGACATCCTTATAGCAAAAGGGGGAAAATGaataat                       >  1:1403795/1‑93 (MQ=255)
                                                                     tCTAACACAATATCAATACGTTCAGCAACGCGCTGGGCATCAAATTCAGACATCAAGGACATCCTTATAGCAAAAGGGGGAAAATGAATAATg                      >  1:1234648/1‑93 (MQ=255)
                                                                     tCTAACACAATATCAATACGTTCAGCAACGCGCTGGGCATCAAATTCAGACATCAAGGACATCCTTATAGCAAAAGGGGGAAAATGAATAATg                      <  2:3023134/93‑1 (MQ=255)
                                                                      cTAACACAATATCAATACGTTCAGCAACGCGCTGGGCATCAAATTCAGACATCAAGGACATCCTTATAGCAAAAGGGGGAAAATGAATAATGc                     <  1:6715990/93‑1 (MQ=255)
                                                                      cTAACACAATATCAATACGTTCAGCAACGCGCTGGGCATCAAATTCAGACATCAAGGACATCCTTATAGCAAAAGGGGGAAAATGAATAATGc                     <  2:4632880/93‑1 (MQ=255)
                                                                      cTAACACAATATCAATACGTTCAGCAACGCGCTGGGCATCAAATTCAGACATCAAGGACATCCTTATAGCAAAAGGGGGAAAATGAATAATGc                     >  2:1804940/1‑93 (MQ=255)
                                                                      cTAACACAATATCAATACGTTCAGCAACGCGCTGGGCATCAAATTCAGACATCAAGGACATCCTTATAGCAAAAGGGGGAAAATGAATAATGc                     <  1:3165441/93‑1 (MQ=255)
                                                                       tAACACAATATCAATACGTTCGGCAACGCGCTGGGCATCAAATTCAGACATCAAGGACATCCTTATAGCAAAAGGGGGAAAATGAATAATGCa                    <  2:2120245/93‑1 (MQ=255)
                                                                       tAACACAATATCAATACGTTCAGCAACGCGCTGGGCATCAAATTCAGACATCAAGGACATCCTTATAGCAAAAGGGGGAAAATGAATAATGCa                    >  1:6074625/1‑93 (MQ=255)
                                                                         acacAATATCAATACGTTCGGCAACGCGCTGGGCATCAAATTCAGACATCAAGGACATCCTTATAGCAAAAGGGGGAAAATGAATAATGCaaa                  >  1:4847752/1‑93 (MQ=255)
                                                                           acaATATCAATACGTTCAGCAACGCGCTGGGCATCAAATTCAGACATCAAGGACATCCTTATAGCAAAAGGGGGAAAATGAATAATGCaaaaa                >  1:7351660/1‑93 (MQ=255)
                                                                           acaATATCAATACGTTCAGCAACGCGCTGGGCATCAAATTCAGACATCAAGGACATCCTTATAGCAAAAGGGGGAAAATGAATAATGCaaaaa                >  2:5876053/1‑93 (MQ=255)
                                                                           acaATATCAATACGTTCAGCAACGCGCTGGGCATCAAATTCAGACATCAAGGACATCCTTATAGCAAAAGGGGGAAAATGAATAATGCaaaaa                <  2:1120873/93‑1 (MQ=255)
                                                                            caATATCAATACGTTCAGCAACGCGCTGGGCATCAAATTCAGACATCAAGGACATCCTTATAGCAAAAGGGGGAAAATGAATAATGCAAAAAt               <  2:5248431/93‑1 (MQ=255)
                                                                            caATATCAATACGTTCAGCAACGCGCTGGGCATCAAATTCAGACATCAAGGACATCCTTATAGCAAAAGGGGGAAAATGAATAATGCAAAAAt               >  2:4100260/1‑93 (MQ=255)
                                                                               tatCAATACGTTCAGCAACGCGCTGGGCATCAAATTCAGACATCAAGGACATCCTTATAGCAAAAGGGGGAAAATGAATAATGCAAAAATTGc            <  2:4912599/93‑1 (MQ=255)
                                                                                 tCAATACGTTCAGCAACGCGCTGGGCATCAAATTCAGACATCAAGGACATCCTTATAGCAAAAGGGGGAAAATGAATAATGCAAAAATTGCCg          >  2:2781477/1‑93 (MQ=255)
                                                                                  cAATACGTTCAGCAACGCGCTGGGCATCAAATTCAGACATCAAGGACATCCTTATAGCAAAAGGGGGAAAATGAATAATGCAAAAATTGCCGc         <  1:1592311/93‑1 (MQ=255)
                                                                                  cAATACGTTCAGCAACGCGCTGGGCATCAAATTCAGACATCAAGGACATCCTTATAGCAAAAGGGGGAAAATGAATAATGCAAAAATTGCCGc         <  1:3552165/93‑1 (MQ=255)
                                                                                   aaTACGTTCAGCAACGCGCTGGGCATCAAATTCAGACATCAAGGACATCCTTATAGCAAAAGGGGGAAAATGAATAATGCAAAAATTGCCGCt        >  2:3735768/1‑93 (MQ=255)
                                                                                   aaTACGTTCAGCAACGCGCTGGGCATCAAATTCAGACATCAAGGACATCCTTATAGCAAAAGGGGGAAAATGAATAATGCAAAAATTGCCGCt        >  2:5783782/1‑93 (MQ=255)
                                                                                     tACGTTCAGCAACGCGCTGGGCATCAAATTCAGACATCAAGGACATCCTTATAGCAAAAGGGGGAAAATGAATAATGCAAAAATTGCCGCTaa      >  2:6873427/1‑93 (MQ=255)
                                                                                       cGTTCAGCAACGCGCTGGGCATCAAATTCAGACATCAAGGACATCCTTATAGCAAAAGGGGGAAAATGAATAATGCAAAAATTGCCGCTaaaa    >  2:1312298/1‑93 (MQ=255)
                                                                                        gTTCGGCAACGCGCTGGGCATCAAATTCAGACATCAAGGACATCCTTATAGCAAAAGGGGGAAAATGAATAATGCAAAAATTGCCGCTAAAag   >  2:6410136/1‑93 (MQ=255)
                                                                                         ttCAGCAACGCGCTGGGCATCAAATTCAGACATCAAGGACATCCTTATAGCAAAAGGGGGAAAATGAATAATGCAAAAATTGCCGCTAAAaga  >  1:1005211/1‑93 (MQ=255)
                                                                                         ttCAGCAACGCGCTGGGCATCAAATTCAGACATCAAGGACATCCTTATAGCAAAAGGGGGAAAATGAATAATGCAAAAATTGCCGCTAAAaga  <  2:1900459/93‑1 (MQ=255)
                                                                                            |                                                                                         
CGCAAGAGTTCTGCTTTAAGTATCTCGAGATTATGGATCGCAGAGTGGTAATCGCCTGCCACAAGGATATCTAACACAATATCAATACGTTCGGCAACGCGCTGGGCATCAAATTCAGACATCAAGGACATCCTTATAGCAAAAGGGGGAAAATGAATAATGCAAAAATTGCCGCTAAAAGA  >  NC_000913/1342512‑1342693

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