Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,360,454 (A)8→7 coding (13/2382 nt) yhcD → putative outer membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,360,4470A.99.3% 1237.0 / 0.3 269coding (6/2382 nt)yhcDputative outer membrane protein
Reads supporting (aligned to +/- strand):  ref base A (2/0);  new base . (150/117);  total (152/117)
Fisher's exact test for biased strand distribution p-value = 5.07e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 299 total aligned reads displayed.

AGTCGAGTTTTACCCAGCCAGCCGATATGGTTATCGCGGTAATAAATATGTCACTGACTTAAAATAACTTTGCCTGGAGCGACAAGGATGTTAAAAAAAACGTTACTGGCCTACACCATTGGTTTTGCCTTTTCTCCCCCAGCCAATGCAGATGGTATAGAGATTGCCGCTGTTGATTTTG  >  NC_000913/3360355‑3360535
                                                                                            |                                                                                        
aGTCGAGTTTTACCCAGCCAGCCGATATGGTTATCGCGGTAATAAATATGTCACTGACTTAAAATAACTTTGCCTGGAGCGACAAGGATGTTa                                                                                          >  1:3987002/1‑93 (MQ=255)
 gTCGAGTTTTACCCAGCCAGCCGATATGGTTATCGCGGTAATAAATATGTCACTGACTTAAAATAACTTTGCCTGGAGCGACAAGGATGTTaa                                                                                         >  1:6274027/1‑93 (MQ=255)
 gTCGAGTTTTACCCAGCCAGCCGATATGGTTATCGCGGTAATAAATAAGTCACTGACTTAAAATAACTTTGCCTGGAGCGACAAGGATGTTaa                                                                                         >  2:4516399/1‑93 (MQ=255)
    gAGTTTTACCCAGCCAGCCGATATGGTTATCGCGGTAATAAATATGTCACTGACTTAAAATAACTTTGCCTGGAGCGACAAGGATGTTaaaaa                                                                                      >  2:357507/1‑93 (MQ=255)
      gTTTTACCCAGCCAGCCGATATGGTTATCGCGGTAATAAATATGTCACTGACTTAAAATAACTTTGCCTGGAGCGACAAGGATGTTaaaaaaa                                                                                    <  1:4118382/93‑1 (MQ=255)
      gTTTTACCCAGCCAGCCGATATGGTTATCGCGGTAATAAATATGTCACTGACTTAAAATAACTTTGCCTGGAGCGACAAGGATGTTaaaaaaa                                                                                    >  1:324645/1‑93 (MQ=255)
       ttttACCCAGCCAGCCGATATGGTTATCGCGGTAATAAATATGTCACTGACTTAAAATAACTTTGCCTGGAGCGACAAGGATGTTaaaaaaac                                                                                   <  1:5073174/93‑2 (MQ=255)
       ttttACCCAGCCAGCCGATATGGTTATCGCGGTAATAAATATGTCACTGACTTAAAATAACTTTGCCTGGAGCGACAAGGATGTTaaaaaaac                                                                                   <  2:1188553/93‑2 (MQ=255)
         ttACCCAGCCAGCCGATATGGTTATCGCGGTAATAAATATGTCACTGACTTAAAATAACTTTGCCTGGAGCGACAAGGATGTTaaaaaaacgt                                                                                 <  2:4634527/93‑4 (MQ=255)
           aCCCAGCCAGCCGATATGGTTATCGCGGTAATAAATATGTCACTGACTTAAAATAACTTTGCCTGGAGCGACAAGGATGTTAAAAAAACGTTa                                                                              >  2:3261792/1‑93 (MQ=255)
             ccagccagCCGATATGGTTATCGCGGTAATAAATATGTCACTGACTTAAAATAACTTTGCCTGGAGCGACAAGGATGTTAAAAAAACGTTACt                                                                            >  2:3062604/1‑93 (MQ=255)
             ccagccagCCGATATGGTTATCGCGGTAATAAATATGTCACTGACTTAAAATAACTTTGCCTGGAGCGACAAGGATGTTAAAAAAACGTTACt                                                                            <  1:5129500/93‑1 (MQ=255)
                 ccagccGATATGGTTATCGCGGTAATAAATATGTCACTGACTTAAAATACCTTTGCCTGGAGCGACAAGGATGTTAAAAAAACGTTACTGGcc                                                                        <  1:2610187/93‑1 (MQ=255)
                 ccagccGATATGGTTATCGCGGTAATAAATATGTCACTGACTTAAAATAACTTTGCCTGGAGCGACAAGGATGTTAAAAAAACGTTACTGGcc                                                                        >  2:4851521/1‑93 (MQ=255)
                   agccGATATGGTTATCGCGGTAATAAATATGTCACTGACTTAAAATAACTTTGCCTGGAGCGACAAGGATGTTAAAAAAACGTTACTGGCCTa                                                                      >  2:3829504/1‑93 (MQ=255)
                    gccGATATGGTTATCGCGGTAATAAATATGTCACTGACTTAAAATAACTTTGCCTGGAGCGACAAGGATGTTAAAAAAACGTTACTGGCCTac                                                                     >  2:2473007/1‑93 (MQ=255)
                     ccGATATGGTTATCGCGGTAATAAATATGTCACTGACTTAAAATAACTTTGCCTGGAGCGACAAGGATGTTAAAAAAACGTTACTGGCCTaca                                                                    <  1:1395762/93‑1 (MQ=255)
                     ccGATATGGTTATCGCGGTAATAAATATGTCACTGACTTAAAATAACTTTGCCTGGAGCGACAAGGATGTTAAAAAAACGTTACTGGCCTaca                                                                    <  1:5772370/93‑1 (MQ=255)
                        agatGGTTATCGCGGTAATAAATATGTCACTGACTTAAAATAACTTTGCCTGGAGCGACAAGGATGTTAAAAAAACGTTACTGGCCTACACCa                                                                 >  1:4202988/3‑93 (MQ=255)
                         tatGGTTATCGCGGTAATAAATATGTCACTGACTTAAAATAACTTTGCCTGGAGCGACAAGGATGTTAAAAAAACGTTACTGGCCTACACCAt                                                                >  2:3787968/1‑93 (MQ=255)
                          atGGTTATCGCGGTAATAAATATGTCACTGACTTAAAATAACTTTGCCTGGAGCGACAAGGATGTTAAAAAAACGTTACTGGCCTACACCAtt                                                               <  1:2922652/93‑1 (MQ=255)
                          atGGTTATCGCGGTAATAAATATGTCACTGACTTAAAATAACTTTGCCTGGAGCGACAAGGATGTTAAAAAAACGTTACTGGCCTACACCAtt                                                               <  2:1951955/93‑1 (MQ=255)
                          atGGTTATCGCGGTAATAAATATGTCACTGACTTAAAATAACTTTGCCTGGAGCGACAAGGATGTTAAAAAAACGTTACTGGCCTACACCAtt                                                               <  2:2171048/93‑1 (MQ=255)
                           tGGTTATCGCGGTAATAAATATGTCACTGACTTAAAATAACTTTGCCTGGAGCGACAAGGATGTTAAAAAAACGTTACTGGCCTACACCATTg                                                              <  2:4259092/93‑1 (MQ=255)
                            ggTTATCGCGGTAATAAATATGTCACTGACTTAAAATAACTTTGCCTGGAGCGACAAGGATGTTAAAAAAACGTTACTGGCCTACACCATTgg                                                             >  2:651051/1‑93 (MQ=255)
                            ggTTATCGCGGTAATAAATATGTCACTGACTTAAAATAACTTTGCCTGGAGCGACAAGGATGTTAAAAAAACGTTACTGGCCTACACCATTgg                                                             <  2:1838234/93‑1 (MQ=255)
                             gTTATCGCGGTAATAAATATGTCACTGACTTAAAATAACTTTGCCTGGAGCGACAAGGATGTTAAAAAAACGTTACTGGCCTACACCATTGGt                                                            <  1:229002/93‑1 (MQ=255)
                              ttATCGCGGTAATAAATATGTCACTGACTTAAAATAACTTTGCCTGGAGCGACAAGGATGTTAAAAAAACGTTACTGGCCTACACCATTGGtt                                                           >  2:5428452/1‑93 (MQ=255)
                               tATCGCGGTAATAAATATGTCACTGACTTAAAATAACTTTGCCTGGAGCGACAAGGATGTTAAAAAAACGTTACTGGCCTACACCATTGGttt                                                          <  1:4981858/93‑1 (MQ=255)
                               tATCGCGGTAATAAATATGTCACTGACTTAAAATAACTTTGCCTGGAGCGACAAGGATGTTAAAAAAACGTTACTGGCCTACACCATTGGttt                                                          <  1:2632192/93‑1 (MQ=255)
                                aTCGCGGTAATAAATATGTCACTGACTTAAAATAACTTTGCCTGGAGCGACAAGGATGTTAAAAAAACGTTACTGGCCTACACCATTGGtttt                                                         >  2:5146885/1‑93 (MQ=255)
                                aTCGCGGTAATAAATATGTCACTGACTTAAAATAACTTTGCCTGGAGCGACAAGGATGTTAAAAAAACGTTACTGGCCTACACCATTGGtttt                                                         <  1:3728414/93‑1 (MQ=255)
                                   gcgGTAATAAATATGTCACTGACTTAAAATAACTTTGCCTGGAGCGACAAGGATGTTAAAAAAACGTTACTGGCCTACACCATTGGTTTTGcc                                                      <  1:3002445/93‑1 (MQ=255)
                                   gcgGTAATAAATATGTCACTGACTTAAAATAACTTTGCCTGGAGCGACAAGGATGTTAAAAAAACGTTACTGGCCTACACCATTGGTTTTGcc                                                      <  2:6140243/93‑1 (MQ=255)
                                     ggTAATAAATATGTCACTGACTTAAAATAACTTTGCCTGGAGCGACAAGGATGTTAAAAAAACGTTACTGGCCTACACCATTGGTTTTGCCtt                                                    <  1:1026863/93‑1 (MQ=255)
                                     ggTAATAAATATGTCACTGACTTAAAATAACTTTGCCTGGAGCGACAAGGATGTTAAAAAAACGTTACTGGCCTACACCATTGGTTTTGCCtt                                                    >  1:4378392/1‑93 (MQ=255)
                                      gTAATAAATATGTCACTGACTTAAAATAACTTTGCCTGGAGCGACAAGGATGTTAAAAAAACGTTACTGGCCTACACCATTGGTTTTGCCttt                                                   <  1:471659/93‑1 (MQ=255)
                                        aataaataTGTCACTGACTTAAAATAACTTTGCCTGGAGCGACAAGGATGTTAAAAAAACGTTACTGGCCTACACCATTGGTTTTGCCTTTtc                                                 >  2:5303303/1‑93 (MQ=255)
                                        aataaataTGTCACTGACTTAAAATAACTTTGCCTGGAGCGACAAGGATGTTAAAAAAACGTTACTGGCCTACACCATTGGTTTTGCCTTTtc                                                 <  1:2537337/93‑1 (MQ=255)
                                        aataaataTGTCACTGACTTAAAATAACTTTGCCTGGAGCGACAAGGATGTTAAAAAAACGTTACTAGCCTACACCATTCGTTTTGCCTTTtc                                                 >  1:5877004/1‑93 (MQ=255)
                                         ataaataGGTCGCTGACTTAAAATAACTTTGCCTGGAGCGACAACGATGTTAAAAAAACGTTACTGGCCTACACCATTGGTTTTGCCTTTtct                                                <  2:3031959/93‑1 (MQ=255)
                                          taaataTGTCACTGACTTAAAATAACTTTGCCTGGAGCGACAAGGATGTTAAAAAAACGTTACTGGCCTACACCATTGGTTTTGCCTTTtctc                                               >  1:5516320/1‑93 (MQ=255)
                                          taaataTGTCACTGACTTAAAATAACTTTGCCTGGAGCGACAAGGATGTTAAAAAAACGTTACTGGCCTACACCATTGGTTTTGCCTTTtctc                                               >  1:1646938/1‑93 (MQ=255)
                                          taaataTGTCACTGACTTAAAATAACTTTGCCTGGAGCGACAAGGATGTTAAAAAAACGTTACTGGCCTACACCATTGGTTTTGCCTTTtctc                                               >  2:1517223/1‑93 (MQ=255)
                                           aaataTGTCACTGACTTAAAATAACTTTGCCTGGAGCGACAAGGATGTTAAAAAAACGTTACTGGCCTACACCATTGGTTTTGCCTTTTCTcc                                              >  2:3925856/1‑93 (MQ=255)
                                            aataTGTCACTGACTTAAAATAACTTTGCCTGGAGCGACAAGGATGTTAAAAAAACGTTACTGGCCTACACCATTGGTTTTGCCTTTTCTccc                                             >  1:4196116/1‑93 (MQ=255)
                                            aataTGTCACTGACTTAAAATAACTTTGCCTGGAGCGACAAGGATGTTAAAAAAACGTTACTGGCCTACACCATTGGTTTTGCCTTTTCTccc                                             >  1:3148370/1‑93 (MQ=255)
                                            aataTGTCACTGACTTAAAATAACTTTGCCTGGAGCGACAAGGATGTTAAAAAAACGTTACTGGCCTACACCATTGGTTTTGCCTTTTCTccc                                             >  1:1948929/1‑93 (MQ=255)
                                             atatGTCACTGACTTAAAATAACTTTGCCTGGAGCGACAAGGATGTTAAAAAAACGTTACTGGCCTACACCATTGGTTTTGCCTTTTCTcccc                                            <  1:6352564/93‑1 (MQ=255)
                                              tatGTCACTGACTTAAAATAACTTTGCCTGGAGCGACAAGGATGTTAAAAAAACGTTACTGGCCTACACCATTGGTTTTGCCTTTTCTccccc                                           >  1:2103297/1‑93 (MQ=255)
                                               atGTCACTGACTTAAAATAACTTTGCCTGGAGCGACAAGGATGTTAAAAAAACGTTACTGGCCTACACCATTGGTTTTGCCTTTTCTCCCCCa                                          <  1:5549571/93‑1 (MQ=255)
                                                  tCACTGACTTAAAATAACTTTGCCTGGAGCGACAAGGATGTTAAAAAAACGTTACTGGCCTACACCATTGGTTTTGCCTTTTCTCCCCCAGcc                                       >  1:6121517/1‑93 (MQ=255)
                                                    aCTGACTTAAAATAACTTTGCCTGGAGCGACAAGGATGTTAAAAAAACGTTACTGGCCTACACCATTGGTTTTGCGTTATCTCCCCCAGCCaa                                     >  2:1147594/1‑93 (MQ=255)
                                                      cgACTTAAAATAACTTTGCCTGGAGCGACAAGGATGTTAAAAAAACGTTACTGGCCTACCCCATTGGTTTTGCCTTTTCTCCCCCAGCCAATg                                   <  1:297891/92‑1 (MQ=255)
                                                       gACTTAAAATAACTTTGCCTGGAGCGACAAGGATGTTAAAAAAACGTTACTGGCCTACACCATTGGTTTTGCCTTTTCTCCCCCAGCCAATGc                                  >  2:113203/1‑93 (MQ=255)
                                                       gACTTAAAATAACTTTGCCTGGAGCGACAAGGATGTTAAAAAAACGTTACTGGCCTACACCATTGGTTTTGCCTTTTCTCCCCCAGCCAATGc                                  >  1:2871045/1‑93 (MQ=255)
                                                       gACTTAAAATAACTTTGCCTGGAGCGACAAGGATGTTAAAAAAACGTTACTGGCCTACACCATTGGTTTTGCCTTTTCTCCCCCAGCCAATGc                                  >  1:6616879/1‑93 (MQ=255)
                                                         cTTAAAATAACTTTGCCTGGAGCGACAAGGATGTTAAAAAAACGTTACTGGCCTACACCATTGGCTTTGCCTTTTCTCCCCCAGCCAATGCAg                                >  1:2733655/1‑93 (MQ=255)
                                                          ttAAAATAACTTTGCCTGGAGCGACAAGGATGTTAAAAAAACGTTACTGGCCTACACCATTGGTTTTGCCTTTTCTCCCCCAGCCAATGCAGa                               >  1:340656/1‑93 (MQ=255)
                                                           tAAAATAACTTTGCCTGGAGCGACAAGGATGTTAAAAAAACGTTACTGGCCTACACCATTGGTTTTGCCTTTTCTCCCCCAGCCAATGCAGAt                              <  1:1968528/93‑1 (MQ=255)
                                                           tAAAATAACTTTGCCTGGAGCGACAAGGATGTTAAAAAAACGTTACTGGCCTACACCATTGGTTTTGCCTTTTCTCCCCCAGCCAATGCAGAt                              <  2:4206278/93‑1 (MQ=255)
                                                           tAAAATAACTTTGCCTGGAGCGACAAGGATGTTAAAAAAACGTTACTGGCCTACACCATTGGTTTTGCCTTTTCTCCCCCAGCCAATGCAGAt                              <  1:950896/93‑1 (MQ=255)
                                                             aaaTAACTTTGCCTGGAGCGACAAGGATGTTAAAAAAACGTTACTGGCCTACACCATTGGTTTTGCCTTTTCTCCCCCAGCCAATGCAGATgg                            <  2:4398297/93‑1 (MQ=255)
                                                             aaaTAACTTTGCCTGGAGCGACAAGGATGTTAAAAAAACGTTACTGGCCTACACCATTGGTTTTGCCTTTTCTCCCCCAGCCAATGCAGATgg                            >  2:5849780/1‑93 (MQ=255)
                                                             aaaTAACTTTGCCTGGAGCGACAAGGATGTTAAAAAAACGTTACTGGCCTACACCATTGGTTTTGCCTTTTCTCCCCCAGCCAATGCAGATgg                            >  2:3318211/1‑93 (MQ=255)
                                                             aaaTAACTTTGCCTGGAGCGACAAGGATGTTAAAAAAACGTTACTGGCCTACACCATTGGTTTTGCCTTTTCTCCCCCAGCCAATGCAGATgg                            >  2:6233032/1‑93 (MQ=255)
                                                             aaaTAACTTTGCCTGGAGCGACAAGGATGTTAAAAAAACGTTACTGGCCTACACCATTGGTTTTGCCTTTTCTCCCCCAGCCAATGCAGATgg                            >  2:2797266/1‑93 (MQ=255)
                                                                tAACTTTGCCTGGAGCGACAAGGATGTTAAAAAAACGTTACTGGCCTACACCATTGGTTTTGCCTTTTCTCCCCCAGCCAATGCAGATGGtat                         >  2:128309/1‑93 (MQ=255)
                                                                tAACTTTGCCTGGAGCGACAAGGATGTTAAAAAAACGTTACTGGCCTACACCATTGGTTTTGCCTTTTCTCCCCCAGCCAATGCAGATGGtat                         >  1:6487768/1‑93 (MQ=255)
                                                                 aaCTTTGCCTGGAGCGACAAGGATGTTAAAAAAACGTTACTGGCCTACACCATTGGTTTTGCCTTTTCTCCCCCAGCCAATGCAGATGGtata                        <  2:2240501/93‑1 (MQ=255)
                                                                  aCTTTGCCTGGAGCGACAAGGATGTTAAAAAAACGTTACTGGCCTACACCATTGGTTTTGCCTTTTCTCCCCCAGCCAATGCAGATGGTATag                       >  1:6205407/1‑93 (MQ=255)
                                                                   cTTTGCCTGGAGCGACAAGGATGTTAAAAAAACGTTACTGGCCTACACCATTGGTTTTGCCTTTTCTCCCCCAGCCAATGCAGATGGTATaga                      >  1:1263607/1‑93 (MQ=255)
                                                                   cTTTGCCTGGAGCGACAAGGATGTTAAAAAAACGTTACTGGCCTACACCATTGGTTTTGCCTTTTCTCCCCCAGCCAATGCAGATGGTATaga                      <  2:2709801/93‑1 (MQ=255)
                                                                    tttGCCTGGAGCGACAAGGATGTTAAAAAAACGTTACTGGCCTACACCATTGGTTTTGCCTTTTCTCCCCCAGCCAATGCAGATGGTATagag                     >  2:953430/1‑93 (MQ=255)
                                                                    tttGCCTGGAGCGACAAGGATGTTAAAAAAACGTTACTGGCCTACACCATTGGTTTTGCCTTTTCTCCCCCAGCCAATGCAGATGGTATagag                     <  1:1320419/93‑1 (MQ=255)
                                                                      tGCCTGGAGCGACAAGGATGTTAAAAAAACGTTACTGGCCTACACCATTGGTTTTGCCTTTTCTCCCCCAGCCAATGCAGATGGTATAGAGAt                   <  2:5537496/93‑1 (MQ=255)
                                                                       gCCTGGAGCGACAAGGATGTTAAAAAAACGTTACTGGCCTACACCATTGGTTTTGCCTTTTCTCCCCCAGCCAATGCAGATGGTATAGAGAtt                  <  1:1584239/93‑1 (MQ=255)
                                                                         cTGGAGCGACAAGGATGTTAAAAAAACGTTACTGGCCTACACCATTGGTTTTGCCTTTTCTCCCCCAGCCAATGCAGATGGTATAGAGATTGc                >  2:2883144/1‑93 (MQ=255)
                                                                         cTGGAGCGACAAGGATGTTAAAAAAACGTTACTGGCCTACACCATTGGTTTTGCCTTTTCTCCCCCAGCCAATGCAGATGGTATAGAGATTGc                >  1:752233/1‑93 (MQ=255)
                                                                         cTGGAGCGACAAGGATGTTAAAAAAACGTTACTGGCCTACACCATTGGTTTTGCCTTTTCTCCCCCAGCCAATGCAGATGGTATAGAGATTGc                >  2:546491/1‑93 (MQ=255)
                                                                            gAGCGACAAGGATGTTAAAAAAACGTTACTGGCCTACACCATTGGTTTTGCCTTTTCTCCCCCAGCCAATGCAGATGGTATAGAGATTGCCGc             <  1:369367/93‑1 (MQ=255)
                                                                            gAGCGACAAGGATGTTAAAAAAACGTTACTGGCCTACACCATTGGTTTTGCCTTTTCTCCCCCAGCCAATGCAGATGGTATAGAGATTGCCGc             >  2:6064929/1‑93 (MQ=255)
                                                                              gCGACAAGGATGTTAAAAAAACGTTACTGGCCTACACCATTGGTTTTGCCTTTTCTCCCCCAGCCAATGCAGATGGTATAGAGATTGCCGCTg           >  2:1701044/1‑93 (MQ=255)
                                                                               cGACAAGGATGTTAAAAAAACTTTACTGGCCTACACCATTGGTTTTGCCTTTTCTCCCCCAGCCAATGCAGATGGTATAGAGATTGCCGCTGt          <  1:3858937/93‑1 (MQ=255)
                                                                                gACAAGGATGTTAAAAAAACGTTACTGGCCTACACCATTGGTTTTGCCTTTTCTCCCCCAGCCAATGGAGATGGTATAGAGATTGCCGCTGtt         <  1:4321706/93‑1 (MQ=255)
                                                                                 aCAAGGATGTTAAAAAAACGTTACTGGCCTACACCATTGGTTTTGCCTTTTCTCCCCCAGCCAATGCAGATGGTATAGAGATTGCCGCTGTTg        >  2:1038456/1‑93 (MQ=255)
                                                                                 aCAAGGATGTTAAAAAAACGTTACTGGCCTACACCATTGGTTTTGCCTTTTCTCCCCCAGCCAATGCAGATGGTATAGAGATTGCCGCTGTTg        >  1:249957/1‑93 (MQ=255)
                                                                                   aaGGATGTTAAAAAAACGTTACTGGCCTACACCATTGGTTTTGCCTTTTCTCCCCCAGCCAATGCAGATGGTATAGAGATTGCCGCTGTTGAt      >  1:2690473/1‑93 (MQ=255)
                                                                                   aaGGATGTTAAAAAAACGTTACTGGCCTACACCATTGGTTTTGCCTTTTCTCCCCCAGCCAATGCAGATGGTATAGAGATTGCCGCTGTTGAt      >  2:838908/1‑93 (MQ=255)
                                                                                   aaGGATGTTAAAAAAACGTTACTGGCCTACACCATTGGTTTTGCCTTTTCTCCCCCAGCCAATGCAGATGGTATAGAGATTGCCGCTGTTGAt      >  2:6686906/1‑93 (MQ=255)
                                                                                   aaGGATGTTAAAAAAACGTTACTGGCCTACACCATTGGTTTTGCCTTTTCTCCCCCAGCCAATGCAGATGGTATAGAGATTGCCGCTGTTGAt      >  1:5418691/1‑93 (MQ=255)
                                                                                    atgATGTTAAAAAAACGTTACTGGCCTACACCATTGGTTTTGCCTTTTCTCCCCCAGCCAATGCAGATGGTATAGAGATTGCCGCTGTTGAtt     <  1:4862930/91‑1 (MQ=255)
                                                                                    aGGATGTTAAAAAAACGTTACTGGCCTACACCATTGGTTTTGCCTTTTCTCCCCCAGCCAATGCAGATGGTATAGAGATTGCCGCTGTTGAtt     <  1:608345/93‑1 (MQ=255)
                                                                                    aGGATGTTAAAAAAACGTTACTGGCCTACACCATTGGTTTTGCCTTTTCTCCCCCAGCCAATGCAGATGGTATAGAGATTGCCGCTGTTGAtt     >  2:5979779/1‑93 (MQ=255)
                                                                                    aGGATGTTAAAAAAACGTTACTGGCCTACACCATTGGTTTTGCCTTTTCTCCCCCAGCCAATGCAGATGGTATAGAGATTGCCGCTGTTGAtt     <  1:3069125/93‑1 (MQ=255)
                                                                                    aGGATGTTAAAAAAACGTAACTGGCCTACACCATTGGTTTTGCCTTTTCTCCCCCAGCCAATGCAGATGGTATAGAGATTGCCGCTGTTGAtt     >  1:805685/1‑93 (MQ=255)
                                                                                     ggATGTTAAAAAAACGTTACTGGCCTACACCATTGGTTTTGCCTTTTCTCCCCCAGCCAATGCAGATGGTATAGAGATTGCCGCTGTTGAttt    <  2:4913695/93‑1 (MQ=255)
                                                                                      gATGTTAAAAAAACGTTACTGGCCTACACCATTGGTTTTGCCTTTTCTCCCCCAGCCAATGCAGATGGTATAGAGATTGCCGCTGTTGAtttt   >  1:523865/1‑93 (MQ=255)
                                                                                      gATGTTAAAAAAACGTTACTGGCCTACACCATTGGTTTTGCCTTTTCTCCCCCAGCCAATGCAGATGGTATAGAGATTGCCGCTGTTGAtttt   <  1:6630018/93‑1 (MQ=255)
                                                                                      gATGTTAAAAAAACGTTACTGGCCTACACCATTGGTTTTGCCTTTTCTCCCCCAGCCAATGCAGATGGTATAGAGATTGCCGCTGTTGAtttt   <  2:2406080/93‑1 (MQ=255)
                                                                                            |                                                                                        
AGTCGAGTTTTACCCAGCCAGCCGATATGGTTATCGCGGTAATAAATATGTCACTGACTTAAAATAACTTTGCCTGGAGCGACAAGGATGTTAAAAAAAACGTTACTGGCCTACACCATTGGTTTTGCCTTTTCTCCCCCAGCCAATGCAGATGGTATAGAGATTGCCGCTGTTGATTTTG  >  NC_000913/3360355‑3360535

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