Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 987,561 G→C 17.8% intergenic (‑579/+24) 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,5610GC17.8% 299.3 / 56.7 118intergenic (‑579/+24)ompF/asnSouter membrane porin 1a (Ia;b;F)/asparaginyl tRNA synthetase
Reads supporting (aligned to +/- strand):  ref base G (40/56);  new base C (14/7);  total (54/64)
Fisher's exact test for biased strand distribution p-value = 5.25e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.24e-01

Only 100 of 127 total aligned reads displayed.

TATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTGCACGCCAGTTACGTAAGCA  >  NC_000913/987462‑987660
                                                                                                   |                                                                                                   
 atTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGACAAAATAAAAAAAGAGCCAGCGGgc                                                                                                    >  2:468069/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:1956812/100‑1 (MQ=255)
       tAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGc                                                                                              <  2:1898058/100‑1 (MQ=255)
       tAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGc                                                                                              <  1:1670716/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)
                      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                                                                  >  1:99127/1‑100 (MQ=255)
                                   cAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGc                                                                  <  1:2114265/100‑1 (MQ=255)
                                   cAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGc                                                                  >  2:189616/1‑100 (MQ=255)
                                      aTAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGtt                                                               >  2:1693807/1‑100 (MQ=255)
                                      aTAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGtt                                                               >  1:833493/1‑100 (MQ=255)
                                            tGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCgg                                                         <  2:1267330/100‑1 (MQ=255)
                                            tGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCgg                                                         <  1:115178/100‑1 (MQ=255)
                                            tGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCgg                                                         <  2:1860737/100‑1 (MQ=255)
                                            tGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCgg                                                         <  2:2723648/100‑1 (MQ=255)
                                            tGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCgg                                                         <  1:2359282/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:2912276/1‑100 (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)
                                                                tAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAatcn                                     <  2:266375/100‑2 (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:2132697/1‑100 (MQ=255)
                                                                              aGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGttt                       <  1:1892470/100‑1 (MQ=255)
                                                                                 tAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTGc                    <  2:432896/100‑1 (MQ=255)
                                                                                 tAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTGc                    >  2:1344267/1‑100 (MQ=255)
                                                                                 tAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTGc                    >  2:2185040/1‑100 (MQ=255)
                                                                                 tAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTGc                    <  1:2148640/100‑1 (MQ=255)
                                                                                 tAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTGc                    >  2:2506927/1‑100 (MQ=255)
                                                                                  aaaaaaaGAGCCAGCGGGCGCTGGCGCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTGCa                   <  1:2717013/100‑1 (MQ=255)
                                                                                   cttcaaGAGCCAGCGCCCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTGCAc                  >  2:767106/5‑100 (MQ=39)
                                                                                   aaaaaaGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTGCAc                  >  1:2193637/1‑100 (MQ=255)
                                                                                    aaaaaGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTGCACg                 >  2:647958/1‑100 (MQ=255)
                                                                                    aaaaaGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTGCACg                 <  2:1154568/100‑1 (MQ=255)
                                                                                    aaaaaGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTGCACg                 <  1:2541468/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             >  1:1558994/1‑100 (MQ=255)
                                                                                        agagCCAGCGCCCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTGCACGCCAg             >  2:435413/1‑100 (MQ=255)
                                                                                        agagCCAGCGCCCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTGCACGCCAg             >  2:844478/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:1031829/1‑100 (MQ=255)
                                                                                           gCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTGCACGCCAGTTa          <  2:738824/100‑1 (MQ=255)
                                                                                            tcAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTGCACGCCAGTTAc         <  2:1288150/99‑1 (MQ=255)
                                                                                             cAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTGCACGCCAGTTAca        >  1:1288145/1‑99 (MQ=255)
                                                                                              aGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTGCACGCCAgatcgga       >  1:212912/1‑94 (MQ=255)
                                                                                                  ggCGCTGGCTCTTGAAGACGCGAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTGCACGCCAGTTACGTAAGc   <  1:2626837/100‑1 (MQ=255)
                                                                                                   gcgcTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTGCACGCCAGTTACGTAAGCa  <  2:56663/100‑1 (MQ=255)
                                                                                                   |                                                                                                   
TATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTGCACGCCAGTTACGTAAGCA  >  NC_000913/987462‑987660

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: