Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,815,269 A→T 5.6% N468I (AAC→ATC)  katE → catalase HPII, heme d‑containing

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,815,2690AT 5.6% 232.3 / 3.5 72N468I (AAC→ATC) katEcatalase HPII, heme d‑containing
Reads supporting (aligned to +/- strand):  ref base A (37/31);  new base T (2/2);  total (39/33)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.01e-01

AATCAGTCGTCTTGGTGGGCCGAATTTCCATGAGATTCCGATTAACCGTCCGACCTGCCCTTACCATAATTTCCAGCGTGACGGCATGCATCGCATGGGGATCGACACTAACCCGGCGAATTACGAACCGAACTCGATTAACGATAACTGGCCGCGCGAAACACCGCCGGGGCCGAAACGCGGCGGTTTTGAATCATACCAGGAGCGCGTGGAAGGCAATAAAGTTCGCGAGCGCAGCCCATCGTTTGGCGAATATTATTCCCATCCGCGTCTGTTCTGGCTAAGTCAGAC  >  NC_000913/1815123‑1815413
                                                                                                                                                  |                                                                                                                                                
aaTCAGTCGTCTTGGTGGGCCGAATTTCCATGAGATTCCGATTAACCGTCCGACCTGCCCTTACCATAATTTCCAGCGTGACGGCATGCATCGCATGGGGATCGACACTAACCCGGCGAATTACGAACCGAACTCGATTAACGATAACt                                                                                                                                                >  2:66602/1‑149 (MQ=255)
aaTCAGTCGTCTTGGTGGGCCGAATTTCCATGAGATTCCGATTAACCGTCCGACCTGCCCTTACCATAATTTCCAGCGTGACGGCATGCATCGCATGGGGATCGACACTAACCCGGCGAATTACGAACCGAACTCGATTAACGATAACt                                                                                                                                                <  7:209460/149‑1 (MQ=255)
    aGTCGTCTTGGTGGGCCGAATTTCCATGAGATTCCGATTAACCGTCCGACCTGCCCTTACCATAATTTCCAGCGTGACGGCATGCATCGCATGGGGATCGACACTACCCCGGCGAATTACGAACCGAACTCGATTAACGATAACTGGcc                                                                                                                                            >  4:304752/1‑149 (MQ=255)
       cgtcTTGGTGGGCCGAATTTCCATGAGATTCCGATTAACCGTCCGACCTGCCCTTACCATAATTTCCAGCGTGACGGCATGCATCGCATGGGGATCGACACTAACCCGGCGAATTACGAACCGAACTCGATTAACGATAACTGGCcgcg                                                                                                                                         >  8:304580/1‑149 (MQ=255)
          cTTGGTGGGCCGAATTTCCATGAGATTCCGATTAACCGTCCGACCTGCCCTTACCATAATTTCCAGCGTGACGGCATGCATCGCATGGGGATCGACACTAACCCGGCGAATTACGAACCGAACTCGATTAACGATAACTg                                                                                                                                               >  8:305856/1‑140 (MQ=255)
          cTTGGTGGGCCGAATTTCCATGAGATTCCGATTAACCGTCCGACCTGCCCTTACCATAATTTCCAGCGTGACGGCATGCATCGCATGGGGATCGACACTAACCCGGCGAATTACGAACCGAACTCGATTAACGATAACTg                                                                                                                                               <  7:305856/140‑1 (MQ=255)
          cTTGGTGGGCCGAATTTCCATGAGATTCCGATTAACCGTCCGACCTGCCCTTACCATAATTTCCAGCGTGACGGCATGCATCGCATGGGGATCGACACTAACCCGGCGAATTACGAACCGAACTCGATTAACGATAACTg                                                                                                                                               <  7:311392/140‑1 (MQ=255)
          cTTGGTGGGCCGAATTTCCATGAGATTCCGATTAACCGTCCGACCTGCCCTTACCATAATTTCCAGCGTGACGGCATGCATCGCATGGGGATCGACACTAACCCGGCGAATTACGAACCGAACTCGATTAACGATAACTg                                                                                                                                               >  8:311392/1‑140 (MQ=255)
           ttGGTGGGCCGAATTTCCATGAGATTCCGATTAACCGTCCGACCTGCAATTACCATAATTTCCAGCGTGACGGCATGCATCGCATGGGGATCGACACTAACCCGTCGAATTACGAACCGAACTCGATTAACGATAACTGGCAGCGCGa                                                                                                                                      >  5:269105/1‑148 (MQ=255)
             ggtggGCCGAATTTCCATGAGATTCCGATTAACCGTCCGACCTGCCCTTACCATAATTTCCAGCGTGACGGCATGCATCGCATGGGGATCGACACTAACCCGGCGAATTACGAACCGAACTCGATTAACGATAACTGGCCGCGCGaaa                                                                                                                                    <  3:363486/148‑1 (MQ=255)
                  gCCGAATTTCCATGAGATTCCGATTAACCGTCCGACCTGCCCTTACCATAATTTCCAGCGTGACGGCATGCATCGCATGGGGATCGACACTAACCCGGCGAATTACGAACCGAACTCGATTAACGATAACTGGCCGCGCGAAACAccgc                                                                                                                              >  5:132903/1‑149 (MQ=255)
                     gAATTTCCATGAGATTCCGATTAACCGTCCGACCTGCCCTTACCATAATTTCCAGCGTGACGGCATGCATCGCATGGGGATCGACACTAACCCGGCGAATTACGAACCGAACTCGATTAACGATAACTGGCCGCGCGAAACACCGCCgg                                                                                                                           <  7:304580/149‑1 (MQ=255)
                            cATGAGATTCCGATTAACCGTCCGACCTGCCCTTACCATAATTTCCAGCGTGACGGCATGCATCGCATGGGGATCGCCACTAACCCGGCGAATTACGAACCGAACTCGATTACCGATAACTGGCCGCGCGAAACACCTCCGGGGCCGaa                                                                                                                    >  8:364178/1‑149 (MQ=255)
                            cATGAGATTCCGATTAACCGTCCGACCTGCCCTTACCATAATTTCCAGCGTGACGGCATGCATCGCATGGGGATCGACACTAACCCGGCGAATTACGAACCGAACTCGATTAACGATAACTGGCCGCGCGAAACACCGCCGGGGCCGa                                                                                                                     >  8:197549/1‑148 (MQ=255)
                             aTGAGATTCCGATTAACCGTCCGACCTGCCCTTACCATAATTTCCAGCGTGACGGCATGCATCGCATGGGGATCGACACTAACCCGGCGAATTACGAACCGAACTCGATTAACGATAACTGGCCGCGCGAAACACCGCCGGGGCCGaaa                                                                                                                   >  1:39896/1‑149 (MQ=255)
                                      cGATTAACCGTCCGACCTGCCCTTACCATAATTTCCAGCGTGACGGCATGCATCGCATGGGGATCGACACTAACCCGGCGAATTACGAACCGAACTCGATTAACGATAACTGGCCGCGCGAAACACCGCCGGGGCCGAAACGcggcgg                                                                                                           <  1:66602/148‑1 (MQ=255)
                                       gATTAACCGTCCGACCTGCCCTTACGATAATTTCCAGCGTGCCGGCATGCATCGCATGGGGATCGACACTAACCCGGCGAATTACGAACCGAACTCGATTAACGATAACTTGCCGCGCGAAACACCGCCGGGGCCTAACAGCGGCGGtt                                                                                                         >  2:238729/1‑149 (MQ=255)
                                            cccGTCCGCCCTGCCCTAACCATAATTTCCAGCGTGACGGCATGGATCGCATGGGGATCGACACTAACCCGGCGAATTACGAACCGAACTCGTTTAACGATAACTGGCCGCGCGAAACACCGCCGGGGCCGAAACGCGGCGGTTTTGa                                                                                                     <  8:135462/147‑1 (MQ=255)
                                            aCCGTCCGACCTTCCATTACCATAATTTCCAGCGTGACGGAATGCATCGCATGGGGATCGACACTAACCCGGCGAATTACGAACCGAACTCGATTAACGATAACTGGCCGCTCGAAACACCTCCGGGGCCGAAACTCGGCGGTTTTGaa                                                                                                    <  3:66808/149‑1 (MQ=255)
                                            aCCGTCCGACCTGCCCTTACCATAATTTCCAGCGTGACGGCATGCATCGCATGGGGATCGACACTAACCCGGCGAATTACGAACCGAACTCGATTAACGATAACTGGCCGCGCGAAACACCGCCGGGGCCGAAACGCGGCGGTTTTGaa                                                                                                    >  1:120289/1‑149 (MQ=255)
                                            aCCGTCCGACCTGCCCTTACCATAATTTCCAGCGTGACGGCATGCATCGCATGGGGATCGACACTAACCCGGCGAATTACGAACCGAACTCGATTAACGATAACTGGCCGCGCGAAACACCGCCGGGGCCGAAACGCGGCGGTTTTGaa                                                                                                    <  1:202050/149‑1 (MQ=255)
                                             ccGTCCGACCTGCCCTTACCATAATTTCCAGCGTGACGGCATGCATCGCATGGGGATCGACACTAACCCGGCGAATTACGAACCGAACTCGATTAACGATAACTGGCCGCGCGAAACACCGCCGGGGCCGAAACGCGGCGGTTTTGAAt                                                                                                   <  1:217791/149‑1 (MQ=255)
                                                  cGACCTGCCCTTACCATAATTTCCAGCGTGACGGCATGCATCGCATGGGGATCGACACTAACCCGGCGAATTACGAACCGAACTCGATTAACGATAACTGGCCGCGCGAAACACCGCCGGGGCCGAAACGCGGCGGTTTTGAATCATa                                                                                               <  7:366674/148‑1 (MQ=255)
                                                    aCCTGCCCATACCATAATTTCCAGAGTGACGGCATGCATCGCATCGGGATCGACACTAACCAGGCGAATTACGAACCGAACTCGATTAACGATAACTGGCCGCGCGAAACACCGCCGGGGCCGAAACGCGCCTGATTTGAATAATAc                                                                                              >  3:271598/1‑147 (MQ=255)
                                                     ccTGCCCTTACCATAATTTCCAGCGTGACGGCATGCATCGCATGGGGATCGACACTAACCCGGCGAATTACGAACCGAACTCGATTAACGATAACTGGCCGCGCGAAACACCGCCGGGGCCGAAACGCGGCGGTTTTGAATCATACCAg                                                                                           <  4:330717/149‑1 (MQ=255)
                                                          ccTTACCATAATTTCCAGCGTGACGGCATGCATCGCATGGGGATCGACACTAACCCGGCGAATTACGAACCGAACTCGATTAACGATAACTGGCCGCGCGAAACACCGCCGGGGCCGAAACGCGGCGGTTTTGAATCATACCAGGAgcg                                                                                      >  7:52287/1‑149 (MQ=255)
                                                            ttACCATAATTTCCAGCGTGACGGCATGCATCGCATGGGGATCGACACTAACCCGGCGAATTACGATCCGAACTCGATTAACGATAACTGGCCGCGCGAAACACCGCCGGGGTCGAAACGCGGCGCTTTTGAATCATACCAGGAGcgcg                                                                                    >  2:377978/1‑149 (MQ=255)
                                                               ccATAATTTCCAGCGTGACGGCATGCATCGCATGGGGATCGACACTAACCCGGCGAATTACGAACCGAACTCGATTAACGATAACTGGCCGCGCGAAACACCGCCGGGGCCGAAACGCGGCGGttt                                                                                                        <  4:48772/126‑1 (MQ=255)
                                                                cATAATTTCCAGCGTGACGGCATGCATCGCATGGGGATCGACACTAACCCGGCGAATTACGAACCGAACTCGATTAACGATATCTGGCCGCGCGAAACACCGCCGGGGCCGAAACGCGGCGGTTTTGAATCATACCAGGAGCGCGTgg                                                                                 >  8:1713/1‑148 (MQ=255)
                                                                cATAATTTCCAGCGTGACGGCATGCATCGCATGGGGATCGACACTAACCCGGCGAATTACGAACCGAACTCGATTAACGATATCTGGCCGCGCGAAACACCGCCGGGGCCGAAACGCGGCGGTTTTGAATCATACCAGGAGCGCGTgg                                                                                 >  8:7149/1‑148 (MQ=255)
                                                                cATAATTTCCAGCGTGACGGCATGCATCGCATGGGGATCGACACTAACCCGGCGAATTACGAACCGAACTCGATTAACGATAACTGGCCGCGCGAAACACCGCCGGGGCCGAAACGCGGCGGTTTTGAATCATACCAGGAGCGCGTgg                                                                                 >  3:234094/1‑148 (MQ=255)
                                                                 aTAATTTCCAGCGTGACGGCATGCATCGCATGGGGATCGACACTAACCCGGCGAATTACGAACCGAACTCGATTAACGATAACTGGCCGCGCGAAACACCGCCGGGGCCGAAACGCGGCGGTTTTGAATCATACCAGGAGCGCGTGGaa                                                                               <  7:197549/149‑1 (MQ=255)
                                                                     tttCCAGCGTGACGGCATGCATCGCATGGGGATCGACACTAACCCGGCGAATTACGAACCGAACTCGATTAACGATAACTGGCCGCGCGAAACACCGCCGGGGCCGAAACGCGGCGGTTTTGAATCATACCAGGAGCGCGTGGAAGGc                                                                            >  7:145543/1‑148 (MQ=255)
                                                                     tttCCAGCGTGACGGCATGCATCGCATGGGGATCGACACTAACCCGGCGAATTACGAACCGAACTCGATTAACGATAACTGGCCGCGCGAAACACCGCCGGGGCCGAAACGCGGCGGTTTTGAATCATACCAGGAGCGCGTGGAAGGCa                                                                           <  5:159053/149‑1 (MQ=255)
                                                                        ccAGCGTGACGGCATGCATCGCATGGGGATCGACACTAACCCGGCGAATTACGAACCGAACTCGATTAACGATAACTGGCCGCGCGAAACACCGCCGGGGCCGAAACGCGGCGGTTTTGAATCATACCAGGAGCGCGTGGAAGGCAAt                                                                         >  2:140614/1‑148 (MQ=255)
                                                                          aGCGTGACGGCATGCATCGCATGGGGATCGACACTAACCCGGCGAATTACGAACCGAACTCGATTAACGATAACTGGCCGCGCGAAACACCGCCGGGGCCGAAACGCGGCGGTTTTGAATCATACCAGGAGCGCGTGGAAGGCAATaaa                                                                      <  2:39896/149‑1 (MQ=255)
                                                                           gCGTGACGGCATGCATCGCATGGGGATCGACACTAACCCGGCGAATTACGAACCGAACTCGATTAACGATAACTGGCCGCGCGAAACACCGCCGGGGCCGAAACGCGGCGGTTTTGAATCATACCAGGAGCGCGTGGAAGGCAATAAAg                                                                     >  6:38695/1‑149 (MQ=255)
                                                                           gCGTGACGGCATGCATCGCATGGGGATCGACACTAACCCGGCGAATTACGAACCGAACTCGATTAACGATAACTGGCCGCGCGAAACACCGCCGGGGCCGAAACGCGGCGGTTTTGAATCATACCAGGAGCGCGTGGAAGGCAATAAAg                                                                     <  5:33898/149‑1 (MQ=255)
                                                                           gCGTGACGGCATGCATCGCATGGGGATCGACACTAACCCGGCGAATTACGAACCGAACTCGATTAACGATAACTGGCCGCGCGAAACACCGCCGGGGCCGAAACGCGGCGGTTTTGAATCATACCAGGAGCGCGTGGAAGGCAATAAAg                                                                     <  5:38695/149‑1 (MQ=255)
                                                                           gCGTGACGGCATGCATCGCATGGGGATCGACACTAACCCGGCGAATTACGAACCGAACTCGATTAACGATAACTGGCCGCGCGAAACACCGCCGGGGCCGAAACGCGGCGGTTTTGAATCATACCAGGAGCGCGTGGAAGGCAATAAAg                                                                     >  6:33898/1‑149 (MQ=255)
                                                                              tGACGGCATGCATCGCATGGGGATCGACACTAACCCGGCGAATTACGAACCGAACTCGATTAACGATATCTGGCCGCGCGAAACACCGCCGGGGCCGAAACGCGGCGGTTTTGAATCATACCAGGAGCGCGTGGAAGGCAATAAAGTTc                                                                  <  7:1713/149‑1 (MQ=255)
                                                                              tGACGGCATGCATCGCATGGGGATCGACACTAACCCGGCGAATTACGAACCGAACTCGATTAACGATATCTGGCCGCGCGAAACACCGCCGGGGCCGAAACGCGGCGGTTTTGAATCATACCAGGAGCGCGTGGAAGGCAATAAAGTTc                                                                  <  7:7149/149‑1 (MQ=255)
                                                                               gccGGCATGCATCGCATGGGGATCGACACTAACCCGGCGAATTACGAACCGAACTCGATTAACGATAACTGGCCGCGCGAAACACCGCCGGGGCCGAAACGCGGCGGTTTTGAATCATACCAGGAGCGCGTGGAAGGCAATAAAGTTcg                                                                 >  2:254038/3‑149 (MQ=255)
                                                                                 cGGCATGCATCGCATGGGGATCGACACTAACCCGGCGAATTACGAACCGAACTCGATTAACGATAACTGGCCGCGCGAAACACCGCCGGGGCCGAAACGCGGCGGTTTTGAATCATACCAg                                                                                           <  8:215012/121‑1 (MQ=255)
                                                                                 cGGCATGCATCGCATGGGGATCGACACTAACCCGGCGAATTACGAACCGAACTCGATTAACGATAACTGGCCGCGCGAAACACCGCCGGGGCCGAAACGCGGCGGTTTTGAATCATACCAg                                                                                           >  7:215012/1‑121 (MQ=255)
                                                                                  ggCATGCATCGCATGGGGATCGACACTAACCCGGCGAATTACGAACCGAACTCGATTAACGATAACTGGCCGCGCGAAACACCGCCGGGGCCGAACCGCGGCGGTTTTGAATCATACCAGGAGCGCGTGGAAGGCAATAAAGTTcgc                                                                >  5:164684/1‑147 (MQ=255)
                                                                                  ggCATGCATCGCATGGGGATCGACACTAACCCGGCGAATTACGAACCGAACTCGATTAACGATAACTGGCCGCGCGAAACACCGCCGGGGCCGAAACGCGGCGGTTTTGAATCATACCAGGAGCGCGTGGAAGGCAATAAAGTTCGCGa                                                              <  2:153165/149‑1 (MQ=255)
                                                                                         atCGCATGGGGATCGACACTAACCCGGCGAATTACGAACCGAACTCGATTAACGATAACTGGCCGCGCGAAACACCGCCGGGGCCGAAACGCGGCGGTTTTGAATCATACCAGGAGCGCGTGGAAGGCAATAAAGTTCGCGAGCGCAGc                                                       >  8:212581/1‑149 (MQ=255)
                                                                                         atCGCATGGGGATCGACACTAACCCGGCGAATTACGAACCGAACTCGATTAACGATAACTGGCCGCGCGAAACACCGCCGGGGCCGAAACGCGGCGGTTTTGAATCATACCAGGAGCGCGTGGAAGGCAATAAAGTTCGCGAGCGCAGc                                                       >  7:197797/1‑149 (MQ=255)
                                                                                                ggggATCGACACTAACCCGGCGAATTACGAACCGAACTCGATTAACGATAACTGGCCGCGCGAAACACCGCCGGGGCCGAAACGCGGCGGTTTTGAATCATACCAGGAGCGCGTGGAAGGCAATAAAGTTCGCGAGCGCAGCCCATCGt                                                >  2:202280/1‑149 (MQ=255)
                                                                                                   gATCGACCCTAACCCCGCGAATTACGAACCGAACTCGATTAACGATAAATGGCCGCTCGAAACACCGCCGAGGCCGAAACGCGGCGGTGTTGAATCATACCATCAGCGCGTGGAAGGCAATAAACTACGCGAGCGAAGCCCATCGTTTg                                             >  8:160117/1‑149 (MQ=255)
                                                                                                   gATCGACACTAACCCGGCGAATTACGAACCGAACTCGATTAACGATAACTGGCCGCGCGAAACACCGCCGGGGCCGAAACGCGGCGGTTTTGAATCATACCAGGAGCGCGTGGAAGGCAATAAAGTTCGCGAGCGCAGCCCATCGTTTg                                             <  7:212581/149‑1 (MQ=255)
                                                                                                   gATCGACACTAACCCGGCGAATTACGAACCGAACTCGATTAACGATAACTGGCCGCGCGAAACACCGCCGGGGCAGAAACGCGGCGGTTTTGAATCATACCAGGAGCGCGTGGAAGGCAATAAAGTTCGCGAGCGCAGCCCATCGTTTg                                             <  6:132903/149‑1 (MQ=255)
                                                                                                        acacTAACCCGGCGAATTACGAACCGAACTCGATTAACGATAACTGGCCGCGCGAAACACCGCCGGGGCCGAAACGCGGCGGTTTTGAATCATACCAGGCGCGCGTGGAAGGCAATAAAGTTCGCGAGCGCAGCCCATCGTTTGGCGaa                                        >  5:55326/1‑149 (MQ=255)
                                                                                                        acacTAACCCGGCGAATTACGAACCGAACTCGATTAACGATAACTGGCCGCGCGAAACACCGCCGGGGCCGAAACGCGGCGGTTTTGAATCATACCAGGAGCGCGTGGAAGGCAATAAAGTTCGCGAGCGCAGCCCATCGTTTGGCGaa                                        <  1:254038/149‑1 (MQ=255)
                                                                                                         cacTAACCCGGCGAATTACGAACCGAACTCGATTAACGATAACTGGCCGCGCGAAACACCGCCGGGGCCGAAACGCGGCGGTTTTGAATCATACCAGGAGCGCGTGGAAGGCAATAAAGTTCGCGAGCGCAGCCCATCGTTTGGCGaa                                        >  4:71783/1‑148 (MQ=255)
                                                                                                           cTAACCCGGCGAATTACGAACCGAACTCGATTAACGATAACTGGCCGCGCGAAACACCGCCGGGGCCGAAACGCGGCGGTTTTGAATCATACCAGGAGCGCGTGGAAGGCAATAAAGTTCGCGAGCGCAGCCCATCGTTTGGCGAatat                                     >  7:262540/1‑149 (MQ=255)
                                                                                                               cccGGCGAATTACGAACCGAACTCGATTAACGATAACTGGCCGCGCGAAACACCGCCGGGGCCGAAACGCGGCGGTTTTGAATCATACCAGGAGCGCGTGGAAGGCAATAAAGTTCGCGAGCGCAGCCCATCGTTTGGCGAATattatt                                 <  1:202280/149‑1 (MQ=255)
                                                                                                                  ggCGAATTACGAACCGAACTCGAATAACGATAACTGGCCGCGCGAAACACCGCCGGGGCCGAAACGCGGCGGTTTTGAAACATACCAGGAGCGCGTGGAAGGCAATAAAGTTCGCGAGCGCAGCACATCGTTTGGCGAATATTATTcc                               >  6:286062/1‑148 (MQ=255)
                                                                                                                    cGAATTACGAACCGAACTCGATTAACGATAACTGGCCGCGCGAAACACCGCCGGGGCCGAAACGCGGCGGTTTTGAATCATACCAGGAGCGCGTGGAAGGCAATAAAGTTCGCGAGCGCAGCCCATCGTTTGGCGAATATTATTCCCAt                            <  1:140614/149‑1 (MQ=255)
                                                                                                                              accgaacTCGATTAACGATAACTGGCCGCGCGAAACACCGCCGGGGCCGAAACGCGGCGGTTTTGAATCATACCAGGAGCGCGTGGAAGGCAATAAAGTTCGCGAGCGCAGCCCATCGTTTGGCGAATATTATTCCCATCCGCGTCTGt                  >  3:75633/1‑149 (MQ=255)
                                                                                                                                   acTCGATTAACGATAACTGGCCGCGCGAAACACCGCCGGGGCCGAAACGCGGCGGTTTTGAATCATACCAGGAGCGCGTGGAAGGCAATAAAGTTCGCGAGCGCAGCCCATCGTTTGGCGAATATTATTCCCATCCGCGTCTGTTCTgg             <  2:120289/149‑1 (MQ=255)
                                                                                                                                     tCGATTAACGATAACTGGCCGCGCGAAACACCGCCGGGGCCGAAACGCGGCGGTTTTGAATCATACCAGGAGCGCGTGGAAGGCAATAAAGTTCGCGAGCGCAGCCCATCGTTTGGCGAATATTATTCCCATCcgc                        >  2:171874/1‑136 (MQ=255)
                                                                                                                                     tCGATTAACGATAACTGGCCGCGCGAAACACCGCCGGGGCCGAAACGCGGCGGTTTTGAATCATACCAGGAGCGCGTGGAAGGCAATAAAGTTCGCGAGCGCAGCCCATCGTTTGGCGAATATTATTCCCATCcgc                        <  1:171874/136‑1 (MQ=255)
                                                                                                                                     tCGATTAACGATAACTGGCCGCGCGAAACACCGCCGGGGCCGAAACGCGGCGGTTTTGAATCATACCAGGAGCGCGTGGAAGGCAATAAAGTTCGCGAGCGCAGCCCATCGTTTGGCGAATATTATTCCCATCCGCGTCTGTTCTGGc            >  4:136756/1‑148 (MQ=255)
                                                                                                                                       gATTAACGATAACTGGCCGCGCGAAACACCGCCGGGGCCGAAACGCGGCGGTTTTGAATCATACCAGGAGCGCGTGGAAGGCAATAAAGTTCGCGAGCGCAGCCCATCGTTTGGCGAATATTATTCCCATCCGCGTCTGTTCTGGCTaa         <  4:75633/149‑1 (MQ=255)
                                                                                                                                       gATTAACGATAACTGGCCGCGCGAAACACCGCCGGGGCCGAAACGCGGCGGTTTTGAATCATACCAGGAGCGCGTGGAAGGCAATAAAGTTCGCGAGCGCAGCCCATCGTTTGGCGAATATTATTCCCATCCGCGTCTGTTCTGGCTaa         <  8:262540/149‑1 (MQ=255)
                                                                                                                                       gATTAACGATAACTGGCCGCGCGAAACACCGCCGGGGCCGAAACGCGGCGGTTTTGAATCATACCAGGAGCGCGTGGAAGGCAATAAAGTTCGCGAGCGCAGCCCATCGTTTGGCGAATATTATTCCCATCCGCGTCTGTTCTGGCTaa         <  4:234094/149‑1 (MQ=255)
                                                                                                                                            aCGATAACTGGCCGCGCGAAACACCGCCGGGGCCGAAACGCGGCGGTTTTGAATCATACAAGGAGCGCGTGGAAGGCAATAAAGTTCGCGAGCGCAGCCCATCGTTTGGCGAATATTATTCCCATCCGCGTCTGTTCTGGCTAAGt       >  5:272649/1‑146 (MQ=255)
                                                                                                                                              gATAACTGGCCGCGCGAAACACCGCCGGGGCCGAAACGCGGCGGTTTTGAATCATACCAGGAGCGCGTGGAAGGCAATAAAGTTcgc                                                                >  7:395274/1‑87 (MQ=255)
                                                                                                                                              gATAACTGGCCGCGCGAAACACCGCCGGGGCCGAAACGCGGCGGTTTTGAATCATACCAGGAGCGCGTGGAAGGCAATAAAGTTcgc                                                                <  8:395274/87‑1 (MQ=255)
                                                                                                                                              gATAACTGGCCGCGCGAAACACCGCCGGGGCCGAAACGCGGCGGTTTTGAATCATACCAGGAGCGCGTGGAAGGCAATAAAGTTCGCGAGCGCAGCCCATCGTTTGGCGAATATTATTCCCATCCGCGTCTGTTCTGGCTAAGTCAGAc  >  1:110370/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                                
AATCAGTCGTCTTGGTGGGCCGAATTTCCATGAGATTCCGATTAACCGTCCGACCTGCCCTTACCATAATTTCCAGCGTGACGGCATGCATCGCATGGGGATCGACACTAACCCGGCGAATTACGAACCGAACTCGATTAACGATAACTGGCCGCGCGAAACACCGCCGGGGCCGAAACGCGGCGGTTTTGAATCATACCAGGAGCGCGTGGAAGGCAATAAAGTTCGCGAGCGCAGCCCATCGTTTGGCGAATATTATTCCCATCCGCGTCTGTTCTGGCTAAGTCAGAC  >  NC_000913/1815123‑1815413

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