Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 235,988 T→C V5V (GTA→GTG rnhA ← RNase HI

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913235,9880TC85.5% 729.9 / 109.1 303V5V (GTA→GTGrnhARNase HI
Reads supporting (aligned to +/- strand):  ref base T (28/16);  new base C (134/125);  total (162/141)
Fisher's exact test for biased strand distribution p-value = 1.91e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.96e-01

Only 100 of 303 total aligned reads displayed.

AAAGGTTTTCTCGCGTCCGCGATAGCGTAAAATAGCGCCGTAACCCCCAGGTCCTGGATTGCCCAGACACGAACCATCGGTGAAAATTTCTACCTGTTTAAGCATCTCTGGTAGACTTCCTGTAATTGAATCGAACTGTAAAACGACAAGTCTGACATAAATGACCGCTATGAGCACTG  >  NC_000913/235898‑236076
                                                                                          |                                                                                        
aaaGGTTTTCTCGCGTCCGCGATAGCGTAAAATAGCGCCGTAACCCCCAGGTCCTGGATTGCCCAGACACGAACCATCGGTGAAAATTTCCAc                                                                                        >  2:1750846/1‑93 (MQ=255)
 aaGGTTTTCTCGCGTCCGCGATAGCGTAAAATAGCGCCGTAACCCCCAGGTCCTGGATTGCCCAGACACGAACCATCGGTGAAAATTTCCAcc                                                                                       <  2:7341875/93‑1 (MQ=255)
  aGGTTTTCTCGCGTCCGCGATAGCGTAAAATAGCGCCGTAACCCCCAGGTCCTGGATTGCCCAGACACGAACCATCGGTGAAAATTTCCACCt                                                                                      <  1:5604221/93‑1 (MQ=255)
    gTTTTCTCGCGTCCGCGATAGCGTAAAATAGCGCCGTAACCCCCAGGTCCTGGATTGCCCAGACACGAACCATCGGTGAAAATTTCCACCTGt                                                                                    <  2:2198548/93‑1 (MQ=255)
     ttttCTCGCGTCCGCGATAGCGTAAAATAGCGCCGTAACCCCCAGGTCCTGGATTGCCCAGACACGAACCATCGGTGAAAATTTCCACCTGtt                                                                                   <  2:503522/93‑1 (MQ=255)
     ttttCTCGCGTCCGCGATAGCGTAAAATAGCGCCGTAACCCCCAGGTCCTGGATTGCCCAGACACGAACCATCGGTGAAAATTTCCACCTGtt                                                                                   <  1:937426/93‑1 (MQ=255)
     ttttCTCGCGTCCGCGATAGCGTAAAATAGCGCCGTAACCCCCAGGTCCTGGATTGCCCAGACACGAACCATCGGTGAAAATTTCCACCTGtt                                                                                   <  1:611715/93‑1 (MQ=255)
     ttttCTCGCGTCCGCGATAGCGTAAAATAGCGCCGTAACCCCCAGGTCCTGGATTGCCCAGACACGAACCATCGGTGAAAATTTCCACCTGtt                                                                                   <  1:5809220/93‑1 (MQ=255)
       ttCTCGCGTCCGCGATAGCGTAAAATAGCGCCGTACCCCCCAGGTCCTGGATTGCCCAGACACGAACCATCGGTGAAAATTTCTACCTGTTTa                                                                                 <  1:6287752/93‑1 (MQ=255)
        tctcGCGTCCGCGATAGCGTAAAATAGCGCCGTAACCCCCAGGTCCTGGATTGCCCAGACACGAACCATCGGTGAAAATTTCCACCTGTTTaa                                                                                <  2:1929864/93‑1 (MQ=255)
        tctcGCGTCCGCGATAGCGTAAAATAGCGCCGTAACCCCCAGGTCCTGGATTGCCCAGACACGAACCATCGGTGAAAATTTCCACCTGTTTaa                                                                                <  2:4031506/93‑1 (MQ=255)
          tcGCGTCCGCGATAGCGTAAAATAGCGCCGTAACCCCCAGGTCCTGGATTGCCCAGACACGAACCATCGGTGAAAATTTCCACCTGTTTAAGc                                                                              >  1:4962244/1‑93 (MQ=255)
          tcGCGTCCGCGATAGCGTAAAATAGCGCCGTAACCCCCAGGTCCTGGATTGCCCAGACACGAACCATCGGTGAAAATTTCCACCTGTTTAAGc                                                                              <  2:7215126/93‑1 (MQ=255)
           cgcgTCCGCGATAGCGTAAAATAGCGCCGTAACCCCCAGGTCCTGGATTGCCCAGACACGAACCATCGGTGAAAATTTCCACCTGTTTAAGCa                                                                             <  2:3588216/93‑1 (MQ=255)
           cgcgTCCGCGATAGCGTAAAATAGCGCCGTAACCCCCAGGTCCTGGATTGCCCAGACACGAACCATCGGTGAAAATTTCCACCTGTTTAAGCa                                                                             <  2:3478122/93‑1 (MQ=255)
            gcgTCCGCGATAGCGTAAAATAGCGCCGTAACCCCCAGGTCCTGGATTGCCCAGACACGAACCATCGGTGAAAATTTCTACCTGTTTAAGCAt                                                                            >  2:5570951/1‑93 (MQ=255)
            gcgTCCGCGATAGCGTAAAATAGCGCCGTAACCCCCAGGTCCTGGATTGCCCAGACACGAACCATCGGTGAAAATTTCCACCTGTTTAAGCAt                                                                            <  1:1966909/93‑1 (MQ=255)
              gTCCGCGATAGCGTAAAATAGCGCCGTAACCCCCAGGTCCTGGATTGCCCAGACACGAACCATCGGTGAAAATTTCCACCTGTTTAAGCAtct                                                                          >  2:5275285/1‑93 (MQ=255)
               tCCGCGATAGCGTAAAATAGCGCCGTAACCCCCAGGTCCTGGATTGCCCAGACACGAACCATCGGTGAAAATTTCCACCTGTTTAAGCAtctc                                                                         >  2:5638379/1‑93 (MQ=255)
               tCCGCGATAGCGTAAAATAGCGCCGTAACCCCCAGGTCCTGGATTGCCCAGACACGAACCATCGGTGAAAATTTCCACCTGTTTAAGCAtctc                                                                         >  1:3811484/1‑93 (MQ=255)
                ccGCGATAGCGTAAAATAGCGCCGTAACCCCCAGGTCCTGGATTGCCCAGACACGAACCATCGGTGAAAATTTCTACCTGTTTAAGCATctct                                                                        <  2:3485247/93‑1 (MQ=255)
                ccGCGATAGCGTAAAATAGCGCCGTAACCCCCAGGTCCTGGATTGCCCAGACACGAACCATCGGTGAAAATTTCTACCTGTTTAAGCATctct                                                                        >  1:861924/1‑93 (MQ=255)
                ccGCGATAGCGTAAAATAGCGCCGTAACCCCCAGGTCCTGGATTGCCCAGACACGAACCATCGGTGAAAATTTCCACCTGTTTAAGCATctct                                                                        <  1:4693037/93‑1 (MQ=255)
                  gcgATAGCGTAAAATAGCGCCGTAACCCCCAGGTCCTGGATTGCCCAGACACGAACCATCGGTGAAAATTTCCACCTGTTTAAGCATCTCTgg                                                                      >  1:6481559/1‑93 (MQ=255)
                    gATAGCGTAAAATAGCGCCGTAACCCCCAGGTCCTGGATTGCCCAGACACGAACCATCGGTGAAAATTTCCACCTGTTTAAGCATCTCTGGTa                                                                    <  2:2727671/93‑1 (MQ=255)
                    gATAGCGTAAAATAGCGCCGTAACCCCCAGGTCCTGGATTGCCCAGACACGAACCATCGGTGAAAATTTCCACCTGTTTAAGCATCTCTGGTa                                                                    >  1:247754/1‑93 (MQ=255)
                     aTAGCGTAAAATAGCGCCGTAACCCCCAGGTCCTGGATTGCCCAGACACGAACCATCGGTGAAAATTTCCACCTGTTTAAGCATCTCTGGTAg                                                                   >  2:5793240/1‑93 (MQ=255)
                      tAGCGTAAAATAGCGCCGTAACCCCCAGGTCCTGGATTGCCCAGACACGAACCATCGGTGAAAATTTCCACCTGTTTAAGCATCTCTGGTAGa                                                                  >  1:5261974/1‑93 (MQ=255)
                        gCGTAAAATAGCGCCGTAACCCCCAGGTCCTGGATTGCCCAGACACGAACCATCGGTGAAAATTTCCACCTGTTTAAGCATCTCTGGTAGACt                                                                >  1:5917841/1‑93 (MQ=255)
                         cGTAAAATAGCGCCGTAACCCCCAGGTCCTGGATTGCCCAGACACGAACCATCGGTGAAAATTTCCACCTGTTTAAGCATCTCTGGTAGACtt                                                               <  2:6124548/93‑1 (MQ=255)
                           tAAAATAGCGCCGTAACCCCCAGGTCCTGGATTGCCCAGACACGAACCATCGGTGAAAATTTCCACCTGTTTAAGCATCTCTGGTAGACTTcc                                                             >  1:5145990/1‑93 (MQ=255)
                           tAAAATAGCGCCGTAACCCCCAGGTCCTGGATTGCCCAGACACGAACCATCGGTGAAAATTTCCACCTGTTTAAGCATCTCTGGTAGACTTcc                                                             <  1:1734387/93‑1 (MQ=255)
                            aaaaTAGCGCCGTAACCCCCAGGTCCTGGATTGCCCAGACACGAACCATCGGTTAAAATTTCCACCTGTTTAAGCATCTCTGGTAGACTTCCt                                                            <  1:7512092/93‑1 (MQ=255)
                             aaaTAGCGCCGTAACCCCCAGGTCCTGGATTGCCCAGACACGAACCATCGGTGAAAATTTCCACCTGTTTAAGCATCTCTGGTAGACTTCCTg                                                           >  2:2907520/1‑93 (MQ=255)
                             aaaTAGCGCCGTAACCCCCAGGTCCTGGATTGCCCAGACACGAACCATCGGTGAAAATTTCCACCTGTTTAAGCATCTCTGGTAGACTTCCTg                                                           >  2:1203319/1‑93 (MQ=255)
                              aaTAGCGCCGTAACCCCCAGGTCCTGGATTGCCCAGACACGAACCATCGGTGAAAATTTCCACCTGTTTAAGCATCTCTGGTAGACTTCCTGt                                                          >  1:6801929/1‑93 (MQ=255)
                                tAGCGCCGTAACCCCCAGGTCCTGGATTGCCCAGACACGAACCATCGGTGAAAATTTCCACCTGTTTAAGCATCTCTGGTAGACTTCCTGTaa                                                        >  2:6311737/1‑93 (MQ=255)
                                 aGCGCCGTAACCCCCAGGTCCTGGATTGCCCAGACACGAACCATCGGTGAAAATTTCCACCTGTTTAAGCATCTCTGGTAGACTTCCTGTAAt                                                       >  2:4681796/1‑93 (MQ=255)
                                 aGCGCCGTAACCCCCAGGTCCTGGATTGCCCAGACACGAACCATCGGTGAAAATTTCCACCTGTTTAAGCATCTCTGGTAGACTTCCTGTAAt                                                       >  1:4511540/1‑93 (MQ=255)
                                   cgcCGTAACCCCCAGGTCCTGGATTGCCCAGACACGAACCATCGGTGAAAATTTCCACCTGTTTAAGCATCTCTGGTAGACTTCCTGTAATTg                                                     <  1:1835170/93‑1 (MQ=255)
                                    gcCGTAACCCCCAGGTCCTGGATTGCCCAGACACGAACCATCGGTGAAAATTTCTACCTGTTTAAGCATCTCTGGTAGACTTCCTGTAATTGa                                                    >  1:2538834/1‑93 (MQ=255)
                                    gcCGTAACCCCCAGGTCCTGGATTGCCCAGACACGAACCATCGGTGAAAATTTCCACCTGTTTAAGCATCTCTGGTAGACTTCCTGTAATTGa                                                    >  2:3875378/1‑93 (MQ=255)
                                      cGTAACCCCCAGGTCCTGGATTGCCCAGACACGAACCATCGGTGAAAATTTCTACCTGTTTAAGCATCTCTGGTAGACTTCCTGTAATTGAAt                                                  <  2:366610/93‑1 (MQ=255)
                                       gTAACCCCCAGGTCCTGGATTGCCCAGACACGAACCATCGGTGAAAATTTCCACCTGTTTAAGCATCTCTGGTAGACTTCCTGTAATTGAATc                                                 <  2:6396220/93‑1 (MQ=255)
                                       gTAACCCCCAGGTCCTGGATTGCCCAGACACGAACCATCGGTGAAAATTTCCACCTGTTTAAGCATCTCTGGTAGACTTCCTGTAATTGAATc                                                 >  2:2309462/1‑93 (MQ=255)
                                        tAACCCCCAGGTCCTGGATTGCCCAGACACGAACCATCGGTGAAAATTTCCACCTGTTTAAGCATCTCTGGTAGACTTCCTGTAATTGAATCg                                                >  1:5552623/1‑93 (MQ=255)
                                         aaCCCCCAGGTCCTGGATTGCCCAGACACGAACCATCGGTGAAAATTTCCACCTGTTTAAGCATCTCTGGTAGACTTCCTGTAATTGAATCGa                                               <  1:5287467/93‑1 (MQ=255)
                                             cccAGGTCCTGGATTGCCCAGACACGAACCATCGGTGAAAATTTCCACCTGTTTAAGCATCTCTGGTAGACTTCCTGTAATTGAATCGAACTg                                           >  2:4756844/1‑93 (MQ=255)
                                             cccAGGTCCTGGATTGCCCAGACACGAACCATCGGTGAAAATTTCCACCTGTTTAAGCATCTCTGGTAGACTTCCTGTAATTGAATCGAACTg                                           >  2:7482726/1‑93 (MQ=255)
                                              ccAGGTCCTGGATTGCCCAGACACGAACCATCGGTGAAAATTTCCACCTGTTTAAGCATCTCTGGTAGACTTCCTGTAATTGAATCGAACTGt                                          >  2:1713153/1‑93 (MQ=255)
                                               cAGGTCCTGGATTGCCCAGACACGAACCATCGGTGAAAATTTCTACCTGTTTAAGCATCTCTGGTAGACTTCCTGTAATTGAATCGAACTGTa                                         >  1:6357797/1‑93 (MQ=255)
                                               cAGGTCCTGGATTGCCCAGACACGAACCATCGGTGAAAATTTCCACCTGTTTAAGCATCTCTGGTAGACTTCCTGTAATTGAATCGAACTGTa                                         <  1:1032006/93‑1 (MQ=255)
                                               cAGGTCCTGGATTGCCCAGACACGAACCATCGGTGAAAATTTCCACCTGTTTAAGCATCTCTGGTAGACTTCCTGTAATTGAATCGAACTGTa                                         <  1:3958716/93‑1 (MQ=255)
                                                aGGTCCTGGATTGCCCAGACACGAACCATCGGTGAAAATTTCCACCTGTTTAAGCATCTCTGGTAGACTTCCTGTAATTGAATCGAACTGTaa                                        >  2:3809958/1‑93 (MQ=255)
                                                aGGTCCTGGATTGCCCAGACACGAACCATCGGTGAAAATTTCCACCTGTTTAAGCATCTCTGGTAGACTTCCTGTAATTGAATCGAACTGTaa                                        <  1:6995823/93‑1 (MQ=255)
                                                 ggTCCTGGATTGCCCAGACACGAACCATCGGTGAAAATTTCCACCTGTTTAAGCATCTCTGGTAGACTTCCTGTAATTGAATCGAACTGTaaa                                       >  2:3766012/1‑93 (MQ=255)
                                                  gTCCTGGATTGCCCAGACACGAACCATCGGTGAAAATTTCCACCTGTTTAAGCATCTCTGGTAGACTTCCTGTAATTGAATCGAACTGTaaaa                                      <  2:6881823/93‑1 (MQ=255)
                                                  gTCCTGGATTGCCCAGACACGAACCATCGGTGAAAATTTCCACCTGTTTAAGCATCTCTGGTAGACTTCCTGTAATTGAATCGAACTGTaaaa                                      <  2:3980671/93‑1 (MQ=255)
                                                   tCCTGGATTGCCCAGACACGAACCATCGGTGAAAATTTCCACCTGTTTAAGCATCTCTGGTAGACTTCCTGTAATTGAATCGAACTGTAAAAc                                     >  1:3751916/1‑93 (MQ=255)
                                                     cTGGATTGCCCAGACACGAACCATCGGTGAAAATTTCCACCTGTTTAAGCATCTCTGGTAGACTTCCTGTAATTGAATCGAACTGTAAAACGa                                   >  1:5336272/1‑93 (MQ=255)
                                                     cTGGATTGCCCAGACACGAACCATCGGTGAAAATTTCCACCTGTTTAAGCATCTCTGGTAGACTTCCTGTAATTGAATCGAACTGTAAAACGa                                   >  2:6557763/1‑93 (MQ=255)
                                                       ggATTGCCCAGACACGAACCATCGGTGAAAATTTCCACCTGTTTAAGCATCTCTGGTAGACTTCCTGTAATTGAATCGAACTGTAAAACGACa                                 <  1:1535691/93‑1 (MQ=255)
                                                        gATTGCCCAGACACGAACCATCGGTGAAAATTTCTACCTGTTTAAGCATCTCTGGTAGACTTCCTGTAATTGAATCGAACTGTAAAACGACaa                                >  1:2170792/1‑93 (MQ=255)
                                                        gATTGCCCAGACACGAACCATCGGTGAAAATTTCCACCTGTTTAAGCATCTCTGGTAGACTTCCTGTAATTGAATCGAACTGTAAAACGACaa                                <  1:7356034/93‑1 (MQ=255)
                                                          ttGCCCAGACACGAACCATCGGTGAAAATTTCTACCTGTTTAAGCATCTCTGGTAGACTTCCTGTAATTGAATCGAACTGTAAAACGACAAGt                              <  2:1060812/93‑1 (MQ=255)
                                                          ttGCCCAGACACGAACCATCGGTGAAAATTTCCACCTGTTTAAGCATCTCTGGTAGACTTCCTGTAATTGAATCGAACTGTAAAACGACAAGt                              >  2:7138696/1‑93 (MQ=255)
                                                           tGCCCAGACACGAACCATCGGTGAAAATTTCCACCTGTTTAAGCATCTCTGGTAGACTTCCTGTAATTGAATCGAACTGTAAAACGACAAGTc                             <  1:362070/93‑1 (MQ=255)
                                                             cccAGACACGAACCATCGGTGAAAATTTCCACCTGTTTAAGCATCTCTGGTAGACTTCCTGTAATTGAATCGAACTGTAAAACGACAAGTCTg                           >  1:1465168/1‑93 (MQ=255)
                                                              ccAGACACGAACCATCGGTGAAAATTTCTACCTGTTTAAGCATCTCTGGTAGACTTCCTGTAATTGAATCGAACTGTAAAACGACAAGTCTGa                          <  1:4901431/93‑1 (MQ=255)
                                                              ccAGACACGAACCATCGGTGAAAATTTCCACCTGTTTAAGCATCTCTGGTAGACTTCCTGTAATTGAATCGAACTGTAAAACGACAAGTCTGa                          >  2:2567129/1‑93 (MQ=255)
                                                              ccAGACACGAACCATCGGTGAAAATTTCCACCTGTTTAAGCATCTCTGGTAGACTTCCTGTAATTGAATCGAACTGTAAAACGACAAGTCTGa                          <  1:81178/93‑1 (MQ=255)
                                                               cAGACACGAACCATCGGTGAAAATTTCCACCTGTTTAAGCATCTCTGGTAGACTTCCTGTAATTGAATCGAACTGTAAAACGACAAGTCTGAc                         <  2:1621493/93‑1 (MQ=255)
                                                                aGACACGAACCATCGGTGAAAATTTCTACCTGTTTAAGCATCTCTGGTAGACTTCCTGTAATTGAATCGAACTGTAAAACGACAAGTCTGACa                        >  1:2333256/1‑93 (MQ=255)
                                                                   cacGAACCATCGGTGAAAATTTCCACCTGTTTAAGCATCTCTGGTAGACTTCCTGTAATTGAATCGAACTGTAAAACGACAAGTCTGACATaa                     >  2:1243528/1‑93 (MQ=255)
                                                                    acGAACCATCGGTGAAAATTTCTACCTGTTTAAGCATCTCTGGTAGACTTCCTGTAATTGAATCGAACTGTAAAACGACAAGTCTGACATaaa                    >  2:207269/1‑93 (MQ=255)
                                                                    acGAACCATCGGTGAAAATTTCCACCTGTTTAAGCATCTCTGGTAGACTTCCTGTAATTGAATCGAACTGTAAAACGACAAGTCTGACATaaa                    <  2:6711559/93‑1 (MQ=255)
                                                                     cGAACCATCGGTGAAAATTTCCACCTGTTTAAGCATCTCTGGTAGACTTCCTGTAATTGAATCGAACTGTAAAACGACAAGTCTGACATAAAt                   <  1:3604852/93‑1 (MQ=255)
                                                                      gAACCATCGGTGAAAATTTCCACCTGTTTAAGCATCTCTGGTAGACTTCCTGTAATTGAATCGAACTGTAAAACGACAAGTCTGACATAAATg                  >  2:890251/1‑93 (MQ=255)
                                                                      gAACCATCGGTGAAAATTTCCACCTGTTTAAGCATCTCTGGTAGACTTCCTGTAATTGAATCGAACTGTAAAACGACAAGTCTGACATAAATg                  >  2:185829/1‑93 (MQ=255)
                                                                        aCCATCGGTGAAAATTTCCACCTGTTTAAGCATCTCTGGTAGACTTCCTGTAATTGAATCGAACTGTAAAACGACAAGTCTGACATAAATGAc                <  2:4135689/93‑1 (MQ=255)
                                                                        aCCATCGGTGAAAATTTCCACCTGTTTAAGCATCTCTGGTAGACTTCCTGTAATTGAATCGAACTGTAAAACGACAAGTCTGACATAAATGAc                >  2:6042502/1‑93 (MQ=255)
                                                                           aTCGGTGAAAATTTCCACCTGTTTAAGCATCTCTGGTAGACTTCCTGTAATTGAATCGAACTGTAAAACGACAAGTCTGACATAAATGACCGc             <  1:3469038/93‑1 (MQ=255)
                                                                            tCGGTGAAAATTTCCACCTGTTTAAGCATCTCTGGTAGACTTCCTGTAATTGAATCGAACTGTAAAACGACAAGTCTGACATAAATGACCGCt            <  2:7069393/93‑1 (MQ=255)
                                                                             cGGTGAAAATTTCCACCTGTTTAAGCATCTCTGGTAGACTTCCTGTAATTGAATCGAACTGTAAAACGACAAGTCTGACATAAATGACCGCTa           <  1:3140126/93‑1 (MQ=255)
                                                                              ggTGAAAATTTCTACCTGTTTAAGCATCTCTGGTAGACTTCCTGTAATTGAATCGAACTGTAAAACGACAAGTCTGACATAAATGACCGCTAt          >  2:4186057/1‑93 (MQ=255)
                                                                              ggTGAAAATTTCCACCTGTTTAAGCATCTCTGGTAGACTTCCTGTAATTGAATCGAACTGTAAAACGACAAGTCTGACATAAATGACCGCTAt          <  1:2770998/93‑1 (MQ=255)
                                                                              ggTGAAAATTTCCACCTGTTTAAGCATCTCTGGTAGACTTCCTGTAATTGAATCGAACTGTAAAACGACAAGTCTGACATAAATGACCGCTAt          <  2:5484786/93‑1 (MQ=255)
                                                                               gTGAAAATTTCCACCTGTTTAAGCATCTCTGGTAGACTTCCTGTAATTGAATCGAACTGTAAAACGACAAGTCTGACATAAATGACCGCTATg         >  2:1435442/1‑93 (MQ=255)
                                                                               gTGAAAATTTCCACCTGTTTAAGCATCTCTGGTAGACTTCCTGTAATTGAATCGAACTGTAAAACGACAAGTCTGACATAAATGACCGCTATg         >  1:2819896/1‑93 (MQ=255)
                                                                               gTGAAAATTTCCACCTGTTTAAGCATCTCTGGTAGACTTCCTGTAATTGAATCGAACTGTAAAACGACAAGTCTGACATAAATGACCGCTATg         >  2:472390/1‑93 (MQ=255)
                                                                                tGAAAATTTCCACCTGTTTAAGCATCTCTGGTAGACTTCCTGTAATTGAATCGAACTGTAAAACGACAAGTCTGACATAAATGACCGCTATGa        >  2:5824249/1‑93 (MQ=255)
                                                                                 gAAAATTTCTACCTGTTTAAGCATCTCTGGTAGACTTCCTGTAATTGAATCGAACTGTAAAACGACAAGTCTGACATAAATGACCGCTATGAg       >  1:509363/1‑93 (MQ=255)
                                                                                 gAAAATTTCCACCTGTTTAAGCATCTCTGGTAGACTTCCTGTAATTGAATCGAACTGTAAAACGACAAGTCTGACATAAATGACCGCTATGAg       >  2:3269374/1‑93 (MQ=255)
                                                                                  aaaaTTTCCACCTGTTTAAGCATCTCTGGTAGACTTCCTGTAATTGAATCGAACTGTAAAACGACAAGTCTGACATAAATGACCGCTATGAGc      <  1:4631150/93‑1 (MQ=255)
                                                                                  aaaaTTTCCACCTGTTTAAGCATCTCTGGTAGACTTCCTGTAATTGAATCGAACTGTAAAACGACAAGTCTGACATAAATGACCGCTATGAGc      <  1:4442090/93‑1 (MQ=255)
                                                                                  aaaaTTTCCACCTGTTTAAGCATCTCTGGTAGACTTCCTGTAATTGAATCGAACTGTAAAACGACAAGTCTGACATAAATGACCGCTATGAGc      <  2:5368975/93‑1 (MQ=255)
                                                                                   aaaTTTCCACCTGTTTAAGCATCTCTGGTAGACTTCCTGTAATTGAATCGAACTGTAAAACGACAAGTCTGACATAAATGACCGCTATGAGCa     <  1:6651550/93‑1 (MQ=255)
                                                                                   aaaTTTCCACCTGTTTAAGCATCTCTGGTAGACTTCCTGTAATTGAATCGAACTGTAAAACGACAAGTCTGACATAAATGACCGCTATGAGCa     <  2:4321752/93‑1 (MQ=255)
                                                                                    aaTTTCCACCTGTTTAAGCATCTCTGGTAGACTTCCTGTAATTGAATCGAACTGTAAAACGACAAGTCTGACATAAATGACCGCTATGAGCAc    <  2:551982/93‑1 (MQ=255)
                                                                                      tttCCACCTGTTTAAGCATCTCTGGTAGACTTCCTGTAATTGAATCGAACTGTAAAACGACAAGTCTGACATAAATGACCGCTATGAGCACTg  >  1:6935369/1‑93 (MQ=255)
                                                                                          |                                                                                        
AAAGGTTTTCTCGCGTCCGCGATAGCGTAAAATAGCGCCGTAACCCCCAGGTCCTGGATTGCCCAGACACGAACCATCGGTGAAAATTTCTACCTGTTTAAGCATCTCTGGTAGACTTCCTGTAATTGAATCGAACTGTAAAACGACAAGTCTGACATAAATGACCGCTATGAGCACTG  >  NC_000913/235898‑236076

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