Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,184,181 (T)5→6 intergenic (+200/‑22) ygiL → / → yi215 putative fimbrial‑like protein/IS2 hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,184,1761.T98.7% 232.3 / ‑3.2 78intergenic (+195/‑27)ygiL/yi215putative fimbrial‑like protein/IS2 hypothetical protein
Reads supporting (aligned to +/- strand):  ref base . (1/0);  new base T (55/22);  total (56/22)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.98e-01

Only 100 of 128 total aligned reads displayed.

AAAATATTTATTGTTCATTATCCGTTATTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGGAGAACA  >  NC_000913/3184084‑3184201
                                                                                             |                         
aaaaTATTTATTGTTCATTATCCGTTATTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAg                            >  1:4435207/1‑93 (MQ=255)
 aaaTATTTATTGTTCATTATCCGTTATTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGt                          <  1:4242220/93‑1 (MQ=255)
  aaTATTTATTGTTCATTATCCGTTATTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGtt                         <  2:7661287/93‑1 (MQ=255)
  aaTATTTATTGTTCATTATCCGTTATTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGtt                         >  2:1663574/1‑93 (MQ=255)
  aaTATTTATTGTTCATTATCCGTTATTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGtt                         >  1:7383879/1‑93 (MQ=255)
  aaTATTTATTGTTCATTATCCGTTATTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGtt                         <  2:4131677/93‑1 (MQ=255)
  aaTATTTATTGTTCATTATCCGTTATTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGtt                         >  2:128917/1‑93 (MQ=255)
  aaTATTTATTGTTCATTATCCGTTATTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGtt                         <  1:2701440/93‑1 (MQ=255)
   atatTTATTGTTCATTATCCGTTATTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGttt                        <  1:4748561/93‑1 (MQ=255)
   atatTTATTGTTCATTATCCGTTATTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGttt                        <  2:6104825/93‑1 (MQ=255)
   atatTTATTGTTCATTATCCGTTATTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGttt                        >  2:1209427/1‑93 (MQ=255)
   atatTTATTGTTCATTATCCGTTATTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGttt                        <  1:62850/93‑1 (MQ=255)
    tatttattGTTCATTATCCGTTATTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGtttt                       <  1:6737534/93‑1 (MQ=255)
    tatttattGTTCATTATCCGTTATTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGtttt                       >  2:6954096/1‑93 (MQ=255)
    tatttattGTTCATTATCCGTTATTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGtttt                       >  2:2135710/1‑93 (MQ=255)
     atttattGTTCATTATCCGTTATTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGttttt                      <  2:5732520/93‑1 (MQ=255)
     atttattGTTCATTATCCGTTATTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGttttt                      >  2:3814411/1‑93 (MQ=255)
     atttattGTTCATTATCCGTTATTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGttttt                      >  2:1823983/1‑93 (MQ=255)
      tttattGTTCATTATCCGTTATTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGtttttt                     >  1:4039923/1‑92 (MQ=255)
      tttattGTTCATTATCCGTTATTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGtttttt                     >  2:1522037/1‑92 (MQ=255)
      ttgattGTTCAGTATCCGGTATTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGtttttt                     <  1:7216032/90‑2 (MQ=14)
       ttattGTTCATTATCCGTTATTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGtttttta                    <  2:4930777/93‑3 (MQ=255)
       ttattGTTCATTATCCGTTATTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGtttttta                    >  2:1355782/1‑91 (MQ=255)
       ttattGTTCATTATCCGTTATTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGtttttta                    <  2:2744712/93‑3 (MQ=255)
       ttattGTTCATTATCCGTTATTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGtttttta                    <  2:5391620/93‑3 (MQ=255)
        tattGTTCATTATCCGTTATTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGttttttag                   <  2:311562/93‑4 (MQ=255)
        tattGTTCATTATCCGTTATTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGttttttag                   >  1:1253766/1‑90 (MQ=255)
        tattGTTCATTATCCGTTATTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGCGATAGTCTTACTACTAGttttttag                   >  1:2629100/1‑90 (MQ=255)
        tattGTTCATTATCCGTTATTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGCGATAGTCTTAATACTAGttttttag                   >  1:2629099/1‑90 (MQ=255)
         attGTTCATTATCCGTTATTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGttttttaga                  >  1:4868630/1‑89 (MQ=255)
         attGTTCATTATCCGTTATTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGttttttaga                  <  2:2912454/93‑5 (MQ=255)
         attGTTCATTATCCGTTATTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGttttttaga                  <  2:1712531/93‑5 (MQ=255)
         attGTTCATTATCCGTTATTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGttttttaga                  >  2:1139315/1‑89 (MQ=255)
         attGTTCATTATCCGTTATTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGttttttaga                  >  2:6592033/1‑89 (MQ=255)
         attGTTCATTATCCGTTATTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGttttttaga                  <  2:4347414/93‑5 (MQ=255)
         attGTTCATTATCCGTTATTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGttttttaga                  >  1:4386941/1‑89 (MQ=255)
         attGTTCATTATCCGTTATTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGttttttaga                  >  1:4868631/1‑89 (MQ=255)
         attGTTCATTATCCGTTATTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGttttttaga                  <  2:5131754/93‑5 (MQ=255)
         attGTTCATTATCCGTTATTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGttttttaga                  <  2:4976109/93‑5 (MQ=255)
         attGTTCATTATCCGTTAATAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGttttttaga                  <  1:3349395/93‑5 (MQ=255)
           tGTTCATTATCCGTTATTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTTAGACt                 <  1:1089921/93‑1 (MQ=255)
           tGTTCATTATCCGTTATTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTTAGACt                 >  2:1792874/1‑93 (MQ=255)
           tGTTCATTATCCGTTATTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTTAGACt                 >  2:1727747/1‑93 (MQ=255)
            gTTCATTATCCGTTATTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTTAGACTa                <  2:4591497/93‑1 (MQ=255)
            gTTCATTATCCGTTATTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTTAGACTa                >  2:5267787/1‑93 (MQ=255)
            gTTCATTATCCGTTATTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTTAGACTa                >  2:1863342/1‑93 (MQ=255)
            gTTCATTATCCGTTATTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTTAGACTa                >  1:7574990/1‑93 (MQ=255)
            gTTCATTATCCGTTATTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTTAGACTa                <  1:5396925/93‑1 (MQ=255)
             ttCATTATCCGTTATTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTTAGACTAg               >  2:363637/1‑93 (MQ=255)
             ttCATTATCCGTTATTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTTAGACTAg               >  1:2107486/1‑93 (MQ=255)
             ttCATTATCCGTTATTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTTAGACTAg               <  2:4168227/93‑1 (MQ=255)
             ttCATTATCCGTTATTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTTAGACTAg               >  1:405119/1‑93 (MQ=255)
             ttCATTATCCGTTATTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTTAGACTAg               >  2:6168780/1‑93 (MQ=255)
              tCATTATCCGTTATTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAGTAAGTGATAGTCTTAATACTAGTTTTTTAGACTAGt              <  1:4289217/93‑1 (MQ=255)
              tCATTATCCGTTATTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTTAGACTAGt              >  1:6057769/1‑93 (MQ=255)
               cATTATCCGTTATTAGACTGGCCCCCTGAATCTCCAGACAACCACTATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTTAGACTAGta             >  1:1675067/1‑92 (MQ=14)
               cATTATCCGTTATTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTTAGACTAGTc             >  2:3295277/1‑93 (MQ=255)
               cATTATCCGTTATTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTTAGACTAGTc             >  1:1675068/1‑93 (MQ=255)
                aTTATCCGTTATTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTTAGACTAGTCa            >  2:2263654/1‑93 (MQ=255)
                 ttATCCGTTATTAGACTGTCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTTAGACTAGTCAt           <  2:2115783/93‑1 (MQ=255)
                 ttATCCGTTATTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTTAGACTAGTCAt           >  2:686442/1‑93 (MQ=255)
                 ttATCCGTTATTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTTAGACTAGTCAt           <  2:4947786/93‑1 (MQ=255)
                  tATCCGTTATTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTTAGACTAGTCAtt          <  2:3142884/93‑1 (MQ=255)
                  tATCCGTTATTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTTAGACTAGTCAtt          <  1:1156159/93‑1 (MQ=255)
                  tATCCGTTATTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTTAGACTAGTCAtt          >  2:1669602/1‑93 (MQ=255)
                  tATCCGTTATTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTTAGACTAGTCAtt          >  1:7908940/1‑93 (MQ=255)
                   aTCCGTTATTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTTAGACTAGTCATTg         >  2:7504943/1‑93 (MQ=255)
                   aTCCGTTATTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTTAGACTAGTCATTg         >  2:2170359/1‑93 (MQ=255)
                   aTCCGTTATTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTTAGACTAGTCATTg         <  1:691596/93‑1 (MQ=255)
                   aTCCGTTATTAGACTGGCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTTAGACTAGTCATTgg        >  1:5001693/1‑93 (MQ=14)
                    tCCGTTATTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTTAGACTAGTCATTgg        >  1:1263165/1‑93 (MQ=255)
                    tCCGTTATTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTTAGACTAGTCATTgg        >  2:1485647/1‑93 (MQ=255)
                    tCCGTTATTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTTAGACTAGTCATTgg        <  1:7189798/93‑1 (MQ=255)
                    tCCGTTATTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTTAGACTAGTCATTgg        >  2:526120/1‑93 (MQ=255)
                    tCCGTTATTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTTAGACTAGTCATTgg        >  2:4712308/1‑93 (MQ=255)
                     ccGTTATTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTTAGACTAGTCATTGga       <  2:3129364/93‑1 (MQ=255)
                     ccGTTATTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTTAGACTAGTCATTGga       <  2:480973/93‑1 (MQ=255)
                      cGTTATTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTTAGACTAGTCATTGgag      >  2:6378837/1‑93 (MQ=255)
                      cGTTATTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTTAGACTAGTCATTGgag      >  1:3908529/1‑93 (MQ=255)
                       gTTATTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTTAGACTAGTCATTGgaga     >  2:2616784/1‑93 (MQ=255)
                       gTTATTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTTAGACTAGTCATTGgaga     <  2:2270618/93‑1 (MQ=255)
                       gTTATTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTTAGACTAGTCATTGgaga     >  1:7545049/1‑93 (MQ=255)
                        ttattaGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTTAGACTAGTCATTGGAGaa    >  2:4733749/1‑93 (MQ=16)
                        ttattaGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTTAGACTAGTCATTGGAGaa    <  2:807903/93‑1 (MQ=16)
                        ttattaGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTTAGACTAGTCATTGGAGaa    >  2:2744614/1‑93 (MQ=16)
                        ttattaGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTTAGACTAGTCATTGGAGaa    >  2:2054386/1‑93 (MQ=16)
                         tattaGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTTAGACTAGTCATTGGAGaa  >  2:4290486/1‑92 (MQ=12)
                         tattaGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTTAGACTAGTCATTGGAGAAc   >  2:4290485/1‑93 (MQ=12)
                         tattaGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTTAGACTAGTCATTGGAGAAc   <  1:4275589/93‑1 (MQ=12)
                         tattaGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTTAGACTAGTCATTGGAGAAc   <  1:5164129/93‑1 (MQ=12)
                          attaGACTGGCCCCCTGAATCTCCGGACAACCAATATCACTTCAATAAGTGATAGTCTTAATACTAGTTTTTTAGACTAGTCATTGGAGAACa  >  2:2929847/1‑93 (MQ=11)
                          attaGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTTAGACTAGTCATTGGAGAACa  >  1:505673/1‑93 (MQ=9)
                          attaGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTTAGACTAGTCATTGGAGAACa  >  1:5116463/1‑93 (MQ=9)
                          attaGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTTAGACTAGTCATTGGAGAACa  <  1:4801033/93‑1 (MQ=9)
                          attaGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTTAGACTAGTCATTGGAGAACa  <  1:4329120/93‑1 (MQ=9)
                          attaGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTTAGACTAGTCATTGGAGAACa  >  1:3734696/1‑93 (MQ=9)
                          attaGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTTAGACTAGTCATTGGAGAACa  >  1:3734695/1‑93 (MQ=9)
                          attaGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTTAGACTAGTCATTGGAGAACa  >  1:1773506/1‑93 (MQ=9)
                          attaGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTTAGACTAGTCATTGGAGAACa  >  1:5476441/1‑93 (MQ=9)
                          attaGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTTAGACTAGTCATTGGAGAACa  >  1:6977534/1‑93 (MQ=9)
                                                                                             |                         
AAAATATTTATTGTTCATTATCCGTTATTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGGAGAACA  >  NC_000913/3184084‑3184201

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