Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,812,482 T→C F198L (TTC→CTC)  katE → catalase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,812,4820TC99.8% 2442.4 / NA 613F198L (TTC→CTC) katEcatalase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  major base C (295/316);  minor base G (1/0);  total (296/317)
Fisher's exact test for biased strand distribution p-value = 4.84e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 623 total aligned reads displayed.

TCTACCGTTCAGGGTGGTGCTGGCTCTGCTGATACCGTGCGTGATATCCGTGGCTTTGCCACCAAGTTCTATACCGAAGAGGGTATTTTTGACTTCGTTGGCAATAACACGCCAATCTTCTTTATCCAGGATGCGCATAAATTCCCCGATTTTGTTCATGCGGTAAAACCAGAACCGCACTGGGCA  >  NC_000913/1812389‑1812574
                                                                                             |                                                                                            
tctACCGTTCAGGGTGGTGCTGGCTCTGCTGATACCGTGCGTGATATCCGTGGCTTTGCCACCAAGTTCTATACCGAAGAGGGTATTTTTGACCTc                                                                                            <  2:10832048/96‑1 (MQ=255)
tctACCGTTCAGGGTGGTGCTGGCTCTGCTGATACCGTGCGTGATATCCGTGGCTTTGCCACCAAGTTCTATACCGAAGAGGGTATTTTTGACCTc                                                                                            >  1:660532/1‑96 (MQ=255)
 ctACCGTTCAGGGTGGTGCTGGCTCTGCTGATACCGTGCGTGATATCCGTGGCTTTGCCACCAAGTTCTATACCGAAGAGGGTATTTTTGACCTCg                                                                                           <  1:10209295/96‑1 (MQ=255)
 ctACCGTTCAGGGTGGTGCTGGCTCTGCTGATACCGTGCGTGATATCCGTGGCTTTGCCACCAAGTTCTATACCGAAGAGGGTATTTTTGACCTCg                                                                                           <  1:3887084/96‑1 (MQ=255)
  tACCGTTCAGGGTGGTGCTGGCTCTGCTGATACCGTGCGTGATATCCGTGGCTTTGCCACCAAGTTCTATACCGAAGAGGGTATTTTTGACCTCGt                                                                                          <  2:3740150/96‑1 (MQ=255)
  tACCGTTCAGGGTGGTGCTGGCTCTGCTGATACCGTGCGTGATATCCGTGGCTTTGCCACCAAGTTCTATACCGAAGAGGGTATTTTTGACCTCGt                                                                                          >  1:3647472/1‑96 (MQ=255)
  tACCGTTCAGGGTGGTGCTAGCTCTGCTGATACCGTGCGTGATATCCGTGGCTTTGCCACCAAGTTCTATACCGAAGAGGGTATTTTTGACCTCGt                                                                                          >  1:7439222/1‑96 (MQ=255)
   aCCGTTCAGGGTGGTGCTGGCTCTGCTGATACCGTGCGTGATATCCGTGGCTTTGCCACCAAGTTCTATACCGAAGAGGGTATTTTTGACCTCGtt                                                                                         >  1:3471011/1‑96 (MQ=255)
    ccGTTCAGGGTGGTGCTGGCTCTGCTGATACCGTGCGTGATATCCGTGGCTTTGCCACCAAGTTCTATACCGAAGAGGGTATTTTTGACCTCGTTg                                                                                        <  2:13545632/96‑1 (MQ=255)
      gTTCAGGGTGGTGCTGGCTCTGCTGATACCGTGCGTGATTCCCGTGGCTTTGCCACCAAGTTCTATACCGAAGAGGGTATTTTTGACCTCGTTGGc                                                                                      <  2:4983693/96‑1 (MQ=255)
      gTTCAGGGTGGTGCTGGCTCTGCTGATACCGTGCGTGATATCCGTGGCTTTGCCACCAAGTTCTATACCGAAGAGGGTATTTTTGACCTCGTTgta                                                                                      <  1:11693768/96‑3 (MQ=255)
      gTTCAGGGTGGTGCTGGCTCTGCTGATACCGTGCGTGATATCCGTGGCTTTGCCACCAAGTTCTATACCGAAGAGGGTATTTTTGACCTCGTTGGc                                                                                      >  2:9840723/1‑96 (MQ=255)
      gTTCAGGGTGGTGCTGGCTCTGCTGATACCGTGCGTGATATCCGTGGCTTTGCCACCAAGTTCTATACCGAAGAGGGTATTTTTGACCTCGTTGGc                                                                                      <  2:3487361/96‑1 (MQ=255)
      gTTCAGGGTGGTGCTGGCTCTGCTGATACCGTGCGTGATATCCGTGGCTTTGCCACCAAGTTCTATACCGAAGAGGGGATTTTTGACCTCGTTGGc                                                                                      >  2:12526397/1‑96 (MQ=255)
       ttCAGGGTGGTGCTGGCTCTGCTGATACCGTGCGTGATATCCGTGGCTTTGCCACCAAGTTCTATACCGAAGAGGGTATTTTTGACCTCGTTGGCa                                                                                     <  2:9133247/96‑1 (MQ=255)
        tCAGGGTGGTGCTGGCTCTGCTGATACCGTGCGTGATATCCGTGGCTTTGCCACCAAGTTCTATACCGAAGAGGGTATTTTTGACCTCGTTGGCaa                                                                                    <  2:13287593/96‑1 (MQ=255)
         cAGGGTGGTGCTGGCTCTGCTGATACCGTGCGTGATATCCGTGGCTTTGCCACCAAGTTCTATACCGAAGAGTGTATTTTTGACCTCGTTGGCAAt                                                                                   >  1:12463402/1‑96 (MQ=255)
         cAGGGTGGTGCTGGCTCTGCTGATACCGTGCGTGATATCCGTGGCTTTGCCACCAAGTTCTATACCGAAGAGGGTATTTTTGACCTCGTTGGCAAt                                                                                   <  1:10950814/96‑1 (MQ=255)
          aGGGTGGTGCTGGCTCTGCTGATACCGTGCGTGATATCCGTGGCTTTGCCACCAAGTTCTATACCGAAGAGGGTATTTTTGACCTCGTTGGCAATa                                                                                  >  2:12999228/1‑96 (MQ=255)
           gggTGGTGCTGGCTCTGCTGATACCGTGCGTGATATCCGTGGCTTTGCCACCAAGTTCTATACCGAAGAGGGTATTTTTGACCTCGTTGGCAATaa                                                                                 >  1:13998070/1‑96 (MQ=255)
           gggTGGTGCTGGCTCTGCTGATACCGTGCGTGATATCCGTGGCTTTGCCACCAAGTTCTATACCGAAGAGGGGATTTTTGACCTCGTTGGCAATaa                                                                                 >  2:6820728/1‑96 (MQ=255)
             gtggtgCTGGCTCTGCTGATACCGTGCGTGATATCCGTGGCTTTGCCACCAAGTTCTATACCGAAGAGGGTATTTTTGACCTCGTTGGCAATAaca                                                                               <  2:9399965/96‑1 (MQ=255)
               ggtgCTGGCTCTGCTGATACCGTGCGTGATATCCGTGGCTTTGCCACCAAGTTCTATACCGAAGAGGGTATTTTTGACCTCGTTGGCAATAACACg                                                                             <  2:9697322/96‑1 (MQ=255)
                gtgCTGGCTCTGCTGATACCGTGCGTGATATCCGTGGCTTTGCCACCAAGTTCTATACCGAAGAGGGTATTTTTGACCTCGTTGGCAATAACACGc                                                                            <  1:12773038/96‑1 (MQ=255)
                 tgCTGGCTCTGCTGATACCGTGCGTGATATCCGTGGCTTTGCCACCAAGTTCTATACCGAAGAGGGTATTTTTGACCTCGTTGGCAATAACACGcc                                                                           >  1:4266537/1‑96 (MQ=255)
                 tgCTGGCTCTGCTGATACCGTGCGTGATATCCGTGGCTTTGCCACCAAGTTCTATACCGAAGAGGGTATTTTTGACCTCGTTGGCAATAACACGcc                                                                           <  2:775372/96‑1 (MQ=255)
                 tgCTGGCTCTGCTGATACCGTGCGTGATATCCGTGGCTTTGCCACCAAGTTCTATACCGAAGAGGGTATTTTTGACCTCGTTGGCAATAACACGcc                                                                           <  1:1460753/96‑1 (MQ=255)
                 tgCTGGCTCTGCTGATACCGTGCGTGATATCCGTGGCTTTGCCACCAAGTTCTATACCGAAGAGGGTATTTTTGACCTCGTTGGCAACAACACGcc                                                                           <  1:12984564/96‑1 (MQ=255)
                 tgCTGGCTCTGCGGATACCGTGCGTGATATCCGTGGCTTTGCCACCAAGTTCTATACCGAAGAGGGGATTTTTGACCTCGTTGGCAATAACACGca                                                                           >  1:4535850/1‑95 (MQ=255)
                    tGGCTCTGCTGATACCGTGCGTGATATCCGTGGCTTTGCCACCAAGTTCTATACCGAAGAGGGGATTTTTGACCTCGTTGGCAATAACACGCCAAt                                                                        >  1:8995353/1‑96 (MQ=255)
                     ggCTCTGCTGATACCGTGCGTGATATCCGTGGCTTTGCCACCAAGTTCTATACCGAAGAGGGTATTTTTGACCTCGTTGGCAATAACACGCCAAtc                                                                       <  1:12641399/96‑1 (MQ=255)
                     ggCTCTGCTGATACCGTGCGTGATATCCGTGGCTTTGCCACCAAGTTCTATACCGAAGAGGGTATTTTTGACCTCGTTGGCAATAACACGCCAAtc                                                                       <  2:3081501/96‑1 (MQ=255)
                      gCTCTGCTGATACCGTGCGTGATATCCGTGGCTTTGCCACCAAGTTCTATACCGAAGAGGGTATTTTTGACCTCGTTGGCAATAACACGCCAAtct                                                                      <  2:2435147/96‑1 (MQ=255)
                       ctctGCTGATACCGTGCGTGATATCCGTGGCTTTGCCACCAAGTTCTATACCGAAGAGGGTATTTTTGACCTCGTTGGCAATAACACGCCAAtctt                                                                     >  1:5111811/1‑96 (MQ=255)
                        tctGCTGATTCCGTGCGTGATATCCGTGGCTTTGCCACCAAGTTCTATACCGAAGAGGGTATTTTTGACCTCGTTGGCAATAACACGCCAAtcttc                                                                    <  1:8720878/96‑1 (MQ=255)
                        tctGCTGATACCGTGCGTGATATCCGTGGCTTTGCCACCAAGTTCTATACCGAAGAGGGTATTTTTGACCTCGTTGGCAATAACACGCCAAtcttc                                                                    <  1:4778983/96‑1 (MQ=255)
                         ctgctgATACCGTGCGTGATATCCGTGGCTTTGCCACCAAGTTCTATACCGAAGAGGGTATTTTTGACCTCGTTGGCAATAACCCGCCAAtcttct                                                                   <  1:7895794/96‑1 (MQ=255)
                         ctgctgATACCGTGCGTGATATCCGTGGCTTTGCCACCAAGTTCTATACCGAAGAGGGTATTTTTGACCTCGTTGGCAATAACACGCCAAtcttct                                                                   <  2:11798530/96‑1 (MQ=255)
                           gctgATTCCGTGCGTGATATCCGTGGCTTTGCCACCAAGTTCTATCCCGAAGAGGGTATTTTTGACCTCGTTGGCAATAACACGCCAATCTTCttt                                                                 <  1:9999654/96‑1 (MQ=255)
                           gctgATACCGTGCGTGATATCCGTGGCTTTGCCACCAAGTTCTATACCGAAGAGGGTATTTTTGACCTCGTTGGCAATAACACGCCAATCTTCttt                                                                 <  2:4225499/96‑1 (MQ=255)
                           gctgATACCGTGCGTGATATCCGTGGCTTTGCCACCAAGTTCTATACCGAAGAGGGTATTTTTGACCTCGTTGGCAATAACACGCCAATCTTCttt                                                                 <  2:3336677/96‑1 (MQ=255)
                            ctgATACCGTGCGTGATATCCGTGGCTTTGCCACCAAGTTCTATACCGAAGAGGGTATTTTTGACCTCGTTGGCAATAACACGCCAATCTTCTTTa                                                                <  2:6576570/96‑1 (MQ=255)
                            ctgATACCGTGCGTGATATCCGTGGCTTTGCCACCAAGTTCTATACCGAAGAGGGTATTTTTGACCTCGTTGGCAATAACACGCCAATCTTCTTTa                                                                <  2:12070148/96‑1 (MQ=255)
                             tgATACCGTGCGTGATATCCGTGGCTTTGCCACCAAGTTCTATACCGAAGAGGGTATTTTTGACCTCGTTGGCAATAACACGCCAATCTTCTTTAt                                                               >  1:5421410/1‑96 (MQ=255)
                             tgATACCGTGCGTGATATCCGTGGCTTTGCCACCAAGTTCTATACCGAAGAGGGTATTTTTGACCTCGTTGGCAATAACACGCCAATCTTCTTTAt                                                               >  1:2285712/1‑96 (MQ=255)
                              gATACCGTGCGTGATATCCGTGGCTTTGCCACCAAGTTCTATACCGAAGAGGGTATTTTTGACCTCGTTGGCAATAACACGCCAATCTTCTTTATc                                                              <  1:6105752/96‑1 (MQ=255)
                              gATACCGTGCGTGATATCCGTGGCTTTGCCACCAAGTTCTATACCGAAGAGGGTATTTTTGACCTCGTTGGCAATAACACGCCAATCTTCTTTATc                                                              <  1:131243/96‑1 (MQ=255)
                               aTACCGTGCGTGATATCCGTGGCTTTGCCACCAAGTTCTATACCGAAGAGGGTATTTTTGACCTCGTTGGCAATAACACGCCAATCTTCTTTATcc                                                             >  1:806914/1‑96 (MQ=255)
                               aTACCGTGCGTGATATCCGTGGCTTTGCCACCAAGTTCTATACCGAAGAGGGTATTTTTGACCTCGTTGGCAATAACACGCCAATCTTCTTTATcc                                                             <  2:1488201/96‑1 (MQ=255)
                                tACCGGGCGGGATATCCGTGGCTTTGCCACCAAGTTCTATACCGAAGAGGGTATTTTTGACCTCGTTGGCAATAACACGCCAATCTTCTTTATCCa                                                            >  1:5650498/1‑96 (MQ=255)
                                 aCCGTGCGTGATATCCGTGGCTTTGCCACCAAGTTCTATACCGAAGAGGGTATTTTTGACCTCGTTGGCAATAACACGCCAATCTTCTTTATCCAg                                                           <  2:13832930/96‑1 (MQ=255)
                                 aCCGTGCGTGATATCCGTGGCTTTGCCACCAAGTTCTATACCGAAGAGGGTATTTTTGACCTCGTTGGCAATAACACGCCAATCTTCTTTATCCAg                                                           >  1:6402802/1‑96 (MQ=255)
                                  ccGTGCGTGATATCCGTGGCTTTGCCACCAAGTTCTATACCGAAGAGGGTATTTTTGACCTCGTTGGCAATAACACGCCAATCTTCTTTATCCAgg                                                          <  2:13107545/96‑1 (MQ=255)
                                       cgtgATATCCGTGGCTTTGCCACCAAGTTCTATACCGAAGAGGGTATTTTTGACCTCGTTGGCAATAACACGCCAATCTTCTTTATCCAGGATgcg                                                     >  1:10040312/1‑96 (MQ=255)
                                       cgtgATATCCGTGGCTTTGCCACCAAGTTCTATACCGAAGAGGGTATTTTTGACCTCGTTGGCAATAACACGCCAATCTTCTTTATCCAGGATgcg                                                     <  1:3248611/96‑1 (MQ=255)
                                       cgtgATATCCGTGGCTTTGCCACCAAGTTCTATACCGAAGAGGGTATTTTTGACCTCGTTGGCAATAACACGCCAATCTTCTTTATCCAGGATgcg                                                     >  1:1145288/1‑96 (MQ=255)
                                       cgtgATATCCGTGGCTTTGCCACCAAGTTCTATACCGAAGAGGGTATTTTTGACCTCGTTGGCAATAACACGCAAATCTTCTTTATCCAGGATgcg                                                     <  1:2537769/96‑1 (MQ=255)
                                        gtgATATCCGTGGCTTTGCCACCAAGTTCTATACCGAAGAGGGCATTTTTGACCTCGTTGGCAATAACACGCCAATCTTCTTTATCCAGGATgcgc                                                    >  1:10753339/1‑96 (MQ=255)
                                         tgATATCCGTGGCTTTGCCACCAAGTTCTATACCGAAGAGGGTATTTTTGACCTCGTTGGCAATAACACGCCAATCTTCTTTATCCAGGATGCGCa                                                   <  2:8916360/96‑1 (MQ=255)
                                           atatCCGTGGCTTTGCCACCAAGTTCTATACCGAAGAGGGTATTTTTGACCTCGTTGGCAATAACACGCCAATCTTCTTTATCCAGGATGCGCATa                                                 >  2:532130/1‑96 (MQ=255)
                                            tatCCGTGGCTTTGCCACCAAGTTCTATACCGAAGAGGGTATTTTTGACCTCGTTGGCAATAACACGCCAATCTTCTTTATCCAGGATGCGCATaa                                                >  1:12231884/1‑96 (MQ=255)
                                               ccGTGGCTTTGCCACCAAGTTCTATACCGAAGAGGGTATTTTTGACCTCGTTGGCAATAACACGCCAATCTTCTTTATCCAGGATGCGCATAAAtt                                             <  1:8340803/96‑1 (MQ=255)
                                                 gTGGCTTTGCCACCAAGTTCTATACCGAAGAGGGTATTTTTGACCTCGTTGGCAATAACACGCCAATCTTCTTTATCCAGGATGCGCATAAATTcc                                           <  1:13675830/96‑1 (MQ=255)
                                                 gTGGCTTTGCCACCAAGTTCTATACCGAAGAGGGTATTTTTGACCTCGTTGGCAATAACACGCCAATCTTCTTTATCCAGGATGCGCATAAATTcc                                           >  1:6859253/1‑96 (MQ=255)
                                                    gCTTTGCCCCCAAGTTCTATACCGAAGAGGGTATTTTTGACCTCGTTGGCAATAACACGCCAATCTTCTTTATCCAGGATGCGCATAAATTCCCCg                                        <  2:585158/96‑1 (MQ=255)
                                                     cTTTGCCACCAAGTTCTATACCGAAGAGGGTATTTTTGACCTCGTTGGCAATAACACGCCAATCTTCTTTATCCAGGATGCGCATAAATTCCCCGa                                       >  2:2109484/1‑96 (MQ=255)
                                                        tGCCACCAAGTTCTATACCGAAGAGGGTATTTTTGACCTCGTTGGCAATAACACGCCAATCTTCTTTATCCAGGATGCGCATAAATTCCCCGAttt                                    <  2:11209608/96‑1 (MQ=255)
                                                          ccaccaAGTTCTATACCGAAGAGGGTATTTTTGACCTCGTTGGCAATAACACGCCAATCTTCTTTATCCAGGATGCGCATAAATTCCCCGATTTTg                                  <  1:582431/96‑1 (MQ=255)
                                                          ccaccaAGTTCTATACCGAAGAGGGTATTTTTGACCTCGTTGGCAATAACACGCCAATCTTCTTTATCCAGGATGCGCATAAATTCCCCGATTTTg                                  >  2:4433587/1‑96 (MQ=255)
                                                               aaGTTCTATACCGAAGAGGGTATTTTTGACCTCGTTGGCAATAACACGCCAATCTTCTTTATCCAGGATGCGCATAAATTCCCCGATTTTGTTCAt                             <  1:1923830/96‑1 (MQ=255)
                                                                   tCTATACCGAAGAGGGTATTTTTGACCTCGTTGGCAATAACACGCCAATCTTCTTTATCCAGGATGCGCATAAATTCCCCGATTTTGTTCATGCgg                         >  2:2770437/1‑96 (MQ=255)
                                                                   tCTATACCGAAGAGGGTATTTTTGACCTCGTTGGCAATAACACGCCAATCTTCTTTATCCAGGATGCGCATAAATTCCCCGATTTTGTTCATGCgg                         >  2:4057463/1‑96 (MQ=255)
                                                                    cTATACCGAAGAGGGTATTTTTGACCTCGTTGGCAATAACACGCCAATCTTCTTTATCCAGGATGCGCATAAATTCCCCGATTTTGTTCATGCGGt                        <  1:14226862/96‑1 (MQ=255)
                                                                     tataCCGAAGAGGGTATTTTTGACCTCGTTGGCAATAACACGCCAATCTTCTTTATCCAGGATGCGCATAAATTCCCCGATTTTGTTCATGCGGTa                       >  2:13409901/1‑96 (MQ=255)
                                                                     tataCCGAAGAGGGTATTTTTGACCTCGTTGGCAATAACACGCCAATCTTCTTTATCCAGGATGCGCATAAATTCCCCGATTTTGTTCATGCGGTa                       <  1:948559/96‑1 (MQ=255)
                                                                     tataCCGAAGAGGGTATTTTTGACCTCGTTGGCAATAACACGCCAATCTTCTTTATCCAGGATGCGCATAAATTCCCCGATTTTGTTCATGCGGTa                       <  2:7438072/96‑1 (MQ=255)
                                                                       taCCGAAGAGGGTATTTTTGACCTCGTTGGCAATAACACGCCAATCTTCTTTATCCAGGATGCGCATAAATTCCCCGATTTTGTTCATGCGGTaaa                     <  2:1828945/96‑1 (MQ=255)
                                                                        aCCGAAGAGGGTATTTTTGACCTCGTTGGCAATAACCCGCCAATCTTCTTTATCCAGGATGCGCATAAATTCCCCGATTTTGTTCATGCGGTaaaa                    <  2:6126246/96‑1 (MQ=255)
                                                                          cGAAGAGGGTATTTTTGACCTCGTTGGCAATAACACGCCAATCTTCTTTATCCAGGATGCGCATAAATTCCCCGATTTTGTTCATGCGGTAAAAcc                  <  2:5650498/96‑1 (MQ=255)
                                                                           aaaGAGGGTATTTTTGACCTCGTTGGCAATAACACGCCAATCTTCTTTATCCAGGATGCGCATAAATTCCCCGCTTTTGTTCATGCGGTAAAACCa                 >  1:11307576/2‑96 (MQ=255)
                                                                            aaGAGGGTATTTTTGACCTCGTTGGCAATAACACGCCAATCTTCTTTATCCAGGATGCGCATAAATTCCCCGATTTTGTTCATGCGGTAAAACCAg                >  2:1274129/1‑96 (MQ=255)
                                                                             agagGGTATTTTTGACCTCGTTGGCAATAACACGCCAATCTTCTTTATCCAGGATGCGCATAAATTCCCCGATTTTGTTCATGCGGTAAAACCAGa               >  1:13383516/1‑96 (MQ=255)
                                                                             agagGGTATTTTTGACCTCGTTGGCAATAACACGCCAATCTTCTTTATCCAGGATGCGCATAAATTCCCCGATTTTGTTCATGCGGTAAAACCAGa               <  2:1145288/96‑1 (MQ=255)
                                                                             agagGGTATTTTTGACCTCGTTGGCAATAACACGCCAATCTTCTTTATCCAGGATGCGCATAAATTCCCCGATTTTGTTCATGCGGTAAAACCAGa               <  1:1754213/96‑1 (MQ=255)
                                                                               agGGTATTTTTGACCTCGTTGGCAATAACACGCCAATCTTCTTTATCCAGGATGCGCATAAATTCCCCGATTTTGTTCATGCGGTAAAACCAGAAc             <  1:7177023/96‑1 (MQ=255)
                                                                               agGGTATTTTTGACCTCGTTGGCAATAACACGCCAATCTTCTTTATCCAGGATGCGCATAAATTCCCCGATTTTGTTCATGCGGTAAAACCAGAAc             >  1:1113983/1‑96 (MQ=255)
                                                                                gggTATTTTTGACCTCGTTGGCAATAACACGCCAATCTTCTTTATCCAGGATGCGCATAAATTCCCCGATTTTGTTCATGCGGTAAAACCAGAAcc            <  2:2550003/96‑1 (MQ=255)
                                                                                 ggTATTTTTGACCTCGTTGGCAATAACACGCCAATCTTCTTTATCCAGGATGCGCATAAATTCCCCGATTTTGTTCATGCGGTAAAACCAGAACCg           >  2:10234463/1‑96 (MQ=255)
                                                                                 ggTATTTTTGACCTCGTTGGCAATAACACGCCAATCTTCTTTATCCAGGATGCGCATAAATTCCCCGATTTTGTTCATGCGGTAAAACCAGAACCg           >  1:9158139/1‑96 (MQ=255)
                                                                                  gTATTTTTGACCTCGTTGGCAATAACACGCCAATCTTCTTTATCCAGGATGCGCATAAATTCCCCGATTTTGTTCATGCGGTAAAACCAGAACCGc          >  2:12338543/1‑96 (MQ=255)
                                                                                  gTATTTTTGACCTCGTTGGCAATAACACGCCAATCTTCTTTATCCAGGATGCGCATAAATTCCCCGATTTTGTTCATGCGGTAAAACCAGAACCGc          >  2:7007428/1‑96 (MQ=255)
                                                                                   tATTTTTGACCTCGTTGGCAATAACACGCCAATCTTCTTTATCCAGGATGCGCATAAATTCCCCGATTTTGTTCATGCGGTAAAACCAGAACCGCa         >  2:14350297/1‑96 (MQ=255)
                                                                                    aTTTTTGACCTCGTTGGCAATAACACGCCAATCTTCTTTATCCAGGATGCGCATAAATTCCCCGATTTTGTTCATGCGGTAAAACCAGAACCGCAc        <  1:10611665/96‑1 (MQ=255)
                                                                                     tttttGACCTCGTTGGCAATAACACGCCAATCTTCTTTATCCAGGATGCGCATAAATTCCCCGATTTTGTTCATGCGGTAAAACCAGAACCGCACt       <  2:10562806/96‑1 (MQ=255)
                                                                                     tttttGACCTCGTTGGCAATAACACGCCAATCTTCTTTATCCAGGATGCGCATAAATTCCCCGATTTTGTTCATGCGGTAAAACCAGAACCGCACt       <  2:806914/96‑1 (MQ=255)
                                                                                     tttttGACCTCGTTGGCAATAACACGCCAATCTTCTTTATCCAGGATGCGCATAAATTCCCCGATTTTGTTCATGCGGTAAAACCAGAACCGCACt       >  2:8476213/1‑96 (MQ=255)
                                                                                       tttGATCTCGTTGGCAATAACACGCCAATCTTCTTTATCCAGGATGCGCATAAATTCCCCGATTTTGTTCATGCGGTAAAACCAGAACCGCACTgg     <  1:2671513/96‑1 (MQ=255)
                                                                                       tttGACCTCGTTGGCAATAACACGCCAATCTTCTTTATCCAGGATGCGCATAAATTCCCCGATTTTGTTCATGCGGTAAAACCAGAACCGCACTgg     <  2:4581991/96‑1 (MQ=255)
                                                                                       tttGACCTCGTTGGCAATAACACGCCAATCTTCTTTATCCAGGATGCGCATAAATTCCCCGATTTTGTTCATGCGGTAAAACCAGAACCGCACTgg     >  1:971826/1‑96 (MQ=255)
                                                                                          gACCTCGTTGGCAATAACACGCCAATCTTCTTTATCCAGGATGCGCATAAATTCCCCGATTTTGTTCATGCGGTAAAACCAGAACCGCACTGGGCa  >  2:10633411/1‑96 (MQ=255)
                                                                                             |                                                                                            
TCTACCGTTCAGGGTGGTGCTGGCTCTGCTGATACCGTGCGTGATATCCGTGGCTTTGCCACCAAGTTCTATACCGAAGAGGGTATTTTTGACTTCGTTGGCAATAACACGCCAATCTTCTTTATCCAGGATGCGCATAAATTCCCCGATTTTGTTCATGCGGTAAAACCAGAACCGCACTGGGCA  >  NC_000913/1812389‑1812574

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 19 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: