Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,060,264 A→C 16.3% intergenic (+34/‑183) typA → / → yihL GTP‑binding protein/putative DNA‑binding transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,060,2640AC16.3% 287.9 / 53.6 124intergenic (+34/‑183)typA/yihLGTP‑binding protein/putative DNA‑binding transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base A (47/56);  new base C (10/10);  total (58/66)
Fisher's exact test for biased strand distribution p-value = 8.08e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.31e-01

Only 100 of 126 total aligned reads displayed.

ATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAATTATCT  >  NC_000913/4060165‑4060363
                                                                                                   |                                                                                                    
aTTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGa                                                                                                      >  2:2420653/1‑100 (MQ=255)
   cttAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATgg                                                                                                   >  2:1754279/3‑100 (MQ=255)
   cGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATgg                                                                                                   >  2:1991059/1‑100 (MQ=255)
   cGTAAACGTCACCGGACGGAAAACGATCGTCGCCGCGCCAACGGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATgg                                                                                                   >  2:1752306/1‑100 (MQ=255)
    gTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGc                                                                                                  <  1:1128008/100‑1 (MQ=255)
    gTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGc                                                                                                  >  1:2028169/1‑100 (MQ=255)
      aaaCGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAg                                                                                                >  2:2204183/1‑100 (MQ=255)
        aCGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCATCGCTGGCAGGt                                                                                             <  1:1938572/100‑1 (MQ=255)
        aCGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCATCGCTGGCAGGt                                                                                             <  2:2593272/100‑1 (MQ=255)
        aCGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCATCGCTGGCAGGt                                                                                             <  2:2081530/100‑1 (MQ=255)
        aCGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCATCGCTGGCAGGt                                                                                             <  2:1232148/100‑1 (MQ=255)
         cGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCATCGCTGGCAGGtt                                                                                            <  2:2515597/100‑1 (MQ=255)
         cGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGtt                                                                                            >  2:1287922/1‑100 (MQ=255)
           tCACCTGACGGAAATCGATCCTCGCCGCGCCAACCGCGCACCGAAAGACGATGAATTTCGTTTTTTAGTCATAAAAAACCTGCCATCGCTGGCAGGtttt                                                                                          <  1:302899/100‑1 (MQ=38)
           tCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGtttt                                                                                          <  2:687504/100‑1 (MQ=255)
           tCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGtttt                                                                                          <  1:93145/100‑1 (MQ=255)
            cACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCATCGCTGGCAGGttttt                                                                                         <  2:1766111/100‑1 (MQ=255)
            cACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGttttt                                                                                         <  2:884978/100‑1 (MQ=255)
            cACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGttttt                                                                                         >  1:1701606/1‑100 (MQ=255)
              ccggACGGAAAACGATCTGCGCCGCGCTAACCGCGCACCGAAAGACGAGTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGtttttt                                                                                        <  1:2002082/96‑1 (MQ=255)
              ccTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCATCGCTGGCAGGtttttta                                                                                       <  2:2329888/100‑2 (MQ=255)
              ccTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGttttttt                                                                                       <  1:1436704/100‑1 (MQ=255)
              ccTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGttttttt                                                                                       >  2:785085/1‑100 (MQ=255)
              ccTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGttttttt                                                                                       >  1:769904/1‑100 (MQ=255)
               cTGACGGAAAACGATTGGCTGCGCGCCAACCGCGCATCGAATGTCGAGTAATTTCGTTTTTTAGTCATAAAATACCTGCCACCGATGGCAGGt                                                                                             <  1:1679169/93‑1 (MQ=255)
                tGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACTATTAATTTCGTTTTTTAGTCATAAAAATCCTGCCATCGCTGGCAGGttttttatg                                                                                     <  1:938164/100‑4 (MQ=38)
                tGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGttttttttt                                                                                     >  1:1474938/1‑100 (MQ=255)
                 gACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCATCGCTGGCAGGttttttatga                                                                                    <  2:279130/100‑5 (MQ=38)
                 gACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTg                                                                                    <  2:1538091/100‑1 (MQ=255)
                 gACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTg                                                                                    <  2:2148418/100‑1 (MQ=255)
                  aCGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGttttttttgcc                                                                                   <  1:1109275/100‑4 (MQ=255)
                  aCGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGc                                                                                   <  2:2570249/100‑1 (MQ=255)
                   cGGAAAACGATCGTCGCCGCGCCAACCGGGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGcc                                                                                  <  1:580201/100‑1 (MQ=255)
                   cGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCATCGCTGGCAGGttttttatgact                                                                                  >  1:642800/1‑94 (MQ=37)
                   cGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGttttttttgcct                                                                                  <  1:62517/100‑5 (MQ=255)
                   cGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGcc                                                                                  >  2:1020810/1‑100 (MQ=255)
                   cGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGc                                                                                   >  1:450113/1‑99 (MQ=255)
                        aaCGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCg                                                                             >  2:521536/1‑100 (MQ=255)
                          cGATCGTCGCCGCGCCAACCGAGCAACGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTATTTGCCTGTCGTg                                                                           >  1:848536/1‑100 (MQ=255)
                            cgatctgccgcGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAg                                                                         <  2:41937/94‑1 (MQ=255)
                            aTCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAACAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGan                                                                         <  1:695636/100‑2 (MQ=255)
                            aTCGTCGCCGCGCCAACCGCGCACAGAAAGACGATTAATTACGCTTTTTAGTCATAAAAAACCTACCAGCGGTGGCAGATTTTTTTTTGCCTgacctgcg                                                                         >  1:848877/1‑93 (MQ=255)
                              cgtcgCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAgcg                                                                       <  2:57503/100‑1 (MQ=255)
                                  gccgcGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAgagatcg                                                                   >  1:41936/1‑94 (MQ=255)
                                   ccgcGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCacac                                                                  <  1:1943547/100‑1 (MQ=255)
                                   ccgcGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCacac                                                                  <  1:2354221/100‑1 (MQ=255)
                                    cgcgTCATCCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAACACCTGCCCGCGATGGCATGTTTTTTTTTGCCTTTCGTGAGCGCACAcc                                                                 <  1:1111794/100‑1 (MQ=255)
                                    cgcgCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACAcc                                                                 >  2:2274634/1‑100 (MQ=255)
                                       gcCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCa                                                                >  1:1564036/1‑98 (MQ=255)
                                        ccAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATg                                                             >  1:431131/1‑100 (MQ=255)
                                        ccAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATg                                                             <  2:2171750/100‑1 (MQ=255)
                                           aCCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGa                                                          <  2:283063/100‑1 (MQ=255)
                                            ccGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAg                                                         <  1:1547630/100‑1 (MQ=255)
                                            ccGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAg                                                         >  1:1954601/1‑100 (MQ=255)
                                             cgcgCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAg                                                         <  2:2156202/100‑1 (MQ=255)
                                               cgcACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAATACCTGCCAGAGATATCTGGTTTTTTTTGATCTGGCGTGAGCGAAGACCAATGCGGGtgt                                                      >  1:763276/1‑98 (MQ=255)
                                                  aCCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTAtctc                                                   <  2:1372899/100‑1 (MQ=255)
                                                  aCCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTAtctc                                                   <  2:1474947/100‑1 (MQ=255)
                                                  aCCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTAtctc                                                   <  2:601259/100‑1 (MQ=255)
                                                    cGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTccc                                                 <  2:6201/100‑1 (MQ=255)
                                                     gAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCt                                                >  2:35651/1‑100 (MQ=255)
                                                      aaaGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCtt                                               <  1:1431728/100‑1 (MQ=255)
                                                      aaaGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCtt                                               <  2:1393418/100‑1 (MQ=255)
                                                         gACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCATCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTtata                                            <  2:697517/100‑1 (MQ=255)
                                                         gACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTtata                                            <  2:355189/100‑1 (MQ=255)
                                                           cGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAg                                          <  2:546108/100‑1 (MQ=255)
                                                           cGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAg                                          <  1:2279836/100‑1 (MQ=255)
                                                            gATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTTCTTTTTTTTAGCGCACACGAATTCGAGTATCTCTTTTATaat                                          >  1:1378885/1‑98 (MQ=255)
                                                             aTTAATTTCGTTTTTTAGTCATAAAACACCTGCCAGCGATGGCCGGTATTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTc                                        >  2:522876/1‑100 (MQ=255)
                                                                 aTTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATCAGTCt                                       >  2:668761/1‑97 (MQ=255)
                                                                 aTTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCtgtg                                    >  2:1443838/1‑100 (MQ=255)
                                                                 aTTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCtgtg                                    <  1:1287916/100‑1 (MQ=255)
                                                                 aTTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTATCGTGCGCGCACACCAATGCGAGTATCTCCCTTATAAGTCtgtg                                    >  1:1537051/1‑100 (MQ=255)
                                                                  tttCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTgg                                   <  2:894111/100‑1 (MQ=255)
                                                                  tttCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTgg                                   <  2:1865928/100‑1 (MQ=255)
                                                                      gTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTAc                               >  2:789113/1‑100 (MQ=255)
                                                                         ccgatctCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTc                            <  2:1897679/94‑1 (MQ=255)
                                                                          tttAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCa                           <  1:1311451/100‑1 (MQ=255)
                                                                          tttAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCa                           >  2:109374/1‑100 (MQ=255)
                                                                          acaggcaaaaaaaaaCCTGCCATCGCTGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCa                           >  2:1649761/10‑100 (MQ=25)
                                                                          acaggcaaaaaaaaaCCTGCCATCGCTGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCa                           >  2:1649712/10‑100 (MQ=25)
                                                                              gcaaaaaaaaaCCTGCCATCGCTGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATag                       >  1:1953298/6‑100 (MQ=38)
                                                                               caaaaaaaaaCCTGCCATCGCTGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATaga                      >  1:1889580/5‑100 (MQ=38)
                                                                                cATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCagatcgg                     >  1:1897669/1‑94 (MQ=255)
                                                                                cATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATagag                     >  2:1682422/1‑100 (MQ=255)
                                                                                cATAAAAAACCTGCCAGCGATGGCAGGTTGTTTTTTGCATGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGACTGTGGTTACGTCAATagag                     >  1:1321226/1‑100 (MQ=255)
                                                                                aaaaaaaaaCCTGCCATCGCTGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATagag                     >  2:225157/4‑100 (MQ=38)
                                                                                 aaaaaaaaCCTGCCATCGCTGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCCCCGTTAGAAGTCTGTGGTTACGTCAATagaga                    >  1:2201145/3‑100 (MQ=37)
                                                                                  tAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATagagag                   >  2:2252423/1‑100 (MQ=255)
                                                                                  aaaaaaaCCTGCCATCGCTGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATagagag                   >  2:2009486/2‑100 (MQ=255)
                                                                                   aaaaaaCCTGCCATCGCTGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGc                  >  1:1146001/1‑100 (MQ=255)
                                                                                    aaaaaCCTGCCAGCGATGGCAGGTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCtt                >  2:1379248/1‑100 (MQ=255)
                                                                                     aaaaCCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGTGTATCTCCCTTCTAAGTCTGTGGTTACGTCAATAGAGAGCt                 >  1:2591942/1‑99 (MQ=255)
                                                                                         ccTGCCAGCGATGGCGGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTc            <  2:2556429/100‑1 (MQ=255)
                                                                                         ccTGCCAGCGATGGCAGGTTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGt             >  2:390955/1‑100 (MQ=255)
                                                                                          cTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATTCGAGTATCTCCCTTCTCAGTCTGTGGTTACGTCAATAGAGAGCTTTGTc            <  1:1203531/99‑1 (MQ=255)
                                                                                          cTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTca           <  1:1443108/100‑1 (MQ=255)
                                                                                               aGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAAtt      <  2:769907/100‑1 (MQ=255)
                                                                                                gCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAATTa     <  2:88330/100‑1 (MQ=255)
                                                                                                   aTGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAATTATCt  <  1:1991048/100‑1 (MQ=255)
                                                                                                   |                                                                                                    
ATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAATTATCT  >  NC_000913/4060165‑4060363

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: