Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 857,347 A→G V313A (GTA→GCA)  ybiU ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913857,3470AG98.6% 751.4 / ‑1.0 216V313A (GTA→GCA) ybiUhypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (2/1);  new base G (114/99);  total (116/100)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.35e-01

Only 100 of 220 total aligned reads displayed.

ACGTCGCAGTGCCACCAGACGGAGTCTCCGGCTTCGAGTTTTGGAATGCTGGTTAACGCCTCAATCAACAGTGGATGCCATTGCTCTGATACCGGCAATACTCTTCCGGGCGCTACGCCGCACAGTTCATCCTCCGGCACATCATCAAGCAGCGGACGTAACAGTACGTACGCCATA  >  NC_000913/857257‑857433
                                                                                          |                                                                                      
aCGTCGCAGTGCCACCAGACGGAGTCTCCGGCTTCGAGTTTTGGAATGCTGGTTAACGCCTCAATCAACAGTGGATGCCATTGCTCTGATGcc                                                                                      <  1:3560435/93‑1 (MQ=255)
aCGTCGCAGTGCCACCAGACGGAGTCTCCGGCTTCGAGTTTTGGAATGCTGGTTAACGCCTCAATCAACAGTGGATGCCATTGCTCTGATGcc                                                                                      <  1:5398124/93‑1 (MQ=255)
aCGTCGCAGTGCCACCAGAAGGAGTCTCCGGCTTCGAGTTTTGGAATGCTGGTTAACGCCTCAATCAACAGTGGATGCCATTGCTCTGATGcc                                                                                      >  1:1933441/1‑93 (MQ=255)
  gTCGCAGTGCCACCAGACGGAGTCTCCGGCTTCGAGTTTTGGAATGCTGGTTAAAGCCTCAATCAACAGCGGATGCAATTGCTCTGATGCCgg                                                                                    >  1:2009644/1‑93 (MQ=255)
   tCGCAGTGCCACCAGACGGAGTCTCCGGCTTCGAGTTTTGGAATGCTGGTTAACGCCTCAATCAACAGTGGATGCCATTGCTCTGATGCCGGc                                                                                   <  2:2822570/93‑1 (MQ=255)
     gCAGTGCCACCAGACGGAGTCTCCGGCTTCGAGTTTTGGAATGCTGGTTAACGCCTCAATCAACAGTGGATGCCATTGCTCTGATGCCGGCaa                                                                                 >  2:1386838/1‑93 (MQ=255)
     gCAGTGCCACCAGACGGAGTCTCCGGCTTCGAGTTTTGGAATGCTGGTTAACGCCTCAATCAACAGTGGATGCCATTGCTCTGATGCCGGCaa                                                                                 <  1:3378827/93‑1 (MQ=255)
     gCAGTGCCACCAGACGGAGTCTCCGGCTTCGAGTTTTGGAATGCTGGTTAACGCCTCAATCAACAGTGGATGCCATTGCTCTGATGCCGGCaa                                                                                 >  1:652347/1‑93 (MQ=255)
      cAGTGCCACCAGACGGAGTCTCCGGCTTCGAGTTTTGGAATGCTGGTTAACGCCTCAATCAACAGTGGATGCCATTGCTCTGATGCCGGCAAt                                                                                <  1:1421644/93‑1 (MQ=255)
       aGTGCCACCAGACGGAGTCTCCGGCTTCGAGTTTTGGAATGCTGGTTAACGCCTCAATCAACAGTGGATGCCATTGCTCTGATGCCGGCAATa                                                                               <  1:1710490/93‑1 (MQ=255)
       aGTGCCACCAGACGGAGTCTCCGGCTTCGAGTTTTGGAATGCTGGTTAACGCCTCAATCAACAGTGGATGCCATTGCTCTGATGCCGGCAATa                                                                               <  1:477438/93‑1 (MQ=255)
         tGCCACCAGACGGAGTCTCCGGCTTCGAGTTTTGGAATGCTGGTTAACGCCTCAATCAACAGTGGATGCCATTGCTCTGATGCCGGCAATAct                                                                             <  2:3915444/93‑1 (MQ=255)
         tGCCACCAGACGGAGTCTCCGGCTTCGAGTTTTGGAATGCTGGTTAACGCCTCAATCAACAGTGGATGCCATTGCTCTGATGCCGGCAATAct                                                                             <  2:6083719/93‑1 (MQ=255)
          gCCACCAGACGGAGTCTCCGGCTTCGAGTTTTGGAATGCTGGTTAACGCCTCAATCAACAGTGGATGCCATTGCTCTGATGCCGGCAATActc                                                                            <  1:669930/93‑1 (MQ=255)
           ccaccaGACGGAGTCTCCGGCTTCGAGTTTTGGAATGCTGGTTAACGCCTCAATCAACAGTGGATGCCATTGCTCTGATGCCGGCAATActct                                                                           >  2:411262/1‑93 (MQ=255)
           ccaccaGACGGAGTCTCCGGCTTCGAGTTTTGGAATGCTGGTTAACGCCTCAATCAACAGTGGATGCCATTGCTCTGATGCCGGCAATActct                                                                           <  2:2611766/93‑1 (MQ=255)
             accaGACGGAGTCTCCGGCTTCGAGTTTTGGAATGCTGGTTAACGCCTCAATCAACAGTGGATGCCATTGCTCTGATGCCGGCAATACTCTTc                                                                         >  1:414469/1‑93 (MQ=255)
              ccaGACGGAGTCTCCGGCTTCGAGTTTTGGAATGCTGGTTAACGCCTCAATCAACAGTGGATGCCATTGCTCTGATGCCGGCAATACTCTTcc                                                                        <  2:1859957/93‑1 (MQ=255)
              ccaGACGGAGTCTCCGGCTTCGAGTTTTGGAATGCTGGTTAACGCCTCAATCAACAGTGGATGCCATTGCTCTGATGCCGGCAATACTCTTcc                                                                        <  2:3377182/93‑1 (MQ=255)
              ccaGACGGAGTCTCCGGCTTCGAGTTTTGGAATGCTGGTTAACGCCTCAATCAACAGTGGATGCCATTGCTCTGATGCCGGCAATACTCTTcc                                                                        >  2:3818944/1‑93 (MQ=255)
               caGACGGAGTCTCCGGCTTCGAGTTTTGGAATGCTGGTTAACGCCTCAATCAACAGTGGATGCCATTGCTCTGATGCCGGCAATACTCTTCCg                                                                       >  2:5736763/1‑93 (MQ=255)
                 gACGGAGTCTCCGGCTTCGAGTTTTGGAATGCTGGTTAACGCCTCAATCAACAGTGGATGCCATTGCTCTGATGCCGGCAATACTCTTCCggg                                                                     >  2:1471818/1‑93 (MQ=255)
                  aCGGAGTCTCCGGCTTCGAGTTTTGGAATGCTGGTTAACGCCTCAATCAACAGTGGATGCCATTGCTCTGATGCCGGCAATACTCTTCCGGgc                                                                    >  2:1189935/1‑93 (MQ=255)
                       gTCTCCGGCTTCGAGTTTTGGAATGCTGGTTAACGCCTCAATCAACAGTGGATGCCATTGCTCTGATGCCGGCAATACTCTTCCGGGCGCTAc                                                               >  2:2281677/1‑93 (MQ=255)
                       gTCTCCGGCTTCGAGTTTTGGAATGCTGGTTAACGCCTCAATCAACAGTGGATGCCATTGCTCTGATGCCGGCAATACTCTTCCGGGCGCTAc                                                               <  2:2759405/93‑1 (MQ=255)
                       gTCTCCGGCTTCGAGTTTTGGAATGCTGGTTAACGCCTCAATCAACAGTGGATGCCATTGCTCTGATGCCGGCAATACTCTTCCGGGCGCTAc                                                               <  1:2379102/93‑1 (MQ=255)
                       gTCTCCGGCTTCGAGTTTTGGAATGCTGGTTAACGCCTCAATCAACAGTGGATGCCATTGCTCTGATGCCGGCAATACTCTTCCGGGCGCTAc                                                               <  2:4604167/93‑1 (MQ=255)
                          tcCGGCTTCGAGTTTTGGAATGCTGGTTAACGCCTCAATCAACAGTGGATGCCATTGCTCTGATGCCGGCAATACTCTTCCGGGCGCTAcgcc                                                            >  2:2095030/1‑93 (MQ=255)
                          tcCGGCTTCGAGTTTTGGAATGCTGGTTAACGCCTCAATCAACAGTGGATGCCATTGCTCTGATGCCGGCAATACTCTTCCGGGCGCTAcgcc                                                            >  2:79779/1‑93 (MQ=255)
                          tcCGGCTTCGAGTTTTGGAATGCTGGTTAACGCCTCAATCAACAGTGGATGCCATTGCTCTGATGCCGGCAATACTCTTCCGGGCGCTAcgcc                                                            >  1:227188/1‑93 (MQ=255)
                             ggCTTCGAGTTTTGGAATGCTGGTTAACGCCTCAATCAACAGTGGATGCCATTGCTCTGATGCCGGCAATACTCTTCCGGGCGCTACGCCGca                                                         >  2:16935/1‑93 (MQ=255)
                             ggCTTCGAGTTTTGGAATGCTGGTTAACGCCTCAATCAACAGTGGATGCCATTGCTCTGATGCCGGCAATACTCTTCCGGGCGCTACGCCGca                                                         >  2:255472/1‑93 (MQ=255)
                              gCTTCGAGTTTTGGAATGCTGGTTAACGCCTCAATCAACAGTGGATGCCATTGCTCTGATGCCGGCAATACTCTTCCGGGCGCTACGCCGcac                                                        <  1:914715/93‑1 (MQ=255)
                              gCTTCGAGTTTTGGAATGCTGGTTAACGCCTCAATCAACAGTGGATGCCATTGCTCTGATGCCGGCAATACTCTTCCGGGCGCTACGCCGcac                                                        <  1:844850/93‑1 (MQ=255)
                              gCTTCGAGTTTTGGAATGCTGGTTAACGCCTCAATCAACAGTGGATGCCATTGCTCTGATGCCGGCAATACTCTTCCGGGCGCTACGCCGcac                                                        >  2:5289110/1‑93 (MQ=255)
                              gCTTCGAGTTTTGGAATGCTGGTTAACGCCTCAATCAACAGTGGATGCCATTGCTCTGATGCCGGCAATACTCTTCCGGGCGCTACGCCGcac                                                        >  2:2782503/1‑93 (MQ=255)
                               cTTCGAGTTTTGGAATGCTGGTTAACGCCTCAATCAACAGTGGATGCCATTGCTCTGATGCCGGCAATACTCTTCCGGGCGCTACGCCGcaca                                                       <  2:6105255/93‑1 (MQ=255)
                                ttCGAGTTTTGGAATGCTGGTTAACGCCTCAATCAACAGTGGATGCCATTGCTCTGATGCCGGCAATACTCTTCCGGGCGCTACGCCGCACAg                                                      >  2:4955783/1‑93 (MQ=255)
                                ttCGAGTTTTGGAATGCTGGTTAACGCCTCAATCAACAGTGGATGCCATTGCTCTGATGCCGGCAATACTCTTCCGGGCGCTACGCCGCACAg                                                      >  1:2161918/1‑93 (MQ=255)
                                ttCGAGTTTTGGAATGCTGGTTAACGCCTCAATCAACAGTGGATGCCATTGCTCTGATGCCGGCAATACTCTTCCGGGCGCTACGCCGCACAg                                                      >  1:244496/1‑93 (MQ=255)
                                   gAGTTTTGGAATGCTGGTTAACGCCTCAATCAACAGTGGATGCCATTGCTCTGATGCCGGCAATACTCTTCCGGGCGCTACGCCGCACAGTTc                                                   >  1:5514812/1‑93 (MQ=255)
                                      ttttGGAATGCTGGTTAACGCCTCAATCAACAGTGGATGCCATTGCTCTGATGCCGGCAATACTCTTCCGGGCGCTACGCCGCACAGTTCAtc                                                <  2:4440692/93‑1 (MQ=255)
                                      ttttGGAATGCTGGTTAACGCCTCAATCAACAGTGGATGCCATTGCTCTGATGCCGGCAATACTCTTCCGGGCGCTACGCCGCACAGTTCAtc                                                >  2:2100749/1‑93 (MQ=255)
                                       tttGGAATGCTGGTTAACGCCTCAATCAACAGTGGATGCCATTGCTCTGATGCCGGCAATACTCTTCCGGGCGCTACGCCGCACAGTTCAtcc                                               <  2:5889534/93‑1 (MQ=255)
                                        ttGGAATGCTGGTTAACGCCTCAATCAACAGTGGATGCCATTGCTCTGATGCCGGCAATACTCTTCCGGGCGCTACGCCGCACAGTTCAtcct                                              >  2:2241894/1‑93 (MQ=255)
                                         tGGAATGCTGGTTAACGCCTCAATCAACAGTGGATGCCATTGCTCTGATGCCGGCAATACTCTTCCGGGCGCTACGCCGCACAGTTCAtcctc                                             <  1:250128/93‑1 (MQ=255)
                                          ggAATGCTGGTTAACGCCTCAATCAACAGTGGATGCCATTGCTCTGATGCCGGCAATACTCTTCCGGGCGCTACGCCGCACAGTTCAtcctcc                                            >  2:1543140/1‑93 (MQ=255)
                                           gAATGCTGGTTAACGCCTCAATCAACAGTGGATGCCATTGCTCTGATGCCTGCAATACTCTTCCGGGCGCTACGCCGCACAGTTCATCCTTCg                                           >  1:5657113/1‑93 (MQ=255)
                                              tGCTGGTTAACGCCTCAATCAACAGTGGATGCCATTGCTCTGATGCCGGCAATACTCTTCCGGGCGCTACGCCGCACAGTTCATCCTCCGGcc                                        >  1:3078503/1‑92 (MQ=255)
                                              tGCTGGTTAACGCCTCAATCAACAGTGGATGCCATTGCTCTGATGCCGGCAATACTCTTCCGGGCGCTACGCCGCACAGTTCATCCTCCGGca                                        >  2:3366373/1‑93 (MQ=255)
                                                 tGGTTAACGCCTCAATCAACAGTGGATGCCATTGCTCTGATGCCGGCAATACTCTTCCGGGCGCTACGCCGCACAGTTCATCCTCCGGCAcat                                     <  1:2534714/93‑1 (MQ=255)
                                                  ggTTAACGCCTCAATCAACAGTGGATGCCATTGCTCTGATGCCGGCAATACTCTTCCGGGCGCTACGCCGCACAGTTCATCCTCCGGCAcatc                                    <  1:717253/93‑1 (MQ=255)
                                                  ggTTAACGCCTCAATCAACAGTGGATGCCATTGCTCTGATGCCGGCAATACTCTTCCGGGCGCTACGCCGCACAGTTCATCCTCCGGCAcatc                                    >  1:3235274/1‑93 (MQ=255)
                                                    ttAACGCCTCAATCAACAGTGGATGCCATTGCTCTGATGCCGGCAATACTCTTCCGGGCGCTACGCCGCACAGTTCATCCTCCGGCAcatcat                                  <  1:2052960/93‑1 (MQ=255)
                                                    ttAACGCCTCAATCAACAGTGGATGCCATTGCTCTGATGCCGGCAATACTCTTCCGGGCGCTACGCCGCACAGTTCATCCTCCGGCAcatcat                                  >  2:5517503/1‑93 (MQ=255)
                                                    ttAACGCCTCAATCAACAGTGGATGCCATTGCTCTGATGCCGGCAATACTCTTCCGGGCGCTACGCCGCACAGTTCATCCTACGGCAcatcat                                  >  1:1119194/1‑93 (MQ=255)
                                                    gtAACGCCTCAATCAACAGTGGATGCCATTGCTCTGATGCCGGCAATTCTCTTCCGGGCGCTACGCCGCACAGTTCATCCTCCGGCAcatcat                                  <  2:1836622/92‑1 (MQ=255)
                                                     tAACGCCTCAATCAACAGTGGATGCCATTGCTCTGATGCCGGCAATACTCTTCCGGGCGCTACGCCGCACAGTTCATCCTCCGGCACatcatc                                 >  1:2929872/1‑93 (MQ=255)
                                                        cGCCTCAATCAACAGTGGATGCCATTGCTCTGATGCCGGCAATACTCTTCCGGGCGCTACGCCGCACAGTTCATCCTCCGGCACATCATCAag                              <  2:2798231/93‑1 (MQ=255)
                                                        cGCCTCAATCAACAGTGGATGCCATTGCTCTGATGCCGGCAATACTCTTCCGGGCGCTACGCCGCACAGTTCATCCTCCGGCACATCATCAag                              <  2:3681398/93‑1 (MQ=255)
                                                        cGCCTCAATCAACAGTGGATGCCATTGCTCTGATGCCGGCAATACTCTTCCGGGCGCTACGCCGCACAGTTCATCCTACGGCACATCATCaac                              >  2:649760/1‑92 (MQ=255)
                                                          ccTCAATCAACAGTGGATGCCATTGCTCTGATGCCGGCAATACTCTTCCGGGCGCTACGCCGCACAGTTCATCCTCCGGCACATCATCAagca                            >  2:2044608/1‑93 (MQ=255)
                                                          ccTCAATCAACAGTGGATGCCATTGCTCTGATGCCGGCAATACTCTTCCGGGCGCTACGCCGCACAGTTCATCCTCCGGCACATCATCAagca                            >  2:553737/1‑93 (MQ=255)
                                                          ccTCAATCAACAGTGGATGCCATTGCTCTGATGCCGGCAATACTCTTCCGGGCGCTACGCCGCACAGTTCATCCTCCGGCACATCATCAagca                            <  1:688448/93‑1 (MQ=255)
                                                          ccTCAATCAACAGTGGATGCCATTGCTCTGATGCCGGCAATACTCTTCCGGGCGCTACGCCGCACAGTTCATCCTCCGGCACATCATCAagca                            <  1:4913553/93‑1 (MQ=255)
                                                          ccTCAATCAACAGTGGATGCCATTGCTCTGATGCCGGCAATACTCTTCCGGGCGCTACGCCGCACAGTTCATCCTCCGGCACATCATCAagaa                            >  1:5159023/1‑91 (MQ=255)
                                                           cTCAATCAACAGTGGATGCCATTGCTCTGATGCCGGCAATACTCTTCCGGGCGCTACGCCGCACAGTTCATCCTCCGGCACATCATCAagcag                           >  2:5195858/1‑93 (MQ=255)
                                                           cTCAATCAACAGTGGATGCCATTGCTCTGATGCCGGCAATACTCTTCCGGGCGCTACGCCGCACAGTTCATCCTCCGGCACATCATCAagcag                           >  1:6025319/1‑93 (MQ=255)
                                                            tcaatcaaCAGTGGATGCCATTGCTCTGATGCCGGCAATACTCTTCCGTGCGCTACGCCGCACAGTTCATCCTCCGGCACATCATCAagcagc                          >  2:2855225/1‑93 (MQ=255)
                                                            tcaatcaaCAGTGGATGCCATTGCTCTGATGCCGGCAATACTCTTCCGGGCGCTACGCCGCACAGTTCATCCTCCGGCACATCATCAagcagc                          <  1:4016655/93‑1 (MQ=255)
                                                            tcaatcaaCAGTGGATGCCATTGCTCTGATGCCGGCAATACTCTTCCGGGCGCTACGCCGCACAGTTCATCCTCCGGCACATCATCAagcagc                          >  1:606205/1‑93 (MQ=255)
                                                            tcaatcaaCAGTGGATGCCATTGCTCTGATACCGGCAATACTCTTCCGGGCGCTACGCCGCACAGTTCATCCTCCGGCACATCATCAagcagc                          <  1:2729561/93‑1 (MQ=255)
                                                               atcaaCAGTGGATGCCATTGCTCTGATGCCGGCAATACTCTTCCGGGCGCTACGCCGCACAGTTCATCCTCCGGCACATCATCAAGCAGCGGa                       <  1:5763517/93‑1 (MQ=255)
                                                               atcaaCAGTGGATGCCATTGCTCTGATGCCGGCAATACTCTTCCGGGCGCTACGCCGCACAGTTCATCCTCCGGCACATCATCAAGCAGCGGa                       <  2:3469422/93‑1 (MQ=255)
                                                               atcaaCAGTGGATGCCATTGCTCTGATGCCGGCAATACTCTTCCGGGCGCTACGCCGCACAGTTCATCCTCCGGCACATCATCAAGCAGCGGa                       <  1:1240878/93‑1 (MQ=255)
                                                               atcaaCAGTGGATGCCATTGCTCTGATGCCGGCAATACTCTTCCGGGCGCTACGCCGCACAGTTCATCCTCCGGCACATCATCAAGCAGCGGa                       >  2:5036418/1‑93 (MQ=255)
                                                               atcaaCAGTGGATGCCATTGCTCTGATGCCGGCAATACTCTTCCGGGCGCTACGCCGCACAGTTCATCCTCCGGCACATCATCAAGCAGCGGa                       >  1:4527894/1‑93 (MQ=255)
                                                                tcaaCAGGGGATGCCATTGCTCTGATGCCGGCAATACTCTTCCGGGCGCTACGCCGCACAGTTCATCCTCCGGCACATCATCAAGCAGCGGAc                      <  1:5983190/93‑1 (MQ=255)
                                                                   aCAGTGGATGCCATTGCTCTGATGCCGGCAATACTCTTCCGGGCGCTACGCCGCACAGTTCATCCTCCGGCACATCATCAAGCAGCGGACGTa                   <  2:3141889/93‑1 (MQ=255)
                                                                   aCAGTGGATGCCATTGCTCTGATGCCGGCAATACTCTTCCGGGCGCTACGCCGCACAGTTCATCCTCCGGCACATCATCAAGCAGCGGACGTa                   >  1:4066515/1‑93 (MQ=255)
                                                                   aCAGTGGATGCCATTGCTCTGATGCCAGCAATACTCTTCCGGGCGCTACGCCGCACAGTTCATCCTCCGGCACATCATCAAGCAGCGGACGTa                   >  2:2631574/1‑93 (MQ=255)
                                                                    cAGTGGATGCCATTGCTCTGATGCCGGCAATACTCTTCCGGGCGCTACGCCGCACAGTTCATCCTCCGGCACATCATCAAGCAGCGGACGTaa                  <  1:4684241/93‑1 (MQ=255)
                                                                      gTGGATGCCATTGCTCTGATGCCGGCAATACTCTTCCGGGCGCTACGCCGCACAGTTCATCCTCCGGCACATCATCAAGCAGCGGACGTAACa                >  1:1062737/1‑93 (MQ=255)
                                                                      gTGGATGCCATTGCTCTGATGCCGGCAATACTCTTCCGGGCGCTACGCCGCACAGTTCATCCTCCGGCACATCATCAAGCAGCGGACGTAACa                <  2:3953090/93‑1 (MQ=255)
                                                                       tGGATGCCATTGCTCTGATGCCGGCAATACTCTTCCGGGCGCTACGCCGCACAGTTCATCCTCCGGCACATCATCAAGCAGCGGACGTAACAg               >  1:2813851/1‑93 (MQ=255)
                                                                       tGGATGCCATTGCTCTGATGCCGGCAATACTCTTCCGGGCGCTACGCCGCACAGTTCATCCTCCGGCACATCATCAAGCAGCGGACGTAACAg               >  1:1477716/1‑93 (MQ=255)
                                                                       gggATGCCATTGCTCTGATGCCGGCAATACTCTTCCGGGCGCTACGCCGCACAGTTCATCCTCCGGCACATCATCAAGCAGCGGACGTAACAg               <  2:3219978/92‑1 (MQ=255)
                                                                         gATGCCATTGCTCTGATGCCGGCAATACTCTTCCGGGCGCTACGCCGCACAGTTCATCCTCCGGCACATCATCAAGCAGCGGACGTAACAgta             <  1:3821249/93‑1 (MQ=255)
                                                                            gCCATTGCTCTGATGCCGGCAATACTCTTCCGGGCGCTACGCCGCACAGTTCATCCTCCGGCACATCATCAAGCAGCGGACGTAACAgtacgt          >  2:4006116/1‑93 (MQ=255)
                                                                            gCCATTGCTCTGATGCCGGCAATACTCTTCCGGGCGCTACGCCGCACAGTTCATCCTCCGGCACATCATCAAGCAGCGGACGTAACAgtacgt          >  2:4376175/1‑93 (MQ=255)
                                                                            gCCATTGCTCTGATGCCGGCAATACTCTTCCGGGCGCTACGCCGCACAGTTCATCCTCCGGCACATCATCAAGCAGCGGACGTAACAgtacgt          <  2:77973/93‑1 (MQ=255)
                                                                             ccATTGCTCTGATGCCGGCAATACTCTTCCGGGCGCTACGCCGCACAGTTCATCCTCCGGCACATCATCAAGCAGCGGACGTAACAgtacgtc         >  2:5650913/1‑92 (MQ=255)
                                                                             ccATTGCTCTGATGCCGGCAATACTCTTCCGGGCGCTACGCCGCACAGTTCATCCTCCGGCACATCATCAAGCAGCGGACGTAACAgtacgta         <  1:165079/93‑1 (MQ=255)
                                                                              cATTGCTCTGATGCCGGCAATTCTCGTCCGGGCGCTACGCCGCACAGTTCATCCTCCGGCACATCATCAAGCAGCGGACGTAACAgtacgtac        <  2:4633840/93‑1 (MQ=255)
                                                                              cATTGCTCTGATGCCGGCAATACTCTTCCGGGCGCTACGCCGCACAGTTCATCCTCCGGCACATCATCAAGCAGCGGACGTAACAgtacgtac        >  1:330748/1‑93 (MQ=255)
                                                                                  gCTCTGATGCCGGCAATACTCTTCCGGGCGCTACGCCGCACAGTTCATCCTCCGGCACATCATCAAGCAGCGGACGTAACAGTACGTACGCCa    >  1:791679/1‑93 (MQ=255)
                                                                                   ctctGATGCCGGCAATACTCTTCCGGGCGCTACGCCGCACAGTTCATCCTCCGGCACATCATCAAGCAGCGGACGTAACAGTACGTACGCCAt   <  1:1829078/93‑1 (MQ=255)
                                                                                   ctctGATGCCGGCAATACTCTTCCGGGCGCTACGCCGCACAGTTCATCCTCCGGCACATCATCAAGCAGCGGACGTAACAGTACGTACGCCAt   >  1:5551067/1‑93 (MQ=255)
                                                                                   ctctGATACCGGCAATACTCTTCCGGGCGCTACGCCGCACAGTTCATCCTCCGGCACATCATCAAGCAGCGGACGTAACAGTACGTACGCCAt   >  2:5344663/1‑93 (MQ=255)
                                                                                    tctGATGCCGGCAATACTCTTCCGGGCGCTACGCCGCACAGTTCATCCTCCGGCACATCATCAAGCAGCGGACGTAACAGTACGTACGCCATa  >  1:3680777/1‑93 (MQ=255)
                                                                                          |                                                                                      
ACGTCGCAGTGCCACCAGACGGAGTCTCCGGCTTCGAGTTTTGGAATGCTGGTTAACGCCTCAATCAACAGTGGATGCCATTGCTCTGATACCGGCAATACTCTTCCGGGCGCTACGCCGCACAGTTCATCCTCCGGCACATCATCAAGCAGCGGACGTAACAGTACGTACGCCATA  >  NC_000913/857257‑857433

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