Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,623,312 A→T 5.3% H159Q (CAT→CAA)  PP_4941 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,623,3120AT 5.3% 292.7 / 6.6 95H159Q (CAT→CAA) PP_4941hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (46/44);  new base T (2/3);  total (48/47)
Fisher's exact test for biased strand distribution p-value = 6.77e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.41e-01

GTCACCACCCAGCTTTCGTCGGCCCGGCTGTACAAACCATACGCGGCCAACTCGGCGTCATCGGCATGTGGCGCGATCACCATCACCCGGCGCCGCTGCAGCTCGACCGAGGGGGTTATCCACAGGCGCAACTGACCCTGCAGCCGGCAATGCCGGCCGCGCAGTCGCAACGTACCGGCCTGCAGTGGCGCAGCAAGACCCGTCAGGTTAAGGAAGCGCACGCCATCGACGC‑CACG‑CT‑CGAACGCTTGCCGATCGGCCTGTTCGCCAGCCAGCAGCTCGACCCGAGGGTCGAGA  >  NC_002947/5623163‑5623456
                                                                                                                                                     |                                                                                                                                                   
gTCACCACCCAGCTTTCGTCGGCCCGGCTGTACAAACCATACGCGGCCAACTCGGCGTCATCGGCATGTGGCGCGATCACCATCACCCGGCGCCGCTGCAGCGCGACCGAGGGGGTTATCCACAGGCGCAACTGACCCTGCAGCCGGCaa                                                                                                                                                     <  1:49555/150‑1 (MQ=255)
 ttaccaccCAGCTTTCGTCGGCACGGCTGTACAAACCATACGCGGCCAACTCGGCGTCATCGGCATGTGGCGCGATCACCATCACCCGGCGCCGCTGCAGCTCGACCGAGGGGGTTATCCACAGGCGCAACTGACCCTGCAGCCGGCAAt                                                                                                                                                    <  2:190879/148‑1 (MQ=255)
  caccacCCAGCTTTCGTCGGCCCGGCTGTACAAACCATACGCGGCCAACTCGGCGTCATCGGCATGTGGCGCGATCACCATCACCCGGCGCCGCTGCAGCTCGACCGAGGGGGTTATCCACAGGCGCAACTGACCCTGCAGCCGGCAATg                                                                                                                                                   <  2:1891849/150‑1 (MQ=255)
        ccAGCTTTCGTCGGCCCGGCTGTACAAACCATACGCGGCCAACTCGGCGTCATCGGCATGTGGCGCGATCACCATCACCCGGCGCCGCTGCAGCTCGACCGAGGGGGTTATCCACAGGCGCAACTGACCCTGCAGCCGGCAATgccggcc                                                                                                                                             >  2:671492/1‑150 (MQ=255)
        ccAGCTTTCGTCGGCCCGGCTGTACAAACCATACGCGGCCAACTCGGCGTCATCGGCATGTGGCGCGATCACCATCACCCGGCGCCGCTGCAGCTCGACCGAGGGGGTTATCCACAGGCGCAACTGACCCTGCAGCCGGCAATgccggcc                                                                                                                                             >  2:652141/1‑150 (MQ=255)
          aGCTTTCGTCGGCCCGGCTGTACAAACCATACGCGGCCAACTCGGCGTCATCGGCATGTGGCGCGATCACCATCACCCGGCGCCGCTGCAGCTCGACCGAGGGGGTTATCCACAGGCGCAACTGACCCTGCAGCCGGCAATGCCGGCcgc                                                                                                                                           <  2:1601130/150‑1 (MQ=255)
           gCTTTCGTCGGCCCGGCTGTACAAACCATACGCGGCCAACTCGGCGTCATCGGCATGTGGCGCGATCACCATCACCCGGCGCCGCTGCAGCTCGACCGAGGGGGTTATCCACAGGCGCAACTGACCCTGCGGCCGGCATTGCCGGCtgca                                                                                                                                          <  2:41020/150‑2 (MQ=255)
             tttCGTCGGCCCGGCTGTACAAACCATACGCGGCCAACTCGGCGTCATCGGCATGTGGCGCGATCACCATCACCCGGCGCCGCTGCAGCTCGACCGAGGGGGTTATCCACAGGCGCAACTGACCCTGCAGCCGGCAATGCCGGCCGCGCa                                                                                                                                        <  1:1652319/150‑1 (MQ=255)
                  tcgGTCCGGCTGTACAAACCATACGCGGCCAACTCGGCGTCATCGGCATGTGGCGCGATCACCATCACCCGGCGCCGCTGCAGCTCGACCGAGGGGGTTATCCACAGGCGCCAAAGTCCCTGCAGCCGGCAATGCCGGCCGCGCAGTCGc                                                                                                                                   <  1:647229/150‑1 (MQ=255)
                  tcgGCCCCGCTGTACAAACCATACCCCGCCAACTCGGCGTCATCGGCATGTGGCGCGATCACCATCACCCGGCGCCGCTGCAGCTCGACCGAGGGGGTTATCCACAGCCGCAACTGACCCTTCAGCCGGCAATGCCGGCCGCGCAGTCGc                                                                                                                                   <  2:147003/150‑1 (MQ=255)
                     gCCCGGCTGTACAAACCATACGCGGCCAACTCGGCGTCATCGGCATGTGGCGCGATCACCATCACCCGGCGCCGCTGCAGCTCGACCGAGGGGGTTATCCACAGGCGCAACTGACCCTGCAGCCGGCAATGCCGGCCGCGCAGTCGCAAc                                                                                                                                >  1:1784661/1‑150 (MQ=255)
                       ccGGCTGTACAAACCATACGCGGCCAACTCGGCGTCATCGGCATGTGGCGCGATCACCATCACCCGGCGCCGCTGCAGCTCGACCGAGGGGGTTATCCACAGGCGCAACTGACCCTGCAGCCGGCAATgccg                                                                                                                                                <  1:1215478/132‑1 (MQ=255)
                       ccGGCTGTACAAACCATACGCGGCCAACTCGGCGTCATCGGCATGTGGCGCGATCACCATCACCCGGCGCCGCTGCAGCTCGACCGAGGGGGTTATCCACAGGCGCAACTGACCCTGCAGCCGGCAATgccg                                                                                                                                                >  2:1215478/1‑132 (MQ=255)
                       ccGGCTGTACAAACCATACGCGGCCAACTCGGCGTCATCGGCATGTGGCGCGATCACCATCACCCGGCGCCGCTGCAGCTCGACCGAGGGGGTTATCCACAGGCGCAACTGACCCTGCAGCCGGCAATGCCGGCCGCGCAGTCGCAACGt                                                                                                                              >  2:1568342/1‑150 (MQ=255)
                       ccGGCTGTACAAACCATACGCGGCCAACTCGGCGTCATCGGCATGTGGCGCGATCACCATCACCCGGCGCCGCTGCAGCTCGACCGAGGGGGTTATCCACAGGCGCAACTGACCCTGCAGCCGGCAATGCCGGCCGCGCAGTCGCAACGt                                                                                                                              >  2:1137348/1‑150 (MQ=255)
                         ggCTGTACAAACCATACGCGGCCAACTCGGCGTCATCGGCATGTGGCGCGATCACCATCACCCGGCGCCGCTGCAGCTCGACCGAGGGGGTTATCCACAGGCGCAACTGACCCTGCAGCCGGCAATGCCGGCCGCGCAGTCGCAACGTAc                                                                                                                            >  2:730360/1‑150 (MQ=255)
                           cTGTACAAACCATACGCGGCCAACTCGGCGTCATCGGCATGTGGCGCGATCACCATCACCCGGCGCCGCTGCAGCTCGACCGAGGGGGTTATCCACAGGCGCAACTGACCCTGCAGCCGGCAATGCCGGCCGCGCAGTCGCAACGTACCg                                                                                                                          <  1:730360/150‑1 (MQ=255)
                           cTGTACAAACCATACGCGGCCAACTCGGCGTCATCGGCATGTGGCGCGATCACCATCACCCGGCGCCGCTGCAGCTCGACCGAGGGGGTTATCCACAGGCGCAACTGACCCTGCAGCCGGCAATGCCGGCCGCGCAGTCGCAACGTACCg                                                                                                                          <  1:1137348/150‑1 (MQ=255)
                             gTACAAACCATACGCGGCCAACTCGGCGTCATCGGCATGTGGCGCGATCACCATCACCCGGCGCCGCTGCAGCTCGACCGAGGGGGTTATCCACAGGCGCAACTGACCCTGCAGCCGGCAATg                                                                                                                                                   <  1:64873/123‑1 (MQ=255)
                             gTACAAACCATACGCGGCCAACTCGGCGTCATCGGCATGTGGCGCGATCACCATCACCCGGCGCCGCTGCAGCTCGACCGAGGGGGTTATCCACAGGCGCAACTGACCCTGCAGCCGGCAATg                                                                                                                                                   >  2:64873/1‑123 (MQ=255)
                             gTACAAACCATACGCGGCCAACTCGGCGTCATCGGCATGTGGCGCGATCACCATCACCCGGCGCCGCTGCAGCTCGACCGAGGGGGTTATCCACAGGCGCAACTGACCCTGCAGCCGGCAATGCCGGCCGCGCAGTCGCAACGTACCGGc                                                                                                                        <  1:1351097/150‑1 (MQ=255)
                                cAAACCATACGCGGCCAACTCGGCGTCATCGGCATGTGGCGCGATCACCATCACCCGGCGCCGCTGCAGCTCGACCGAGGGGGTTATCCACAGGCGCAACTGACCCTGCAGCCGGCAATGCCGGCCGCGCAGTCGCAACGTACCGGCCTg                                                                                                                     <  2:422192/150‑1 (MQ=255)
                                 aaaCCATACGCGGCCAACTCGGCGTCATCGGCATGTGGCGCGATCTCCATCACCCGGCGCCGCTGCAGCTCGACCGAGGGGGTTATCCACAGGCGCAACTGACCCTGCAGCCGGCAATGCCGGCCGCGCAGTCGCAACGTACCGGCCTGc                                                                                                                    >  1:1427262/1‑150 (MQ=255)
                                 aaaCCATACGCGGCCAACTCGGCGTCATCGGCATGTGGCGCGATCACCATCACCCGGCGCCGCTGCAGCTCGACCGAGGGGGTTATCCACAGGCGCAACTGACCCTGCAGCCGGCAATGCCGGCCGCGCAGTCGCAACGTACCGGCCTGc                                                                                                                    >  2:1655858/1‑150 (MQ=255)
                                 aaaCCATACGCGGCCAACTCGGCGTCATCGGCATGTGGCGCGATCACCATCACCCGGCGCCGCTGCAGCTCGACCGAGGGGGTTATCCACAGGCGCAACTGACCCTGCAGCCGGCAATGCCGGCCGCGCAGTCGCAACGTACCGGCCTGc                                                                                                                    >  1:1427406/1‑150 (MQ=255)
                                 aaaCCATACGCGGCCAACTCGGCGTCATCGGCATGTGGCGCGATCACCATCACCCGGCGCCGCTGCAGCTCGACCGAGGGGGTTATCCACAGGCGCAACTGACCCTGCAGCCGGCAATGCCGGCCGCGCAGTCGCAACGTACCGGCCTGc                                                                                                                    >  2:356299/1‑150 (MQ=255)
                                  aaCCATACGCGGCCAACTCGGCGTCATCGGCATGTGGCGCGATCACCATCACCCGGCGCCGCTGCAGCTCGACCGAGGGGGTTATCCACAGGCGCAACTGACCCTGCAGCCGGCAATGCCGGCCGCGCAGTCGCAACGTACCGGCCTGCa                                                                                                                   >  1:1591679/1‑150 (MQ=255)
                                   aCCATACGCGGCCAACTCGGCGTCATCGGCATGTGGCGCGATCACCATCACCCGGCGCCGCTGCAGCTCGACCGAGGGGGTTATCCACAGGCGCAACTGACCCTGCAGCCGGCAATGCCGGCCGCGCAGTCGCAACGTACCGGCCTGCAg                                                                                                                  >  1:1706631/1‑150 (MQ=255)
                                     cATACGCGGCCAACTCGGCGTCATCGGCATGTGGCGCGATCACCATCACCCGGCGCCGCTGCAGCTCGACCGAGGGGGTTATCCACAGGCGCAACTGACCCTGCAGCCGGCAATGCCGGCCGCGCAGTCGCAACGTACCGGCCTGCAGTg                                                                                                                >  2:998633/1‑150 (MQ=255)
                                      aTACGCGGCCAACTCGGCGTCATCGGCATGTGGCGCGATCACCATCACCCGGCGCCGCTGCAGCTCGACCGAGGGGGTTATCCACAGGCGCAACTGACCCTGCAGCCGGCAATgcc                                                                                                                                                 <  2:300330/116‑1 (MQ=255)
                                      aTACGCGGCCAACTCGGCGTCATCGGCATGTGGCGCGATCACCATCACCCGGCGCCGCTGCAGCTCGACCGAGGGGGTTATCCACAGGCGCAACTGACCCTGCAGCCGGCAATgcc                                                                                                                                                 >  1:300330/1‑116 (MQ=255)
                                      aTACGCGGCCAACTCGGCGTCATCGGCATGTGGCGCGATCACCATCACCCGGCGCCGCTGCAGCTCGACCGAGGGGGTTATCCACAGGCGCAACTGACCCTGCAGCCGGCAATGCCGGCCGCGCAGTCGCAACGTACCGGcc                                                                                                                       <  2:1604011/142‑1 (MQ=255)
                                      aTACGCGGCCAACTCGGCGTCATCGGCATGTGGCGCGATCACCATCACCCGGCGCCGCTGCAGCTCGACCGAGGGGGTTATCCACAGGCGCAACTGACCCTGCAGCCGGCAATGCCGGCCGCGCAGTCGCAACGTACCGGcc                                                                                                                       >  1:1604011/1‑142 (MQ=255)
                                            ggCCAACTCGGCGTCATCGGCATGTGGCGCGATCACCATCACCCGGCGCCGCTGCAGCTCGACCGAGGGGGTTATCCACAGGCGCAACTGACCCTGCAGCCGGCAATGCCGGCCGCGCAGTCGCAACGTACCGGCCTGCAGTGGCgcagc                                                                                                         >  1:874649/1‑150 (MQ=255)
                                                accTCGGCGTCATCGGCATGTGGCGCGATCACCATCACCCGGCGCCGCTGCAGCTCGACCGAGGGGCTTATCCACAGGCGCAACTGACCCTGCAGCCGGCAATGCCGGCCGCGCAGTCGCAACGTACCGGCCTGCAGTGGCGCAGCAAGa                                                                                                     <  1:207735/148‑1 (MQ=255)
                                                aaCTCGGCGTCATCGGCATGTGGCGCGATCACCATCACCCGGCGCCGCTGCAGCTCGACCGAGGGGGTTATCCACAGGCGCAACTGACCCTGCAGCCGGCAATGCCGGCCGCGCAGTCGCAACGTACCGGCCTGCAGTGGCGCAGCAAGa                                                                                                     >  1:994246/1‑150 (MQ=255)
                                                aaCTCGGCGTCATCGGCATGTGGCGCGATCACCATCACCCGGCGCCGCCGCAGCTCGACCGAGGGGGTTATCCACAGGCGCAACTGACCCTGCGGCCGGCATTgccggctgcag                                                                                                                                         <  2:772194/114‑6 (MQ=255)
                                                aaCTCGGCGTCATCGGCATGTGGCGCGATCACCATCACCCGGCGCCGCCGCAGCTCGACCGAGGGGGTTATCCACAGGCGCAACTGACCCTGCGGCCGGCATTgccggctgcag                                                                                                                                         >  1:772194/1‑109 (MQ=255)
                                                 aCTCGGCGTCATCGGCATGTGGCGCGATCACCATCACCCGGCGCCGCTGCAGCTCGACCGAGGGGGTTATCCACAGGCGCAACTGACCCTGCAGCCGGCAATgc                                                                                                                                                  <  1:1258218/104‑1 (MQ=255)
                                                 aCTCGGCGTCATCGGCATGTGGCGCGATCACCATCACCCGGCGCCGCTGCAGCTCGACCGAGGGGGTTATCCACAGGCGCAACTGACCCTGCAGCCGGCAATgc                                                                                                                                                  >  2:1258218/1‑104 (MQ=255)
                                                   tCGGCGTCATCGGCATGTGGCGCCATCACCATCACCCGGCGCCGCTGCAGCTCGACCGAGGGGGTTATCCACAGGCGCAACTGACCCTGCAGCCGGCAATGCCGGCCGCGCAGTCGCAACGTACCGGCCTGCAGTGGCGCAGCAAGAccc                                                                                                  <  1:392799/150‑1 (MQ=255)
                                                        gTCATCGGCATGTGGCGCGATCACCATCACCCGGCGCCGCTGCAGCTCGACCGAGGGGGTTATCCACAGGCGCAACTGACCCTGCAGCCGGCAATGCCGGCCGCGCAGTCGCAACGTACCGGCCTGCAGTgg                                                                                                               >  1:1048757/1‑132 (MQ=255)
                                                        gTCATCGGCATGTGGCGCGATCACCATCACCCGGCGCCGCTGCAGCTCGACCGAGGGGGTTATCCACAGGCGCAACTGACCCTGCAGCCGGCAATGCCGGCCGCGCAGTCGCAACGTACCGGCCTGCAGTgg                                                                                                               <  2:1048757/132‑1 (MQ=255)
                                                          cATCGGCATGTGGCGCGATCACCATCACCCGGCGCCGCTGCAGCTCGACCGAGGGGGTTATCCACAGGCGCAACTGACCCTGCAGCCGGCAATGCCGGCCGCGCAGTCGCAACGTACCGGCCTGCAGTGGCGCAGCAAGACCCGTCAGGt                                                                                           >  2:1243046/1‑150 (MQ=255)
                                                             cGGCATGTGGCGCGATCACCATCACCCGGCGCCGCTGCAGCTCGACCGAGGGGGTTATCCACAGGCGCAACTGACCCTGCAGCCGGCAATGCCGGCCGCGCAGTCGCAACGTACCGGCCTGCAGTGGCGCAGCAAGACCCGTCAGGTTaa                                                                                        <  2:1706631/150‑1 (MQ=255)
                                                             cGGCATGTGGCGCGATCACCATCACCCGGCGCCGCTGCAGCTCGACCGAGGGGGTTATCCACAGGCGCAACTGACCCTGCAGCCGGCAATGCCGGCCGCGCAGTCGCAACGTACCGGCCTGCAGTGGCGCAGCAAGACCCGTCAGGTTaa                                                                                        <  1:998633/150‑1 (MQ=255)
                                                             cGGCATGTGGCGCGATCACCATCACCCGGCGCCGCTGCAGCTCGACCGAGGGGGTTATCCACAGGCGCAACTGACCCTGCAGCCGGCAATGCCGGCCGCGCAGTCGCAACGTACCGGCCTGCAGTGGCGCAGCAAGACCCGTCAGGTTaa                                                                                        <  2:1591679/150‑1 (MQ=255)
                                                                  tgtgGCGCGATCACCATCACCCGGCGCCGCTGCAGCTCGACCGAGGGGGTTATCCACAGGCGCAACTGACCCTGCAGCCGGCAATGCCGGCCGCGCAGTCGCAACGTACCGGCCTGCAGTGGCGCAGCAAGACCCGTCAGGTTAAGGAAg                                                                                   <  2:1784661/150‑1 (MQ=255)
                                                                      gcgcGATCACCATCACCCGGCGCCGCTGCAGCTCGACCGAGGGGGTTATCCACAGGCGCAACTGACCCTGCAGCCGGCAATGCCGGCCGCGCAGTCGCAACGTACCGGCCTGCAGTGGCGCAGCAAGACCCGTCAGGTTAAGGAAGCGCa                                                                               >  2:1536193/1‑150 (MQ=255)
                                                                          gATCACCATCACCCGGCGCCGCTGCAGCTCGACCGAGGGGGTTATCCACAGGCGCAACTGACCCTGCAGCCGGCAATgcc                                                                                                                                                 >  1:747243/1‑80 (MQ=255)
                                                                          gATCACCATCACCCGGCGCCGCTGCAGCTCGACCGAGGGGGTTATCCACAGGCGCAACTGACCCTGCAGCCGGCAATgcc                                                                                                                                                 <  2:747243/80‑1 (MQ=255)
                                                                            tcaccatcaccCGGCGCCGCTGCAGCTCGACCGAGGGGGTTATCCACAGGCGCAACTGACCCTGCAGCCGGCAATGCCGGCCGCGCAGTCGCAACGTACCGGCCTGCAGTGGCGCAGCAAGACCCGTCAGGTTAAGGAAGCGCACGCCAt                                                                         >  2:1511894/1‑150 (MQ=255)
                                                                              accatcaccTGGCGCCGCTACAGCTCGACCGAGGGGGTTATCCACAGGCGCCACTGACCCTGCAGCCGGCAATGCCGGCCGCGCAGTCGCAACGTACCGGCCTGCAGTGGCGCAGCAAGACCCGTCAGGTTAAGGAAGCGCACGCCATCg                                                                       >  2:707090/1‑150 (MQ=255)
                                                                                catgaccCGGCGCCGCTGCAGCTCGACCGAGGGGGTTATCCCCAGGCGCAACTGACCCTGCAGCCGGCAATGGCGGCCGCGCAGTCGCAACGTAGTGGCCTGCAGTGGCGCAGCAAGACCAGACAGGTTAAGGAAGCGCACGGCGTCggc                                                                     >  2:770018/1‑148 (MQ=255)
                                                                                catcaccCGGCGCCGCTGCAGCTCGACCGAGGGGGTTATCCACAGGCGCAACTGACCCTGCAGCCGGCAATGCCGGCCGCGCAGTCGCAACGTACCGGCCTGCAGTGGCGCAGCAAGACCCGTCAGGTTAAGGAAGCGCACGCCATCGAc                                                                     >  1:1582091/1‑150 (MQ=255)
                                                                                 atcaccCGGCGCCGCTGCAGCTCGACCGAGGGGGTTATCCACAGGCGCAACTGACCCTGCAGCCGGCAATGCCGGCCGCGCAGTCGCAACGTACCGGCCTGCAGTGGCGCAGCAAGACCCGTCAGGTTAAGGAAGCGCACGCCATCGACg                                                                    <  2:546716/150‑1 (MQ=255)
                                                                                  tcaccCGGCGCCGCTGCAGCTCGACCGAGGGGGTTATCCACAGGCGCAACTGACCCTGCAGCCGGCAATGCCGGCCGCGCAGTCGCAACGTACCGGCCTGCAGTGGCGCAGCAAGACCCGTCATGTTAAGGAAGCGCACGCCATCGACGc                                                                   <  1:356299/150‑1 (MQ=255)
                                                                                  ccaccCGGCGCCGCTGCAGCTCTACCGAGGGGGTTATCCACAGGCGCAACTGACCCTGCAGCCGGCTATACAGGTCTCGCAGTCGCAACGTACCGGCCTGCAGTGGCGCAGCAAGACCCGTCAGGTTAAGGAAGCGCACGCCATCGACGc                                                                   <  1:128831/149‑1 (MQ=255)
                                                                                    accCGGCGCCGCTGCAGCTCGACCGAGGGGGTTATCCACAGGCGCAACTGACCCTGCGGCCGGCATTGCCGGCtgcag                                                                                                                                         <  1:604281/78‑3 (MQ=38)
                                                                                    accCGGCGCCGCTGCAGCTCGACCGAGGGGGTTATCCACAGGCGCAACTGACCCTGCAGCCGGCAATGCCGGCCGCGCAGTCGCAACGTACCGGCCTGCAGTGGCGCAGCAAGACCCGTCAGGTTAAGGAAGCGCACGCCATCGACGC‑Ca                                                                >  2:1464423/1‑150 (MQ=255)
                                                                                     cccGGCGCCGCTGCAGCTCGACCGAGGGGGTTATCCACAGGCGCAACTGACCCTGCAGCCGGCAATGCCGGCCGCGCAGt                                                                                                                                      <  2:1212168/80‑1 (MQ=255)
                                                                                     cccGGCGCCGCTGCAGCTCGACCGAGGGGGTTATCCACAGGCGCAACTGACCCTGCAGCCGGCAATGCCGGCCGCGCAGt                                                                                                                                      >  1:1212168/1‑80 (MQ=255)
                                                                                      ccGGCGCCGCTGCAGCTCGACCGAGGGGGTTATCCACAGGCGCAACTGACCCTGCAGCCGGCAATGCCGGCCGCGCAGTCGCAACGTACCGGCCTGCAGTGGCGCAGCAAGACCCGTCAGGTTAAGGAAGCGCACGCCATCGACGC‑CACg                                                              >  2:255525/1‑150 (MQ=255)
                                                                                      ccGGCGCCGCTGCAGCTCGACCGAGGGGGTTATCCACAGGCGCAACTGACCCTGCAGCCGGCAATGCCGGCCGCGCAGTCGCAACGTACCGGCCTGCAGTGGCGCAGCAAGACCCGTCAGGTTAAGGAAGCGCACGCCATCGACGC‑CACg                                                              >  2:1880917/1‑150 (MQ=255)
                                                                                            ccgcTGCAGCTCGACCGAGGGGGTTATCCACAGGCGCAACTGACCCTGCAGCCGGCAATGCCGGCCGCGCAGTCGCAACGTACCGGCCTGCAGTGGCGCAGCAAGACCCGTCAGGTTAAGGAAGCGCACGCCATCGACGC‑CACG‑CT‑CGaa                                                      <  1:1536193/150‑1 (MQ=255)
                                                                                              gcTGCAGCTCGACCGAGGGGGTTATCCACAGGCGCAACTGACCCTGCAGCCGGCAATGCCGGCCGCGCAGTc                                                                                                                                     <  1:325689/72‑1 (MQ=255)
                                                                                              gcTGCAGCTCGACCGAGGGGGTTATCCACAGGCGCAACTGACCCTGCAGCCGGCAATGCCGGCCGCGCAGTc                                                                                                                                     >  2:325689/1‑72 (MQ=255)
                                                                                                tGCAGCTCGACCGAGGGGGTTATCCACAGGCGCAACTGACCCTGCAGCCGGCAATGCCGGCCGCGCAGTCGCAACGTACCGGCCTGCAGTGGCGCAGCAAGACCCGTCAGGTTAAGGAAGCGCACGCCATCGACGC‑CACG‑CT‑CGAACGCt                                                  <  1:1568342/150‑1 (MQ=255)
                                                                                                   aGCTCGACCGAGGGGGTTATCCACAGGCGCAACTGACCCTGCAGCCGGCAATGCCGGCCGCGCAGTCGCAACGTACCGGCCTGCAGTGGCGCAGCAAGACCCGTCAGGTTAAGGAAGCGCACGCCATCGACGC‑CACG‑CT‑CGAACGCTTGc                                               <  1:671492/150‑1 (MQ=255)
                                                                                                     cTCGACCGAGGGGGTTATCCACAGGCGCAACTGACCCTGCAGCCGGCAATGCCGGCCGCGCAGTCGCAACGTACCGGCCTGCAGTGGCGCAGCAAGACCCGTCAGGTTAAGGAAGCGCACGCCATCGACGC‑CACG‑c                                                            <  1:1387991/136‑1 (MQ=255)
                                                                                                     cTCGACCGAGGGGGTTATCCACAGGCGCAACTGACCCTGCAGCCGGCAATGCCGGCCGCGCAGTCGCAACGTACCGGCCTGCAGTGGCGCAGCAAGACCCGTCAGGTTAAGGAAGCGCACGCCATCGACGC‑CACG‑c                                                            >  2:1387991/1‑136 (MQ=255)
                                                                                                     cTCGACCGAGGGGGTTATCCACAGGCGCAACTGACCCTGCAGCCGGCAATGCCGGCCGCGCAGTCGCAACGTACCGGCCTGCAGTGGCGCAGCAAGACCCGTCAGGTTAAGGAAGCGCACGCCATCGACGC‑CACG‑CT‑CGAACGCTTGCCg                                             >  2:259744/1‑150 (MQ=255)
                                                                                                             aGGGGGTTATCCACAGGCGCAACTGACCCTGCAGCCGGCAATGCCGGCCGCGCAGTCGCAACGTACCGGCCTGCAGTGGCGCAGCAAGACCCGTCAGGTTaa                                                                                        <  1:625371/102‑1 (MQ=255)
                                                                                                             aGGGGGTTATCCACAGGCGCAACTGACCCTGCAGCCGGCAATGCCGGCCGCGCAGTCGCAACGTACCGGCCTGCAGTGGCGCAGCAAGACCCGTCAGGTTaa                                                                                        >  2:625371/1‑102 (MQ=255)
                                                                                                                 ggTTATCCACAGGCGCAACTGACCCTGCAGCCGGCAATGCCTGCCGCGCAGTCGCAACGTACCGGCCTGCAGt                                                                                                                 <  1:47138/73‑1 (MQ=255)
                                                                                                                 ggTTATCCACAGGCGCAACTGACCCTGCAGCCGGCAATGCCTGCCGCGCAGTCGCAACGTACCGGCCTGCAGt                                                                                                                 >  2:47138/1‑73 (MQ=255)
                                                                                                                   ttATCCACAGGCGCAACTGACCCTGCAGCCGGCAATGCCGGCCGCGCAGTCGCAACGTACCGGCCTGCAGTGGCGCAGCAAGACCCGTCAGGTTAAGGAAGCGCACGCCATCGACGC‑CACG‑CT‑CGAACGCTTGCCGATCGGCCTGTTCgc                               <  1:1880917/150‑1 (MQ=255)
                                                                                                                      tCCACATGCGCAACTGACCCTGCAGCCGGCAATGCCGGCCGCGCAGTCGCAACGTACCGGCCTGCAGTGGCGCAGCAAGACCCGTCAGGTTAAGGAAGCGCACGCCATCGACGC‑CACG‑CT‑CGAACGCTTGCCGATCGGCCTGTTCgccag                            >  2:1485596/1‑150 (MQ=255)
                                                                                                                       ccACAGGCGCAACTGACCCTGCAGCCGGCAATGCCGGCCGCGCAGTCGCAACGTACCGGCCTGCAGTGGCGCAGCAAGACCCGTCAGGTTAAGGAAGCGCACGCCATCGACGC‑CACG‑CT‑CGAACGCTTGCCGATCGGCCTGTTCgccagc                           <  1:1243046/150‑1 (MQ=255)
                                                                                                                        cacaGGCGCAACTGACCCTGCAGCCGGCAATGCCGGCCGCGCAGTCGCAACGTACCGGCCTGCAGTGGCGCAGCAAGACCCGTCAGGTTAAGGAAGCGCACGCCATCGACGC‑CACG‑CT‑CGAACGCTTGCCGATCGGCCTGTTCgccagcc                          <  1:255525/150‑1 (MQ=255)
                                                                                                                         acaGGCGCAACTGACCCTGCAGCCGGCAATGCCGGCCACGCAGTCGCAACGTACCGGCCTGCAGTGGCGCAGCAAGACCCGTCAGGTTAAGGAAGCGCACGCCATCGACGC‑CACG‑CT‑CGaa                                                      >  2:1470725/1‑121 (MQ=255)
                                                                                                                         acaGGCGCAACTGACCCTGCAGCCGGCAATGCCGGCCACGCAGTCGCAACGTACCGGCCTGCAGTGGCGCAGCAAGACCCGTCAGGTTAAGGAAGCGCACGCCATCGACGC‑CACG‑CT‑CGaa                                                      <  1:1470725/121‑1 (MQ=255)
                                                                                                                            gccgcAACTGACCCGGCACCCGCCAATGCCGGCCGCGCAGTCGCAACGTACCGGCCTGCAGTGGCGCAGCAAGACCCGTCAGGTTAAGGAAGCGCACGCCATCGACGC‑CACG‑CT‑CGAACGCTTGCCGATCGGCCTGTTCGCCAGCcagca                      <  1:770018/148‑1 (MQ=255)
                                                                                                                             gcgcAACTGACCCTGCAGCCGGCAATGCCGGCCGCGCAGTCGCAACGTACCGGCCTGCAGTGGCGCAGCAAGACCCGTCAGGTTAAGGAAGCGCACGCCATCGACGC‑CACG‑CT‑CGAACGCTTGCCGATCGGCCTGTTCGCCAGCcagcag                     >  2:927371/1‑150 (MQ=255)
                                                                                                                                   cTGACCCTGCAGCCGGCAATGCCGGCCGCGCAGTCGCAACGTACCGGCCTGCAGTGGCGCAGCAAGACCCGTCAGGTTAAGGAAGCGCACGCCATCGACGC‑CACG‑CT‑CGAACGCTTGCCGATCGGCCTGTTCGCCAGCCAGCAGCTCGAc               <  1:418349/150‑1 (MQ=255)
                                                                                                                                    tGACCCTGCAGCCGGCAATGCCGGCCGCGCAGTCGCAACGTACCGGCCTGCAGTGGCGCAGCAAGACCCGTCAGGTTAAGGAAGCGCACGCCATCGACGC‑CACG‑CT‑CGAACGCTTGCCGATCGGCCTGTTCGCCAGCCAGCAGCTCGAcc              <  1:1655858/150‑1 (MQ=255)
                                                                                                                                       cccTGCAGCCGGCAATGCCGGCCGCGCAGTCGCAACGTACCGGCCTGCAGTGGCGCAGCAAGACCCGTCAGGTTAAGGAAGCGCACGCCATCGACGC‑CACG‑CT‑CGAACGCTTGCCGATCGGCCTGTTCGCCAGCCAGCAGCTCGACCCGa           >  2:491768/1‑150 (MQ=255)
                                                                                                                                        ccTGCAGCCGGCAATGCCGGCCGCGCAGTCGCAACGTACCGGCCTGCAGTGGCGCAGCAAGACCCGTCAGGTTAAGGAAGCGCACGCCATCGACGC‑CACG‑CT‑CGAACGCTTGCCGATCGGCCTGTTCGCCAGCCAGCAGCTc                  <  2:33252/142‑1 (MQ=255)
                                                                                                                                        ccTGCAGCCGGCAATGCCGGCCGCGCAGTCGCAACGTACCGGCCTGCAGTGGCGCAGCAAGACCCGTCAGGTTAAGGAAGCGCACGCCATCGACGC‑CACG‑CT‑CGAACGCTTGCCGATCGGCCTGTTCGCCAGCCAGCAGCTc                  >  1:33252/1‑142 (MQ=255)
                                                                                                                                        ccTGCAGCCGGCAATGCCGGCCGCGCAGTCGCAACGTACCGGCCTGCAGTGGCGCAGCAAGACCCGTCAGGTTAAGGAAGCGCACGCCATCGACGC‑CACG‑CT‑CGAACGCTTGCCGATCGGCCTGTTCGCCAGCCAGCAGCTCGACCCGAg          >  2:1693039/1‑150 (MQ=255)
                                                                                                                                          tGCAGCCGGCAATGCCGGCCGCGCAGTCGCAACGTACCGGCCTGCAGTGGCGCAGCAAGACCCGTCAGGTTAAGGAAGCGCACGCCATCGACGC‑CACG‑CT‑CGAACGCTTGCCGATCGGCCTGTTCGCCAGCCAGCAGCTCGACCCGAggg        <  2:994246/150‑1 (MQ=255)
                                                                                                                                            cggCCGGCATTGCCGGCCGCGCAGTCGCAACGTACCGGCCTGCAGTGGCGCAGCAAGACCCGTCAGGTTAAGGAAGCGCACGCCATCGACGC‑CACG‑CT‑CGAACGCTTGCCGATCGGCCTGTTCGCCAGCCAGCAGCGCGACCCGAGGGTc      >  1:1172680/3‑150 (MQ=255)
                                                                                                                                              gCCGGCAATGCCGGCCGCGCAGTCGCAACGTACCGGCCTGCAGTGGCGCAGCAAGACCCGTCAGGTTAAGGAAGCGCACGCCATCGACGC‑CACG‑CT‑CGAACGCTTGCCGATCGGCCTGTTCGCCAGCCAGCAGCTCGACCCGAGGGTCga    >  2:161851/1‑150 (MQ=255)
                                                                                                                                              gCCGGCAATGCCGGCCGCGCAGTCGCAACGTACCGGCCTGCAGTGGCGCAGCAAGACCCGTCAGGTTAAGGAAGCGCACGCCATCGACGC‑CACG‑CT‑CGAACGCTTGCCGATCGGCCTGTTCGCCAGCCAGCAGCTCGACCCGAGGGTCga    <  1:927371/150‑1 (MQ=255)
                                                                                                                                               ccGGCAATGCCGGCCGCGCAGTCGCAACGTACCGGCCTGCAGTGGCGCAGCAAGACCCGTCAGGTTAAGGAAGCGCACGCCATCGACGC‑CAAG‑CT‑CGAACGCTTGCCGATCGGCCTGTTCGCCAGCCAGCAGCTCGACCCGAGGGTCgag   >  1:1403236/1‑150 (MQ=255)
                                                                                                                                                cGGCAATGCCGGCCGCGCAGTCGCAACGTCAGGGCCTGCAGTGGCGCAGCAAGACCCGTCAGCTTAAGGATGAGATCGGAA‑‑GA‑GCACACGTCTGCACTCCCGTGCCGATCGGCCTGTTCGCCAGCCAGCAGCTCGACCCGAGGGCCgaga  >  1:443604/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                   
GTCACCACCCAGCTTTCGTCGGCCCGGCTGTACAAACCATACGCGGCCAACTCGGCGTCATCGGCATGTGGCGCGATCACCATCACCCGGCGCCGCTGCAGCTCGACCGAGGGGGTTATCCACAGGCGCAACTGACCCTGCAGCCGGCAATGCCGGCCGCGCAGTCGCAACGTACCGGCCTGCAGTGGCGCAGCAAGACCCGTCAGGTTAAGGAAGCGCACGCCATCGACGC‑CACG‑CT‑CGAACGCTTGCCGATCGGCCTGTTCGCCAGCCAGCAGCTCGACCCGAGGGTCGAGA  >  NC_002947/5623163‑5623456

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: ‑