Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,470,558 G→A 55.9% 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,5580GA55.9% 34.2 / 126.3 104pseudogene (42/3497 nt)ydbApseudogene, autotransporter homolog; interrupted by IS2 and IS30
Reads supporting (aligned to +/- strand):  ref base G (20/26);  new base A (32/26);  total (52/52)
Fisher's exact test for biased strand distribution p-value = 3.24e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 106 total aligned reads displayed.

AAAACACCGAAAGAGATAATTGAAAGGGGTGTTGCATTGACAGATTGAATCTACACTGTTGAAATTCAGAATTTAGGTTTTGCTTTTGTTACTGGTGAAAATACAACAGGTATAAATAGTGGCACGATCTCGTTATTACAAAATGGTAAAGATCCGGCACCGTCTCCCATTGTTTTACTGGCTACTAACGGAGGG  >  NC_000913/1470462‑1470656
                                                                                                |                                                                                                  
aaaaCACCGAAAGAGATAATTGAAAGGGGTGTTGCATTGACAGATTGAATCTACACTGTTGAAATTCAGAATTTAGGTTTTGCTTTTGTTACTGGTaaaa                                                                                                 >  2:1416597/1‑100 (MQ=255)
aaaaCACCGAAAGAGATAATTGAAAGGGGTGTTGCATTGACAGATTGAATCTACACTGTTGAAATTCAGAATTTAGGTTTTGCTTTTGTTACTGGTaaaa                                                                                                 >  2:1434195/1‑100 (MQ=255)
 aaaCACCGAAAGAGATAATTGAAAGGGGTATTGCATTGACAGATTGAATCTACACTGTTGAAATTCAGAATTTAGGTTTTGCTTTTGTTACTGGTGaaaa                                                                                                <  1:1528385/100‑1 (MQ=255)
    cacCGAAAGAGATAATTGAAAGGGGTGTTGCATTGACAGATTGAATCTACACTGTTGAAATTCAGAATTTAGGTTTTGCTTTTGTTACTGGTAAAAATac                                                                                             <  1:1106110/100‑1 (MQ=255)
     acCGAAAGAGATAATTGAAAGGGGTGTTGCATTGACAGATTGAATCTACACTGTTGAAATTCAGAATTTAGGTTTTGCTTTTGTTACTGGTGAAAATaca                                                                                            <  1:1557338/100‑1 (MQ=255)
     acCGAAAGAGATAAGTGAAAGGGGTGTTGCATTGACAGATTGAATCTACACTGTTGAAATTCAGAATTTAGGTTTTGCTTTTGTTACTGGTGAAAATaca                                                                                            <  2:2571505/100‑1 (MQ=255)
      ccGAAAGAGATAATTGAAAGGGGTGTTGCATTGACAGATTGAATCTACACTGTTGAAATTCAGAATTTAGGTTTTGCTTTTGTTACTGGTGAAAATacaa                                                                                           >  1:2510149/1‑100 (MQ=255)
       cGAAAGAGATAATTGAAAGGGGTGTTGCATTGACAGATTGAATCTACACTGTTGAAATTCAGAATTTAGGTTTTGCTTTTGTTACTGGTGAAAATacaac                                                                                          <  2:1302940/100‑1 (MQ=255)
        gAAAGAGATAATTGAAAGGGGTGTTGCATTGACAGATTGAATCTACACTGTTGAAATTCAGAATTTAGGTTTTGCTTTTGTTACTGGTGAAAATacaaca                                                                                         >  2:1053889/1‑100 (MQ=255)
         aaaGAGATAATTGAAAGGGGTGTTGCATTGACAGATTGAATCTACACTGTTGAAATTCAGAATTTAGGTTTTGCTTTTGTTACTGGTAAAAATACAACAg                                                                                        >  1:256207/1‑100 (MQ=255)
          aaGAGATAATTGAAAGGGGTGTTGCAGTGACAGATTGAATCTACACTGTTGAAATTCAGAATTTAGGTTTTGCTTTTGTTACTGGTAAAAATACAACAgg                                                                                       <  1:2052275/100‑1 (MQ=255)
              gaTAATTGAAAGGGGTGTTGCATTGACAGATTGAATCTACACTGTTGAAATTCAGAATTTAGGTTTTGCTTTTGTTACTGGTGAAAATACAACAGGtata                                                                                   >  1:333710/1‑100 (MQ=255)
               aTAATTGAAAGGGGTGTTGCATTGACAGATTGAATCTACACTGTTGAAATTCAGAATTTAGGTTTTGCTTTTGTTACTGGTAAAAATACAACAGGTATaa                                                                                  <  2:831074/100‑1 (MQ=255)
                  aTTGAAAGGGGTGTTGCATTGACAGATTGAATCTACACTGTTGAAATTCAGAATTTAGGTTTTGCTTTTGTTACTGGTAAAAATACAACAGGTATAAATa                                                                               >  1:1050101/1‑100 (MQ=255)
                   ttGAAAGGGGTGTTGCATTGACAGATTGAATCTACACTGTTGAAATTCAGAATTTAGGTTTTTTTTTTTTTTCTGGTAAAAATACAACAGGTATAAATAg                                                                              >  1:119795/1‑100 (MQ=37)
                   ttGAAAGGGGTGTTGCATTGACAGATTGAATCTACACTGTTGAAATTCAGAATTTAGGTTTTGCTTTTGTTACTGGTAAAAATACAACAGGTATAAATAg                                                                              >  2:2084766/1‑100 (MQ=255)
                     gAAAGGGGTGTTGCATTGACAGATTGAATCTACACTGTTGAAATTCAGAATTTAGGTTTTGCTTTTGTTACTGGTAAAAATACAACAGGTATAAATAg                                                                              <  2:1935267/98‑1 (MQ=255)
                     gAAAGGGGTGTTGCATTGACAGATTGAATCTACACTGTTGAAATTCAGAATTTAGGTTTTGCTTTTGTTACTGGTAAAAATACAACAGGTATAAATAGTg                                                                            >  2:1722044/1‑100 (MQ=255)
                     gAAAGGGGTGTTGCATTGACAGATTGAATCTACACTGTTGAAATTCAGAATTTAGGTTTTGCTTTTGTAACTGGTAAAAATACAACTGGTATAAATAGTg                                                                            <  2:2279196/100‑1 (MQ=255)
                      aaaGGGGTGTTGCATTGACAGATTGAATCTACACTGTTGAAATTCAGAATTTAGGTTTTGCTTTTGTTACTGGTAAAAATACAACAGGTATAAATAGTgg                                                                           <  1:2265075/100‑1 (MQ=255)
                        aGGGGTGTTGCATTGACAGATTGAATCTACACTGTTGAAATTCAGAATTTAGGTTTTGCTTTTGTTACTGGTAAAAATACAACAGGTATAAATAGTGGCa                                                                         <  2:1676683/100‑1 (MQ=255)
                        aGGGGTGTTGCATTGACAGATTGAATCTACACTGTTGAAATTCAGAATTTAGGTTTTGCTTTTGTTACTGGTAAAAATACAACAGGTATAAATAGTGGCa                                                                         <  2:950206/100‑1 (MQ=255)
                          gggTGTTGCATTGACAGATTGAATCTACACTGTTGAAATTCAGAATTTAGGTTTTGCTTTTGTTACTGGTGAAAATACAACAGGTATAAATAGTGGCACg                                                                       >  1:1576417/1‑100 (MQ=255)
                           ggTGTTGCATTGACAGATTGAATCTACACTGTTGAAATTCAGAATTTAGGTTTTGCTTTTGTTACTGGTGAAAATACAACAGGTATAAATAGTGGCACGa                                                                      >  1:3120246/1‑100 (MQ=255)
                             tgtTGCATTGACAGATTGAATCTACACTGTTGAAATTCAGAATTTAGGTTTTGCTTTTGTTACTGGTAAAAATACAACAGGTATAAATAGTGGCACGAtc                                                                    <  1:1416591/100‑1 (MQ=255)
                              gtTGCATTGACAGATTGAATCTACACTGTTGAAATTCAGAATTTAGGTTTTGCTTTTGTTACTGGTGAAAATACAACAGGTATAAATAGTGGCACGAtc                                                                    <  1:1072945/99‑1 (MQ=255)
                              gtTGCATTGACAGATTGAATCTACACTGTTGAAATTCAGAATTTAGGTTTTGCTTTTGTTACTGGTGAAAATACAACAGGTATAAATAGTGGCACGAtc                                                                    >  2:3079255/1‑99 (MQ=255)
                              gtTGCATTGACAGATTGAATCTACACTGTTGAAATTCAGAATTTAGGTTTTGCTTTTGTTACTGGTAAAAATACAACAGGTATAAATAGTGGCACGAtct                                                                   >  1:1158816/1‑100 (MQ=255)
                              gtTGCATTGACAGATTGAATCTACACTGTTGAAATTCAGAATTTAGGTTTTGCTTTTGTTACTGGTAAAAATACAACAGGTATAAATAGTGGCACGAtct                                                                   >  1:2339836/1‑100 (MQ=255)
                                 acATTGACAGATTGAATCTACACTGTTGAAATTCAGAATTTAGGTTTTGCTTTTGTTACTGGTAAAAATACAACAGGTATAAATAGTGGCACGATCTCGt                                                                >  2:68504/2‑100 (MQ=255)
                                  cATTGACAGATTGAATCTACACTGTTGAAATTCAGAATTTAGGTTTTGCTTTTGTTACTGGTGAAAATACAACAGGTATAAATAGTGGCACGATCTCGtt                                                               >  1:804013/1‑100 (MQ=255)
                                  cATTGACAGATTGAATCTACACTGTTGAAATTCAGAATTTAGGTTTTGCTTTTGTTACTGGTAAAAATACAACAGGTATAAATAGTGGCACGATCTCGtt                                                               <  1:2129336/100‑1 (MQ=255)
                                  cATTGACAAATTGAATCTACACTGTTGAAATTCAGAATTTAGGTTTTGCTTTTGTTACTGGTAAAAATACAACAGGTATAAATAGTGGCACGATCTCGt                                                                <  1:907620/99‑1 (MQ=255)
                                   aTTGACAGATTGAATCTACACTGTTGAAATTCAGAATTTAGGTTTTGCTTTTGTTACTGGTAAAAATACAACAGGTATAAATAGTGGCACGATCTCGtta                                                              >  1:373900/1‑100 (MQ=255)
                                   aTTGACAGATTGAATCTACACTGTTGAAATTCAGAATTTAGGTTTTGCTTTTGTTACTGGTAAAAATACAACAGGTATAAATAGTGGCACGATCTCGtta                                                              >  2:1108273/1‑100 (MQ=255)
                                    ttGACAGATTGAATCTACACTGTTGAAATTCAGAATTTAGGTTTTGCTTTTGTTACTGGTAAAAATACAACAGGTATAAATAGTGGCACGATCTCGttat                                                             >  1:2537361/1‑100 (MQ=255)
                                    ttGACAGATTGAATCTACACTGTTGAAATTCAGAATTTAGGTTTTGCTTTTGTTACTGGTAAAAATACAACAGGTATAAATAGTGGCACGATCTCGttat                                                             <  2:1158820/100‑1 (MQ=255)
                                     tGACAGATTGAATCTACACTGTTGAAATTCAGAATTTAGGTTTTGCTTTTGTTACTGGTGAAAATACAACAGGTATAAATAGTGGCACGATCTCGttatt                                                            <  2:2096257/100‑1 (MQ=255)
                                      gACAGATTGAATCTACACTGTTGAAATTCAGAATTTAGGTTTTGCTTTTGTTACTGGTGAAAATACAACAGGTATAAATAGTGGCACGATCTCGttatta                                                           <  1:2636872/100‑1 (MQ=255)
                                      gACAGATTGAATCTACACTGTTGAAATTCAGAATTTAGGTTTTGCTTTTGTTACTGGTGAAAATACAACAGGTATAAATAGTGGCACGATCTCGttatta                                                           >  1:1723615/1‑100 (MQ=255)
                                        cAGATTGAATCTACACTGTTGAAATTCAGAATTTAGGTTTTGCTTTTGTTACTGGTAAAAATACAACAGGTATAAATAGTGGCACGATCTCGTTATTACa                                                         >  2:3251431/1‑100 (MQ=255)
                                        cAGATTGAATCTACACTGTTGAAATTCAGAATTTAGGTTTTGCTTTTGTTACTGGTAAAAATACAACAGGTATAAATAGTGGCACGATCTCGTTATTACa                                                         <  1:450729/100‑1 (MQ=255)
                                        cAGATTGAATCTACACTGTTGAAATTCAGAATTTAGGTTTTGCTTTTGTTACTGGTAAAAATACAACAGGTATAAATAGTGGCACGATCTCGTTATTACa                                                         >  1:1059815/1‑100 (MQ=255)
                                          gATTGAAGCTAGACTGTTGAAATTCAGAATTTAGGTTTTGCTTTTGTTACTGGTAAAAATACAACCGGTATAAATAGTGGCACGATCTCGTTATTACaaa                                                       >  1:1775339/1‑100 (MQ=255)
                                           aTTGAATCTACACTGTTGAAATTCAGAATTTAGGTTTTGCTTTTGTTACTGGTGAAAATACAACAGGTATAAATAGTGGCACGATCTCGTTATTACaaaa                                                      >  1:900473/1‑100 (MQ=255)
                                           aTTGAATCTACACTGTTGAAATTCAGAATTTAGGTTTTGCTTTTGTTACTGGTAAAAATACAACAGGTATAAATAGTGGCACGATCTCGTTATTACaaaa                                                      >  2:235380/1‑100 (MQ=255)
                                            ttGAATCTACACTGTTGAAATTCAGAATTTAGGTTTTGCTTTTGTTACTGGTAAAAATACAACAGGTATAAATAGTGGCACGATCTCGTTATTACAAAAt                                                     <  1:726576/100‑1 (MQ=255)
                                             tGAATCTACACTGTTGAAATTCAGAATTTAGGTTTTGCTTTTGTTACTGGTAAAAATACAACAGGTATAAATAGTGGCACGATCTCGTTATTACAAAATg                                                    >  1:1102753/1‑100 (MQ=255)
                                                 tCTACACTGTTGAAATTCAGAATTTAGGTTTTGCTTTTGTTACTGGTAAAAATACAACAGGTATAAATAGTGGCACGATCTCGTTATTACAAAATGGTaa                                                <  1:1938919/100‑1 (MQ=255)
                                                   tACACTGTTGAAATTCAGAATTTAGGTTTTGCTTTTGTTACTGGTGAAAATACAACAGGTATAAATAGTGGCACGATCTCGTTATTACAAAATGGTAAAg                                              <  2:1870977/100‑1 (MQ=255)
                                                   tACACTGTTGAAATTCAGAATTTAGGTTTTGCTTTTGTTACTGGTAAAAATACAACAGGTATAAATAGTGGCACGATCTCGTTATTACAAAATGGTaaa                                               >  1:2351143/1‑99 (MQ=255)
                                                    acacTGTTGAACATCAGAATTTAGTCTTTGCTTTTGTTACTGGCAAACATACAACTGGTATAAATGGTGGCATGATCTCTTTATTAAAACCTGGTacata                                             >  1:3074063/1‑96 (MQ=255)
                                                      agtagtGAAATTCAGAATTTAGGTTTCTCTTTTGTTACTGGTGAAAATACAACAGGTATAAATAGTGGCATGATCTCGTTATTACAAAATGGTaaa                                               <  1:3258932/91‑1 (MQ=255)
                                                      acTGTTGAAATTCAGAATTTAGGTTTTGCTTTTGTTACTGGTAAAAATACAACAGGTATAAATAGTGGCACGATCTCGTTATTACAAAATGGTAAAGATc                                           >  1:1796033/1‑100 (MQ=255)
                                                        tGTTGAAATTCAGAATTTAGGTTTTGCTTTTGTTACTGGTGAAAATACAACAGGTATAAATAGTGGCACGATCTCGTTATTACAAAATGGTAAAGATCCg                                         <  2:1520420/100‑1 (MQ=255)
                                                        tGTTGAAATTCAGAATTTAGGTTTTGCTTTTGTTACTGGTGAAAATACAACAGGTATAAATAGTGGCACGATCTCGTTATTACAAAATGGTAAAGATCCg                                         <  1:2240120/100‑1 (MQ=255)
                                                         gTTGAAATTCAGAATTTAGGTTTTTTTTTTTTTTTTTGTAAAAATACAACAGGTATAAATAGTGGCACGATCTCGTTATTACAAAATGGTAAAGATCCgg                                        >  1:1563892/1‑100 (MQ=255)
                                                         gTTGAAATTCAGAATTTAGGTTTTGCTTTTGTTACTGGTGAAAATACAACAGGTATAAAGAGTGGCACGATCTCGTTATTACAAAATGGTAAAGATCCgg                                        <  2:1323395/100‑1 (MQ=255)
                                                          ttGAAATTCAGAATTTAGGTTTTGCTTTTGTTACTGGTGAAAATACAACAGGTATAAATAGTGGCACGATCTCGTTATTACAAAATGGTAAAGATCCGGc                                       <  1:2816473/100‑1 (MQ=255)
                                                           tGAAATTCAGAATTTAGGTTTTGCTTTTGTTACTGGTAAAAATACAACAGGTATAAATAGTGGCACGATCTCGTTATTACAAAATGGTAAAGATCCGGCa                                      >  1:354094/1‑100 (MQ=255)
                                                            gAAATTCAGAATTTAGGTTTTGCTTTTGTTACTGGTGAAAATACAACAGGTATAAATAGTGGCACGATCTCGTTATTACGAAATGGTAAAGATCCGGCAc                                     <  1:306170/100‑1 (MQ=255)
                                                            gAAATTCAGAATTTAGGTTTTGCTTTTGTTACTGGTGAAAATACAACAGGTATAAATAGTGGCACGATCTCGTTATTACAAAATGGTAAAGATCCGGCAc                                     >  2:1167146/1‑100 (MQ=255)
                                                            gAAATTCAGAATTTAGGTTTTGCTTTTGTTACTGGTAAAAATACAACAGGTATAAATAGTGGCACGATCTCGTTATTACAAAATGGTAAAGATCCGGCAc                                     <  1:1619652/100‑1 (MQ=255)
                                                             aaaTTCAGAATTTAGGTTTTGCTTTTGTTACTGGTGAAAATACAACAGGTATAAATAGTGGCACGATCTCGTTATTACAAAATGGTAAAGATCCGGCa                                      <  1:2457695/98‑1 (MQ=255)
                                                              aaTTCAGAATTTAGGTTTTGCTTTTGTTACTGGTGAAAATACAACAGGTATAAATAGTGGCACGATCTCGTTATTACAAAATGGTAAAGATCCGGCACCg                                   <  2:2510168/100‑1 (MQ=255)
                                                                ttCAGAATTTAGGTTTTGCTTTTGTTACTGGTGAAAATACAACAGGTATAAATAGTGGCACGATCTCGTTATTACAAAATGGTAAAGATCCGGCACCGtc                                 >  2:2846435/1‑100 (MQ=255)
                                                                ttCAGAATTTAGGTTTTGCTTTTGTTACTGGTGAAAATACAACAGGTATAAATAGTGGCACGATCTCGTTATTACAAAATGGTAAAGATCCGGCACCGtc                                 >  1:619373/1‑100 (MQ=255)
                                                                ttCAGAATTTAGGTTTTGCTTTTGTTACTGGTGAAAATACAACAGGTATAAATAGTGGCACGATCTCGTTATTACAAAATGGTAAAGATCCGGCACCGtc                                 >  1:2949498/1‑100 (MQ=255)
                                                                 tCAGAATTTAGGTTTTGCTTTTGTTACTGGTGAAAATACAACAGGTATAAATAGTGGCACGATCTCGTTATTACAAAATGGTAAAGATCCGGCACCGtct                                >  1:2199193/1‑100 (MQ=255)
                                                                   aGAATTTAGGTTTTGCTTTTGTTACTGGTGAAAATACAACAGGTATAAATAGTGGCACGATCTCGTTATTACAAAATGGTAAAGATCCGGCACCGTCTcc                              >  2:1871902/1‑100 (MQ=255)
                                                                   aGAATTTAGGTTTTGCTTTTGTTACTGGTGAAAATACAACAGGTATAAATAGTGGCACGATCTCGTTATTACAAAATGGTAAAGATCCGGCACCGTCTcc                              >  2:1871940/1‑100 (MQ=255)
                                                                   aGAATTTAGGTTTTGCTTGTGTTACTGGTAAAAATACAACAGGTATAAATAGTGGCACGATCTCGTTATTACAAAATGGTAAAGATCCGGCACCGTCTcc                              <  2:3230159/100‑1 (MQ=255)
                                                                       tttAGGTTTTGCTTTTGTTACTGGTAAAAATACAACAGGTATAAATAGTGGCACGATCTCGTTATTACAAAATGGTAAAGATCCGGCACCGTCTCCCAtt                          >  2:1204674/1‑100 (MQ=255)
                                                                        ttAGGTTTTGCTTTTGTTACTGGTAAAAATACAACAGGTATAAATAGTGGCACGATCTCGTTATTACAAAATGGTAAAGATCCGGCACCGTCTCCCATTg                         >  2:162465/1‑100 (MQ=255)
                                                                        tgaGGTTTTGCTTTTGTTACTGGTGAAAATACAACAGGTATAAATAGTGGCACGATCTCGTTATTACAAAATGGTAAAGATCCGGCACCGTCTCCCATTg                         <  2:1825026/98‑1 (MQ=255)
                                                                         tcggTTTTGCTTTTGGTCATGGTAAAAATACAACAGGTATAATGAGTGGCACTATCTCGTTATAACAAGGTGGTAAAGATCCGGCACCGTCTCCCATTGt                        <  1:2794019/98‑1 (MQ=255)
                                                                             ttttGCTTTTGTTACTGGTAAAAATACAACAGGTATAAATAGTGGCACGATCTCGTTATTACAAAATGGTAAAGATCCGGCACCGTCTCCCATTGTTTTa                    <  2:342675/100‑1 (MQ=255)
                                                                               ttGCTTTTGTTACTGGTGAAAATACAACAGGTATAAATAGTGGCACGATCTCGTTATTACAAAATGGTAAAGATCCGGCACCGTCTCCCATTGTTTTACt                  <  1:3079229/100‑1 (MQ=255)
                                                                                tGCTTTTGTTACTGGTGAAAATACAACAGGTATAAATAGTGGCACGATCTCGTTATTACAAAATGGTAAAGATCCGGCACCGTCTCCCATTGTTTTACTg                 <  2:1425332/100‑1 (MQ=255)
                                                                                 gCTTTTGTTACTGGTGAAAATACAACAGGTATAAATAGTGGCACGATCTCGTTATTACAAAATGGTAAAGGTCCGGCACCGTCTCCCATTGTTTTACTgg                <  2:2839240/100‑1 (MQ=255)
                                                                                 gCTTTTGTTACTGGTGAAAATACAACAGGTATAAATAGTGGCACGATCTCGTTATTACAAAATGGTAAAGATCCGGCACCGTCTCCCATTGTTTTACTgg                <  2:2839936/100‑1 (MQ=255)
                                                                                  cTTTTGTTACTGGTGAAAATACAACAGGTATAAATAGTGGCACGATCTCGTTATTACAAAATGGTAAAGATCCGGCACCGTCTCCCATTGTTTTACTGGc               <  2:2259104/100‑1 (MQ=255)
                                                                                  cTTTTGTTACTGGTGAAAATACAACAGGTATAAATAGTGGCACGATCTCGTTATTACAAAATGGTAAAGATCCGGCACCGTCTCCCATTGTTTTACTGGc               >  1:3226646/1‑100 (MQ=255)
                                                                                  cTTTTGTTACTGGTGAAAATACAACAGGTATAAATAGTGGCACGATCTCGTTATTACAAAATGGTAAAGATCCGGCACCGTCTCCCATTGTTTTACTGGc               <  1:90603/100‑1 (MQ=255)
                                                                                  cTTTTGTTACTGGTAAAAATACAACAGGTATAAATAGTGGCACGATCTCGTTATTACAAAATGGTAAAGATCCGGCACCGTCTCCCATTGTTTTACTGGc               <  2:1856251/100‑1 (MQ=255)
                                                                                  cTTTTGTTACTGGTAAAAATACAACAGGTATAAATAGTGGCACGATCTCGTTATTACAAAATGGTAAAGATCCGGCACCGTCTCCCATTGTTTTACTGGc               >  2:1680760/1‑100 (MQ=255)
                                                                                   ttttGTTACTGGTAAAAATACAACAGGTATAAATAGTGGCACGATCTCGTTATTACAAAATGGTAAAGATCCGGCACCGTCTCCCATTGTTTTACTGGct              <  2:3150076/100‑1 (MQ=255)
                                                                                     ttGTTACTGGTAAAAATACAACAGGTATAAATAGTGGCACGATCTCGTTATTACAAAATGGTAAAGATCCGGCACCGTCTCCCATTGTTTTACTGGctac            <  1:1049456/100‑1 (MQ=255)
                                                                                       gTTACTGGTGAAAATACAACAGGTATAAATAGTGGCACGATCTCGTTATTACAAAATGGTAAAGATCCGGCACCGTCTCCCATTGTTTTACTGGctacta          <  2:1745176/100‑1 (MQ=255)
                                                                                       gTTACTGGTAAAAATACAACAGGTATAAATAGTGGCACGATCTCGTTATTACAAAATGGTAAAGATCCGGCACCGTCTCCCATTGTTTTACTGGctacta          >  2:96553/1‑100 (MQ=255)
                                                                                        ttACTGGTGAAAATACAACAGGTATAAATAGTGGCACGATCTCGTTATTACAAAATGGTAAAGATCCGGCACCGTCTCCCATTGTTTTACTGGCTACTaa         <  1:2962087/100‑1 (MQ=255)
                                                                                         tACTGGTAAAAATACAACAGGTATAAATAGTGGCACGATCTCGTTATTACAAAATGGTCAAGATCCGGCACCGTCTCCCATTGTTTTACTGGCTACTAAc        >  1:1042392/1‑100 (MQ=255)
                                                                                           cTGGTGAAAATACAACAGGTATAAATAGTGGCACGATCTCGTTATTACAAAATGGTAAAGATCCGGCACCGTCTCCCATTGTTTTACTGGCTACTAACgg      >  1:2117380/1‑100 (MQ=255)
                                                                                           cTGGTAAAAATACAACAGGTATAAATAGTGGCACGATCTCGTTATTACAAAATGGTAAAGATCCGGCACCGTCTCCCATTGTTTTACTGGCTACTAACgg      >  1:1389104/1‑100 (MQ=255)
                                                                                            tGGTAAAAATACAACAGGTATAAATAGTGGCACGATCTCGTTATTACAAAATGGTAAAGATCCGGCACCGTCTCCCATTGTTTTACTGGCTACTAACgga     <  1:2879021/100‑1 (MQ=255)
                                                                                            tGGTAAAAATACAACAGGTATAAATAGTGGCACGATCTCGTTATTACAAAATGGTAAAGATCCGGCACCGTCTCCCATTGTTTTACTGGCTACTAACgga     <  2:732559/100‑1 (MQ=255)
                                                                                             ggTAAAAATATAACAGGTATAAATAGTGGCACGATCTCGTTATTACAAAATGGTAAAGATCCAGCACCGTCTCCCATTGTTTTACTGGCTACAAAc        >  1:1290138/1‑96 (MQ=255)
                                                                                             ggTAAAAATACAACAGGTATAAATAGTGGCACGATCTCGTTATTACAAAATGGTAAAGATCCGGCACCGTCTCCCATTGTTTTACTGGCTACTAACggag    >  2:1971324/1‑100 (MQ=255)
                                                                                             ggTAAAAATACAACAGGTATAAATAGTGGCACGATCTCGTTATTACAAAATGGTAAAGATCCGGCACCGTCTCCCATTGTTTTACTGGCTACTAACggag    >  2:3319782/1‑100 (MQ=255)
                                                                                               tGAAAATACAACAGGTATAAATAGTGGCACGATCTCGTTATTACAAAATGGTAAAGATCCGGCACCGTCTCCCATTGTTTTACTGGCTACTAACggaggg  <  1:887765/100‑1 (MQ=255)
                                                                                                |                                                                                                  
AAAACACCGAAAGAGATAATTGAAAGGGGTGTTGCATTGACAGATTGAATCTACACTGTTGAAATTCAGAATTTAGGTTTTGCTTTTGTTACTGGTGAAAATACAACAGGTATAAATAGTGGCACGATCTCGTTATTACAAAATGGTAAAGATCCGGCACCGTCTCCCATTGTTTTACTGGCTACTAACGGAGGG  >  NC_000913/1470462‑1470656

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: