Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 565,068 G→A 5.6% G182D (GGC→GAC)  katA → catalase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947565,0680GA 5.6% 245.2 / 2.1 82G182D (GGC→GAC) katAcatalase
Reads supporting (aligned to +/- strand):  ref base G (37/40);  new base A (3/2);  total (40/42)
Fisher's exact test for biased strand distribution p-value = 6.72e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TTCGCGACCCACTTAAGTTCCCGGACTTTATCCACACCCAGAAGCGCCACCCGCAGTCCAACCTGAAGAATGCTCAGATGATGTGGGACTTCTGGTCGCACTCTCCAGAAGCGCTGCATCAGGTCACCATCCTGTTTTCCGATCGTGGCATTCCGGACGGCTACCGTCACATGCATGGCTTTGGCAGCCACACCTACAGCCTGATCAATGCAAAGGGTGAGCGCACTTGGGTCAAGTGGCACTTCAAGACCCAGCAAGGCATCAAGAACCTCACCCCGGCAGACGCTGCACGCCTG  >  NC_002947/564921‑565216
                                                                                                                                                   |                                                                                                                                                    
ttCGCGACCCACTTAAGTTCCCGGACTTTATCCACACCCAGAAGCGCCACCCGCAGTCCAACCTGAAGAATGCTCAGATGATGTGGGACTTCTGGTCGCACTCTCCAGAAGCGCTGCATCAGGTCACCATCCTGTTTTCCGATCGTGGCa                                                                                                                                                    <  1:893937/150‑1 (MQ=255)
    cgACCCACTTAAGTTCCCGGACTTTATCCACACCCAGAAGCGCCACCCGCAGTCCAACCTGAAGAATGCTCAGATGATGTGGGACTTCTGGTCGCACTCTCCAGAAGCGCTGCATCAGGTCACCATCCTGTTTTCCGATCGTGGCATTcc                                                                                                                                                <  1:1873582/150‑1 (MQ=255)
    cgACCCACTTAAGTTCCCGGACTTTATCCACACCCAGAAGCGCCACCCGCAGTCCAACCTGAAGAATGCTCAGATGATGTGGGACTTCTGGTCGCACTCTCCAGAAGCGCTGCATCAGGTCACCATCCTGTTTTCCGATCGTGGCATTcc                                                                                                                                                <  1:341433/150‑1 (MQ=255)
        ccACTTAAGTTCCCGGACTTTATCCACACCCAGAAGCGCCGCCCGCAGTCCCACCGGAAGAATGCTCAGATGATGTGGGACTTCTGGTCGCACTCTCCAGAAGCGCTGCATCAGGTCACCATCCTGTTTTCCGATCGTGGCATTCCGGAc                                                                                                                                            <  1:623806/150‑1 (MQ=255)
          aCTTAAGTTCCCGGACTTTATCCACACCCAGAAGCGCCACCCGCAGTCCAACCTGAAGAATGCTCAGATGATGTGGGACTTCTGGTCGCACTCTCCAGAAGCGCTGCATCAGGTCACCATCCTGTTTTCCGATCGTGGCATTCCGGACgg                                                                                                                                          >  1:302569/1‑150 (MQ=255)
                 ttCCCGGACTTTATCCACACCCAGAAGCGCCACCCGCAGTCCAACCTGAAGAATGCTCAGATGATGTGGGACTTCTGGTCGCACTCTCCAGAAGCGCTGCATCAGGTCACCATCCTGTTTTCCGATCGTGGCATTCCGGACGGCTACCGt                                                                                                                                   >  2:98351/1‑150 (MQ=255)
                    ccGGACTTTATCCACACCCAGAAGCGCCACCCGCAGTCCAACCTGAAGAATGCTCAGATGATGTGGGACTTCTGGTCGCACTCTCCAGAAGCGCTGCATCAGGTCACCATCCTGTTTTCCGATCGTGGCATTCCGGACGGCTACCGTcac                                                                                                                                >  2:1100074/1‑150 (MQ=255)
                    ccGGACTTTATCCACACCCAGAAGCGCCACCCGCAGTCCAACCTGAAGAATGCTCAGATGATGTGGGACTTCTGGTCGCACTCTCCAGAAGCGCTGCATCAGGTCACCATCCTGTTTTCCGATCGTGGCATTCCGGACGGCTACCGTcac                                                                                                                                >  1:939030/1‑150 (MQ=255)
                     cGGACTTTATCCACACCCAGAAGCGCCACCCGCAGTCCAACCTGAAGAATGCTCAGATGATGTGGGACTTCTGGTCGCACTCTCCAGAAGCGCTGCATCAGGTCACCATCCTGTTTTCCGATCGTGACATTc                                                                                                                                                 <  1:444230/132‑1 (MQ=255)
                     cGGACTTTATCCACACCCAGAAGCGCCACCCGCAGTCCAACCTGAAGAATGCTCAGATGATGTGGGACTTCTGGTCGCACTCTCCAGAAGCGCTGCATCAGGTCACCATCCTGTTTTCCGATCGTGACATTc                                                                                                                                                 >  2:444230/1‑132 (MQ=255)
                         cTTTCTCCCCACCCAGAAGCGCCCCCCGCAGTCCAACCTGAAGAATGCTCAGATGATGTGGGACTCCTGGTCGCACTCTCCAGAAGCGCTGCATCAGGTCACCATCCTGTTTTCCGATCGTGGCATTCCGGACGGCTACCGTCAcatgca                                                                                                                           <  1:574215/150‑1 (MQ=255)
                          tttATCCACACCCAGAAGCGCCACCCGCAGTCCAACCTGAAGAATGCTCAGATGATGTGGGACTTCTGGTCGCACTCTCCAGAAGCGCTGCATCAGGTCACCATCCTGTTTTCCGATCGTGGCATTCCGGACGGCTACCGTCAcatgcat                                                                                                                          >  1:831867/1‑150 (MQ=255)
                                cacaCCCAGAAGCGCCACCCGCAGTCCAACCTGAAGAATGCTCAGATGATGTGGGACTTCTGGTCGCACTCTCCAGAAGCGCTGCATCAGGTCACCATCCTGTTTTCCGATCGTGGCATTCCGGACGGCTACCGTCACATGCATGGCttt                                                                                                                    <  1:1263689/150‑1 (MQ=255)
                                cacaCCCAGAAGCGCCACCCGCAGTCCAACCTGAAGAATGCTCAGATGATGTGGGACTTCTGGTCGCACTCTCCAGAAGCGCTGCATCAGGTCACCATCCTGTTTTCCGATCGTGGCATTCCGGACGGCTACCGTCACATGCATGGCttt                                                                                                                    <  1:1334693/150‑1 (MQ=255)
                                 acacCCAGAAGCGCCACCCGCAGTCCAACCTGTAGACTGCTCTCTTCCTGTGTGTCTTCTCGTCTCACGCTCCAGAAGCGCTGCATCAGGTCACCATCCTGTTTTCCGATCGTGGCATTCCGGACGGCTACCGTCACATGCATGGCTTTg                                                                                                                   <  2:1876549/150‑1 (MQ=255)
                                 acacCCAGAAGCGCCACCCGCAGTCCAACCTGAAGAATGCTCAGATGATGTGGGACTTCTGGTCGCACTCTCCAGAAGCGCTGCATCAGGTCACCATCCTGTTTTCCGATCGTGGCATTCCGGACGGCTACCGTCACATGCATGGCTTTg                                                                                                                   >  1:1876549/1‑150 (MQ=255)
                                  cacCCAGAAGCGCCACCCGCAGTCCAACCTGAAGAATGCTCAGATGATGTGGGACTTCTGGTCGCACTCTCCAGAAGCGCTGCATCAGGTCACCATCCTGTTTTCCGATCGTGGCATTCCGGACGGCTACCGTCACATGCATGGCTTTgg                                                                                                                  <  2:1896718/150‑1 (MQ=255)
                                      cAGAAGCGCCACCCGCAGTCCAACCTGAAGAATGCTCAGATGATGTGGGACTTCTGGTCGCACTCTCCAGAAGCGCTGCATCAGGTCACCATCCTGTTTTCCGATCGTGGCATTCCGGACGGCTACCGTCACATGCATGGCTTTGGCAGc                                                                                                              <  1:733565/150‑1 (MQ=255)
                                      cAGAAGCGCCACCCGCAGTCCAACCTGAAGAATGCTCAGATGATGTGGGACTTCTGGTCGCACTCTCCAGAAGCGCTGCATCAGGTCACCATCCTGTTTTCCGATCGTGACATTCCGGACGGCTACCGTCACATGCATGGCTTTGGCAGc                                                                                                              >  1:1118367/1‑150 (MQ=255)
                                        gAAGCGCCACCCGCAGTCCAACCTGAAGAATGCTCAGATGATGTGGGACTTCTGGTCGCACTCTCCAGAAGCGCTGCATCAGGTCACCATCCTGTTTTCCGATCGTGGCATTCCGGACGGCTACCGTCACATGCATGGCTTTGGCAGCca                                                                                                            >  2:444030/1‑150 (MQ=255)
                                           gcgcCACCCGCAGTCCAACCTGAAGAATGCTCAGATGATGTGGGACTTCTGGTCGCACTCTCCAGAAGCGCTGCATCAGGTCACCATCCTGTTTTCCGATCGTGGCATTCCGGACGGCTACCGTCACATGCATGGCTTTGGCAGCCacac                                                                                                         <  1:1019102/150‑1 (MQ=255)
                                                  ccGCAGTCCAACCTGAAGAATGCTCAGATGATGTGGGACTTCTGGTCGCACTCTCCAGAAGCGCTGCATCAGGTCACCATCCTGTTTTCCGATCGTGGCATTCCGGACGGCTACCGTCACATGCATGGCTTTGGCAGCCACACCTACAGc                                                                                                  >  2:806576/1‑150 (MQ=255)
                                                      aGTCCAACCTGAAGAATGCTCAGATGATGTGGGACTTCTGGTCGCACTCTCCAGAAGCGCTGCATCAGGTCACCATCCTGTTTTCCGATCGTGGCATTCCGGACGGCTACCGTCACATGCATGGCTTTGGCAGCCACACCTACAGCCTGa                                                                                              <  1:569413/150‑1 (MQ=255)
                                                          cAACCTGAAGAATGCTCAGATGATGTGGGACTTCTGGTCGCACTCTCCAGAAGCGCTGCATCAGGTCACCATCCTGTTTTCCGATCGTGGCATTCCGGACGGCTACCGTCACATGCATGGCTTTGGCAGCCACACCTACAGCCTGATCaa                                                                                          >  2:1578601/1‑150 (MQ=255)
                                                               tGAAGAGTGCTCAGATGATGTGGGACTTCTGGTCGCACTCTCCAGAAGCGCTGCATCAGGTCACCATCCTGTTTTCCGATCGTGGCATTCCGGACGGCTACCGTCACATGCATGGCTTTGGCAGCCACACCTACAGCCTGATCAATGCaa                                                                                     >  2:1735295/1‑150 (MQ=255)
                                                                gaagaaTGCTCAGATGATGTGGGACTTCTGGTCGCACTCTCCAGAAGCGCTGCATCAGGTCACCATCCTGTTTTCCGATCGTGGCATTCCGGACGGCTACCGTCACATGCATGGCTTTGGCAGCCACACCTACAGCCTGATCAATGCaaa                                                                                    <  1:806576/150‑1 (MQ=255)
                                                                 aagaaTGCTCAGATGATGTGGGACTTCTGGTCGCACTCTCCAGAAGCGCTGCATCAGGTCACCATCCTGTTTTCCGATCGTGGCATTCCGGACGGCTACCGTCACATGCATGGCTTTGGCAGCCACACCTACAGCCTGATCAATGCAAAg                                                                                   <  1:468420/150‑1 (MQ=255)
                                                                    aaTGCTCAGATGATGTGGGACTTCTGGTCGCACTCTCCAGAAGCGCTGCATCAGGTCACCATCCTGTTTTCCGATCGTGGCATTCCGGACGGCTACCGTCACATGCATGGCTTTGGCAGCCACACCTACAGCCTGATCAATGCAAAGGGt                                                                                <  2:1513047/150‑1 (MQ=255)
                                                                       gCTCAGATGATGTGGGACTTCTGGTCGCACTCTCCAGAAGCGCTGCATCAGGTCACCATCCTGTTTTCCGATCGTGGCATTCCGGACGGCTACCGTCACATGCATGGCTTTGGCAGCCACACCTACAGCCTGATCAATGCAAAGGGTGAg                                                                             <  1:98351/150‑1 (MQ=255)
                                                                         tCAGATGATGTGGGACTTCTGGTCGCACTCTCCAGAAGCGCTGCATCAGGTCACCATCCTGTTTTCCGATCGTGGCATTCCGGACGGCTACCGTCACATGCATGGCTTTGGCAGCCACACCTACAGCCTGATCAATGCAAAGGGTGAgcg                                                                           >  1:1056424/1‑150 (MQ=255)
                                                                           aGATGATGTGGGACTTCTGGTCGCACTCTCCAGAAGCGCTGCATCAGGTCACCATCCTGTTTTCCGATCGTGGCATTCCGGACGGCTACCGt                                                                                                                                   >  1:1536180/1‑92 (MQ=255)
                                                                           aGATGATGTGGGACTTCTGGTCGCACTCTCCAGAAGCGCTGCATCAGGTCACCATCCTGTTTTCCGATCGTGGCATTCCGGACGGCTACCGt                                                                                                                                   <  2:1536180/92‑1 (MQ=255)
                                                                            gatgatGTGGGACTTCTGGTCGCACTCTCCAGAAGCGCTGCATCAGGTCACCATCCTGTTTTCCGATCGTGGCATTCCGGACGGCTACCGTCAcatgcat                                                                                                                          <  2:868260/100‑1 (MQ=255)
                                                                            gatgatGTGGGACTTCTGGTCGCACTCTCCAGAAGCGCTGCATCAGGTCACCATCCTGTTTTCCGATCGTGGCATTCCGGACGGCTACCGTCAcatgcat                                                                                                                          >  1:868260/1‑100 (MQ=255)
                                                                            gatgatGTGGGACTTCTGGTCGCACTCTCCAGAAGCGCTGCATCAGGTCACCATCCTGTTTTCCGATCGTGGCATTCCGGACGGCTACCGTCACATGCATTGCTTTGGCAGCCACACCTCCAGCCTGATCAATGCAACGGGTGAgcgccc                                                                        >  2:838188/1‑148 (MQ=255)
                                                                             atgatgTGGGACTTCTGGTCGCACTCTCCAGAAGCGCTGCATCAGGTCACCATCCTGTTTTCCGATCGTGGCATTCCGGACGGCTACCGTCACATGCATGGCTTTGGCAGCCACACCTACAGCCTGATCAATGCAAAGGGTGAGCGCACt                                                                       <  1:1896596/150‑1 (MQ=255)
                                                                              tgatgTGGGACTTCTGGTCGCACTCTCCAGAAGCGCTGCATCAGCTAACCATCTTGTTTTCCGATCGTGGCATTCCGGACGGCTACCGTCACATGCATGGCTTTGGCAGCCACACCTACAGCCTGATCAATGCAAAGGGTGAGCGCACtt                                                                      <  1:1100074/150‑1 (MQ=255)
                                                                               gatgTGGGACTTCTGGTCGCACTCTCCAGAAGCGCTGCATCAGGTCACCATCCTGTTTTCCGATCGTGGCATTCCGGACGGCTACCGTCACATGCATGGCTTTGGCAGCCACACCTACAGCCTGATCAATGCAAAGGGTGAGCGCACTTg                                                                     <  2:831867/150‑1 (MQ=255)
                                                                                atgTGGGACTTCTGGTCGCACTCTCCAGAAGCGCTGCATCAGGTCACCATCCTGTTTTCCGATCGTGGCATTCCGGACGGCTACCGTCACATGCATGGCTTTGGCAGCCACACCTACAGCCTGATCAATGCAAAGGGTGAGCGCACTTgg                                                                    <  2:1817870/150‑1 (MQ=255)
                                                                                 tgtgGGACTTCTGGTCGCACTCTCCAGAAGCGCTGCATCAGGTCACCATCCTGTTTTCCGATCGTGGCATTCCGGACGGCTACCGTCACATGCATGGCTTTGGCAGCCACACCTACAGCCTGATCAATGCAAAGGGTGAGCGCACTTggg                                                                   <  2:22320/150‑1 (MQ=255)
                                                                                 tgtgGGACTTCTGGTCGCACTCTCCAGAAGCGCTGCATCAGGTCACCATCCTGTTTTCCGATCGTGGCATTCCGGACGGCTACCGTCACATGCATGGCTTTGGCAGCCACACCTACAGCCTGATCAATGCAAAGGGTGAGCGCACTTggg                                                                   >  2:331499/1‑150 (MQ=255)
                                                                                   tgGGACTTCTGGTCGCACTCTCCAGAAGCGCTGCATCAGGTCACCATCCTGTTTTCCGATCGTGGCATTCCGGACGGCTACCGTCACATGCATGGCTTTGGCAGCCACACCTACAGCCTGAt                                                                                             <  1:867427/122‑1 (MQ=255)
                                                                                   tgGGACTTCTGGTCGCACTCTCCAGAAGCGCTGCATCAGGTCACCATCCTGTTTTCCGATCGTGGCATTCCGGACGGCTACCGTCACATGCATGGCTTTGGCAGCCACACCTACAGCCTGAt                                                                                             >  2:867427/1‑122 (MQ=255)
                                                                                    gggACTTCTGGTCGCACTCTCCAGAGGCGCTGCATCTGGTCACCATCCTGGTGTCCGAGCGTGGCATTCGGGACGGCTACCGTCACATGCATGGCTTTGGCAGCCACACCTACAGCCTGATCAATGCAAAGGGTGAGCGCACTTGGGTCa                                                                >  1:1106424/1‑150 (MQ=255)
                                                                                       aCTTCTGGTCGCACTCTCGAGAAGCGCTGCATCAGGTCACCATCCTGTTTTCCGATCGTGGCATTCCGGACGGCTACCGTCACATGCATGGCTTTGGCAGCCACACCTACAGCCTGATCAATGCAAAGGGTGAGCGCACTTGGGTCAAGt                                                             >  2:1142875/1‑150 (MQ=255)
                                                                                         ttCTGGTCGCACTCTCCAGAAGCGCTGCATCAGGTCACCATCCTGTTTTCCGATCGTGGCATTCCGGACGGCTACCGTCACATGCATGGCTTTg                                                                                                                   <  1:88058/94‑1 (MQ=255)
                                                                                         ttCTGGTCGCACTCTCCAGAAGCGCTGCATCAGGTCACCATCCTGTTTTCCGATCGTGGCATTCCGGACGGCTACCGTCACATGCATGGCTTTg                                                                                                                   >  2:88058/1‑94 (MQ=255)
                                                                                              gTCGCACTCTCCAGAAGCGCTGCACCAGGTCACCATCCTGTTTTCCGATCGTGGCATTCCGGACCACTACCGTCACATGCATGGCTTTGGCAGCCACACCTACAGCCTGATCAATGCAAAGGGTGAGCGCACTTGGGTCAAGTGTCACtt                                                      >  2:950495/1‑150 (MQ=255)
                                                                                                 gCACTCTCCAGAAGCGCTGCATCAGGTCACCATCCTGTTTTCCGATCGTGGCATTCCGGACGGCTACCGTCACATGCATGGCTTTGGCa                                                                                                                <  2:1361535/89‑1 (MQ=255)
                                                                                                 gCACTCTCCAGAAGCGCTGCATCAGGTCACCATCCTGTTTTCCGATCGTGGCATTCCGGACGGCTACCGTCACATGCATGGCTTTGGCa                                                                                                                >  1:1361535/1‑89 (MQ=255)
                                                                                                  cACTCTCCAGAAGCGCTGCATCAGGTCACCATCCTGTTTTCCGATCGTGGCATTCCGGACGGCTACCGTCACATGCATGGCTTTGGCAGCCACACCTACAGCCTGATCAATGCAAAGGGTGAGCGCACTTGGGTCAAGTGGCACTTCAAg                                                  <  2:1056424/150‑1 (MQ=255)
                                                                                                       tcCAGAAGCGCTGCATCAGGTCACCATCCTGTTTTCCGATCGTGGCATTCCGGACGGCTACCGTCACATGCATGGCTTTGGCAGCCACACCTACAGCCTGATCAATGCAAAGGGTGAGCGCACTTGGGTCAAGTGGCACTTCAAGACCCa                                             >  2:257523/1‑150 (MQ=255)
                                                                                                       tcCAGAAGCGCTGCATCAGGTCACCATCCTGTTTTCCGATCGTGGCATTCCGGACGGCTACCGTCACATGCATGGCTTTGGCAGCCACACCTACAGCCTGATCAATGCAAAGGGTGAGCGCACTTGGGTCAAGTGGCACTTCAAGACCCa                                             >  2:1633549/1‑150 (MQ=255)
                                                                                                         cAGAAGCGCTGCATCAGGTCACCATCCTGTTTTCCGATCGTGGCATTCCGGACGGCTACCGTCACATGCATGGCTTTGGCAGCCACACCTACa                                                                                                    <  1:1519972/93‑1 (MQ=255)
                                                                                                         cAGAAGCGCTGCATCAGGTCACCATCCTGTTTTCCGATCGTGGCATTCCGGACGGCTACCGTCACATGCATGGCTTTGGCAGCCACACCTACa                                                                                                    >  2:1519972/1‑93 (MQ=255)
                                                                                                           gAAGCGCTGCATCAGGTCACCATCCTGTTTTCCGATCGTGGCATTCCGGACGGCTACCGTCACATGCATGGCTTTGGCAGCCACACCTACAGCCTGATCAATGCAAAGGGTGAGCGCACTTGGGt                                                                  <  1:431217/125‑1 (MQ=255)
                                                                                                           gAAGCGCTGCATCAGGTCACCATCCTGTTTTCCGATCGTGGCATTCCGGACGGCTACCGTCACATGCATGGCTTTGGCAGCCACACCTACAGCCTGATCAATGCAAAGGGTGAGCGCACTTGGGt                                                                  >  2:431217/1‑125 (MQ=255)
                                                                                                              gcgcTGCATCAGGTCACCATCCTGTTTTCCGATCGTGGCATTCCGGACGGCTACCTTCACATGCATGGCTTTGGCAGCCACACCTACAGCCTGATCAATGCAAAgg                                                                                  <  1:1220325/106‑1 (MQ=255)
                                                                                                              gcgcTGCATCAGGTCACCATCCTGTTTTCCGATCGTGGCATTCCGGACGGCTACCTTCACATGCATGGCTTTGGCAGCCACACCTACAGCCTGATCAATGCAAAgg                                                                                  >  2:1220325/1‑106 (MQ=255)
                                                                                                                gcTGCATCAGGTCACCATCCTGTTTTCCGATCGTGGCATTCCGGACGGCTACCGTCACATGCATGGCTTTGGCAGCCACACCTACAGCCTGATCAATGCAAAGGGTGAGCGCACTTGGGTCAAGTGGCACTTCAAGACCCAGCAAGGCAt                                    >  1:1461887/1‑150 (MQ=255)
                                                                                                                gcTGCATCAGGTCACCATCCTGTTTTCCGATCGTGGCATCCCGGACGGCTACTGTACCATGAATGGCTGTGGGAGCCACTGCTACAGCCTTATCAATGCAAAGGGTTATCGCACTTGGGTCAAGTGGCACTTCAAGACCCAGCAAGGCAt                                    <  2:1189713/150‑1 (MQ=255)
                                                                                                                gcTGCATCAGGTCACCATCCTGTTTTCCGATCGTGACATTCCGGACGGCTACCGTCACATGCATGGCTTTGGCAGCCACACCTACAGCCTGATCAATGCAAAGGGTGAGCGCACTTGGGTCAAGTGGCACTTc                                                     <  1:27728/133‑1 (MQ=255)
                                                                                                                gcTGCATCAGGTCACCATCCTGTTTTCCGATCGTGACATTCCGGACGGCTACCGTCACATGCATGGCTTTGGCAGCCACACCTACAGCCTGATCAATGCAAAGGGTGAGCGCACTTGGGTCAAGTGGCACTTc                                                     >  2:27728/1‑133 (MQ=255)
                                                                                                                 cTGCATCAGGTCACCATCCTGTTTTCCGATCGTGGCATTCCGGACGGCTACCGTCACATGCATGGCTTTGGCAGCCACACCTACAGc                                                                                                  <  1:453698/87‑1 (MQ=255)
                                                                                                                 cTGCATCAGGTCACCATCCTGTTTTCCGATCGTGGCATTCCGGACGGCTACCGTCACATGCATGGCTTTGGCAGCCACACCTACAGc                                                                                                  >  2:453698/1‑87 (MQ=255)
                                                                                                                    cATCAGGTCACCATCCTGTTTTCCGATCGTGGCATTCCGGACGGCTACCGTCACATGCATGGCTTTGGCAGCCACACCTACAGCCTGATCAATGCAAAGGGTGAGCGCACTTGGGTCAAGTGGCACTTCAAGACCCAGCAAGGCATCAAg                                <  1:1633549/150‑1 (MQ=255)
                                                                                                                     aTCAGGTCACCATCCTGTTTTCCGATCGTGGCATTCCGGACGGCTACCGTCACATGCATGGCTTTGGCAGCCACACCTACAGCCTGATCAATGCAAAGGGTGAGCGCACTTGGGTCAAGTGGCACTTCAAGACCCAGCAAGGCATCAAGa                               >  2:1408733/1‑150 (MQ=255)
                                                                                                                      tCAGGTCACCATCCTGTTTTCCGATCGTGGCATTCCGGACGGCTACCGTCACATGCATGGCTTTGGCAGCCACACCTACAGCCTGATCAATGCAAAGGGTGAGCGCACTTGGGTCAAGTGGCACTTCAAGACCCAGCAAGGCATCAAGaa                              <  2:535270/150‑1 (MQ=255)
                                                                                                                          gTCACCATCCTGTTTTCCGATCGTGGCATTCCGGACGGCTACCGTCACATGCATGGCTTTGGCAGCCACACCTACAGCCTGATCAATGCAAAGGGTGAGCGCACTTGGGTCAAGTGGCACTTCAAGACCCAGCAAGGCATCAAGAACCTc                          >  2:1721311/1‑150 (MQ=255)
                                                                                                                          gTCACCATCCTGTTTTCCGATCGTGGCATTCCGGACGGCTACCGTCACATGCATGGCTTTGGCAGCCACACCTACAGCCTGATCAATGCAAAGGGTGAGCGCACTTGGGTCAAGTGGCACTGCAAGTCCCAGCAAGGCATCAAGAACCTc                          >  2:1772363/1‑150 (MQ=255)
                                                                                                                           tCACCATCCTGTTTTCCGATCGTGGCATTCCGGACGGCTACCGTCACATGCATGGCTTTGGCAGCCACACCTACAGCCTGATCAATGCAAAGGGTGAGCGCACTTGGGTCAAGTGGCACTTCAAGACCCAGCAAGGCATCAAGAACCTCa                         <  1:257523/150‑1 (MQ=255)
                                                                                                                            cACCATCCTGTTTTCCGATCGTGGCATTCCGGACGGCTACCGTCACATGCATGGCTTTGGCAGCCACACCTACAGCCTGATCAATGCAAAGGGTGAGCGCACTTGGGTCAAGTGGCACTTCAAGACCCAGCAAGGCATCAAGAACCTCAc                        <  2:302569/150‑1 (MQ=255)
                                                                                                                               cATCCTGTTTTCCGATCGTGGCATTCCGGACGGCTACCGTCACATGCATGGCTTTGGCAGCCACACCTACAGCCTGATCAATGCAAAGGGTGAGCGCACTTGGGTCAAGTGGCACTTCAAGa                                                 >  2:169971/1‑122 (MQ=255)
                                                                                                                               cATCCTGTTTTCCGATCGTGGCATTCCGGACGGCTACCGTCACATGCATGGCTTTGGCAGCCACACCTACAGCCTGATCAATGCAAAGGGTGAGCGCACTTGGGTCAAGTGGCACTTCAAGa                                                 <  1:169971/122‑1 (MQ=255)
                                                                                                                               cATCCTGTTTTCCGATCGTGGCATTCCGGACGGCTACCGTCACATGCATGGCTTTGGCAGCCACACCTACAGCCTGATCAATGCAAAGGGTGAGCGCACTTGGGTCAAGTGGCACTTCAAGACCCAGCAAGGCATCAAGaa                              >  1:348320/1‑141 (MQ=255)
                                                                                                                               cATCCTGTTTTCCGATCGTGGCATTCCGGACGGCTACCGTCACATGCATGGCTTTGGCAGCCACACCTACAGCCTGATCAATGCAAAGGGTGAGCGCACTTGGGTCAAGTGGCACTTCAAGACCCAGCAAGGCATCAAGaa                              <  2:348320/141‑1 (MQ=255)
                                                                                                                               cATCCTGTTTTCCGATCGTGGCATTCCGGACGGCTACCGTCACATGCATGGCTTTGGCAGCCACACCTACAGCCTGATCAATGCAAAGGGTGAGCGCACTTGGGTCAAGTGGCACTTCAAGACCCAGCAAGGCATCAAGAACCTCAcccc                     >  2:414040/1‑150 (MQ=255)
                                                                                                                                   cTGTTTTCCGATCGTGGCATTCCGGACGGCTACCGTCACATGCATGGCTTTGGCAGCCACACCTACAGCCTGATCAATg                                                                                        <  2:860907/79‑1 (MQ=255)
                                                                                                                                   cTGTTTTCCGATCGTGGCATTCCGGACGGCTACCGTCACATGCATGGCTTTGGCAGCCACACCTACAGCCTGATCAATg                                                                                        >  1:860907/1‑79 (MQ=255)
                                                                                                                                         tCCGATCGTGGCATTCCGGACGGCTACCGTCACATGCATGGCTTTGGCAGCCACACCTACAGCCTGATCAATGCAAAGGGTGAGCGCACTTGGGTCAAGTGGCACTTCAAGACCCAGCAAGGCATCAAGAACCTCACCCCGGCAGACGCt           <  1:1578601/150‑1 (MQ=255)
                                                                                                                                           cGATCGTGGCATTCCGGACGGCTACCGTCACATGCATGGCTTTGGCAGCCACACCTACAGCCTGATCAATGCAAAGGGTGAGCGCACTTGGGTCAAGTGGCACTTCAAGACCCAGCAAGGCATCAAGAACCTCACCCCGGCAGACGCTGc         >  2:885572/1‑150 (MQ=255)
                                                                                                                                                 tGGCATTCCGGACGGCTACCGTCACATGCATGGCTTTGGCAGCCACACCTACAGCCTGATCAATGCAAAGGGTGAGCGCACTTGGGTCAAGTGGCACTTCAAGACCCAGCAAGGCATCAAGAACCTCACCCCGGCAGACGCTGCACGCCt   >  1:506623/1‑150 (MQ=255)
                                                                                                                                                  ggCATTCCGGACGGCTCCCGTCACATGCATGGCTTTGGCAGCCACACCTACAGCCTGATCAATGCAAAGGGTGAGCGCACTTGGGTCAAGTGGCACTTCAAGACCCAGCAAGGCATCAAGAACCTCACCCCGGCAGACGCTGCACGCCTg  <  1:885572/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
TTCGCGACCCACTTAAGTTCCCGGACTTTATCCACACCCAGAAGCGCCACCCGCAGTCCAACCTGAAGAATGCTCAGATGATGTGGGACTTCTGGTCGCACTCTCCAGAAGCGCTGCATCAGGTCACCATCCTGTTTTCCGATCGTGGCATTCCGGACGGCTACCGTCACATGCATGGCTTTGGCAGCCACACCTACAGCCTGATCAATGCAAAGGGTGAGCGCACTTGGGTCAAGTGGCACTTCAAGACCCAGCAAGGCATCAAGAACCTCACCCCGGCAGACGCTGCACGCCTG  >  NC_002947/564921‑565216

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: