Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,186,925 A→G T504A (ACT→GCT)  yqiG → putative membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,186,9250AG83.0% 648.7 / 129.4 294T504A (ACT→GCT) yqiGputative membrane protein
Reads supporting (aligned to +/- strand):  ref base A (31/19);  new base G (139/105);  total (170/124)
Fisher's exact test for biased strand distribution p-value = 5.34e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.38e-01

Only 100 of 311 total aligned reads displayed.

TATCAGAATGAAAAAGAGAGTTATATTGTCACGTTTAACCAGTATCTTGAATCATTAAGGTTAAATACATACGTAAGTTTGGCTCGTAATACTTACTGGGACGCCAGCAGTAATGTGAATTATTCATTATCACTTAGCCGCGATTTTGATATCGGGCCATTAAAAAACGTCTCCACTTCTC  >  NC_000913/3186835‑3187015
                                                                                          |                                                                                          
tATCAGAATGAAAAAGAGAGTTATATTGTCACGTTTAACCAGTATCTTGAATCATTAAGGTTAAATACATACGTAAGTTTGGCTCGTAAtgct                                                                                          <  2:6358485/93‑1 (MQ=255)
tATCAGAATGAAAAAGAGAGTTATATTGTCACGTTTAACCAGTATCTTGAATCATTAAGGTTAAATACATACGTAAGTTTGGCTCGTAAtgct                                                                                          >  2:4144848/1‑93 (MQ=255)
  tCAGAATGAAAAAGAGAGTTATATTGTCACGTTTAACCAGTATCTTGAATCATTAAGGTTAAATACATACGTAAGTTTGGCTCGTAAtgctta                                                                                        <  2:6844523/93‑1 (MQ=255)
   cAGAATGAAAAAGAGAGTTATATTGTCACGTTTAACCAGTATCTTGAATCATTAAGGTTAAATACATACGTAAGTTTGGCTCGTAAtgcttac                                                                                       <  2:4827880/93‑1 (MQ=255)
    aGAATGAAAAAGAGAGTTATATTGTCACGTTTAACCAGTATCTTGAATCATTAAGGTTAAATACATACGTAAGTTTGGCTCGTAAtgcttact                                                                                      >  2:5648958/1‑93 (MQ=255)
    aGAATGAAAAAGAGAGTTATATTGTCACGTTTAACCAGTATCTTGAATCATTAAGGTTAAATACATACGTAAGTTTGGCTCGTAAtgcttact                                                                                      >  1:3716401/1‑93 (MQ=255)
    aGAATGAAAAAGAGAGTTATATTGTCACGTTTAACCAGTATCTTGAATCATTAAGGTTAAATACATACGTAAGTTTGGCTCGTAAtgcttact                                                                                      >  1:3802150/1‑93 (MQ=255)
     gAATGAAAAAGAGAGTTATATTGTCACGTTTAACCAGTATCTTGAATCATTAAGGTTAAATACATACGTAAGTTTGGCTCGTAATGCTTACTg                                                                                     >  2:5265028/1‑93 (MQ=255)
     gAATGAAAAAGAGAGTTATATTGTCACGTTTAACCAGTATCTTGAATCATTAAGGTTAAATACATACGTAAGTTTGGCTCGTAATGCTTACTg                                                                                     >  1:4765721/1‑93 (MQ=255)
        tGAAAAAGAGAGTTATATTGTCACGTTTAACCAGTATCTTGAATCATTAAGGTTAAATACATACGTAAGTTTGGCTCGTAATGCTTACTGGGa                                                                                  <  2:862272/93‑1 (MQ=255)
        tGAAAAAGAGAGTTATATTGTCACGTTTAACCAGTATCTTGAATCATTAAGGTTAAATACATACGTAAGTTTGGCTCGTAATACTTACTGGGa                                                                                  <  2:2810791/93‑1 (MQ=255)
            aaaGAGAGTTATATTGTCACGTTTAACCAGTATCTTGAATCATTAAGGTTAAATACATACGTAAGTTTGGCTCGTAATGCTTACTGGGACGcc                                                                              <  2:2483202/93‑1 (MQ=255)
               gagagTTATATTGTCACGTTTAACCAGTATCTTGAATCATTAAGGTTAAATACATACGTAAGTTTGGCTCGTAATGCTTACTGGGACGCcagc                                                                           <  1:16707/93‑1 (MQ=255)
               gagagTTATATTGTCACGTTTAACCAGTATCTTGAATCATTAAGGTTAAATACATACGTAAGTTTGGCTCGTAATGCTTACTGGGACGCcagc                                                                           >  1:1000350/1‑93 (MQ=255)
               gagagTTATATTGTCACGTTTAACCAGTATCTTGAATCATTAAGGTTAAATACATACGTAAGTTTGGCTCGTAATGCTTACTGGGACGCcagc                                                                           >  1:2893414/1‑93 (MQ=255)
                agagTTATATTGTCACGTTTAACCAGTATCTTGAATCATTAAGGTTAAATACATACGTAAGTTTGGCTCGTAATGCTTACTGGGACGCcagca                                                                          >  2:5784318/1‑93 (MQ=255)
                agagTTATATTGTCACGTTTAACCAGTATCTTGAATCATTAAGGTTAAATACATACGTAAGTTTGGCTCGTAATACTTACTGGGACGCcagca                                                                          <  2:7297752/93‑1 (MQ=255)
                   gTTATATTGTCACGTTTAACCAGTATCTTGAATCATTAAGGTTAAATACATACGTAAGTTTGGCTCGTAATGCTTACTGGGACGCCAGCAGTa                                                                       >  1:3632247/1‑93 (MQ=255)
                    ttATATTGTCACGTTTAACCAGTATCTTGAATCATTAAGGTTAAATACATACGTAAGTTTGGCTCGTAATGCTTACTGGGACGCCAGCAGTaa                                                                      <  2:1848105/93‑1 (MQ=255)
                    ttATATTGTCACGTTTAACCAGTATCTTGAATCATTAAGGTTAAATACATACGTAAGTTTGGCTCGTAATGCTTACTGGGACGCCAGCAGTaa                                                                      >  2:4956537/1‑93 (MQ=255)
                     tataTTGTCACGTTTAACCAGTATCTTGAATCATTAAGGTTAAATACATACGTAAGTTTGGCTCGTAATGCTTACTGGGACGCCAGCAGTAAt                                                                     <  1:1429425/93‑1 (MQ=255)
                     tataTTGTCACGTTTAACCAGTATCTTGAATCATTAAGGTTAAATACATACGTAAGTTTGGCTCGTAATGCTTACTGGGACGCCAGCAGTAAt                                                                     >  1:2084085/1‑93 (MQ=255)
                     tataTTGTCACGTTTAACCAGTATCTTGAATCATTAAGGTTAAATACATACGTAAGTTTGGCTCGTAATGCTTACTGGGACGCCAGCAGTAAt                                                                     >  2:2639815/1‑93 (MQ=255)
                       tatTGTCACGTTTAACCAGTATCTTGAATCATTAAGGTTAAATACATACGTAAGTTTGGCTCGTAATGCTTACTGGGACGCCAGCAGTAAtgt                                                                   >  1:5563203/1‑93 (MQ=255)
                       tatTGTCACGTTTAACCAGTATCTTGAATCATTAAGGTTAAATACATACGTAAGTTTGGCTCGTAATGCTTACTGGGACGCCAGCAGTAAtgt                                                                   <  1:4438851/93‑1 (MQ=255)
                       tatTGTCACGTTTAACCAGTATCTTGAATCATTAAGGTTAAATACATACGTAAGTTTGGCTCGTAATGCTTACTGGGACGCCAGCAGTAAtgt                                                                   <  1:518012/93‑1 (MQ=255)
                       tatTGTCACGTTTAACCAGTATCTTGAATCATTAAGGTTAAATACATACGTAAGTTTGGCTCGTAATGCTTACTGGGACGCCAGCAGTAAtgt                                                                   <  2:4338645/93‑1 (MQ=255)
                        atTGTCACGTTTAACCAGTATCTTGAATCATTAAGGTTAAATACATACGTAAGTTTGGCTCGTAATGCTTACTGGGACGCCAGCAGTAAtgtg                                                                  >  2:6552858/1‑93 (MQ=255)
                         ttGTCACGTTTAACCAGTATCTTGAATCATTAAGGTTAAATACATACGTAAGTTTGGCTCGTAATACTTACTGGGACGCCAGCAGTAATGTGa                                                                 >  2:3354576/1‑93 (MQ=255)
                          tGTCACGTTTAACCAGTATCTTGAATCATTAAGGTTAAATACATACGTAAGTTTGGCTCGTAATGCTTACTGGGACGCCAGCAGTAATGTGaa                                                                >  1:7050058/1‑93 (MQ=255)
                           gTCACGTTTAACCAGTATCTTGAATCATTAAGGTTAAATACATACGTAAGTTTGGCTCGTAATGCTTACTGGGACGCCAGCAGTAATGTGAat                                                               <  2:4919063/93‑1 (MQ=255)
                           gTCACGTTTAACCAGTATCTTGAATCATTAAGGTTAAATACATACGTAAGTTTGGCTCGTAATGCTTACTGGGACGCCAGCAGTAATGTGAat                                                               <  2:2729043/93‑1 (MQ=255)
                           gTCACGTTTAACCAGTATCTTGAATCATTAAGGTTAAATACATACGTAAGTTTGGCTCGTAATGCTTACTGGGACGCCAGCAGTAATGTGAat                                                               <  1:4102942/93‑1 (MQ=255)
                               cGTTTAACCAGTATCTTGAATCATTAAGGTTAAATACATACGTAAGTTTGGCTCGTAATGCTTACTGGGACGCCAGCAGTAATGTGAattatt                                                           >  2:6785463/1‑93 (MQ=255)
                                gTTTAACCAGTATCTTGAATCATTAAGGTTAAATACATACGTAAGTTTGGCTCGTAATGCTTACTGGGACGCCAGCAGTAATGTGAATTATTc                                                          <  1:6966462/93‑1 (MQ=255)
                                gTTTAACCAGTATCTTGAATCATTAAGGTTAAATACATACGTAAGTTTGGCTCGTAATGCTTACTGGGACGCCAGCAGTAATGTGAATTATTc                                                          <  2:5060267/93‑1 (MQ=255)
                                 tttAACCAGTATCTTGAATCATTAAGGTTAAATACATACGTAAGTTTGGCTCGTAATGCTTACTGGGACGCCAGCAGTCATGTGAATTATTCa                                                         <  2:7485580/93‑1 (MQ=255)
                                 tttAACCAGTATCTTGAATCATTAAGGTTAAATACATACGTAAGTTTGGCTCGTAATGCTTACTGGGACGCCAGCAGTAATGTGAATTATTCa                                                         <  2:1937280/93‑1 (MQ=255)
                                 tttAACCAGTATCTTGAATCATTAAGGTTAAATACATACGTAAGTTTGGCTCGTAATGCTTACTGGGACGCCAGCAGTAATGTGAATTATTCa                                                         >  2:3886401/1‑93 (MQ=255)
                                  ttAACCAGTATCTTGAATCATTAAGGTTAAATACATACGTAAGTTTGGCTCGTAATGCTTACTGGGACGCCAGCAGTAATGTGAATTATTCAt                                                        <  1:5458271/93‑1 (MQ=255)
                                    aaCCAGTATCTTGAATCATTAAGGTTAAATACATACGTAAGTTTGGCTCGTAATGCTTACTGGGACGCCAGCAGTAATGTGAATTATTCATTa                                                      >  1:2440498/1‑93 (MQ=255)
                                     aCCAGTATCTTGAATCATTAAGGTTAAATACATACGTAAGTTTGGCTCGTAATGCTTACTGGGACGCCAGCAGTAATGTGAATTATTCATTAt                                                     <  2:1761495/93‑1 (MQ=255)
                                      ccAGTATCTTGAATCATTAAGGTTAAATACATACGTAAGTTTGGCTCGTAATGCTTACTGGGACGCCAGCAGTAATGTGAATTATTCATTATc                                                    <  1:5311540/93‑1 (MQ=255)
                                      ccAGTATCTTGAATCATTAAGGTTAAATACATACGTAAGTTTGGCTCGTAATACTTACTGGGACGCCAGCAGTAATGTGAATTATTCATTATc                                                    >  1:3020449/1‑93 (MQ=255)
                                       cAGTATCTTGAATCATTAAGGTTAAATACATACGTAAGTTTGGCTCGTAATGCTTACTGGGACGCCAGCAGTAATGTGAATTATTCATTATCa                                                   <  1:4619060/93‑1 (MQ=255)
                                       cAGTATCTTGAATCATTAAGGTTAAATACATACGTAAGTTTGGCTCGTAATGCTTACTGGGACGCCAGCAGTAATGTGAATTATTCATTATCa                                                   <  1:4177543/93‑1 (MQ=255)
                                       cAGTATCTTGAATCATTAAGGTTAAATACATACGTAAGTTTGGCTCGTAATGCTTACTGGGACGCCAGCAGTAATGTGAATTATTCATTATCa                                                   <  2:2395724/93‑1 (MQ=255)
                                        aGTATCTTGAATCATTAAGGTTAAATACATACGTAAGTTTGGCTCGTAATGCTTACTGGGACGCCAGCAGTAATGTGAATTATTCATTATCAc                                                  >  1:7162520/1‑93 (MQ=255)
                                        aGTATCTTGAATCATTAAGGTTAAATACATACGTAAGTTTGGCTCGTAATGCTTACTGGGACGCCAGCAGTAATGTGAATTATTCATTATCAc                                                  >  2:6040899/1‑93 (MQ=255)
                                          tATCTTGAATCATTAAGGTTAAATACATACGTAAGTTTGGCTCGTAATGCTTACTGGGACGCCAGCAGTAATGTGAATTATTCATTATCACtt                                                <  1:6498083/93‑1 (MQ=255)
                                          tATCTTGAATCATTAAGGTTAAATACATACGTAAGTTTGGCTCGTAATGCTTACTGGGACGCCAGCAGTAATGTGAATTATTCATTATCACtt                                                >  1:7422783/1‑93 (MQ=255)
                                          tATCTTGAATCATTAAGGTTAAATACATACGTAAGTTTGGCTCGTAATACTTACTGGGACGCCAGCAGTAATGTGAATTATTCATTATCACtt                                                >  2:2139759/1‑93 (MQ=255)
                                           aTCTTGAATCATTAAGGTTAAATACATACGTAAGTTTGGCTCGTAATACTTACTGGGACGCCAGCAGTAATGTGAATTATTCATTATCACTTa                                               >  2:1397551/1‑93 (MQ=255)
                                            tCTTGAATCATTAAGGTTAAATACATACGTAAGTTTGGCTCGTAATGCTTACTGGGACGCCAGCAGTAATGTGAATTATTCATTATCACTTAg                                              >  2:1550834/1‑93 (MQ=255)
                                            tCTTGAATCATTAAGGTTAAATACATACGTAAGTTTGGCTCGTAATGCTTACTGGGACGCCAGCAGTAATGTGAATTATTCATTATCACTTAg                                              >  1:6325187/1‑93 (MQ=255)
                                            tCTTGAATCATTAAGGTTAAATACATACGTAAGTTTGGCTCGTAATACTTACTGGGACGCCAGCAGTAATGTGAATTATTCATTATCACTTAg                                              <  1:1241146/93‑1 (MQ=255)
                                              ttGAATCATTAAGGTTAAATACATACGTAAGTTTGGCTCGTAATGCTTACTGGGACGCCAGCAGTAATGTGAATTATTCATTATCACTTAGcc                                            >  2:552424/1‑93 (MQ=255)
                                              ttGAATCATTAAGGTTAAATACATACGTAAGTTTGGCTCGTAATGCTTACTGGGACGCCAGCAGTAATGTGAATTATTCATTATCACTTAGcc                                            <  2:1240597/93‑1 (MQ=255)
                                              ttGAATCATTAAGGTTAAATACATACGTAAGTTTGGCTCGTAATGCTTACTGGGACGCCAGCAGTAATGTGAATTATTCATTATCACTTAGcc                                            <  1:5680389/93‑1 (MQ=255)
                                               tGAATCATTAAGGTTAAATACATACGTAAGTTTGGCTCGTAATGCTTACTGGGACGCCAGCAGTAATGTGAATTATTCATTATCACTTAGCcg                                           >  1:57116/1‑93 (MQ=255)
                                               tGAATCATTAAGGTTAAATACATACGTAAGTTTGGCTCGTAATACTTACTGGGACGCCAGCAGTAATGTGAATTATTCATTATCACTTAGCcg                                           >  2:6198318/1‑93 (MQ=255)
                                                gAATCATTAAGGTTAAATACATACGTAAGTTTGGCTCGTAATGCTTACTGGGACGCCAGCAGTAATGTGAATTATTCATTATCACTTAGCcgc                                          >  2:5479930/1‑93 (MQ=255)
                                                 aaTCATTAAGGTTAAATACATACGTAAGTTTGGCTCGTAATACTTACTGGGACGCCAGCAGTAATGTGAATTATTCATTATCACTTAGCcgcg                                         >  2:2986563/1‑93 (MQ=255)
                                                  aTCATTAAGGTTAAATACATACGTAAGTTTGGCTCGTAATGCTTACTGGGACGCCAGCAGTAATGTGAATTATTCATTATCACTTAGCCGCGa                                        <  2:253701/93‑1 (MQ=255)
                                                     aTTAAGGTTAAATACATACGTAAGTTTGGCTCGTAATGCTTACTGGGACGCCAGCAGTAATGTGAATTATTCATTATCACTTAGCCGCGAttt                                     >  1:223590/1‑93 (MQ=255)
                                                     aTTAAGGTTAAATACATACGTAAGTTTGGCTCGTAATGCTTACTGGGACGCCAGCAGTAATGTGAATTATTCATTATCACTTAGCCGCGAttt                                     >  2:356313/1‑93 (MQ=255)
                                                     aTTAAGGTTAAATACATACGTAAGTTTGGCTCGTAATGCTTACTGGGACGCCAGCAGTAATGTGAATTATTCATTATCACTTAGCCGCGAttt                                     >  2:3801371/1‑93 (MQ=255)
                                                      ttAAGGTTAAATACATACGTAAGTTTGGCTCGTAATGCTTACTGGGACGCCAGCAGTAATGTGAATTATTCATTATCACTTAGCCGCGAtttt                                    <  1:6004283/93‑1 (MQ=255)
                                                      ttAAGGTTAAATACATACGTAAGTTTGGCTCGTAATGCTTACTGGGACGCCAGCAGTAATGTGAATTATTCATTATCACTTAGCCGCGAtttt                                    >  2:5338192/1‑93 (MQ=255)
                                                       tAAGGTTAAATACATACGTAAGTTTGGCTCGTAATACTTACTGGGACGCCAGCAGTAATGTGAATTATTCATTATCACTTAGCCGCGATTTTg                                   <  2:7172991/93‑1 (MQ=255)
                                                        aaGGTTAAATACATACGTAAGTTTGGCTCGTAATGCTTACTGGGACGCCAGCAGTAATGTGAATTATTCATTATCACTTAGCCGCGATTTTGa                                  <  2:5147214/93‑1 (MQ=255)
                                                          ggTTAAATACATACGTAAGTTTGGCTCGTAATGCTTACTGGGACGCCAGCAGTAATGTGAATTATTCATTATCACTTAGCCGCGATTTTGata                                <  1:4848633/93‑1 (MQ=255)
                                                          ggTTAAATACATACGTAAGTTTGGCTCGTAATGCTTACTGGGACGCCAGCAGTAATGTGAATTATTCATTATCACTTAGCCGCGATTTTGata                                <  1:1839204/93‑1 (MQ=255)
                                                          ggTTAAATACATACGTAAGTTTGGCTCGTAATACTTACTGGGACGCCAGCAGTAATGTGAATTATTCATTATCACTTAGCCGCGATTTTGata                                >  2:5826455/1‑93 (MQ=255)
                                                           gTTAAATACATACGTAAGTTTGGCTCGTAATGCTTACTGGGACGCCAGCAGTAATGTGAATTATTCATTATCACTTAGCCGCGATTTTGatat                               <  1:4258034/93‑1 (MQ=255)
                                                             tAAATACATACGTAAGTTTGGCTCGTAATGCTTACTGGGACGCCAGCAGTAATGTGAATTATTCATTATCACTTAGCCGCGATTTTGATATCg                             <  1:1086046/93‑1 (MQ=255)
                                                               aaTACATACGTAAGTTTGGCTCGTAATGCTTACTGGGACGCCAGCAGTAATGTGAATTATTCATTATCACTTAGCCGCGATTTTGATATCggg                           >  2:230034/1‑93 (MQ=255)
                                                               aaTACATACGTAAGTTTGGCTCGTAATGCTTACTGGGACGCCAGCAGTAATGTGAATTATTCATTATCACTTAGCCGCGATTTTGATATCggg                           >  2:4248272/1‑93 (MQ=255)
                                                               aaTACATACGTAAGTTTGGCTCGTAATACTTACTGGGACGCCAGCAGTAATGTGAATTATTCATTATCACTTAGCCGCGATTTTGATATCggg                           >  2:3446550/1‑93 (MQ=255)
                                                                 tacatacGTAAGTTTGGCTCGTAATGCTTACTGGGACGCCAGCAGTAATGTGAATTATTCATTATCACTTAGCCGCGATTTTGATATCGGGcc                         <  1:3356202/93‑1 (MQ=255)
                                                                 tacatacGTAAGTTTGGCTCGTAATGCTTACTGGGACGCCAGCAGTAATGTGAATTATTCATTATCACTTAGCCGCGATTTTGATATCGGGcc                         <  1:991195/93‑1 (MQ=255)
                                                                  acatacGTAAGTTTGGCTCGTAATGCTTACTGGGACGCCAGCAGTAATGTGAATTATTCATTATCACTTAGCCGCGATTTTGATATCGGGCCa                        <  2:4543180/93‑1 (MQ=255)
                                                                     tacGTAAGTTTGGCTCGTAATGCTTACTGGGACGCCAGCAGTAATGTGAATTATTCATTATCACTTAGCCGCGATTTTGATATCGGGCCATTa                     >  2:6930325/1‑93 (MQ=255)
                                                                         tAAGTTTGGCTCGTAATGCTTACTGGGACGCCAGCAGTAATGTGAATTATTCATTATCACTTAGCCGCGATTTTGATATCGGGCCATTaaaaa                 <  1:6135624/93‑1 (MQ=255)
                                                                           aGTTTGGCTCGTAATGCTTACTGGGACGCCAGCAGTAATGTGAATTATTCATTATCACTTAGCCGCGATTTTGATATCGGGCCATTAAAAAAc               <  2:1295542/93‑1 (MQ=255)
                                                                           aGTTTGGCTCGTAATGCTTACTGGGACGCCAGCAGTAATGTGAATTATTCATTATCACTTAGCCGCGATTTTGATATCGGGCCATTAAAAAAc               <  1:2767213/93‑1 (MQ=255)
                                                                           aGTTTGGCTCGTAATGCTTACTGGGACGCCAGCAGTAATGTGAATTATTCATTATCACTTAGCCGCGATTTTGATATCGGGCCATTAAAAAAc               >  1:1914017/1‑93 (MQ=255)
                                                                           aGTTTGGCTCGTAATGCTTACTGGGACGCCAGCAGTAATGTGAATTATTCATTATCACTTAGCCGCGATTTTGATATCGGGCCATTAAAAAAc               >  2:3658215/1‑93 (MQ=255)
                                                                           aGTTTGGCTCGTAATGCTTACTGGGACGCCAGCAGTAATGTGAATTATTCATTATCACTTAGCCGCGATTTTGAAATCGGGCCATTAAAAAAc               <  1:4914085/93‑1 (MQ=255)
                                                                            gTTTGGCTCGTAATGCTTACTGGGACGCCAGCAGTAATGTGAATTATTCATTATCACTTAGCCGCGATTTTGATATCGGGCCATTAAAAAACg              >  1:3424237/1‑93 (MQ=255)
                                                                            gTTTGGCTCGTAATGCTTACTGGGACGCCAGCAGTAATGTGAATTATTCATTATCACTTAGCCGCGATTTTGATATCGGGCCATTAAAAAACg              >  2:6669754/1‑93 (MQ=255)
                                                                              ttGGCTCGTAATGCTTACTGGGACGCCAGCAGTAATGTGAATTATTCATTATCACTTAGCCGCGATTTTGATATCGGGCCATTAAAAAACGtc            <  1:2347187/93‑1 (MQ=255)
                                                                               tGGCTCGTAATGCTTACTGGGACGCCAGCAGTAATGTGAATTATTCATTATCACTTAGCCGCGATTTTGATATCGGGCCATTAAAAAACGtct           <  2:4689042/93‑1 (MQ=255)
                                                                                ggCTCGTAATGCTTACTGGGACGCCAGCAGTAATGTGAATTATTCATTATCACTTAGCCGCGATTTTGATATCGGGCCATTAAAAAACGtctc          <  2:4022222/93‑1 (MQ=255)
                                                                                ggCTCGTAATGCTTACTGGGACGCCAGCAGTAATGTGAATTATTCATTATCACTTAGCCGCGATTTTGATATCGGGCCATTAAAAAACGtctc          <  1:722702/93‑1 (MQ=255)
                                                                                ggCTCGTAATACTTACTGGGACGCCAGCAGTAATGTGAATTATTCATTATCACTTAGCCGCGATTTTGATATCGGGCCATTAAAAAACGtctc          >  1:138818/1‑93 (MQ=255)
                                                                                 gCTCGTAATGCTTACTGGGACGCCAGCAGTAATGTGAATTATTCATTATCACTTAGCCGCGATTTTGATATCGGGCCATTAAAAAACGTCTcc         >  1:4484617/1‑93 (MQ=255)
                                                                                 gCTCGTAATGCTTACTGGGACGCCAGCAGTAATGTGAATTATTCATTATCACTTAGCCGCGATTTTGATATCGGGCCATTAAAAAACGTCTcc         >  2:4508798/1‑93 (MQ=255)
                                                                                  cTCGTAATGCTTACTGGGACGCCAGCAGTAATGTGAATTATTCATTATCACTTAGCCGCGATTTTGATATCGGGCCATTAAAAAACGTCTCCa        <  1:3269510/93‑1 (MQ=255)
                                                                                    cGTAATGCTTACTGGGACGCCAGCAGTAATGTGAATTATTCATTATCACTTAGCCGCGATTTTGATATCGGGCCATTAAAAAACGTCTCCACt      <  1:1988807/93‑1 (MQ=255)
                                                                                          |                                                                                          
TATCAGAATGAAAAAGAGAGTTATATTGTCACGTTTAACCAGTATCTTGAATCATTAAGGTTAAATACATACGTAAGTTTGGCTCGTAATACTTACTGGGACGCCAGCAGTAATGTGAATTATTCATTATCACTTAGCCGCGATTTTGATATCGGGCCATTAAAAAACGTCTCCACTTCTC  >  NC_000913/3186835‑3187015

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