Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 987,560 G→C 18.1% intergenic (‑578/+25) ompF ← / ← asnS outer membrane porin 1a (Ia;b;F)/asparaginyl tRNA synthetase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913987,5600GC18.1% 304.8 / 58.5 120intergenic (‑578/+25)ompF/asnSouter membrane porin 1a (Ia;b;F)/asparaginyl tRNA synthetase
Reads supporting (aligned to +/- strand):  ref base G (41/57);  new base C (14/8);  total (55/65)
Fisher's exact test for biased strand distribution p-value = 9.62e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.20e-01

Only 100 of 127 total aligned reads displayed.

ATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTGCACGCCAGTTACGTAAGC  >  NC_000913/987461‑987659
                                                                                                   |                                                                                                   
atatTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCgg                                                                                                     >  2:1618537/1‑100 (MQ=255)
    tgcggaaCAGGATGATCGTAATGGCGGGGTTGCAAATAGATTGAATGACTTAACAACCAGATAACTTGAAAAAAAAAAAAAAAAAAGAGCCAGCggg                                                                                                    <  1:233240/92‑1 (MQ=255)
    tCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGgcgc                                                                                                 <  2:2743918/100‑1 (MQ=255)
     ccGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCt                                                                                                <  2:189659/100‑1 (MQ=255)
     ccGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCt                                                                                                >  1:1553295/1‑100 (MQ=255)
        tAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGc                                                                                             <  2:2409156/100‑1 (MQ=255)
        tAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGc                                                                                             <  2:1807838/100‑1 (MQ=255)
        tAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGc                                                                                             <  2:1898058/100‑1 (MQ=255)
        tAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGc                                                                                             <  1:1670716/100‑1 (MQ=255)
        tAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGCCCGCTGGc                                                                                             <  2:2355101/100‑1 (MQ=255)
         aaaaGGATGATCGTTATAGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAAATAAAAAAAGAGCCAGCGGGCGCTGGct                                                                                            <  1:1007508/97‑1 (MQ=255)
         aaCAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGct                                                                                            >  1:195473/1‑100 (MQ=255)
          aCAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGctc                                                                                           >  1:2534060/1‑100 (MQ=255)
           cAGGATGATCGTTATGGCGGGGTTGCAAAGAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGctct                                                                                          <  1:2202588/100‑1 (MQ=255)
             ggATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTg                                                                                        <  1:1726733/100‑1 (MQ=255)
              gatgatCGTTATGGGTGAGTTGCAAAAAGATAAATTTACGTAACAAGCTTATAACTTGAAAAAAAAAAAAAAAAAAGAGCCAGCGCCCGCTGGctct                                                                                          <  1:223908/97‑1 (MQ=255)
              gatgatCGTTATGGCGGGGTTGCAAATAGATTGTTTGACTGAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGa                                                                                       <  1:535618/100‑1 (MQ=255)
              gatgatCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGCCCGCTGGCTCtttt                                                                                       <  1:1820099/100‑3 (MQ=255)
              gatgatCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGCCCGCTGGCTCtttt                                                                                       <  2:575084/100‑3 (MQ=255)
               atgatCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGaa                                                                                      <  1:315508/100‑1 (MQ=255)
                tgatCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAg                                                                                     <  1:2361314/100‑1 (MQ=255)
                 gatCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGa                                                                                    >  2:2558208/1‑100 (MQ=255)
                  atCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGCCCGCTGGCTCttttttta                                                                                   <  2:750436/100‑7 (MQ=38)
                   tCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGCCCGCTGGCTCtttttttat                                                                                  <  1:2690098/100‑8 (MQ=38)
                    cGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGc                                                                                 >  2:548689/1‑100 (MQ=255)
                       tATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCt                                                                                <  1:2594850/98‑1 (MQ=255)
                        aTGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTaaa                                                                             >  1:1154563/1‑100 (MQ=255)
                        aTGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTaaa                                                                             >  2:1509031/1‑100 (MQ=255)
                         tGGCGGGGTTGCAAATAGATTACTTGACTTAACACCCAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAACGGGCGCTGGCTCTTGAAGACGCTAAAt                                                                            >  1:1026152/1‑100 (MQ=255)
                          ggCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATATCTTGCCAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAAtt                                                                           <  1:1919235/100‑1 (MQ=255)
                           gCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAAtn                                                                           <  1:1372160/99‑2 (MQ=255)
                              gggTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGaa                                                                       >  1:2660077/1‑100 (MQ=255)
                              gggTTGAAAATAGATTGCTTTACTTACCAAACAGATAACTTGACACAAAAGATAAAAAAAGAGTCAGCGGTCGCTGGCTCTTTAAGACGCTAAATTAg                                                                         <  1:2449762/98‑1 (MQ=255)
                               ggTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAg                                                                      >  2:71995/1‑100 (MQ=255)
                                    cAAATAGATTGCTTGACTTAACACACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGc                                                                 >  2:2470142/1‑100 (MQ=255)
                                    cAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGc                                                                 >  2:1892479/1‑100 (MQ=255)
                                    cAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGc                                                                 >  2:189616/1‑100 (MQ=255)
                                    cAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGc                                                                 <  1:2114265/100‑1 (MQ=255)
                                    cAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGc                                                                 >  1:99127/1‑100 (MQ=255)
                                       aTAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGtt                                                              >  1:833493/1‑100 (MQ=255)
                                             tGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCgg                                                        <  2:1267330/100‑1 (MQ=255)
                                             tGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCgg                                                        <  1:2359282/100‑1 (MQ=255)
                                             tGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCgg                                                        <  2:2723648/100‑1 (MQ=255)
                                             tGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCgg                                                        <  1:115178/100‑1 (MQ=255)
                                              gCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTCAAGACGCTAAATTAGAAGCTGGCGTTACGCGGa                                                       <  1:1744766/100‑1 (MQ=255)
                                               cTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAg                                                      >  2:2597932/1‑100 (MQ=255)
                                                ttGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGt                                                     <  2:2173704/100‑1 (MQ=255)
                                                 tGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTa                                                    <  2:1832904/100‑1 (MQ=255)
                                                 tGAATTAACAAACCTAGAACTTGTCAGGAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTTAATTAGAAGCTGGCGTTACGCGTAg                                                      <  1:2327714/98‑1 (MQ=255)
                                                   aCTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGt                                                     >  1:2651898/1‑97 (MQ=255)
                                                   aCTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACg                                                  <  1:2777049/100‑1 (MQ=255)
                                                    cTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGt                                                 <  2:2547107/100‑1 (MQ=255)
                                                      tAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTgg                                               <  2:28838/100‑1 (MQ=255)
                                                         caaacaGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGaa                                            >  1:1108442/1‑100 (MQ=255)
                                                              aGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGaa                                       >  1:738821/1‑100 (MQ=255)
                                                              aGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGaa                                       <  1:1407519/100‑1 (MQ=255)
                                                               gATAACTTGACAGAAAAGATAAAAAAAGAGCCATCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAat                                      <  2:622539/100‑1 (MQ=255)
                                                                 tAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAatca                                    <  2:1054821/100‑1 (MQ=255)
                                                                 tAACTGGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAatc                                     <  1:1095051/99‑1 (MQ=255)
                                                                  aagttaaCAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGCCTCTTGAAGACGCCAAATAAGAAGCCGGGCGTACGCGGACTACGTGGGAACGGAatca                                    <  1:1246852/93‑1 (MQ=255)
                                                                  aaCTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAatca                                    <  1:1469875/99‑1 (MQ=255)
                                                                   aCTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAatcaca                                  <  1:1930251/100‑1 (MQ=255)
                                                                   aCTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAatcaca                                  >  2:678942/1‑100 (MQ=255)
                                                                   aCTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGAGAACGGAatcaca                                  >  1:320008/1‑100 (MQ=255)
                                                                     ttGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTCGCGTTACGCGGAGTACGTGGGAACGGAatcacatc                                >  1:715582/1‑100 (MQ=255)
                                                                         cAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCCGAGTACGTGGGAACGGAATCACATCacgt                            >  2:1290909/1‑100 (MQ=255)
                                                                               aGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGttt                      <  1:1892470/100‑1 (MQ=255)
                                                                               aGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGttt                      >  1:2132697/1‑100 (MQ=255)
                                                                                  tAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTGc                   >  2:1344267/1‑100 (MQ=255)
                                                                                  tAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTGc                   <  2:432896/100‑1 (MQ=255)
                                                                                  tAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTGc                   <  1:2148640/100‑1 (MQ=255)
                                                                                   aaaaaaaGAGCCAGCGGGCGCTGGCGCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTGCa                  <  1:2717013/100‑1 (MQ=255)
                                                                                    cttcaaGAGCCAGCGCCCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTGCAc                 >  2:767106/5‑100 (MQ=39)
                                                                                    aaaaaaGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTGCAc                 >  1:2193637/1‑100 (MQ=255)
                                                                                     aaaaaGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTGCACg                <  1:2541468/100‑1 (MQ=255)
                                                                                     aaaaaGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTGCACg                >  2:647958/1‑100 (MQ=255)
                                                                                     aaaaaGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTGCACg                <  2:1154568/100‑1 (MQ=255)
                                                                                      tcaaGAGCCAGCGCCCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTGCACGc               >  2:2010508/3‑100 (MQ=255)
                                                                                      aaaaGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTGCACGc               >  1:1827385/1‑100 (MQ=255)
                                                                                       tccgatctaGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTGCACGcc              <  2:212913/92‑1 (MQ=255)
                                                                                       caaGAGCCAGCGCCCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTGCACGc               >  2:820305/2‑99 (MQ=255)
                                                                                       caaGAGCCAGCGCCCGCTGGCTATTGAAGACGCTAAATTAGAAGCTCGCGTTACGCGGATTACGTGGGAACGGAATCACATCACGTACGTTTTGCACGcc              >  1:2816763/2‑100 (MQ=255)
                                                                                         agagCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTGCACGCCAg            <  1:384977/100‑1 (MQ=255)
                                                                                         agagCCAGCGCCCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTGCACGCCAg            >  2:435413/1‑100 (MQ=255)
                                                                                         agagCCAGCGCCCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTGCACGCCAg            >  2:844478/1‑100 (MQ=255)
                                                                                         agagCCAGCGCCCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTGCACGCCAg            >  1:1558994/1‑100 (MQ=255)
                                                                                         agagCCAGCGCCCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTGCACGCCAg            >  2:2696498/1‑100 (MQ=255)
                                                                                         agagCCAGCGCCCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTGCACGCCAg            >  2:2512560/1‑100 (MQ=255)
                                                                                         agagCCAGCGCCCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTGCACGCCAg            >  2:1727889/1‑100 (MQ=255)
                                                                                         agagCCAGCGCCCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTGCACGCCAg            >  2:1344769/1‑100 (MQ=255)
                                                                                           agCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTGCACGCCAGtt          <  2:177939/100‑1 (MQ=255)
                                                                                            gCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTGCACGCCAGTTa         <  1:1509025/100‑1 (MQ=255)
                                                                                            gCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTGCACGCCAGTTa         <  2:738824/100‑1 (MQ=255)
                                                                                            gCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTGCACGCCAGTTa         >  2:1031829/1‑100 (MQ=255)
                                                                                             tcAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTGCACGCCAGTTAc        <  2:1288150/99‑1 (MQ=255)
                                                                                              cAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTGCACGCCAGTTAca       >  1:1288145/1‑99 (MQ=255)
                                                                                               aGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTGCACGCCAgatcgga      >  1:212912/1‑94 (MQ=255)
                                                                                                   ggCGCTGGCTGTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTGCACGCCAGTTACGTAAGc  <  2:19829/100‑1 (MQ=255)
                                                                                                   ggCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTGCACGCCAGTTACGTAAGc  >  2:429166/1‑100 (MQ=255)
                                                                                                   ggCGCTGGCTCTTGAAGACGCGAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTGCACGCCAGTTACGTAAGc  <  1:2626837/100‑1 (MQ=255)
                                                                                                   |                                                                                                   
ATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTGCACGCCAGTTACGTAAGC  >  NC_000913/987461‑987659

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