Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,570,594 G→A A878V (GCG→GTG)  pqqL ← putative peptidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,570,5940GA99.5% 665.1 / ‑5.3 189A878V (GCG→GTG) pqqLputative peptidase
Reads supporting (aligned to +/- strand):  ref base G (1/0);  new base A (111/77);  total (112/77)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.69e-01

Only 100 of 190 total aligned reads displayed.

TTCTCTACCGTCATTTGTTTCAACAATTGCTCCTGCTCAGTCCATGCTGCAGGATCGTCATATTGAATAAGACTATTTACAATAGTGTTCGCTAATTGTTGAACGCTACGCTGTTGGATATCGAGGCTGCGCTGAACGTTTTGCTGGTATTCATTCAGTTCTTGCTCACTGATCCCTTTA  >  NC_000913/1570504‑1570683
                                                                                          |                                                                                         
ttCTCTACCGTCATTTGTTTCAACAATTGCTCCTGCTCAGTCCATGCTGCAGGATCGTCATATTGAATAAGACTATTTACAATAGTGTTCACt                                                                                         <  1:128444/93‑1 (MQ=255)
ttCTCTACCGTCATTTGTTTCAACAATTGCTCCTGCTCAGTCCATGCTGCAGGATCGTCATATTGAATAAGACTATTTACAATAGTGTTCACt                                                                                         <  2:2319633/93‑1 (MQ=255)
ttCTCTACCGTCATTTGTTTCAACAATTGCTCCTGCTCAGTCCATGCTGCAGGATCGTCATATTGAATAAGACTATTTACAATAGTGTTCACt                                                                                         >  2:498711/1‑93 (MQ=255)
ttCTCTACCGTCATTTGTTTCAACAATTGCTCCTGCTCAGTCCATGCTGCAGGATCGTCATATTGAATAAGACTATTTACAATAGTGTTCACt                                                                                         >  2:579695/1‑93 (MQ=255)
 tctcTACCGTCATTTGTTTCAACAATTGCTCCTGCTCAGTCCATGCTGCAGGATCGTCATATTGAATAAGACTATTTACAATAGTGTTCACTa                                                                                        >  2:2358462/1‑93 (MQ=255)
   tctACCGTCATTTGTTTCAACAATTGCTCCTGCTCAGTCCATGCTGCAGGATCGTCATATTGAATAAGACTATTTACAATAGTGTTCACTAAt                                                                                      <  1:4197396/93‑1 (MQ=255)
   tctACCGTCATTTGTTTCAACAATTGCTCCTGCTCAGTCCATGCTGCAGGATCGTCATATTGAATAAGACTATTTACAATAGTGTTCACTAAt                                                                                      >  2:1632590/1‑93 (MQ=255)
      aCCGTCATTTGTTTCAACAATTGCTCCTGCTCAGTCCATGCTGCAGGATCGTCATATTGAATAAGACTATTTACAATAGTGTTCACTAAttgt                                                                                   >  2:327832/1‑93 (MQ=255)
      aCCGTCATTTCTTTCAACAAATGCTCCTGCTCAGTCCATGCTCCAGGATCGTCATATTAAATAAGACTATTTACAATAGTGTTCACTAAttgt                                                                                   <  1:3500406/93‑1 (MQ=255)
       ccGTCATTTGTTTCAACAATTGCTCCTGCTCAGTCCATGCTGCAGGATCGTCATATTGAATAAGACTATTTACAATAGTGTTCACTAAttgtt                                                                                  >  1:2699646/1‑93 (MQ=255)
       ccGTCATTTGTTTCAACAATTGCTCCTGCTCAGTCCATGCTGCAGGATCGTCATATTGAATAAGACTATTTACAATAGTGTTCACTAAttgtt                                                                                  >  1:2153489/1‑93 (MQ=255)
         gTCATTTGTTTCAACAATTGCTCCTGCTCAGTCCATGCTGCAGGATCGTCATATTGAATAAGACTATTTACAATAGTGTTCACTAATTGTTGa                                                                                <  2:1318057/93‑1 (MQ=255)
         gTCATTTGTTTCAACAATTGCTCCTGCTCAGTCCATGCTGCAGGATCGTCATATTGAATAAGACTATTTACAATAGTGTTCACTAATTGTTGa                                                                                >  1:670928/1‑93 (MQ=255)
         gTCATTTGTTTCAACAATTGCTCCTGCTCAGTCCATGCTGCAGGATCGTCATATTGAATAAGACTATTTACAATAGTGTTCACTAATTGTTGa                                                                                >  2:3098565/1‑93 (MQ=255)
         gTCATTTGTTTCAACAATTGCTCCTGCTCAGTCCATGCTGCAGGATCGTCATATTGAATAAGACTATTTACAATAGTGTTCACTAATTGTTGa                                                                                <  1:3815253/93‑1 (MQ=255)
          tCATTTGTTTCAACAATTGCTCCTGCTCAGTCCATGCTGCAGGATCGTCATATTGAATAAGACTATTTACAATAGTGTTCACTAATTGTTGaa                                                                               >  1:2343182/1‑93 (MQ=255)
           cATTTGTTTCAACAATTGCTCCTGCTCAGTCCATGCTGCAGGATCGTCATATTGAATAAGACTATTTACAATAGTGTTCACTAATTGTTGAac                                                                              <  1:4505908/93‑1 (MQ=255)
            aTTTGTTTCAACAATTGCTCCTGCTCAGTCCATGCTGCAGGATCGTCATATTGAATAAGACTATTTACAATAGTGTTCACTAATTGTTGAacg                                                                             >  1:2403181/1‑93 (MQ=255)
            aTTTGTTTCAACAATTGCTCCTGCTCAGTCCATGCTGCAGGATCGTCATATTGAATAAGACTATTTACAATAGTGTTCACTAATTGTTGAacg                                                                             >  2:3752587/1‑93 (MQ=255)
            aTTTGTTTCAACAATTGCTCCTGCTCAGTCCATGCTGCAGGATCGTCATATTGAATAAGACTATTTACAATAGTGTTCACTAATTGTTGAacg                                                                             >  1:4265997/1‑93 (MQ=255)
            aTTTGTTTCAACAATTGCTCCTGCTCAGTCCATGCTGCAGGATCGTCATATTGAATAAGACTATTTACAATAGTGTTCACTAATTGTTGAacg                                                                             <  1:3161159/93‑1 (MQ=255)
             tttGTTTCAACAATTGCTCCTGCTCAGTCCATGCTGCAGGATCGTCATATTGAATAAGACTATTTACAATAGTGTTCACTAATTGTTGAacgc                                                                            >  1:353866/1‑93 (MQ=255)
             tttGTTTCAACAATTGCTCCTGCTCAGTCCATGCTGCAGGATCGTCATATTGAATAAGACTATTTACAATAGTGTTCACTAATTGTTGAacgc                                                                            >  2:503347/1‑93 (MQ=255)
               tGTTTCAACAATTGCTCCTGCTCAGTCCATGCTGCAGGATCGTCATATTGAATAAGACTATTTACAATAGTGTTCACTAATTGTTGAacgcta                                                                          >  2:2571591/1‑93 (MQ=255)
               tGTTTCAACAATTGCTCCTGCTCAGTCCATGCTGCAGGATCGTCATATTGAATAAGACTATTTACAATAGTGTTCACTAATTGTTGAacgcta                                                                          <  1:3967834/93‑1 (MQ=255)
                ttttCAACAATTGCTCCTGCTCAGTCCATGCTGCAGGATCGTCATATTGAATAAGACTATTTACAATAGTGTTCACTAATTGTTGAacgctac                                                                         <  1:1644191/92‑1 (MQ=255)
                gTTTCAACAATTGCTCCTGCTCAGTCCATGCTGCAGGATCGTCATATTGAATAAGACTATTTACAATAGTGTTCACTAATTGTTGAacgctac                                                                         <  1:3624877/93‑1 (MQ=255)
                  ttCAACAATTGCTCCTGCTCAGTCCATGCTGCAGGATCGTCATATTGAATAAGACTATTTACAATAGTGTTCACTAATTGTTGAacgctacgc                                                                       <  2:477171/93‑1 (MQ=255)
                   tCAACAATTGCTCCTGCTCAGTCCATGCTGCAGGATCGTCATATTGAATAAGACTATTTACAATAGTGTTCACTAATTGTTGAacgctacgct                                                                      <  1:4428564/93‑1 (MQ=255)
                   tCAACAATTGCTCCTGCTCAGTCCATGCTGCAGGATCGTCATATTGAATAAGACTATTTACAATAGTGTTCACTAATTGTTGAacgctacgct                                                                      >  1:4352/1‑93 (MQ=255)
                   tCAACAATTGCTCCTGCTCAGTCCATGCTGCAGGATCGTCATATTGAATAAGACTATTTACAATAGTGTTCACTAATTGTTGAacgctacgct                                                                      >  1:312731/1‑93 (MQ=255)
                     aacaaTTGCTCCTGCTCAGTCCATGCTGCAGGATCGTCATATTGAATAAGACTATTTACAATAGTGTTCACTAATTGTTGAACGCTACGCTGt                                                                    <  2:1553163/93‑1 (MQ=255)
                       caaTTGCTCCTGCTCAGTCCATGCTGCAGGATCGTCATATTGAATAAGACTATTTACAATAGTGTTCACTAATTGTTGAACGCTACGCTGTTg                                                                  <  1:2954020/93‑1 (MQ=255)
                       caaTTGCTCCTGCTCAGTCCATGCTGCAGGATCGTCATATTGAATAAGACTATTTACAATAGTGTTCACTAATTGTTGAACGCTACGCTGTTg                                                                  >  1:1882756/1‑93 (MQ=255)
                            tcTCCTGCTCAGTCCATGCTGCAGGATCGTCATATTGAATAAGACTATTTACAATAGTGTTCACTAATTGTTGAACGCTACGCTGTTGGatat                                                             <  1:3896927/92‑1 (MQ=255)
                              tCCTGCTCAGTCCATGCTGCAGGATCGTCATATTGAATAAGACTATTTACAATAGTGTTCACTAATTGTTGAACGCTACGCTGTTGGATATCg                                                           <  1:2481262/93‑1 (MQ=255)
                              tCCTGCTCAGTCCATGCTGCAGGATCGTCATATTGAATAAGACTATTTACAATAGTGTTCACTAATTGTTGAACGCTACGCTGTTGGATATCg                                                           >  1:965484/1‑93 (MQ=255)
                                 tGCTCAGTCCATGCTGCAGGATCGTCATATTGAATAAGACTATTTACAATAGTGTTCACTAATTGTTGAACGCTACGCTGTTGGATATCGAgg                                                        >  1:1758803/1‑93 (MQ=255)
                                 tGCTCAGTCCATGCTGCAGGATCGTCATATTGAATAAGACTATTTACAATAGTGTTCACTAATTGTTGAACGCTACGCTGTTGGATATCGAgg                                                        >  1:474065/1‑93 (MQ=255)
                                 cgCTCAGTCCATTCTGCAGGATCGTCATATTGAATAAGACTATTTACAATAGTGTTCACTAATTGTTGAACGCTACGCTGTTGGATATCGAgg                                                        <  1:1670966/92‑1 (MQ=255)
                                    tCAGTCCATGCTGCAGGATCGTCATATTGAATAAGACTATTTACAATAGTGTTCACTAATTGTTGAACGCTACGCTGTTGGATATCGCGGCTg                                                     >  1:742788/1‑93 (MQ=255)
                                    tCAGTCCATGCTGCAGGATCGTCATATTGAATAAGACTATTTACAATAGTGTTCACTAATTGTTGAACGCTACGCTGTTGGATATCGAGGCTg                                                     <  2:561360/93‑1 (MQ=255)
                                       gTCCATGCTGCAGGATCGTCATATTGAATAAGACTATTTACAATAGTGTTCACTAATTGTTGAACGCTACGCTGTTGGATATCGAGGCTgcgc                                                  >  1:2770292/1‑93 (MQ=255)
                                           aTGCTGCAGGATCGTCATATTGAATAAGACTATTTACAATAGTGTTCACTAATTGTTGAACGCTACGCTGTTGGATATCGAGGCTGCGCTGaa                                              >  1:3581873/1‑93 (MQ=255)
                                             gctgcAGGATCGTCATATTGAATAAGACTATTTACAATAGTGTTCACTAATTGTTGAACGCTACGCTGTTGGATATCGAGGCTGCGCTGAACg                                            <  1:144781/93‑1 (MQ=255)
                                             gctgcAGGATCGTCATATTGAATAAGACTATTTACAATAGTGTTCACTAATTGTTGAACGCTACGCTGTTGGATATCGAGGCTGCGCTGAACg                                            >  2:2547067/1‑93 (MQ=255)
                                              ctgcAGGATCGTCATATTGAATAAGACTATTTACAATAGTGTTCACTAATTGTTGAACGCTACGCTGTTGGATATCGAGGCTGCGCTGAACGt                                           <  2:2855761/93‑1 (MQ=255)
                                               tgcAGGATCGTCATATTGAATAAGACTATTTACAATAGTGTTCACTAATTGTTGAACGCTACGCTGTTGGATATCGAGGCTGCGCTGAACGtt                                          >  2:1085824/1‑93 (MQ=255)
                                                gcAGGATCGTCATATTGAATAAGACTATTTACAATAGTGTTCACTAATTGTTGAACGCTACGCTGTTGGATATCGAGGCTGCGCTGAACGttt                                         >  1:2046849/1‑93 (MQ=255)
                                                 cAGGATCGTCATATTGAATAAGACTATTTACAATAGTGTTCACTAATTGTTGCACGCTACGCTGTTGGATATCGAGGCTGCGCTGAACGtttt                                        >  2:4503767/1‑93 (MQ=255)
                                                  aGGATCGTCATATTGAATAAGACTATTTACAATAGTGTTCACTAATTGTTGAACGCTACGCTGTTGGATATCGAGGCTGCGCTGAACGTTTTg                                       >  2:2471314/1‑93 (MQ=255)
                                                   ggATCGTCATATTGAATAAGACTATTTACAATAGTGTTCACTAATTGTTGAACGCTACGCTGTTGGATATCGAGGCTGCGCTGAACGTTTTGc                                      >  2:220713/1‑93 (MQ=255)
                                                   ggATCGTCATATTGAATAAGACTATTTACAATAGTGTTCACTAATTGTTGAACGCTACGCTGTTGGATATCGAGGCTGCGCTGAACGTTTTGc                                      >  2:1751000/1‑93 (MQ=255)
                                                   ggATCGTCATATTGAATAAGACTATTTACAATAGTGTTCACTAATTGTTGAACGCTACGCTGTTGGATATCGAGGCTGCGCTGAACGTTTTGc                                      >  2:3442188/1‑93 (MQ=255)
                                                         tCATATTGAATAAGACTATTTACAATAGTGTTCACTAATTGTTGAACGCTACGCTGTTGGATATCGAGGCTGCGCTGAACGTTTTGCTGGTat                                >  2:4663302/1‑93 (MQ=255)
                                                          cATATTGAATAAGACTATTTACAATAGTGTTCACTAATTGTTGAACGCTACGCTGTTGGATATCGAGGCTGCGCTGAACGTTTTGCTGGTatt                               <  1:2313460/93‑1 (MQ=255)
                                                           atatTGAATAAGACTATTTACAATAGTGTTCACTAATTGTTGAACGCTACGCTGTTGGATATCGAGGCTGCGCTGAACGTTTTGCTGGTattc                              >  1:2139117/1‑93 (MQ=255)
                                                           atatTGAATAAGACTATTTACAATAGTGTTCACTAATTGTTGAACGCTACGCTGTTGGATATCGAGGCTGCGCTGAACGTTTTGCTGGTattc                              >  2:608525/1‑93 (MQ=255)
                                                            tatTGAATAAGACTATTTACAATAGTGTTCACTAATTGTTGAACGCTACGCTGTTGGATATCGAGGCTGCGCTGAACGTTTTGCTGGTattca                             <  2:4042350/93‑1 (MQ=255)
                                                             atTGAATAAGACTATTTACAATAGTGTTCACTAATTGTTGAACGCTACGCTGTTGGATATCGAGGCTGCGCTGAACGTTTTGCTGGTattcat                            >  1:908486/1‑93 (MQ=255)
                                                             atTGAATAAGACTATTTACAATAGTGTTCACTAATTGTTGAACGCTACGCTGTTGGATATCGAGGCTGCGCTGAACGTTTTGCTGGTattcat                            >  2:2716448/1‑93 (MQ=255)
                                                               tGAATAAGACTATTTACAATAGTGTTCACTAATTGTTGAACGCTACGCTGTTGGATATCGAGGCTGCGCTGAACGTTTTGCTGGTattcattc                          <  2:1008356/93‑1 (MQ=255)
                                                               tGAATAAGACTATTTACAATAGTGTTCACTAATTGTTGAACGCTACGCTGTTGGATATCGAGGCTGCGCTGAACGTTTTGCTGGTattcattc                          <  2:3189836/93‑1 (MQ=255)
                                                               tGAATAAGACTATTTACAATAGTGTTCACTAATTGTTGAACGCTACGCTGTTGGATATCGAGGCTGCGCTGAACGTTTTGCTGGTattcattc                          <  2:148023/93‑1 (MQ=255)
                                                                gAATAAGATTATATAAAATAGTGTTAACTAATTGTTGAACTCTACGCTGTTGGATATCGAGGCTGCGCTGAACGTTTTGCTGGTAttcattca                         <  1:861784/93‑1 (MQ=255)
                                                                gAATAAGACTATTTACAATAGTGTTCACTAATTGTTGAACGCTACGCTGTTGGATATCGAGGCTGCGCTGAACGTTGTGCTGGTAttcattca                         >  1:4220226/1‑93 (MQ=255)
                                                                 aaTAAGACTATTTACAATAGTGTTCACTAATTGTTGAACGCTACGCTGTTGGATATCGAGGCTGCGCTGAACGTTTTGCTGGTATTCATTCAg                        <  2:4736543/93‑1 (MQ=255)
                                                                 aaTAAGACTATTTACAATAGTGTTCACTAATTGTTGAACGCTACGCTGTTGGATATCGAGGCTGCGCTGAACGTTTTGCTGGTATTCATTCAg                        >  2:2025073/1‑93 (MQ=255)
                                                                  aTAAGACTATTTACAATAGTGTTCACTAATTGTTGAACGCTACGCTGTTGGATATCGAGGCTGCGCTGAACGTTTTGCTGGTATTCATTCAGt                       >  1:2835039/1‑93 (MQ=255)
                                                                   tAAGACTATTTACAATAGTGTTCACTAATTGTTGAACGCTACGCTGTTGGATATCGAGGCTGCGCTGAACGTTTTGCTGGTATTCATTCAGtt                      <  1:2714028/93‑1 (MQ=255)
                                                                     aGACTATTTACAATAGTGTTCACTAATTGTTGAACGCTACGCTGTTGGATATCGAGGCTGCGCTGAACGTTTTGCTGGTATTCATTCAGTTCt                    <  1:4303494/93‑1 (MQ=255)
                                                                     aGACTATTTACAATAGTGTTCACTAATTGTTGAACGCTACGCTGTTGGATATCGAGGCTGCGCTGAACGTTTTGCTGGTATTCATTCAGTTCt                    <  1:4693943/93‑1 (MQ=255)
                                                                     aGACTATTTACAATAGTGTTCACTAATTGTTGAACGCTACGCTGTTGGATATCGAGGCTGCGCTGAACGTTTTGCTGGTATTCATTCAGTTCt                    <  1:3364884/93‑1 (MQ=255)
                                                                      gACTATTTACAATAGTGTTCACTAATTGTTGAACGCTACGCTGTTGGATATCGAGGCTGCGCTGAACGTTTTGCTGGTATTCATTCAGTTCtt                   <  2:1949435/93‑1 (MQ=255)
                                                                       aCTATTTACAATAGTGTTCACTAATTGTTGAACGCTACGCTGTTGGATATCGAGGCTGCGCTGAACGTTTTGCTGGTATTCATTCAGTTCTTg                  >  1:3045228/1‑93 (MQ=255)
                                                                         tATTTACAATAGTGTTCACTAATTGTTGAACGCTACGCTGTTGGATATCGAGGCTGCGCTGAACGTTTTGCTGGTATTCATTCAGTTCTTGCt                <  2:4292385/93‑1 (MQ=255)
                                                                         tATTTACAATAGTGTTCACTAATTGTTGAACGCTACGCTGTTGGATATCGAGGCTGCGCTGAACGTTTTGCTGGTATTCATTCAGTTCTTGCt                >  2:623189/1‑93 (MQ=255)
                                                                         tATTTACAATAGTGTTCACTAATTGTTGAACGCTACGCTGTTGGATATCGAGGCTGCGCTGAACGTTTTGCTGGTATTCATTCAGTTCTTGCt                >  1:1928513/1‑93 (MQ=255)
                                                                          aTTTACAATAGTGTTCACTAATTGTTGAACGCTACGCTGTTGGATATCGAGGCTGCGCTGAACGTTTTGCTGGTATTCATTCAGTTCTTGCTc               >  1:1290284/1‑93 (MQ=255)
                                                                           tttACAATAGTGTTCACTAATTGTTGAACGCTACGCTGTTGGATATCGAGGCTGCGCTGAACGTTTTGCTGGTATTCATTCAGTTCTTGCTCa              >  2:831030/1‑93 (MQ=255)
                                                                           tttACAATAGTGTTCACTAATTGTTGAACGCTACGCTGTTGGATATCGAGGCTGCGCTGAACGTTTTGCTGGTATTCATTCAGTTCTTGCTCa              >  1:2588879/1‑93 (MQ=255)
                                                                           tttACAATAGTGTTCACTAATTGTTGAACGCTACGCTGTTGGATATCGAGGCTGCGCTGAACGTTTTGCTGGTATTCATTCAGTTCTTGCTCa              >  1:3340954/1‑93 (MQ=255)
                                                                           tttACAATAGTGTTCACTAATTGTTGAACGCTACGCTGTTGGATATCGAGGCTGCGCTGAACGTTTTGCTGGTATTCATTCAGTTCTTGCTCa              >  1:718239/1‑93 (MQ=255)
                                                                            ttACAATAGTGTTCACTAATTGTTGAACGCTACGCTGTTGGATATCGAGGCTGCGCTGAACGTTTTGCTGGTATTCATTCAGTTCTTGCTCAc             <  2:4152298/93‑1 (MQ=255)
                                                                             tACAATAGTGTTCACTAATTGTTGAACGCTACGCTGTTGGATATCGAGGCTGCGCTGAGCGTTTTGCTGGTATTCATTCAGTTCTTGCTCACt            <  1:3926404/93‑1 (MQ=255)
                                                                             tACAATAGTGTTCACTAATTGTTGAACGCTACGCTGTTGGATATCGAGGCTGCGCTGAACGTTTTGCTGGTATTCATTCAGTTCTTGCTCACt            >  2:1207613/1‑93 (MQ=255)
                                                                              aCAATAGTGTTCACTAATTGTTGAACGCTACGCTGTTGGATATCGAGGCTGCGCTGAACGTTTTGCTGGTATTCATTCAGTTCTTGCTCACTg           >  1:4455492/1‑93 (MQ=255)
                                                                              aCAATAGTGTTCACTAATTGTTGAACGCTACGCTGTTGGATATCGAGGCTGCGCTGAACGTTTTGCTGGTATTCATTCAGTTCTTGCTCACTg           >  2:1474553/1‑93 (MQ=255)
                                                                              aCAATAGTGTTCACTAATTGTTGAACGCTACGCTGTTGGATATCGAGGCTGCGCTGAACGTTTTGCTGGTATTCATTCAGTTCTTGCTCACTg           >  2:2705037/1‑93 (MQ=255)
                                                                               cAATAGTGTTCACTAATTGTTGAACGCTACGCTGTTGGATATCGAGGCTGCGCTGAACGTTTTGCTGGTATTCATTCAGTTCTTGCTCACTGa          <  2:305237/93‑1 (MQ=255)
                                                                               cAATAGTGTTCACTAATTGTTGAACGCTACGCTGTTGGATATCGAGGCTGCGCTGAACGTTTTGCTGGTATTCATTCAGTTCTTGCTCACTGa          >  1:4125423/1‑93 (MQ=255)
                                                                               cAATAGTGTTCACTAATTGTTGAACGCTACGCTGTTGGATATCGAGGCTGCGCTGAACGTTTTGCTGGTATTCATTCAGTTCTTGCTCACTGa          <  2:682869/93‑1 (MQ=255)
                                                                               cAATAGTGTTCACTAATTGTTGAACGCTACGCTGTTGGATATCGAGGCTGCGCTGAACGATTTGCTGGTATTCATTCAGTTCTTGCTCACTGa          >  2:2946807/1‑93 (MQ=255)
                                                                                aaTAGTGTTCACTAATTGTTGAACGCTACGCTGTTGGATATCGAGGCTGCGCTGAACGTTTTGCTGGTATTCATTCAGTTCTTGCTCACTGAt         <  1:4619767/93‑1 (MQ=255)
                                                                                   aGTGTTCACTAATTGTTGAACGCTACGCTGTTGGATATCGAGGCTGCGCTGAACGTTTTGCTGGTATTCATTCAGTTCTTGCTCACTGATccc      <  1:3715002/93‑1 (MQ=255)
                                                                                    gtgtTCACTAATTGTTGAACGCTACGCTGTTGGATATCGAGGCTGCGCTGAACGTTTTGCTGGTATTCATTCAGTTCTTGCTCACTGATCCCt     >  2:4371458/1‑93 (MQ=255)
                                                                                    gtgtTCACTAATTGTTGAACGCTACGCTGTTGGATATCGAGGCTGCGCTGAACGATTTGCTGGTATTCATTCAGTTCTTGCTCACTGATCCCt     >  1:1710713/1‑93 (MQ=255)
                                                                                     tgtTCACTAATTGTTGAACGCTACGCTGTTGGATATCGAGGCTGCGCTGAACGTTTTGCTGGTATTCATTCAGTTCTTGCTCACTGATCCCtt    <  1:2750643/93‑1 (MQ=255)
                                                                                      gtTCACTAATTGTTGAACGCTACGCTGTTGGATATCGAGGCTGCGCTGAACGTTTTGCTGGTATTCATTCAGTTCTTGCTCACTGATCCCttt   <  2:160079/93‑1 (MQ=255)
                                                                                      gtTCACTAATTGTTGAACGCTACGCTGTTGGATATCGAGGCTGCGCTGAACGTTTTGCTGGTATTCATTCAGTTCTTGCTCACTGATCCCttt   >  1:2446021/1‑93 (MQ=255)
                                                                                          |                                                                                         
TTCTCTACCGTCATTTGTTTCAACAATTGCTCCTGCTCAGTCCATGCTGCAGGATCGTCATATTGAATAAGACTATTTACAATAGTGTTCGCTAATTGTTGAACGCTACGCTGTTGGATATCGAGGCTGCGCTGAACGTTTTGCTGGTATTCATTCAGTTCTTGCTCACTGATCCCTTTA  >  NC_000913/1570504‑1570683

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: