Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,470,558 G→A 20.5% pseudogene (42/3497 nt) ydbA → pseudogene, autotransporter homolog; interrupted by IS2 and IS30

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,470,5580GA20.5% 226.5 / 52.6 113pseudogene (42/3497 nt)ydbApseudogene, autotransporter homolog; interrupted by IS2 and IS30
Reads supporting (aligned to +/- strand):  ref base G (51/38);  new base A (9/14);  total (60/53)
Fisher's exact test for biased strand distribution p-value = 1.60e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.01e-02

Only 100 of 115 total aligned reads displayed.

CAAAACACCGAAAGAGATAATTGAAAGGGGTGTTGCATTGACAGATTGAATCTACACTGTTGAAATTCAGAATTTAGGTTTTGCTTTTGTTACTGGTGAAAATACAACAGGTATAAATAGTGGCACGATCTCGTTATTACAAAATGGTAAAGATCCGGCACCGTCTCCCATTGTTTTACTGGCTACTAACGGAGGGA  >  NC_000913/1470461‑1470657
                                                                                                 |                                                                                                   
cAAAACACCGAAAGAGATAATTGAAAGGGGTGTTGCATTGACAGATTGAATCTACACTGTTGAAATTCAGAATTTAGGTTTTGCTTTTGTTACTGGTGaa                                                                                                   >  2:1708683/1‑100 (MQ=255)
 aaaaCACCGAAAGAGATAATTGAAAGGGGTGTTGCATTGACAGATTGAATCTACACTGTTGAAATTCAGAATTTAGGTTTTGCTTTTGTTACTGGTGaaa                                                                                                  >  2:1193368/1‑100 (MQ=255)
 aaaaCACCGAAAGAGATAATTGAAAGGGGTGTTGCATTGACAGATTGAATCTACACTGTTGAAATTCAGAATTTAGGTTTTGCTTTTGTTACTGGTGaaa                                                                                                  >  1:1465620/1‑100 (MQ=255)
    acacCGAAAGAGATAATTGAAAGGGGTGTTGCATTGACAGATTGAATCTACACTGTTGAAATTCAGAATTTAGGTTTTGCTTTTGTTACTGGTAAAAATa                                                                                               <  2:3834092/100‑1 (MQ=255)
     cacCGAAAGAGATAATTGAAAGGGGCGTTGCATTGACAGATTGAATCTACACTGTTGAAATTCAGAATTTAGGTTTTGCTTTTGTTACTGGTGAAAATac                                                                                              >  1:3220477/1‑100 (MQ=255)
      acCGAAAGAGATTATTTAATTTAATGATGCAGTGGCAGATTAGATCTACACTGTTGAAATTCAGAATTTAGGTTTTGCTTTTGTTACTGGTGAAAATaca                                                                                             <  2:3994229/100‑1 (MQ=255)
      acCGAAAGAGATAATTGAAAGGGGTGTTGCATTGACAGATTGAATCTACACTGTTGAAATTCAGAATTTAGGTTTTGCTTTTGTTACTGGTGAAAATaca                                                                                             <  2:1550525/100‑1 (MQ=255)
      acCGAAAGAGATAATTGAAAGGGGTGTTGCATTGACAGATTGAATCTACACTGTTGAAATTCAGAATTTAGGTTTTGCTTTTGTTACTGGTGAAAATaca                                                                                             <  1:360693/100‑1 (MQ=255)
      acCGAAAGAGATAATTGAAAGGGGTGTTGCATTGACAGATTGAATCTACACTGTTGAAATTCAGAATTTAGGTTTTGCTTTTGTTACTGGTAAAAATaca                                                                                             <  2:485732/100‑1 (MQ=255)
       ccGAAAGAGATAATTGAAAGGGGTGTTGCATTGACAGATTGAATCTACACTGTTGAAATTCAGAATTTAGGTTTTGCTTTTGTTACTGGTGAAAATacaa                                                                                            >  1:175167/1‑100 (MQ=255)
         gAAAGAGATAATTGAAAGGGGTGTTGCATTGACAGATTGAATCTACACTGTTGAAATTCAGAATTTAGGTTTTGCTTTTGTTACTGGTGAAAATacaaca                                                                                          >  1:1814682/1‑100 (MQ=255)
          aaaGAGATAATTGAAAGGGGTGTTGCATTGACAGATTGAATCTACACTGTTGAAATTCAGAATTTAGGTTTTGCTTTTGTTACTGGTGAAAATACAACAg                                                                                         <  2:1690524/100‑1 (MQ=255)
             gagaTAATTGAAAGGGGTGTTGCATTGACAGATTGAATCTACACTGTTGAAATTCAGAATTTAGGTTTTGCTTTTGTTACTGGTGAAAATACAACAGGta                                                                                      <  1:3746783/100‑1 (MQ=255)
               gaTAATTGAAAGGGGTGTTGCATTGACAGATTGAATCTACACTGTTGAAATTCAGAATTTAGGTTTTGCTTTTGTTACTGGTGAAAATACAACAGGtata                                                                                    >  2:3155330/1‑100 (MQ=255)
               gaTAATTGAAAGGGGTGTTGCATTGACAGATTGAATCTACACTGTTGAAATTCAGAATTTAGGTTTTGCTTTTGTTACTGGTGAAAATACAACAGGtata                                                                                    >  2:1940944/1‑100 (MQ=255)
               gaTAATTGAAAGGGGTGTTGCATTGACAGATTGAATCTACACTGTTGAAATTCAGAATTTAGGTTTTGCTTTTGTTACTGGTGAAAATACAACAGGtata                                                                                    >  1:775447/1‑100 (MQ=255)
                aTAATTGAAAGGGGTGTTGCATTGACAGATTGAATCTACACTGTTGAAATTCAGAATTTAGGTTTTGCTTTTGTTACTGGTAAAAATACAACAGGTATaa                                                                                   >  2:3376190/1‑100 (MQ=255)
                 tAATTGAAAGGGGTGTTGCATTGACAGATTGAATCTACACTGTTGAAATTCAGAATTTAGGTTTTGCTTTTGTTACTGGTGAAAATACAACAGGTATaaa                                                                                  >  1:1904602/1‑100 (MQ=255)
                 tAATTGAAAGGGGTGTTGCATTGACAGATTGAATCTACACTGTTGAAATTCAGAATTTAGGTTTTGCTTTTGTTACTGGTGAAAATACAACAGGTATaaa                                                                                  >  1:1017267/1‑100 (MQ=255)
                    ttGAAAGGGGTGTTGCATGGACAGATTGAATGTACACTGTTGAAATTCAGAAATTAGGTTTTGCTTTTCATACTGGTGAAAATACAACAGGTATAAATAg                                                                               >  2:3958648/1‑100 (MQ=38)
                     tGAAAGGGGTGTTGCATTGACAGATTGAATCTACACTGTTGAAATTCAGAATTTAGGTTTTGATTTTGTTACTGGTAAAAATACAACAGGTATAAATAGt                                                                              <  2:2204214/100‑1 (MQ=255)
                     tGAAAGGAGTGTTGCAGTGACAAATTTAATCTACACTGTTGAAATTCAGAATTTAGGTTTTGCTTTTGTTACTGGTAAAAATACAACAGGTATAAATAGt                                                                              <  2:2801468/100‑1 (MQ=255)
                        aaGGGGTGTTGCATTGACAGATTGAATCTACACTGTTGAAATTCAGAATTTAGGTTTTGCTTTTGTTACTGGTGAAAATACAACAGGTATAAATAGTGGc                                                                           >  1:2184457/1‑100 (MQ=255)
                          ggggTGTTGCATTGACAGATTGAATCTACACTGTTGAAATTCAGAATTTAGGTTTTGCTTTTGTTACTGGTGAAAATACAACAGGTATAAATAGTGacac                                                                         >  2:845643/1‑100 (MQ=255)
                          ggggTGTTGCATTGACAGATTGAATCTACACTGTTGAAATTCAGAATTTAGGTTTTGCTTTTGTTACTGGTGAAAATACAACAGGTATAAATAGTGGCAc                                                                         >  1:3392629/1‑100 (MQ=255)
                           gggTGTTGCATTGACAGATTGAATCTACACTGTTGAAATTCTGAATTTAGGTTTTGCTTTTGTTACTGGTGAAAATACAACAGGTATAAATAGTGGCACg                                                                        <  1:1019181/100‑1 (MQ=255)
                             gtgtTGCATTGACAGATTGGATCTACACTGTTGAAATTCAGAATTTAGGTTTTGCTTTTGTTACTGGTGAAAATACAACAGGTATAAATAGTg                                                                             <  2:3911234/93‑1 (MQ=255)
                               gtTGCATTGACAGATTGAATCTACACTGTTGAAATTCAGAATTTAGGTTTTGCTTTTGTTACTGGTGAAAATACAACAGGTATAAATAGTGGCACGAtct                                                                    >  2:2262224/1‑100 (MQ=255)
                                 tGCATTGACAGATTGAATCTACACTGTTGAAATTCAGAATTTAGGTTTTGCTTTTGTTACTGGTGAAAATACAACAGGTATAAATAGTGGCACGAtctc                                                                   >  2:2758479/1‑99 (MQ=255)
                                 tGCATTGACAGATTGAATCTACACTGTTGAAATTCAGAATTTAGGTTTTGCTTTTGTTACTGGTGAAAATACAACAGGTATAAATAGTGGCACGATCTCg                                                                  <  1:1346369/100‑1 (MQ=255)
                                 tGCATTGACAGATTGAATCTACACTGTTGAAATTCAGAATTTAGGTTTTGCTTTTGTTACTGGTAAAAATACAACAGGTATAAATAGTGGCACGATCTCg                                                                  >  1:338490/1‑100 (MQ=255)
                                  ncATTGACAGATTGAATCTACACTGTTGAAATTCAGAATTTAGGTTTTGCTTTTGTTACTGGTGAAAATACAACAGGTATAAATAGTGGCACGATCTCGt                                                                 >  1:3221892/2‑100 (MQ=255)
                                   cATTGACAGATTGAATCTACACTGTTGAAATTCAGAATTTAGGTTTTGCTTTTGTTACTGGTGAAAATACAACAGGTATAAATAGTGGCACGAtctc                                                                   >  2:1957542/1‑97 (MQ=255)
                                   cATTGACAGATTGAATCTACACTGTTGAAATTCAGAATTTAGGTTTTGCTTTTGTTACTGGTGAAAATACAACAGGTATAAATAGTGGCACGATCTCGtn                                                                <  2:3221872/100‑2 (MQ=255)
                                    aTTGACAGATTGAATCTACACTGTTGAAATTCAGAATTTAGGTTTTGCTTTTGTTACTGGTGAAAATACAACAGGTATAAATAGTGGCACGATCTCGtta                                                               <  1:510583/100‑1 (MQ=255)
                                     ttGACAGATTGAATCTACACTGTTGAAATTCAGAATTTAGGTTTTGCTTTTGTTACTGGTGAAAATACAACAGGTATAAATAGTGGCACGATCTCGttat                                                              <  2:4090561/100‑1 (MQ=255)
                                        aCAGATTGAATCTACACTGTTGAAATTCAGAATTTAGGTTTTGCTTTTGTTACTGGTGAAAATACAACAGGTATAAATAGTGGCACGATCTCGTTATTAc                                                           <  1:1023887/100‑1 (MQ=255)
                                        aCAGATTGAATCTACACTGTTGAAATTCAGAATTTAGGTTTTGCTTTTGTTACTGGTGAAAATACAACAGGTATAAATAGTGGCACGATCTCGTTATTAc                                                           >  1:2631522/1‑100 (MQ=255)
                                         cAGATTGAATCTACACTGTTGAAATTCAGAATTTAGGTTTTGCTTTTGTTACTGGTGAAAATACAACAGGTATAAATAGTGGCACGATCTCGTTATTACa                                                          >  2:2750258/1‑100 (MQ=255)
                                         cAGATTGAATCTACACTGTTGAAATTCAGAATTTAGGTTTTGCTTTTGTTACTGGTGAAAATACAACAGGTATAAATAGTGGCACGATCTCGTTATTACa                                                          >  2:673916/1‑100 (MQ=255)
                                           gATTGAATCTACACTGTTGAAATTCAGAATTTAGGTTTTGCTTTTGTTACTGGTGAAAATACAACAGGTATAAATAGTGGCACGATCTCGTTATTACaaa                                                        >  2:879558/1‑100 (MQ=255)
                                            aTTGAATCTACACTGTTGAAATTCAGAATTTAGGTTTTGCTTTTGTTACTGGTGAAAATACAACAGGTATAAATAGTGGCACGATCTCGTTATTACaaaa                                                       >  1:2562039/1‑100 (MQ=255)
                                            aTTGAATCTACACTGTTGAAATTCAGAATTTAGGTTTTGCTTTTGTTACTGGTAAAAATACAACAGGTATAAATAGTGGCACGATCTCGTTATTACaaaa                                                       <  1:885390/100‑1 (MQ=255)
                                               gAATCTACACTGTTGAAATTCAGAATTTAGGTTTTGCTTTTGTTACTGGTGAAAATACAACAGGTATAAATAGTGGCACGATCTCGTTATTACAAAATgg                                                    <  2:3772013/100‑1 (MQ=255)
                                               gAATCTACACTGTTGAAATTCAGAATTTAGGTTTTGCTTTTGTTACTGGTGAAAATACAACAGGTATAAATAGTGGCACGATCTCGTTATTACAAAATg                                                     <  2:330360/99‑1 (MQ=255)
                                                aaTCTACACTGTTGAAATTCAGAATTTAGGTTTTGCTTTTGTTACTGGTGAAAATACAACAGGTATAAATAGTGGCACGATCTCGTTATTACAAAATGGt                                                   <  2:2616108/100‑1 (MQ=255)
                                                aaTCTACACTGTTGAAATTCAGAATTTAGGTTTTGCTTTTGTTACTGGTGAAAATACAACAGGTATAAATAGTGGCACGATCTCGTTATTACAAAATGGt                                                   <  2:484445/100‑1 (MQ=255)
                                                 aTCTACACTGTTGAAATTCAGAATTTAGGTTTTGCTTTTGTTACTGGTGAAAATACAACAGGTATAAATAGTGGCACGATCTCGTTATTACAAAATGGTa                                                  <  2:3111819/100‑1 (MQ=255)
                                                   cTACACTGTTGAAATTCAGAATTTAGGTTTTGCTTTTGTTACTGGTGAAAATACAACAGGTATAAATAGTGGCACGATCTCGTTATTACAAAATGGTaaa                                                <  1:214859/100‑1 (MQ=255)
                                                    tACACTGTTGAAATTCAGAATTTAGGTTTTGCTTTTGTTACTGGTGAAAATACAACAGGTATAAATAGTGGCACGATCTCGTTATTACAAAATGGTAAAg                                               >  2:1142739/1‑100 (MQ=255)
                                                      cacTGTTGAAATTCAGAATTTAGGTTTTGCTTTTGTTACTGGTGAAAATACAACAGGTATAAATAGTGGCACGATCTCGTTATTACAAAATGGTAAAGAt                                             >  2:3004215/1‑100 (MQ=255)
                                                      cacTGTTGAAATTCAGAATTTAGGTTTTGCTTTTGTTACTGGTGAAAATACAACAGGTATAAATAGTGGCACGATCTCGTTATTACAAAATGGTAAAGAt                                             >  2:3004084/1‑100 (MQ=255)
                                                       acTGTTGAAATTCAGAATTTAGGTTTTGCTTTTGTTACTGGTAAAAATACAACAGGTATAAATAGTGGCACGATCTCGTTATTACAAAATGGTAAAGATc                                            >  2:1196129/1‑100 (MQ=255)
                                                        cTGTTGAAATTCAGAATTTAGGTTTTGCTTTTGTTACTGGTGAAAATACAACAGGTATAAATAGTGGCACGATCTCGTTATTACAAAATGGTAAAGATcc                                           >  2:4085160/1‑100 (MQ=255)
                                                         tGTTGAAATTCAGAATTTAGGTTTTGCTTTTGTTACTGGTGAAAATACAACAGGTATAAATAGTGGCACGATCTCGTTATTACAAAATGGTAAAGATCCg                                          <  2:4002188/100‑1 (MQ=255)
                                                         aacctaaaTTCTGAATTTAGGTTTTGCTTTTGTTACTGGTGAAAATACAACAGGTATAAATAGTGGCACGATCTCGTTATTACAAAATGGTAAAGATCCg                                          >  1:597249/6‑100 (MQ=255)
                                                          gTTGAAATTCAGAATTTAGGTTTTGCTTTTGTTACTGGTGAAAATACAACAGGTATAAATAGTGGCACGATCTCGTTATTACAAAATGGTAAAGATCCgg                                         >  2:409178/1‑100 (MQ=255)
                                                          gTTGAAATTCAGAATTTAGGTTTTGCTTTTGTTACTGGTGAAAATACAACAGGTATAAATAGTGGCACGATCTCGTTATTACAAAATGGTAAAGATCCg                                          <  1:418535/99‑1 (MQ=255)
                                                           ttGAAATTCAGAATTTAGGTTTTGCTTTTGTTACTGGTGAAAATACAACAGGTATAAATAGTGGCACGATCTCGTTATTACAAAATGGTAAAGATCCGGc                                        >  2:1039964/1‑100 (MQ=255)
                                                               aaTTCAGAATTTAGGTTTTGCTTTTGTTACTGGTAAAAATACAACAGGTATAAATAGTGGCACGATCTCGTTATTACAAAATGGTAAAGATCCGGCACCg                                    >  2:3119823/1‑100 (MQ=255)
                                                                aTTCAGAATTTAGGTTTTGCTTTTGTTACTGGTAAAAATACAACAGGTATAAATAGTGGCACGATCTCGTTATTACAAAATGGTAAAGATCCGGCAcc                                     >  2:630895/1‑98 (MQ=255)
                                                                 ttCAGAATTTAGGTTTTGCTTTTGTTACTGGTGAAAATACAACAGGTATAAATAGTGGCACGATCTCGTTATTACAAAATGGTAAAGATCCGGCACCGtc                                  >  1:247599/1‑100 (MQ=255)
                                                                   cAGAATTTAGGTTTTGCTTTTGTTACTGGTGAAAATACAACAGGTATAAATAGTGGCACGATCTCGTTATTACAAAATGGTAAAGATCCGGCACCGtctc                                <  1:122803/100‑1 (MQ=255)
                                                                    aGAATTTAGGTTTTGCTTTTGTTACTGGTGAAAATACAACAGGTATAAATAGTGGCACGATCTCGTTATTACAAAATGGTAAAGATCCGGCACCGTCTcc                               <  2:2308837/100‑1 (MQ=255)
                                                                    aGAATTTAGGTTTTGCTTTTGTTACTGGTGAAAATACAACAGGTATAAATAGTGGCACGATCTCGTTATTACAAAATGGTAAAGATCCGGCACCGTCTcc                               <  2:3289622/100‑1 (MQ=255)
                                                                       aTTTAGGTTTTGCTTTTGTTACTGGTGAAAATACAACAGGTATAAATAGTGGCACGATCTCGTTATTACAAAATGGTAAAGATCCGGCACCGTCTCCCAt                            <  1:1039970/100‑1 (MQ=255)
                                                                         ttAGGTTTTGCTTTTGTTACTGGTGAAAATACAACAGGTATAAATAGTGGCACGATCTGGTTATTACAAAATGGTAAAGATCCGGCACCGTCTCCCATTg                          <  2:2412407/100‑1 (MQ=255)
                                                                            ggTTTTGCTTTTGTTACTGGTAAAAATACAACAGGTATAAATAGTGGCACGATCTCGTTATTACAAAATGGTAAAGATCCGGCACCGTCTCCCATTGttt                       <  1:3912440/100‑1 (MQ=255)
                                                                              ttttGCTTTTGTTACTGGTGAAAATACAACAGGTATAAATAGTGGCACGATCTCGTTATTACAAAATGGTAAAGATCCGGCACCGTCTCCCATTGTTTTa                     >  1:3480659/1‑100 (MQ=255)
                                                                               tttGCTTTTGTTACTGGTGAAAATACAACAGGTATAAATAGTGGCACGATCTCGTTATTACAAAATGGTAAAGATCCGGCACCGTCTCCCATTGTTTTAc                    >  2:1097439/1‑100 (MQ=255)
                                                                                 tGCTTTTGTTACTGGTGAAAATACAACAGGTATAAATAGTGGCACGATCTCGTTATTACAAAATGGTAAAGATCCGGCACCGTCTCCCATTGTTTTACTg                  >  2:3733769/1‑100 (MQ=255)
                                                                                 tGCTTTTGTTACTGGTGAAAATACAACAGGTATAAATAGTGGCACGATCTCGTTATTACAAAATGGTAAAGATCCGGCACCGTCTCCCATTGTTTTACTg                  <  1:1176969/100‑1 (MQ=255)
                                                                                 tGCTTTTGTTACTGGTAAAAATACAACAGGTATAAATAGTGGCACGATCTCGTTATTACAAAATGGTAAAGATCCGGCACCGTCTCCCATTGTTTTACTg                  <  1:31453/100‑1 (MQ=255)
                                                                                  gCTTTTGTTACTGGTGAAAATACAACAGGTATAAATAGTGGCACGATCTCGTTATTACAAAATGGTAAAGATCCGGCACCGTCTCCCATTGTTTTACTgg                 >  1:3069668/1‑100 (MQ=255)
                                                                                  gCTTTTGTTACTGGTAAAAATACAACAGGTATAAATAGTGGCACGATCTCGTTATTACAAAATGGTAAAGATCCGGCACCGTCTCCCATTGTTTTACTgg                 <  1:3862197/100‑1 (MQ=255)
                                                                                   cTTTTGTTACTGGTGAAAATACAACAGGTATAAATAGTGGCACGATCTCGTTATTACAAAATGGTAATGATCCGGCACCGTCTCCCATTGTTTTACTGGc                <  1:3926196/100‑1 (MQ=255)
                                                                                   cTTTTGTTACTGGTGAAAATACAACAGGTATAAATAGTGGCACGATCTCGTTATTACAAAATGGTAAAGATCCGGCACCGTCTCCCATTGTTTTACTGGc                <  2:1704198/100‑1 (MQ=255)
                                                                                    ttttGTTACTGGTGAAAATACAACAGGTATAAATAGTGGCACGATCTCGTTATTACAAAATGGTAAAGATCCGGCACCGTCTCCCATTGTTTTACTGGct               >  2:1741763/1‑100 (MQ=255)
                                                                                    ttttGTTACTAGTAAAAATACAACAGGTATAAATAGTGGCACGATCTCGTTATTACAAAATGGTAAAGATCCGGCACCGTCTCCCATTGTTTTACTGGct               <  2:3476787/100‑1 (MQ=255)
                                                                                     tttGTTACTGGTGAAAATACAACAGGTATAAATAGTGGCACGATCTCGTTATTACAAAATGGTAAAGATCCGGCACCGTCTCCCATTGTTTTACTGGcta              <  2:56818/100‑1 (MQ=255)
                                                                                      ttGTTACTGGTTAAAATACAACAGGTATAAATAGTGGCACCATCTCGTTATTACAAAATGGTAAAGATCCGGCACCGTCTCCCATTGTTTTACTGGc                <  2:3626005/97‑1 (MQ=255)
                                                                                      ttGTTACTGGTGAAAATACAACAGGTATAAATAGTGGCACGATCTCGTTATTACAAAATGGTAAAGATCCGGCACCGTCTCCCATTGTTTTACTGGctac             <  1:322893/100‑1 (MQ=255)
                                                                                      ttGTTACTGGTAAAAATACAACAGGTATAAATAGTGGCACGATCTCGTTATTACAAAATGGTAAAGATCCGGCACCGTCTCCCATTGTTTTACTGGctac             <  1:2836555/100‑1 (MQ=255)
                                                                                       tGTTACTGGTAAAAATACAACAGGTATAAATAGTGGCACGATCTCGTTATTACAAAATGGTAAAGATCCGGCACCGTCTCCCATTGTTTTACTGGctacn            <  2:1748515/100‑2 (MQ=255)
                                                                                        gTTACTGGTGAAAATACAACAGGTATAAATAGTGGCACGATCTCGTTATTACAAAATGGTAAAGATCCGGCACCGTCTCCCATTGTTTTAcccggcggga           >  2:232380/1‑91 (MQ=255)
                                                                                         ttACTGGTAAAAATACAACAGGTATAAATAGTGGCACGATCTCGTTATTACAAAATGGTAAAGATCCGGCACCGTCTCCCATTGTTTTACTGGCTACTaa          >  1:2563762/1‑100 (MQ=255)
                                                                                          tACTGGTGAAAATACAACAGGTATAAATAGTGGCACGATCTCGTTATTACAAAATGGTAAAGATCCGGCACCGTCTCCCATTGTTTTACTGGCTACTAAc         >  2:2567647/1‑100 (MQ=255)
                                                                                          tACTGGTAAAAATACAACAGGTATAAATAGTGGCACGATCTCGTTATTACAAAATGGTAAAGATCCGGCACCGTCTCCCATTGTTTTACTGGCTACTAAc         >  1:2652525/1‑100 (MQ=255)
                                                                                          tACTGGTAAAAATACAACAGGTATAAATAGTGGCACGATCTCGTTATTACAAAATGGTAAAGATCCGGCACCGTCTCCCATTGTTTTACTGGCTACTAAc         >  1:2675406/1‑100 (MQ=255)
                                                                                           aCTGGTGAAAATACAACAGGTATAAATAGTGGCACGATCTCGTTATTACAAAATGGTAAAGATCCGGCACCGTCTCCCATTGTTTTACTGGCTACTAACg        <  2:928766/100‑1 (MQ=255)
                                                                                            cTGGTGAAAATACAACAGGTATAAATAGTGGCACGATCTCGTTATTACAAAATGGTAAAGATCCGGCACCGTCTCCCATTGTTTTACTGGCTACTAACgg       >  2:1440116/1‑100 (MQ=255)
                                                                                             tGGTGAAAATACAACAGGTATAAATAGTGGCACGATCTCGTTATTACAAAATGGTAAAGATCCGGCACCGTCTCCCATTGTTTTACTGGCTACTAACgga      <  1:409181/100‑1 (MQ=255)
                                                                                              ggTGAAAATACAACAGGTATAAATAGTGGCACGATCTCGTTATTACAAAATGGTAAAGATCCGGCACCGTCTCCCATTGTTTTACTGGCTACTAACggag     >  2:934660/1‑100 (MQ=255)
                                                                                              ggTGAAAATACAACAGGTATAAATAGTGGCACGATCTCGTTATTACAAAATGGTAAAGATCCGGCACCGTCTCCCATTGTTTTACTGGCTACTAACggag     >  1:3884842/1‑100 (MQ=255)
                                                                                              ggTGAAAATACAACAGGTATAAATAGTGGCACGATCTCGTTATTACAAAATGGTAAAGATCCGGCACCGTCTCCCATTGTTTTACTGGCTACTAACg        >  1:2873207/1‑97 (MQ=255)
                                                                                               gTGAAAATACAACAGGTATAAATAGTGGCACGATCTCGTTATTACAAAATGGTAAAGATCCGGCACCGTCTCCCATTGTTTTACTGGCTACTAACggagg    >  1:1251704/1‑100 (MQ=255)
                                                                                                tGAAAATACAACAGGTATAAATAGTGGCACGATCTCGTTATTACAAAATGGTAAAGATCCGGCACCGTCTCCCATTGTTTTACTGGCTACTAACggaggg   <  2:890792/100‑1 (MQ=255)
                                                                                                tGAAAATACAACAGGTATAAATAGTGGCACGATCTCGTTATTACAAAATGGTAAAGATCCGGCACCGTCTCCCATTGTTTTACTAGctac             >  2:3052508/1‑90 (MQ=255)
                                                                                                 gAAAATACAACAGGTATAAATAGTGGCACGATCTCGTTATTACAAAATGGTAAAGATCCGGCACCGTCTCCCATTGTTTTACTGGCTACTAACggaggga  <  2:247598/100‑1 (MQ=255)
                                                                                                 gAAAATACAACAGGTATAAATAGTGGCACGATCTCGTTATTACAAAATGGTAAAGATCCGGCACCGTCTCCCATTGTTTTACTGGCTACTAACggaggga  <  2:3034858/100‑1 (MQ=255)
                                                                                                 |                                                                                                   
CAAAACACCGAAAGAGATAATTGAAAGGGGTGTTGCATTGACAGATTGAATCTACACTGTTGAAATTCAGAATTTAGGTTTTGCTTTTGTTACTGGTGAAAATACAACAGGTATAAATAGTGGCACGATCTCGTTATTACAAAATGGTAAAGATCCGGCACCGTCTCCCATTGTTTTACTGGCTACTAACGGAGGGA  >  NC_000913/1470461‑1470657

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: