Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 477,992 A→G intergenic (‑145/+13) ylaB ← / ← ylaC hypothetical protein/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913477,9920AG99.4% 1564.0 / ‑4.0 477intergenic (‑145/+13)ylaB/ylaChypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (2/1);  new base G (244/230);  total (246/231)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 492 total aligned reads displayed.

TGGAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCTTTACGGACGATCATTTTTT  >  NC_000913/477902‑478081
                                                                                          |                                                                                         
     aTAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTGTGTAAcc                                                                                    >  2:1744864/1‑93 (MQ=255)
     aTAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTGTGTAAcc                                                                                    <  1:3490099/93‑1 (MQ=255)
      ttgggaggTAGAGCTGGTTTTTCGTGGAGTGATGGGTGGAAACGGTCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTGTGTAAccc                                                                                   <  1:7717740/90‑1 (MQ=255)
      tAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTGTGTAAccc                                                                                   >  1:2797229/1‑93 (MQ=255)
        aggaggTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTGTGTAACCCTa                                                                                 >  2:9627325/1‑93 (MQ=255)
        aggaggTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTa                                                                                 >  2:1933255/1‑93 (MQ=255)
         ggaggTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTGTGTAACCCTaa                                                                                <  1:7928534/93‑1 (MQ=255)
         ggaggTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTGTGTAACCCTaa                                                                                <  2:3903025/93‑1 (MQ=255)
          gaggTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTGTGTAACCCTAAc                                                                               <  2:6936667/93‑1 (MQ=255)
          gaggTAGAGCTGGTATTTCATGAAGTGATGGGGGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTGGGTAACCCTAAc                                                                               >  2:3542509/1‑93 (MQ=255)
           aggTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTGTGTAACCCTAACt                                                                              >  2:9027727/1‑93 (MQ=255)
           aggTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTGTGTAACCCTAACt                                                                              >  2:2587031/1‑93 (MQ=255)
            ggTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTGTGTAACCCTAACtt                                                                             >  1:1809248/1‑93 (MQ=255)
             gTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTGTGTAACCCTAACTTa                                                                            <  1:980334/93‑1 (MQ=255)
               agagCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTGTGTAACCCTAACTTAAg                                                                          <  2:2086764/93‑1 (MQ=255)
                  gCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTGTGTAACCCTAACTTAAGATg                                                                       >  1:138588/1‑93 (MQ=255)
                   cTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTGTGTAACCCTAACTTAAGATGt                                                                      <  2:4169406/93‑1 (MQ=255)
                   cTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTGTGTAACCCTAACTTAAGATGt                                                                      >  1:3788954/1‑93 (MQ=255)
                    tGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTGTGTAACCCTAACTTAAGATGtt                                                                     >  1:7335313/1‑93 (MQ=255)
                    tGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTGTGTAACCCTAACTTAAGATGtt                                                                     <  1:5534814/93‑1 (MQ=255)
                    tGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTGTGTAACCCTAACTTAAGATGtt                                                                     >  1:6403885/1‑93 (MQ=255)
                    tGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTGTGTAACCCTAACTTAAGATGtt                                                                     >  1:3646383/1‑93 (MQ=255)
                    tGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTGTGTAACCCTAACTTAAGATGtt                                                                     >  2:489082/1‑93 (MQ=255)
                     gttATTTCATGAAGTAATGGGTGAAAACGATCATTGTCATTTGCCCGGGCGCAGACCTTTACGCCTGTTGTGTACCCCTAACTTAAGATGTTg                                                                    <  2:626828/91‑1 (MQ=255)
                     ggTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTGTGTAACCCTAACTTAAGATGTTg                                                                    <  2:2809318/93‑1 (MQ=255)
                      gTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTGTGTAACCCTAACTTAAGATGTTGa                                                                   >  1:4057300/1‑93 (MQ=255)
                      gTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTGTGTAACCCTAACTTAAGATGTTGa                                                                   >  1:7465142/1‑93 (MQ=255)
                        aTTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTGTGTAACCCTAACTTAAGATGTTGAtt                                                                 <  2:1785261/93‑1 (MQ=255)
                        aTTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTGTGTAACCCTAACTTAAGATGTTGAtt                                                                 <  2:5195059/93‑1 (MQ=255)
                          ttCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTGTGTAACCCTAACTTAAGATGTTGATTCg                                                               <  2:609521/93‑1 (MQ=255)
                            cATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTGTGTAACCCTAACTTAAGATGTTGATTCGgc                                                             >  2:869680/1‑93 (MQ=255)
                             aTGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTGTGTAACCCTAACTTAAGATGTTGATTCGgcg                                                            >  1:10042526/1‑93 (MQ=255)
                              tGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTGTGTAACCCTAACTTAAGATGTTGATTCGgcgc                                                           >  1:4544266/1‑93 (MQ=255)
                                aaGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTGTGTAACCCTAACTTAAGATGTTGATTCGGCgcga                                                         <  2:1639577/93‑1 (MQ=255)
                                 aGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTGTGTAACCCTAACTTAAGATGTTGATTCGGCgcgag                                                        <  1:9524530/93‑1 (MQ=255)
                                  gTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTGTGTAACCCTAACTTAAGATGTTGATTCGGCgcgagc                                                       <  1:10004586/93‑1 (MQ=255)
                                  gTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTGTGTAACCATAACTTAAGATGTTGATTCGGCgcgagc                                                       <  2:3087585/93‑1 (MQ=255)
                                   tGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTGTGTAACCCTAACTTAAGATGTTGATTCGGCgcgagcg                                                      <  1:513667/93‑1 (MQ=255)
                                    gATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTGTGTAACCCTAACTTAAGATGTTGATTCGGCgcgagcga                                                     <  2:4563393/93‑1 (MQ=255)
                                    gATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTGTGTAACCCTAACTTAAGATGTTGATTCGGCgcgagcga                                                     <  2:8488714/93‑1 (MQ=255)
                                       gggTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTGTGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAggg                                                  >  2:4908379/1‑93 (MQ=255)
                                       gggTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTGTGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAggg                                                  >  1:5951380/1‑93 (MQ=255)
                                       gggTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTGTGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAggg                                                  <  2:4069374/93‑1 (MQ=255)
                                         gTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTGTGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTa                                                >  2:3210418/1‑93 (MQ=255)
                                           ggAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTGTGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTaaa                                              <  2:5715347/93‑1 (MQ=255)
                                             aaaCGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTGTGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTaaaaa                                            >  1:2884131/1‑93 (MQ=255)
                                              aaCGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTGTGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAat                                           <  2:7515744/93‑1 (MQ=255)
                                              aaCGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTGTGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAat                                           >  2:2252937/1‑93 (MQ=255)
                                              aaCGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTGTGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAat                                           >  1:4834017/1‑93 (MQ=255)
                                              aaCGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTGTGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAat                                           >  2:8887713/1‑93 (MQ=255)
                                               aCGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTGTGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAata                                          <  1:427570/93‑1 (MQ=255)
                                                 gATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTGTGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATc                                        >  1:5447129/1‑93 (MQ=255)
                                                     aTTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTGTGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTaa                                    >  1:8229016/1‑93 (MQ=255)
                                                      ttGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTGTGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTaaa                                   <  1:631533/93‑1 (MQ=255)
                                                       tGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTGTGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTaaaa                                  <  2:6584286/93‑1 (MQ=255)
                                                       tGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTGTGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTaaaa                                  <  2:7024661/93‑1 (MQ=255)
                                                          cATTTGACCAGGCGCAGAGCTTTACGCCTGGTGTGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGa                               >  2:968044/1‑93 (MQ=255)
                                                          cATTTGACCAGGCGCAGAGCTTTACGCCTGGTGTGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGa                               <  2:113479/93‑1 (MQ=255)
                                                          cATTTGACCAGGCGCAGAGCTTTACGCCTGGTGTGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGa                               >  2:8200178/1‑93 (MQ=255)
                                                           aTTTGACCAGGCGCAGAGCTTTACGCCTGGTGTGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGAc                              >  1:2523777/1‑93 (MQ=255)
                                                           aTTTGACCAGGCGCAGAGCTTTACGCCTGGTGTGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGAc                              >  2:3630635/1‑93 (MQ=255)
                                                            tttGACCAGGCGCAGAGCTTTACGCCTGGTGTGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACa                             <  2:5402010/93‑1 (MQ=255)
                                                            tttGACCAGGCGCAGAGCTTTACGCCTGGCGTGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACa                             >  1:1517756/1‑93 (MQ=255)
                                                              tGACCAGGCGCAGAGCTTTACGCCTGGTGTGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGt                           <  1:2007837/93‑1 (MQ=255)
                                                               gACCAGGCGCAGAGCTTTACGCCTGGTGTGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGtt                          <  2:7986101/93‑1 (MQ=255)
                                                               gACCAGGCGCAGAGCTTTACGCCTGGTGTGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGtt                          >  1:8463718/1‑93 (MQ=255)
                                                                aCCAGGCGCAGAGCTTTACGCCTGGTGTGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTc                         <  2:4217755/93‑1 (MQ=255)
                                                                aCCAGGCGCAGAGCTTTACGCCTGGTGTGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTc                         <  2:2497957/93‑1 (MQ=255)
                                                                 ccAGGCGCAGAGCTTTACGCCTGGTGTGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCa                        <  2:5669839/93‑1 (MQ=255)
                                                                  cAGGCGCAGAGCTTTACGCCTGGTGTGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCAc                       <  1:2337189/93‑1 (MQ=255)
                                                                  cAGGCGCAGAGCTTTACGCCTGGTGTGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCAc                       <  2:5490382/93‑1 (MQ=255)
                                                                   aGGCGCAGAGCTTTACGCCTGGTGTGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCAcc                      <  2:7790661/93‑1 (MQ=255)
                                                                    ggCGCAGAGCTTTACGCCTGGTGTGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCt                     >  2:7365445/1‑93 (MQ=255)
                                                                      cgcAGAGCTTTACGCCTGGTGTGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTCAAAAGACAGTTCACCttt                   >  2:2941453/1‑93 (MQ=255)
                                                                       gcAGAGCTTTACGCTTGTTGTGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCTTTa                  <  1:7210462/93‑1 (MQ=255)
                                                                       gcAGAGATTTACACCAGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCTTTa                  <  2:8955552/93‑1 (MQ=255)
                                                                        cAGAGCTTTACGCCTGGTGTGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCTTTAc                 >  1:7046228/1‑93 (MQ=255)
                                                                        cAGAGCTTTACGCCTGGTGTGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCTTTAc                 >  1:6678380/1‑93 (MQ=255)
                                                                        cAGAGCTTTACGCCTGGTGTGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCTTTAc                 >  1:5331761/1‑93 (MQ=255)
                                                                        cAGAGCTTTACGCCTGGTGTGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCTTTAc                 >  2:3806941/1‑93 (MQ=255)
                                                                         agagCTTTACGCCTGGTGTGTAACCCTCACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCTTTACg                >  1:3508549/1‑93 (MQ=255)
                                                                         agagCTTTACGCCTGGTGTGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCTTTACg                >  2:9305482/1‑93 (MQ=255)
                                                                          gagCTTTACGCCTGGTGTGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCTTTACgg               >  2:808927/1‑93 (MQ=255)
                                                                          gagCTTTACGCCTGGTGTGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCTTTACgg               >  2:6754834/1‑93 (MQ=255)
                                                                              tttACGCCTGGTGTGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCTTTACGGACGa           >  1:1190408/1‑93 (MQ=255)
                                                                              tttACGCCTGGTGTGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCTTTACGGACGa           <  1:9920194/93‑1 (MQ=255)
                                                                               ttACGCCTGGTGTGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCTTTATGGACGAt          >  2:129770/1‑93 (MQ=255)
                                                                               ttACGCCTGGTGTGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCTTTACGGACGAt          >  2:4767278/1‑93 (MQ=255)
                                                                               ttACGCCTGGTGTGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCTTTACGGACGAt          >  1:6124682/1‑93 (MQ=255)
                                                                                tACGCCTGGTGTGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCTTTACGGACGATc         >  2:9478340/1‑93 (MQ=255)
                                                                                 aCGCCTGGTGTGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCTTTACGGACGATCa        <  1:4653169/93‑1 (MQ=255)
                                                                                   gCCTGGTGTGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCTTTACgatat            >  1:8750956/1‑83 (MQ=255)
                                                                                    ccTGGTGTGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCTTTACGGACGATCAttt     >  1:9111676/1‑93 (MQ=255)
                                                                                      tGGTGTGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCTTTACGGACGATCAttttt   >  1:2601780/1‑93 (MQ=255)
                                                                                      tGGTGTGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCTTTACGGACGATCAttttt   <  2:8383937/93‑1 (MQ=255)
                                                                                      tGGTGTGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCTTTACGGACGATCAttttt   >  1:4986270/1‑93 (MQ=255)
                                                                                      tGGTGTGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCTTTACGGACGATCAttttt   >  1:3096118/1‑93 (MQ=255)
                                                                                      tGGTGTGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCTTTACGGACGATCAttttt   >  1:3188473/1‑93 (MQ=255)
                                                                                      tGGTGTGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCTTTACGGACGATCAttttt   >  2:7216058/1‑93 (MQ=255)
                                                                                       ggTGTGTAACCCTAACTTAAGATGTTGATTCGGTGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCTTTACGGACGATCAtttttt  <  1:9675770/93‑1 (MQ=255)
                                                                                          |                                                                                         
TGGAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGGTGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCTTTACGGACGATCATTTTTT  >  NC_000913/477902‑478081

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