Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,534,711 T→G intergenic (‑36/‑80) yjhZ ← / → yjhR pseudogene, rimK paralog, C‑terminal fragment/pseudogene, helicase family; putative frameshift suppressor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,534,7110TG98.2% 766.3 / ‑2.9 219intergenic (‑36/‑80)yjhZ/yjhRpseudogene, rimK paralog, C‑terminal fragment/pseudogene, helicase family; putative frameshift suppressor
Reads supporting (aligned to +/- strand):  ref base T (2/2);  new base G (107/108);  total (109/110)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 221 total aligned reads displayed.

CCGAGCCAAAACACATCGCATTGCTGACGTAAATAGATATCACCAATCAAAACCTGCGTCATTTGTTGACAATGGCCAATTGAATATCGCGAGCAAGTGTAAGAACTCCCAGCAAGACCCTGTAAATAATTCTGTGTAACTGCCCACGATTAAAGGTGGCCGCTCAGGCGGTCACCGTATTCGATAATAAATG  >  NC_000913/4534614‑4534806
                                                                                                 |                                                                                               
ccGAGCCAAAACACATCGCATTGCTGACGTAAATAGATATCACCAATCAAAACCTGCGTCATTTGTTGACAATGGCCAATTGAATATCGCGAGCAAgggt                                                                                               <  1:5346310/100‑1 (MQ=255)
 cGAGCCAAAACACATCGCATTGCTGACGTAAATAGATATCACCAATCAAAACCTGCGTCATTTGTTGACAATGGCCAATTGAATATCGCGAGCAAGGGTa                                                                                              >  2:27927/1‑100 (MQ=255)
 cGAGCCAAAACACATCGCATTGCTGACGTAAATAGATATCACCAATCAAAACCTGCGTCATTTGTTGACAATGGCCAATTGAATATCGCGAGCAAGGGTa                                                                                              >  2:1241148/1‑100 (MQ=255)
  gAGCCAAAACACATCGCATTGCTGACGTAAATAGATATCACCAATCAAAACCTGCGTCATTTGTTGACAATGGCCAATTGAATATCGCGAGCAAGGGTaa                                                                                             <  1:1025356/100‑1 (MQ=255)
  gAGCCAAAACACATCGCATTGCTGACGTAAATAGATATCACCAATCAAAACCTGCGTCATTTGTTGACAATGGCCAATTGAATATCGCGAGCAAGGGTaa                                                                                             <  1:4728590/100‑1 (MQ=255)
  gAGCCAAAACACATCGCATTGCTGACGTAAATAGATATCACCAATCAAAACCTGCGTCATTTGTTGACAATGGCCAATTGAATATCGCGAGCAAGGGTaa                                                                                             <  1:2882056/100‑1 (MQ=255)
    gCCAAAACACATCGCATTGCTGACGTAAATAGATATCACCAATCAAAACCTGCGTCATTTGTTGACAATGGCCAATTGAATATCGCGAGCAAGGGTAAGa                                                                                           >  1:4810020/1‑100 (MQ=255)
      cAAAACACATCGCATTGCTGACGTAAATAGATATCACCAATCAAAACCTGCGTCATTTGTTGACAATGGCCAATTGAATATCGCGAGCAAGGGTAAGAAc                                                                                         <  1:4252740/100‑1 (MQ=255)
      cAAAACACATCGCATTGCTGACGTAAATAGATATCACCAATCAAAACCTGCGTCATTTGTTGACAATGGCCAATTGAATATCGCGAGCAAGGGTAAGAAc                                                                                         <  1:2116152/100‑1 (MQ=255)
      cAAAACACATCGCATTGCTGACGTAAATAGATATCACCAATCAAAACCTGCGTCATTTGTTGACAATGGCCAATTGAATATCGCGAGCAAGGGTAAGAAc                                                                                         >  2:1615838/1‑100 (MQ=255)
       aaaaCACATCGCATTGCTGACGTAAATAGATATCACCAATCAAAACCTGCGTCATTTGTTGACAATGGCCAATTGAATATCGCGAGCAAGGGTAAGaaag                                                                                        >  2:4772961/1‑98 (MQ=255)
        aaaCACATCGCATTGCTGACGTAAATAGATATCACCAATCAAAACCTGCGTCATTTGTTGACAATGGCCAATTGAATATCGCGAGCAAGGGTAAGAACTc                                                                                       <  1:2659784/100‑1 (MQ=255)
        aaaCACATCGCATTGCTGACGTAAATAGATATCACCAATCAAAACCGGCGTCATTTGTTGACAATGGCCAATTGAATATCGCGAGCAAGGGTAAGAACTc                                                                                       <  2:3843853/100‑1 (MQ=255)
          acacATCGCATTGCTGACGTAAATAGATATCACCAATCAAAACCTGCGTCATTTGTTGACAATGGCCAATTGAATATCGCGAGCAAGGGTAAGAACTccc                                                                                     <  1:1746791/100‑1 (MQ=255)
           cacaTCGCATTGCTGACGTAAATAGATATCACCAATCAAAACCTGCGTCATTTGTTGACAATGGCCAATTGAATATCGCGAGCAAGGGTAAGAACTCCCa                                                                                    >  1:40946/1‑100 (MQ=255)
           cacaTCGCATTGCTGACGTAAATAGATATCACCAATCAAAACCTGCGTCATTTGTTGACAATGGCCAATTGAATATCGCGAGCAAGGGTAAGAACTCCCa                                                                                    >  2:4514753/1‑100 (MQ=255)
           cacaTCGCATTGCTGACGTAAATAGATATCACCAATCAAAACCTGCGTCATTTGTTGACAATGGCCAATTGAATATCGCGAGCAAGGGTAAGAACTCCCa                                                                                    <  1:1703356/100‑1 (MQ=255)
           cacaTCGCATTGCTGACGTAAATAGATATCACCAATCAAAACCTGCGTCATTTGTTGACAATGGCCAATTGAATATCGCGAGCAAGGGTAAGAACTCCCa                                                                                    >  1:1306855/1‑100 (MQ=255)
            acaTCGCATTGCTGACGTAAATAGATATCACCAATCAAAACCTGCGTCATTTGTTGACAATGGCCAATTGAATATCGCGAGCAAGGGTAAGAACTCCCAg                                                                                   <  2:4722144/100‑1 (MQ=255)
              aTCGCATTGCTGACGTAAATAGATATCACCAATCAAAACCTGCGTCATTTGTTGACAATGGCCAATTGAATATCGCGAGCAAGGGTAAGAACTCCCAGCa                                                                                 <  2:4245388/100‑1 (MQ=255)
                  cATTGCTGACGTAAATAGATATCACCAATCAAAACCTGCGTCATTTGTTGACAATGGCCAATTGAATATCGCGAGCAAGGGTAAGAACTCCCAGCAAGAc                                                                             >  2:2030132/1‑100 (MQ=255)
                     tGCTGACGTAAATAGATATCACCAATCAAAACCTGCGTCATTTGTTGACAATGGCCAATTGAATATCGCGAGCAAGGGTAAGAACTCCCAGCAAGACCCt                                                                          <  2:4384375/100‑1 (MQ=255)
                       cTGACGTAAATAGATATCACCAATCAAAACCTGCGTCATTTGTTGACAATGGCCAATTGAATATCGCGAGCAAGGGTAAGAACTCCCAGCAAGACCCTGt                                                                        >  1:1221886/1‑100 (MQ=255)
                        tGACGTAAATAGATATCACCAATCAAAACCTGCGTCATTTGTTGACAATGGCCAATTGAATATCGCGAGCAAGGGTAAGAACTCCCAGCAAGACCCTGTa                                                                       >  2:2508515/1‑100 (MQ=255)
                        tGACGTAAATAGATATCACCAATCAAAACCTGCGTCATCTGTTGACAATGGCCAATTGAATATCGCGAGCAAGGGTAAGAACTCCCATCAAGACCCTGTa                                                                       <  2:5295143/100‑1 (MQ=255)
                         gACGTAAATAGATATCACCAATCAAAACCTGCGTCATTTGTTGACAATGGCCAATTGAATATCGCGAGCAAGGGTAAGAACTCCCAGCAAGACCCTGTaa                                                                      >  2:3133352/1‑100 (MQ=255)
                          aCGTAAATAGATATCACCAATCAAAACCTGCGTCATTTGTTGACAATGGCCAATTGAATATCGCGAGCAAGGGTAAGAACTCCCAGCAAGACCCTGTaaa                                                                     <  1:852534/100‑1 (MQ=255)
                           cGTAAATAGATATCACCAATCAAAACCTGCGTCATTTGTTGACAATGGCCAATTGAATATCGCGAGCAAGGGTAAGAACTCCCAGCAAGACCCTGTAaat                                                                    <  1:3889051/100‑1 (MQ=255)
                           cGTAAATAGATATCACCAATCAAAACCTGCGTCATTTGTTGACAATGGCCAATTGAATATCGCGAGCAAGGGTAAGAACTCCCAGCAAGACCCTGTAaat                                                                    >  1:4781892/1‑100 (MQ=255)
                           cGTAAATAGATATCACCAATCAAAACCTGCGTCATTTGTTGACAATGGCCAATTGAATATCGCGAGCAAGGGTAAGAACTCCCAGCAAGACCCTGTAaat                                                                    >  1:4497351/1‑100 (MQ=255)
                            gTAAATAGATATCACCAATCAAAACCTGCGTCATTTGTTGACAATGGCCATTTGAATATCGCGAGCAAGGGTAAGTACTCCCAGCAAGACCCTGCAcata                                                                   >  1:5196091/1‑100 (MQ=255)
                            gTAAATAGATATCACCAATCAAAACCTGCGTCATTTGTTGACAATGGCCAATTGAATATCGCGAGCAAGGGTAAGAACTCCCAGCAAGACCCTGTAaata                                                                   >  1:2398198/1‑100 (MQ=255)
                            gTAAATAGATATCACCAATCAAAACCTGCGTCATTTGTTGACAATGGCCAATTGAATATCGCGAGCAAGGGTAAGAACTCCCAGCAAGACCCTGTAaata                                                                   >  2:1685/1‑100 (MQ=255)
                            gTAAATAGATATCACCAATCAAAACCTGCGTCATTTGTTGACAATGGCCAATTGAATATCGCGAGCAAGGGTAAGAACTCCCAGCAAGACCCTGTAaata                                                                   <  2:5221247/100‑1 (MQ=255)
                             tAAATAGATATCACCAATCAAAACCTGCGTCATTTGTTGACAATGGCCAATTGAATATCGCGAGCAAGGGTAAGAACTCCCAGCAAGACCCTGTAaataa                                                                  >  2:4920907/1‑100 (MQ=255)
                             tAAATAGATATCACCAATCAAAACCTGCGTCATTTGTTGACAATGGCCAATTGAATATCGCGAGCAAGGGTAAGAACTCCCAGCAAGACCCTGTAaataa                                                                  >  1:199288/1‑100 (MQ=255)
                              aaaTAGATATCACCAATCAAAACCTGCGTCATTTGTTGACAATGGCCAATTGAATATCGCGAGCAAGGGTAAGAACTCCCAGCAAGACCCTGTAaataan                                                                 <  1:4416786/100‑2 (MQ=255)
                               aaTAGATATCACCAATCAAAACCTGCGTCATTTGTTGACAATGGCCAATTGAATATCGCGAGCAAGGGTAAGAACTCCCAGCAAGACCCTGTAAATAAtt                                                                >  1:52291/1‑100 (MQ=255)
                                aTAGATATCACCAATCAAAACCTGCGTCATTTGTTGACAATGGCCAATTGAATATCGCGAGCAAGGGTAAGAACTCCCAGCAAGACCCTGTAAATAATTc                                                               >  1:1422545/1‑100 (MQ=255)
                                aTAGATATCACCAATCAAAACCTGCGTCATTTGTTGACAATGGCCAATTGAATATCGCGAGCAAGGGTAAGAACTCCCAGCAAGACCCTGTAAATAATTc                                                               <  1:2727385/100‑1 (MQ=255)
                                 tAGATATCACCAATCAAAACCTGCGTCATTTGTTGACAATGGCCAATTGAATATCGCGAGCAAGGGTAAGAACTCCCAGCAAGACCCTGTAAATAATTCt                                                              >  2:2645711/1‑100 (MQ=255)
                                  aGATATCACCAATCAAAACCTGCGTCATTTGTTGACAATGGCCAATTGAATATCGCGAGCAAGGGTAAGAACTCCCAGCAAGACCCTGTAAATAATTCtg                                                             <  1:3996630/100‑1 (MQ=255)
                                  aGATATCACCAATCAAAACCTGCGTCATTTGTTGACAATGGCCAATTGAATATCGCGAGCAAGGGTAAGAACTCCCAGCAAGACCCTGTAAATAATTCtg                                                             <  2:1174425/100‑1 (MQ=255)
                                  aGATATCACCAATCAAAACCTGCGTCATTTGTTGACAATGGCCAATTGAATACCGCGAGCAAGGGTAAGAACTCCCAGCAAGACCCTGTAAATAATTCtg                                                             <  1:1618290/100‑1 (MQ=255)
                                   gATATCACCAATCAAAACCTGCGTCATTTGTTGACAATGGCCAATTGAATATCGCGAGCAAGGGTAAGAACTCCCAGCAAGACCCTGTAAATAATTCtgt                                                            <  2:3680020/100‑1 (MQ=255)
                                    atatCACCAATCAAAACCTGCGTCATTTGTTGACAATGGCCAATTGAATATCGCGAGCAAGGGTAAGAACTCCCAGCAAGACCCTGTAAATAATTCtgtg                                                           <  2:2374166/100‑1 (MQ=255)
                                    atatCACCAATCAAAACCTGCGTCATTTGTTGACAATGGCCAATTGAATATCGCGAGCAAGGGTAAGAACTCCCAGCAAGACCCTGTAAATAATTCtgtg                                                           <  2:197375/100‑1 (MQ=255)
                                    atatCACCAATCAAAACCTGCGTCATTTGTTGACAATGGCCAATTGAATATCGCGAGCAAGGGTAAGAACTCCCAGCAAGACCCTGTAAATAATTCtgtg                                                           >  2:2908913/1‑100 (MQ=255)
                                    atatCACCAATCAAAACCTGCGTCATTTGTTGACAATGGCCAATTGAATATCGCGAGCAAGGGTAAGAACTCCCAGCAAGACCCTGTAAATAATTCtgtg                                                           <  1:799930/100‑1 (MQ=255)
                                     tatCACCAATCAAAACCTGCGTCATTTGTTGACAATGGCCAATTGAATATCGCGAGCAAGGGTAAGAACTCCCAGCAAGACCCTGTAAATAATTCTgtgt                                                          <  1:2439561/100‑1 (MQ=255)
                                     tatCACCAATCAAAACCTGCGTCATTTGTTGACAATGGCCAATTGAATATCGCGAGCAAGGGTAAGAACTCCCAGCAAGACCCTGTAAATAATTCTgtgt                                                          >  1:3677713/1‑100 (MQ=255)
                                      atCACCAATCAAAACCTGCGTCATTTGTTGACAATGGCCAATTGAATATCGCGAGCAAGGGTAAGAACTCCCAGCAAGACCCTGTAAATAATTCTGTGTa                                                         <  1:676244/100‑1 (MQ=255)
                                      atCACCAATCAAAACCTGCGTCATTTGTTGACAATGGCCAATTGAATATCGCGAGCAAGGGTAAGAACTCCCAGCAAGACCCTGTAAATAATTCTGTGTa                                                         >  1:3307636/1‑100 (MQ=255)
                                        cACCAATCAAAACCTGCGTCATTTGTTGACAATGGCCAATTGAATATCGCGAGCAAGTGTAAGAACTCCCAGCAAGACCCTGTAAATAATTCTGTGTAAc                                                       <  1:2528890/100‑1 (MQ=255)
                                        cACCAATCAAAACCTGCGTCATTTGTTGACAATGGCCAATTGAATATCGCGAGCAAGGGTAAGAACTCCCAGCAAGACCCTGTAAATAATTCTGTGTAAc                                                       >  2:3273516/1‑100 (MQ=255)
                                        cACCAATCAAAACCTGCGTCATTTGTTGACAATGGCCAATTGAATATCGCGAGCAAGGGTAAGAACTCCCAGCAAGACCCTGTAAATAATTCTGTGTAAc                                                       >  2:4127346/1‑100 (MQ=255)
                                        cACCAATCAAAACCTGCGTCATTTGTTGACAATGGCCAATTGAATATCGCGAGCAAGGGTAAGAACTCCCAGCAAGACCCTGTAAATAATTCTGTGTAAc                                                       >  2:1548207/1‑100 (MQ=255)
                                        cACCAATCAAAACCTGCGTCATTTGTTGACAATGGCCAATTGAATATCGCGAGCAAGGGTAAGAACTCCCAGCAAGACCCTGTAAATAATTCTGTGTAAc                                                       <  1:882649/100‑1 (MQ=255)
                                         aCCAATCAAAACCTGCGTCATTTGTTGACAATGGCCAATTGAATATCGCGAGCAAGGGTAAGAACTCCCAGCAAGACCCTGTAAATAATTCTGTGTAACt                                                      >  2:2442738/1‑100 (MQ=255)
                                          ccAATCAAAACCTGCGTCATTTGTTGACAATGGCCAATTGAATATCGCGAGCCAGGGGAAGAACTCCCAGCAAGACCCTGTACATAATTCTGTGTAACTg                                                     >  2:4776206/1‑100 (MQ=255)
                                          ccAATCAAAACCTGCGTCATTTGTTGACAATGGCCAATTGAATATCGCGAGCAAGGGTAAGAACTCCCAGCAAGACCCTGTAAATAATTCTGTGTAACTg                                                     >  1:5450630/1‑100 (MQ=255)
                                           cAATCAAAACCTGCGTCATTTGTTGACAATGGCCAATTGAATATCGCGAGCAAGGGTAAGAACTCCCAGCAAGACCCTGTAAATAATTCTGTGTAACTGc                                                    >  1:3007747/1‑100 (MQ=255)
                                            aaTCAAAACCTGCGTCATTTGTTGACAATGGCCAATTGAATATCGCGAGCAAGGGTAAGAACTCCCAGCAAGACCCTGTAAATAATTCTGTGTAACTGcc                                                   <  2:2122225/100‑1 (MQ=255)
                                             aTCAAAAGCTGCGTCATTTGTTGACAATGGCGAATTGAGTATCGCGGGCAGGGGTAAGCACTCCCAGCAAGACCCGGTGTATAATTCTGTGTAACTGctc                                                  <  2:656041/100‑3 (MQ=255)
                                               cAAAACCTGCGTCATTTGTTGACAATGGCCAATTGAATATCGCGAGCAAGGGTAAGAACTCCCAGCAAGACCCTGTAAATAATTCTGTGTAACTGCCCAc                                                >  1:5503894/1‑100 (MQ=255)
                                                aaaaCCTGCGTCATTTGTTGACAATGGCCAATTGAATATCGCGAGCAAGGGTAAGAACTCCCAGCAAGACCCTGTAAATAATTCTGTGTAACTGCCCACg                                               >  2:685061/1‑100 (MQ=255)
                                                aaaaCCTGCGTCATTTGTTGACAATGGCCAATTGAATATCGCGAGCAAGGGTAAGAACTCCCAGCAAGACCCTGTAAATAATTCTGTGTAACTGCCCACg                                               >  1:3170862/1‑100 (MQ=255)
                                                aaaaCCTGCGTCATTTGTTGACAATGGCCAATTGAATATCGCGAGCAAGGGTAAGAACTCCCAGCAAGACCCTGTAAATAATTCTGTGTAACTGCCCACg                                               >  1:4660072/1‑100 (MQ=255)
                                                aaaaCCTGCGTCATTTGTTGACAATGGCCAATTGAATATCGCGAGCAAGGGTAAGAACTCCCAGCAAGACCCTGTAAATAATTCTGTGTAACTGCCCACg                                               <  2:146520/100‑1 (MQ=255)
                                                 aaaCCTGCGTCATTTGTTGACAATGGCCAATTGAATATCGCGAGCAAGGGTAAGAACTCCCAGCAAGACCCTGTAAATAATTCTGTGTAACTGCCCACGa                                              >  1:2796796/1‑100 (MQ=255)
                                                     cTGCGTCATTTGTTGACACTGTCCAATTGAATATCGCGAGCAAGGGTAAGAACTCCCAGCAAGACCCAGTAACTAATTCTGTGTAACTGCCCACGATTaa                                          >  2:553267/1‑100 (MQ=255)
                                                      tGCGTCATTTGTTGACAATGGCCAATTGAATATCGCGAGCAAGGGTTAGAACTCCCAGCAAGACCATGTAAATAATTCTGTGTAACTGCCCACGATTaaa                                         <  2:5081064/100‑1 (MQ=255)
                                                           cATTTGTTGACAATGGCCAATTGAATATCGCGAGCAAGGGTAAGAACTCCCAGCAAGACCCTGTAAATAATTCTGTGTAACTGCCCACGATTAAAGGTgg                                    >  2:4829065/1‑100 (MQ=255)
                                                              ttgttgACAATGGCCAATTGAATATCGCGAGCAAGGGTAAGAACTCCCAGCAAGACCCTGTAAATAATTCTGTGTAACTGCCCACGATTAAAGGTGGCCg                                 >  2:747837/1‑100 (MQ=255)
                                                                gttgACAATGGCCAATTGAATATGGCGAGCAAGGGTAAGAACTCCCAGCAAGACGCTGTAAATAATTCTATGTAACTGCCCACGCTTAAAGGTGGCCGCt                               <  1:4844648/100‑1 (MQ=255)
                                                                 ttgACAATGGCCAATTGAATATCGCGAGCAAGGGTAAGAACTCCCAGCAAGACCCTGTAAATAATTCTGTGTAACTGCCCACGATTAAAGGTGGCCGCTc                              >  2:5147843/1‑100 (MQ=255)
                                                                  tgACAATGGGCAATTGAATATCGCGAGCAAGGGGAAGAACTCCCAGCAAGACCCTGAAAATAATACGGTGTAACTGACCACGATTAAAGGTGGCAGCTCa                             >  2:1165083/1‑100 (MQ=37)
                                                                  tgACAATGGCCAATTGAATATCGCGAGCAAGGGTAAGAACTCCCAGCAAGACCCTGTAAATAATTCTGTGTAACTGCCCACGATTAAAGGTGGCCGCTCa                             <  1:1645349/100‑1 (MQ=255)
                                                                  tgACAATGGCCAATTGAATATCGCGAGCAAGGGTAAGAACTCCCAGCAAGACCCTGTAAATAATTCTGTGTAACTGCCCACGATTAAAGGTGGCCGCTCa                             <  2:1289571/100‑1 (MQ=255)
                                                                       aTGGCCAATTGAATATCGCGAGCAAGGGTAAGAACTCCCAGCAAGACCCTGTAAATAATTCTGTGTAACTGCCCACGATTAAAGGTGGCCGCTCAGGCgg                        >  2:2139419/1‑100 (MQ=255)
                                                                       aTGGCCAATTGAATATCGCGAGCAAGGGTAAGAACTCCCAGCAAGACCCTGTAAATAATTCTGTGTAACTGCCCACGATTAAAGGTGGCCGCTCAGGCgg                        <  1:4172215/100‑1 (MQ=255)
                                                                        tGGCCAATTGAATATCGCGAGCAAGGGTAAGAACGCCCAGCAAGACCCTGTAAATAATTCTGTGTAACTGCCCACGATTAAAGGTGGCCGCTCAGGCGGt                       <  1:2376965/100‑1 (MQ=255)
                                                                         ggCCAATTGAATATCGCGAGCAAGGGTAAGAACTCCCAGCAAGACCCTGTAAATAATTCTGTGTAACTGCCCACGATTAAAGGTGGCCGCTCAGGCGGTc                      <  2:3511578/100‑1 (MQ=255)
                                                                         ggCCAATTGAATATCGCGAGCAAGGGTAAGAACTCCCAGCAAGACCCTGTAAATAATTCTGTGTAACTGCCCACGATTAAAGGTGGCCGCTCAGGCGGTc                      >  1:4995155/1‑100 (MQ=255)
                                                                         ggCCAATTGAATATCGCGAGCAAGGGTAAGAACTCCCAGCAAGACCCTGTAAATAATTCTGTGTAACTGCCCACGATTAAAGGTGGCCGCTCAGGCGGTc                      >  1:3617244/1‑100 (MQ=255)
                                                                          gCCAATTGAATATCGCGAGCAAGGGTAAGAACTCCCAGCAAGACCCTGTAAATAATTCTGTGTAACTGCCCACGATTAAAGGTGGCCGCTCAGGCGGTCa                     >  1:4403781/1‑100 (MQ=255)
                                                                           ccAATTGAATATCGCGAGCAAGGGTAAGAACGCCCAGCAAGACCCTGTAAATAATTCTGTGTAACTGCCCACGATTAAAGGTGGCCGCTCAGGCGGTCAc                    <  1:3120166/100‑1 (MQ=255)
                                                                             aaTTGAATATCGCGAGCAAGGGTAAGAACTCCCAGCAAGACCCTGTAAATAATTCTGTGTAACTGCCCACGATTAAAGGTGGCCGCTCAGGCGGTCACCg                  <  2:3367145/100‑1 (MQ=255)
                                                                              aTTGAATATCGCGAGCAAGGGTAAGAACTCCCAGCAAGACCCTGTAAATAATTCTGTGTAACTGCCCACGATTAAAGGTGGCCGCTCAGGCGGTCACCGt                 <  2:3928374/100‑1 (MQ=255)
                                                                              aTTGAATATCGCGAGCAAGGGTAAGAACTCCCAGCAAGACCCTGTAAATAATTCTGTGTAACTGCCCACGATTAAAGGTGGCCGCTCAGGCGGTCACCGt                 >  1:3525080/1‑100 (MQ=255)
                                                                                   atatCGCGAGCAAGGGTAAGAACTCCCAGCAAGACCCTGTAAATAATTCTGTGTAACTGCCCACGATTAAAGGTGGCCGCTCAGGCGGTCACCGTATTCg            <  1:338317/100‑1 (MQ=255)
                                                                                    tatCGCGAGCAAGGGTAAGAACTCCCAGCAAGACCCTGTAAATAATTCTGTGTAACTGCCCACGATTAAAGGTGGCCGCTCAGGCGGTCACCGTATTCGa           >  1:1053434/1‑100 (MQ=255)
                                                                                     atCGCGAGCAAGGGTAAGAACTCCCAGCAAGACCCTGTAAATAATTCTGTGTAACTGCCCACGATTAAAGGTGGCCGCTCAGGCGGTCACCGTATTCGat          >  1:95623/1‑100 (MQ=255)
                                                                                     atCGCGAGCAAGGGTAAGAACTCCCAGCAAGACCCTGTAAATAATTCTGTGTAACTGCCCACGATTAAAGGTGGCCGCTCAGGCGGTCACCGTATTCGat          >  2:4066620/1‑100 (MQ=255)
                                                                                     atCGCGAGCAAGGGTAAGAACGCCCAGCAAGACCCTGTAAATAACTCTGTGTAACTGCCCACGATTAAAGGTGGCCGCTCAGGCGGTCACCGTATTCGat          <  2:2541217/100‑1 (MQ=255)
                                                                                      tCGCGAGCAAGGGTAAGAACTCCCAGCAAGACCCTGTAAATAATTCTGTGTAACTGCCCACGATTAAAGGTGGCCGCTCAGGCGGTCACCGTATTCGata         >  1:3840000/1‑100 (MQ=255)
                                                                                       cgcgAGCAAGGGTAAGAACTCCCAGCAAGACCCTGTAAATAATTCTGTGTAACTGCCCACGATTAAAGGTGGCCGCTCAGGCGGTCACCGTATTCGataa        >  1:322939/1‑100 (MQ=255)
                                                                                         cgAGCAAGGGTAAGGACTCCCAGCAAGACCCTGTAAATAATTCTGTGTAACTGCCCACGATTAAAGGTGGCCGCTCAGGCGGTCACCGTATTCGataata      >  2:2228298/1‑100 (MQ=255)
                                                                                         cgAGCAAGGGTAAGAACTCCCAGCAAGACCCTGTAAATAATTCTGTGTAACTGCCCACGATTAAAGGTGGCCGCTCAGGCGGTCACCGTATTCGataata      <  2:5377537/100‑1 (MQ=255)
                                                                                            gCAAGTGTAAGAACTCCCAGCAAGACCCTGTAAATAATTCTGTGTAACTGCCCACGATTAAAGGTGGCCGCTCAGGCGGTCACCGTATTCGATAATAAAt   <  2:3302806/100‑1 (MQ=255)
                                                                                                 |                                                                                               
CCGAGCCAAAACACATCGCATTGCTGACGTAAATAGATATCACCAATCAAAACCTGCGTCATTTGTTGACAATGGCCAATTGAATATCGCGAGCAAGTGTAAGAACTCCCAGCAAGACCCTGTAAATAATTCTGTGTAACTGCCCACGATTAAAGGTGGCCGCTCAGGCGGTCACCGTATTCGATAATAAATG  >  NC_000913/4534614‑4534806

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