Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,949,113 T→C 6.6% S141G (AGC→GGC)  PP_3485 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,949,1130TC 6.6% 268.4 / 6.5 91S141G (AGC→GGC) PP_3485hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (42/43);  new base C (2/4);  total (44/47)
Fisher's exact test for biased strand distribution p-value = 6.78e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.00e-01

TGCCCGATGAGCGTTGCCCGAGCCAAGGCCAGGCCGTGTTCGCGTATGACCAATGGATTATCACTTTGGTTCTGCGCTCCTTGCAGTCTGTCCTGGTTGGCATCAAAAAATGCCTCCACTT‑CGCCCTCGCTCGGCGC‑GTCGACTGCGCTGAGC‑AGGGCGTACACTTGCTGGGCCGCCATCTCATGCCGGGTGTATTC‑AGTGTACTGCGCCCTGTCGAAGCCCGCCTCTGCCAGTCGCCTTTCGAATAACGCCGGGCTGCCAAATGCCGCTTCGATCTCGCCAACCTGCGCCGACACT  >  NC_002947/3948965‑3949261
                                                                                                                                                      |                                                                                                                                                      
tcgccGATGAGCGTTGCCCGAGCCAAGGCCTGGCCGTGTTCGCGTATGACCAATTGTTTATCAGTTTGGTTCTGCTGTCCTTGCAGTCTGTCCTTGTTGGCATCAAAAAATTCCTCGACTT‑CGCCCTCGCTCGGCGC‑GTCGACTGCGCt                                                                                                                                                        <  5:356631/146‑1 (MQ=255)
tGCCCGATGAGCGTTGCCCGAGCCAAGGCCAGGCCGTGTTCGCGTATGACCAATGGATTATCACTTTGGTTCTGCGCTCCTTGCAGTCTGTCCTGGTTGGCATCAAAAAATGCCTCCACTT‑CGCCCTCGCTCGGCGC‑GTCGACTGCGCt                                                                                                                                                        >  1:9559/1‑149 (MQ=255)
     gATGAGCGTTGCCCGAGCCAAGGCCAGGCCGTGTTCGCGTATGACCAATGGATTATCACTTTGGTTCTGCGCTCCTTGCAGTCTGTCCTGGTTGGCATCAAAAAATGCCTCCACTT‑CGCCCTCGCTCGGCGC‑GTCGACTGCGCTGAg                                                                                                                                                     >  2:53611/1‑147 (MQ=255)
       tGAGCGTTGCCCGAGGCAAGGCCAGGCCGTGGTAGCGAATGACCAATGCATTATCACTTTGGTTCTGCGCTCCTTGCTGTCTGTCCTGGTTGGCATCAAAAAATGCCTCCACTT‑CGCCCT‑GCTCAGCGCAGTCGAC‑GCGCCGAGcgagg                                                                                                                                                <  7:211468/149‑5 (MQ=37)
       tGAGCGTTGCCCGAGCCAAGGCCAGGCCGTGTTCGCGTATGACCAATGGATTATCACTTTGGTTCTGCGCTCCTTGCAGTCTGTCCTGGTTGGCATCAAAAAATGCCTCCACTT‑CGCCCTCGCTCGGCGC‑GTCGACTGCGCTGAGC‑Agg                                                                                                                                                <  4:38650/149‑1 (MQ=255)
       tGAGCGTTGCCCGAGCCAAGGCCAGGCCGTGTTCGCGTATGACCAATGGATTATCACTTTGGTTCTGCGCTCCTTGCAGTCTGTCCTGGTTGGCATCAAAAAATGCCTCCACTT‑CGCCCT‑GCTCAGCGCAGTCGAC‑GCGCCGAGcgagg                                                                                                                                                <  8:346036/149‑5 (MQ=38)
        gatcGTTGGCCGAGCCAAGTCCAGGCCGTGTTCGCTTATGACCAATGGTTTATCACTTTGGTTCTGCTCTCCTTGCAGTCTGTCCTGGTTTGCATCAAAAAATGCCTCCACTT‑CGCCCTCGCTCGGCGC‑GTCGACTGCTCTGAGC‑Agg                                                                                                                                                <  8:177066/145‑1 (MQ=255)
         aGCGTTGCCCGAGCCAAGGCCAGGCCGTGTTCGCGTATGACCAATGGATTATCACTTTGGTTCTGCGCTCCTTGCAGTCTGTCCTGGTTGGCATCAAAAAATGCCTCCACTT‑CGCCCT‑GCTCAGCGCAGTCGAC‑GCGCCGAGcgaggg                                                                                                                                               <  1:239362/148‑6 (MQ=38)
         aGCGTTGCCCGAGCCAAGGCCAGGCCGTGTTCGCGTATGACCAATGGATTATCACTTTGGTTCTGCGCTCCTTGCAGTCTGTCCTGGTTGGCATCAAAAAATGCCTCCACTT‑CGCCCT‑GCTCAGCGCAGTCGAC‑GCGCCGAGcgaggg                                                                                                                                               <  5:69851/148‑6 (MQ=38)
            gTTGCCCGAGCCAAGGCCAGGCCGTGTTCGCGTATGACCAATGGATTATCACTTTGGTTCTGCGCTCCTTGCAGTCTGTCCTGGTTGGCATCAAAAAATGCCTCCACTT‑CGCCCTCGCTCGGCGC‑GTCGACTGCGCTGAGC‑AGGGCGTa                                                                                                                                           >  4:309667/1‑149 (MQ=255)
            gTTGCCCGAGCCAAGGCCAGGCCGTGTTCGCGTATGACCAATGGATTATCACTTTGGTTCTGCGCTCCTTGCAGTCTGTCCTGGTTGGCATCAAAAAATGCCTCCACTT‑CGCCCTCGCTCGGCGC‑GTCGACTGCGCTGAGC‑AGGGCGTa                                                                                                                                           <  5:238897/149‑1 (MQ=255)
             ttGCCCGAGCCAAGGCCAGGCCGTGTTCGCGTATGACCAATGGATTATCACTTTGGTTCTGCGCTCCTTGCAGTCTGTCCTGGTTGGCATCAAAAAATGCCTCCACTT‑CGCCCTCGCTCGGCGC‑GTCGACTGCGCTGAGC‑AGGGCGTac                                                                                                                                          <  8:274042/149‑1 (MQ=255)
               gCCCGAGCCAAGGCCAGGCCGTGTTCGCGTATGACCAATGGATTATCACTTTGGTTCTGCGCTCCTTGCAGTCTGTCCTGGTTGGCATCAAAAAATGCCTCCACTT‑CGCCCTCGCTCGGCGC‑GTCGACTGCGCTGAGC‑AGGGCGTaca                                                                                                                                         <  3:309667/148‑1 (MQ=255)
                 ccGAGCCAATGCCAGGCCGTGTTCGCGTATGACCAATGGATTATCACTTTGGTTCTGCGCTCCTTGCAGTCTGTCCTGGTTGGCATCAAAAAATGCCTCCACTT‑CGCCCTCGCTCGGCGC‑GTCGACTGCTCTGAGC‑AGGGCGTACACtt                                                                                                                                      <  3:327923/149‑1 (MQ=255)
                  cGAGCCAAGGCCAGGCCGTGTTCGCGTATGACCAATGGATTATCACTTTGGTTCTGCGCTCCTTGCAGTCTGTCCTGGTTGGCATCAAAAAATGCCTCCACTT‑CGCCCTCGCTCGGCGC‑GTCGACTGCGCTGAGC‑AGGGCGTACACTTg                                                                                                                                     <  2:149197/149‑1 (MQ=255)
                     gCCAAGGCCAGGCCGTGTTCGCGTATGACCAATGGATTATCACTTTGGTTCTGCGCTCCTTGCAGTCTGTCCTGGTTGGCATCAAAAAATGCCTCCACTT‑CGCCCTCGCTCGGCGC‑GTCGACTGCGCTGAGC‑AGGGCGTACACTTGCTg                                                                                                                                  <  2:129124/149‑1 (MQ=255)
                      ccAAGGCCAGGCCGTGTTCGCGTATGACCAATGGATTATCACTTTGGTTCTGCGCTCCTTGCAGTCTGTCCTGGTTGGCATCAAGAAATGCCTCCACTT‑CGCCCTCGCTCGGCGC‑GTCGACTGCGCTGAGC‑AGGGCGTACACTTGCTgg                                                                                                                                 >  8:243970/1‑149 (MQ=255)
                      ccAAGGCCAGGCCGTGTTCGCGTATGACCAATGGATTATCACTTTGGTTCTGCGCTCCTTGCAGTCTGTCCTGGTTGGCATCAAAAAATGCCTCCACTT‑CGCCCTCGCTCGGCGC‑GTCGACTGCGCTGAGC‑AGGGCGTACACTTGCTgg                                                                                                                                 >  7:235764/1‑149 (MQ=255)
                      ccAAGGCCAGGCCGTGTTCGCGTATGACCAATGGATTATCACTTTGGTTCTGCGCTCCTTGCAGTCTGTCCTGGTTGGCATCAAAAAATGCCTCCACTT‑CGCCCTCGCTCGGCGC‑GTCGACTGCGCTGAGC‑AGGGCGTACACTTGCTgg                                                                                                                                 >  5:391968/1‑149 (MQ=255)
                      ccAAGGCCAGGCCGTGTTCGCGTATGACCAATGGATTATCACTTTGGTTCTGCGCTCCTTGCAGTCTGTCCTGGTTGGCATCAAAAAATGCCTCCACTT‑CGCCCTCGCTCGGCGC‑GTCGACTGCGCTGAGC‑AGGGCGTACACTTGCTgg                                                                                                                                 >  3:256370/1‑149 (MQ=255)
                       cAAGGCCAGGCCGTGTTCGCGTATGACCAATGGATTATCACTTTGGTTCTGCGCTCCTTGCAGTCTGTCCTGGTTGGCATCAAAAAATGCCTCCACTT‑CGCCCTCGCTCGGCGC‑GTCGACTGCGCTGAGC‑AGGGCGTACACTTGCTggg                                                                                                                                <  2:314704/149‑1 (MQ=255)
                         aggccaggccGTGTTCGCGTATGACCAATGGATTATCACTTTGGTTCTGCGCTCCTTGCAGTCTGTCCTGGTTGGCATCAAAAAATGCCTCCACTT‑CGCCCTCGCTCGGCGC‑GTCGACTGCGCTGAGC‑AGGGCGTACACTTGCTGGgcc                                                                                                                              <  6:359469/149‑1 (MQ=255)
                           gccaggccGTGTTCGCGTATGACCAATGGATTATCACTTTGGTTCTGCGCTCCTTGCAGTCTGTCCTGGTTGGCATCAAAAAATGCCTCCACTT‑CGCCCTCGCTCGGCGC‑GTCGACTGCGCTGAGC‑AGGGCGTACACTTGCTGGgccgc                                                                                                                            >  2:404153/1‑149 (MQ=255)
                             caggccGTGTTCGCGTATGACCAATGGATTATCACTTTGGTTCTGCGCTCCTTGCAGTCTGTCCTGGTTGGCATCAAAAAATGCCTCCACTT‑CGCCCTCGCTCGGCGC‑GTCGACTGCGCTGAGC‑AGGGCGTACACTTGCTGGGCCGCCa                                                                                                                          >  3:107302/1‑149 (MQ=255)
                               ggccGTGTTCGCGTATGACCAATGGATTATCACTTTGGTTCTGCGCTCCTTGCAGTCTGTCCTGGTTGGCATCAAAAAATGCCTCCACTT‑CGCCCTCGCTCGGCGC‑GTCGACTGCGCTGAGC‑AGGGCGTACACTTGCTGGGCCGCCAtc                                                                                                                        <  8:235764/149‑1 (MQ=255)
                                gccGTGTTCGCGTATGACCAATGGATTATCACTTTGGTTCTGCGCTCCTTGCAGTCTGTCCTGGTTGGCATCAAAAAATGCCTCCACTT‑CGCCCTCGCTCGGCGC‑GTCGACTGCGCTGAGC‑AGGGCGTACACTTGCTGGGCCGCCAtct                                                                                                                       >  5:180027/1‑149 (MQ=255)
                                 ccGTGTTCGCGTATGACCAATGGATTATCACTTTGGTTCTGCGCTCCTTGCAGTCTGTCCTGGTTGGCATCAAAAAATGCCTCCACTT‑CGCCCTCGCTCGGCGC‑GTCGACTGCGCTGAGC‑AGGGCGTACACTTGCTGGGCCGCCAtct                                                                                                                       >  7:41235/1‑148 (MQ=255)
                                        cgcgTATGACCAATGGATTATCACTTTGGTTCTGCGCTCCTTGCAGTCTGTCCTGGTTGGCATCAAAAAATGCCTCCACTT‑CGCCCTCGCTCGGCGC‑GTCGACTGCGCTGAGC‑AGGGCGTACACtt                                                                                                                                      >  6:323791/1‑126 (MQ=255)
                                        cgcgTATGACCAATGGATTATCACTTTGGTTCTGCGCTCCTTGCAGTCTGTCCTGGTTGGCATCAAAAAATGCCTCCACTT‑CGCCCTCGCTCGGCGC‑GTCGACTGCGCTGAGC‑AGGGCGTACACtt                                                                                                                                      <  5:323791/126‑1 (MQ=255)
                                          cgTATGACCAATTGATTATCACTTTGGTTCTGCGCTCCTTGCAGTCTGTCCTGGTTGGCATCAAAAAATGCCTCCACTT‑CGCCCTCGCTCGGCGC‑GTCGACTGCGCTGAGC‑AGGGCGTACACTTGCTGGGCCGCCATCTCATGCCGGgt                                                                                                             <  2:167531/149‑1 (MQ=255)
                                           gTATGACCAATGGATTATCACTTTGGTTCTGCGCTCCTTGCAGTCTGTCCTGGTTGGCATCAAAAAATGCCTCCACTT‑CGCCCTCGCTCGGCGC‑GTCGACTGCGCTGAGC‑AGGGCGTACACTTGCTGGGCCGCCATCTCATGCCGGgtg                                                                                                            <  4:271718/149‑1 (MQ=255)
                                              tGACCAATGGATTATCACTTTGGTTCTGCGCTCCTTGCAGTCTGTCCTGGTTGGCATCAAAAAATGCCTCCACTT‑CGCCCTCGCTCGGCGC‑GTCGACTGCGCTGAGC‑AGGGCGTACACTTGCTGGGCCGCCATCTCATGCCGGGTGTAt                                                                                                         <  2:64393/149‑1 (MQ=255)
                                                aCCAATGGATTATCACTTTGGTTCTGCGCTCCTTGCAGTCTGTCCTGGTTGGCATCAAAAAATGCCTCCACTT‑CGCCCTCGCTCGGCGC‑GTCGACTGCGCTGAGC‑AGGGCGTACACTTGCTGGGCCGCCATCTCATGCCGGGTGTAtt                                                                                                        <  6:214512/148‑1 (MQ=255)
                                                aCCAATGGATTATCACTTTGGTTCTGCGCTCCTTGCAGTCTGTCCTGGTTGGCATCAAAAAATGCCTCCACTT‑CGCCCTCGCTCGGCGC‑GTCGACTGCGCTGAGC‑AGGGCGTACACTTGCTGGGCCGCCATCTCATGCCGGGTGTAtt                                                                                                        <  4:256370/148‑1 (MQ=255)
                                                  cAATGGATTATCACTTTGGTTCTGCGCTCCTTGCAGTCTGTCCTGGTTGGCATCAAAAAATGCCTCCACTT‑CGCCCTCGCTCGGCGC‑GTCGACTGCGCTGAGC‑AGGGCGTACACTTGCTGGGCCGCCATCTCATGCCGGgt                                                                                                             <  4:30655/141‑1 (MQ=255)
                                                  cAATGGATTATCACTTTGGTTCTGCGCTCCTTGCAGTCTGTCCTGGTTGGCATCAAAAAATGCCTCCACTT‑CGCCCTCGCTCGGCGC‑GTCGACTGCGCTGAGC‑AGGGCGTACACTTGCTGGGCCGCCATCTCATGCCGGgt                                                                                                             >  3:30655/1‑141 (MQ=255)
                                                  cAATGGATTATCACTTTGGTTCTGCGCTCCTTGCAGTCTGTCCTGGTTGGCATCAAAAAATGCCTCCACTT‑CGCCCTCGCTCGGCGC‑GTCGACTGCGCTGAGC‑AGGGCGTACACTTGCTGGGCCGCCATCTCATGCCGGGTGTATTC‑Ag                                                                                                    >  2:313355/1‑149 (MQ=255)
                                                  cAATGGATTATCACTTTGGTTCTGCGCTCCTTGCAGTCTGTCCTGGTTGGCATCAAAAAATGCCTCCACTT‑CGCACTCGCTCGGCGC‑GTCGACTGCGCTGAGC‑AGGGCGTACACTTGCTGGGCCGCCATCTCATGCCGGGTGTATTA‑Ag                                                                                                    >  5:365527/1‑149 (MQ=255)
                                                          tATCACTTTGGTTCTGCGCTCCTTGCAGTCTGTCCTGGTTGGCATCAAAAAATGCCTCCACTT‑CGCCCTCGCTCGGCGC‑GTCGACTGCGCTGAGC‑AGGGCGTACACTTGCTGGGCCGCCATCTCATGCCGGGTGTATTC‑AGTGTACTgc                                                                                            >  2:319788/1‑149 (MQ=255)
                                                            tCACTTTGGTTCTGCGCTCCTTGCAGTCTGTCCTGGTTGGCATCAAAAAATGCCTCCACTT‑CGCCCTCGCTCGGCGC‑GTCGACTGCGCTGAGC‑AGGGCGTACACTTGCTGGGCCGCCATCTCAt                                                                                                                    <  7:20825/124‑1 (MQ=255)
                                                            tCACTTTGGTTCTGCGCTCCTTGCAGTCTGTCCTGGTTGGCATCAAAAAATGCCTCCACTT‑CGCCCTCGCTCGGCGC‑GTCGACTGCGCTGAGC‑AGGGCGTACACTTGCTGGGCCGCCATCTCAt                                                                                                                    >  8:20825/1‑124 (MQ=255)
                                                             cACTTTGGTTCTGCGCTCCTTGCAGTCTGTCCTGGTTGGCATCAAAAAATGCCTCCACTT‑CGCCCTCGCTCGGCGC‑GTCGACTGCGCTGAGC‑AGGGCGTACACTTGc                                                                                                                                    >  7:54081/1‑107 (MQ=255)
                                                             cACTTTGGTTCTGCGCTCCTTGCAGTCTGTCCTGGTTGGCATCAAAAAATGCCTCCACTT‑CGCCCTCGCTCGGCGC‑GTCGACTGCGCTGAGC‑AGGGCGTACACTTGc                                                                                                                                    <  8:54081/107‑1 (MQ=255)
                                                              aCTTTGGTTCTGCGCTCCTTGCAGTCTGTCCTGGTTGGCATCAAAAAATGCCTCCACTT‑CGCCCTCGCTCGGCGC‑GTCGACTGCGCTGAGC‑AGGGCGTACACTTGCTGGGCCGCCATCTCATGCCGGGTGTATTC‑AGTGTACTGCGcc                                                                                         <  8:41235/148‑1 (MQ=255)
                                                               cTTTGGTTCTGCGCTCCTTGCAGTCTGTCCTGGTTGGCATCAAAAAATGCCTCCACTTACGCCCTCGCTCGGCGC‑GTCGACTGCGCTGAGC‑AGGGCGTACACTTGCTGGGCCGCCATCTCATGCCGGGTGTATTC‑AGTGTACTGCGccc                                                                                        >  1:319616/1‑149 (MQ=255)
                                                               cTTTGGTTCTGCGCTCCTTGCAGTCTGTCCTGGTTGGCATCAAAAAATGCCTCCACTT‑CGCCCTCGCTCGGCGC‑GTCGACTGCGCTGAGC‑AGGGCGTACACTTGCTGGGCCGCCATCTCATGCCGGGTGTATTC‑AGTGTACTGCGCCCt                                                                                       <  2:9559/149‑1 (MQ=255)
                                                                    gTTCTGCGCTCCTTGCAGTCTGTCCTGGTTGGCATCAAAAAATGCCTCCACTT‑CGCCCTCGCTCGGCGC‑GTCGACTGCGCTGAGC‑AGGGCGTACACTTGCTGGGCCGCCATCTCATGCCGGGTGTATTC‑AGTGTACTGCGCCCTGTCg                                                                                   <  1:53611/148‑1 (MQ=255)
                                                                         gcgcTCCTTGCAGTCTGTCCTGGTTGGCATCAAAAAATGCCTCCACTT‑CGCCCTCGCTCGGCGC‑GTCGACTGCGCTGAGC‑AGGGCGTACACTTGCTGGGCCGCCATCTCATGCCGGGTGTATTC‑AGTGTACTGCGCCCTGTCGAAGccc                                                                             >  3:75893/1‑149 (MQ=255)
                                                                           gcTCCTTGCAGTCTGTCCTGGTTGGCATCAAAAAATGCCTCCACTT‑CGCCCTCGCTCGGCGC‑GTCGACTGCGCTGAGC‑AGGGCGTACACTTGCTGGGCCGCCATCTCATGCCGGGTGTATTC‑AGTGTACTGCGCCCTGTCGAAGCCCGc                                                                           >  8:144800/1‑149 (MQ=255)
                                                                               cTTGCAGTCTGTCCTGGTTGGCATCAAAAAATGCCTCCACTT‑CGCCCTCGCTCGGCGC‑GTCGACTGCGCTGAGC‑AGGGCGTACACTTGCTGGGCCGCCATCTCATGCCGGGTGTATTC‑AGTGTACTGCGCCCTGTCGAAGCCCGCctc                                                                        <  3:295052/148‑1 (MQ=255)
                                                                                ttGCAGTCTGTCCTGGTTGGCATCAAAAAATGCCTCCACTT‑CGCCCTCGCTCGGCGC‑GTCGACTGCGCTGAGC‑AGGGCGTACACTTGCTGGGCCGCCATCTCATGCCGGGTGTATTC‑AGTGTACTGCGCCCTGTCGAAGCCCGCCTCTg                                                                      >  5:179818/1‑149 (MQ=255)
                                                                                  gccgTCGGTCCTGGTGGGCATCAAAAAATGCCTCCACTG‑CGCCCTCGCTCGGCGC‑GTCGACTGCGCTGAGC‑AGGGCGTACACTTGCTGGGCCGCCATCTCATGCCGTGTGTATTC‑AGTGTACTGCGCCCTGTCGAAGCCCGCCTCTGc                                                                     <  8:158721/145‑1 (MQ=255)
                                                                                  gCAGTCTGTCCTGGTTGGCATCAAAAAATGCCTCCACTT‑CGCCCTCGCTCGGCGC‑GTCGACTGCGCTGAGC‑AGGGCGTACACTTGCTGGGCCGCCATCTCATGCCGGGTGTATTC‑Agt                                                                                                   >  7:322563/1‑118 (MQ=255)
                                                                                  gCAGTCTGTCCTGGTTGGCATCAAAAAATGCCTCCACTT‑CGCCCTCGCTCGGCGC‑GTCGACTGCGCTGAGC‑AGGGCGTACACTTGCTGGGCCGCCATCTCATGCCGGGTGTATTC‑Agt                                                                                                   <  8:322563/118‑1 (MQ=255)
                                                                                  gCAGTCTGTCCTGGTTGGCATCAAAAAATGCCTCCACTT‑CGCCCTCGCTCGGCGC‑GTCGACTGCGCTGAGC‑AGGGCGTACACTTGCTGGGCCGCCATCTCATGCCGGGTGTATTC‑AGTGTACTGCGCCCTGTCGAAGCCCGCCTCTGc                                                                     <  7:144800/148‑1 (MQ=255)
                                                                                   cAGTCTGTCCTGGTTGGCATCAAAAAATGCCTCCACTT‑CGCCCTCGCTCGGCGC‑GTCGACTGCGCTGAGC‑AGGGCGTACACTTGCTGGGCCGCCATCTCATGCCGGGTGTATTC‑AGTGTACTGCGCCCTGTCGAAGCCCGCCTCTGcc                                                                    <  1:404153/148‑1 (MQ=255)
                                                                                          tCCTGGTTGGCATCAAAAAATGCCTCCACTT‑CGCCCTCGCTCGGCGC‑GTCGACTGCGCTGAGC‑AGGGCGTACACTTGCTGGGCCGCCATCTCATGCCGGGTGTATTC‑AGTGTACTGCGCCCTGTCGAAGCCCGCCTCTGCCAGTCGcc                                                             <  4:75893/148‑1 (MQ=255)
                                                                                                  ggCATCAAAAAATGCCTCCACTT‑CGCCCTCGCTCGGCGC‑GTCGACTGCGCTGAGA‑AGCGCGTCCAATTGCTGGGCCGCCATCTCATCCCGGGTGTATTC‑AGTGTACTGCTCCCTGTCGAAGCCCGCCTCTGCCAGTCGCCTTTCGAAt                                                     >  1:268122/1‑148 (MQ=255)
                                                                                                   gCATCAAAAAATGCCTCCACTT‑CGCCCTCGCTCGGCGC‑GTCGACTGCGCTGAGC‑AGGGCGTACACTTGCTGGGCCGCCATCTCATGCCGGGTGTATTC‑AGTGTACTGCGCCCTGTCGAAGCCCGCCTCTGCCAGTCGCCTTTCGAATaa                                                   <  2:319616/149‑1 (MQ=255)
                                                                                                   gCATCAAAAAATGCATGCACTT‑GGCCCTCGCTCGGGGG‑GGGGACTGCGCTTAGC‑AGGGCTTACACTTGCTGGGCAGCCATCTCATGCCTGGTGTAATC‑AGTGTACTGCGCCCTGTCGAAGCCCGCCTCTGCCAGTCGCCTTTCGAATa                                                    <  5:349401/148‑1 (MQ=255)
                                                                                                      tCAAAAAATGCCTCCACTT‑CGCCCTCGCTCGGCGC‑GTCGACTGCGCTGAGC‑AGGGCTTACACTTGCTGGGCCGCCATCTCATGCCGGGTGTATTC‑AGTGTACTGCGCCCTGTCGAAGCCCGCCTCTGCCAGTCGCCTTTCGAATAACGc                                                >  4:217011/1‑149 (MQ=255)
                                                                                                      tCAAAAAATGCCTCCACTT‑CGCCCTCGCTCGGCGC‑GTCGACTGCGCTGAGC‑AGGGCGTACACTTGCTGGGCCGCCATCTCATGCCGGGTGTATTC‑AGTGTACTGCGCCCTGTCGAAGCCCGCCTCTGCCAGTCGCCTTTCGAATAACGc                                                >  8:174697/1‑149 (MQ=255)
                                                                                                        aaaaaaTGCCTCCACTT‑CGCCCTCGCTCGGCGC‑GTCGACTGCGCTGAGC‑AGGGCGTACACTTGCTGGGCCGCCATCTCATGCCGGGTGTATTC‑AGTGTACTGCGCCCTGTCGAAGCCCGCCTCTGCCAGTCGCCTTTCGAATAACGcc                                               <  6:179818/148‑1 (MQ=255)
                                                                                                          aaaaTGCCTCCACTT‑CGCCCTCGCTCGGCGC‑GTCGACTGCGCTGAGC‑AGGGCGTACACTTGCTGGGCCGCCATCTCATGCCGGGTGTATTC‑AGTGTACTGCGCCCTGTCGAAGCCCGCCTCTGCCAGTCGCCTTTCGAATAACGCCggg                                            >  2:379822/1‑149 (MQ=255)
                                                                                                          aaaaTGCCTCCACTT‑CGCCCTCGCTCGGCGC‑GTCGACTGCGCTGAGC‑AGGGCGTACACTTGCTGGGCCGCCATCTCATGCCGGGTGTATTC‑AGTGTACTGCGCCCTGTCGAAGCCCGCCTCTGCCAGTCGCCTTTCGAATAACGCCggg                                            <  6:391968/149‑1 (MQ=255)
                                                                                                             aTGCCTCCAGTT‑CGCCCTCGCTCGGCGG‑GTCGACTGCGCTGAGC‑AGGGCGTACACTTGCTGGGCCGCCATCTCATGCCGGGTGTTTTCTAGTGTACTGCGCCCTGTCGAAGCCCGCCTCTGCCAGTCGCCTTTCGAATAACGCCGGGCt                                          <  1:319788/149‑1 (MQ=255)
                                                                                                              tGCCTCCACTT‑CGCCCTCGCTCGGCGC‑GTCGACTGCGCTGAGC‑AGGGCGTACACTTGCTGGGCCGCCATCTCATGCCGGGTGTATTC‑AGTGTAc                                                                                               >  5:155023/1‑94 (MQ=255)
                                                                                                              tGCCTCCACTT‑CGCCCTCGCTCGGCGC‑GTCGACTGCGCTGAGC‑AGGGCGTACACTTGCTGGGCCGCCATCTCATGCCGGGTGTATTC‑AGTGTAc                                                                                               <  6:155023/94‑1 (MQ=255)
                                                                                                                ccTCCACTT‑CGCCCTCGCTCGGCGC‑GTCGACTGCGCTGAGC‑AGGGCGTACACTTGCTGGGCCGCCATCTCATGCCGGGTGTATTC‑AGTGTACTGCGCCCTGTCGAAGCCCGCCTCTGCCAGTCGCCTTTCGAATAACGCCGGGCTGcc                                       >  3:103096/1‑148 (MQ=255)
                                                                                                                  tCCACTT‑CGCCCTCGCTCGGCGC‑GTCGACTGCGCTGAGC‑AGGGCGTACACTTGCTGGGCCGCCATCTCATGCCGGGTGTATTC‑a                                                                                                     >  4:147284/1‑84 (MQ=255)
                                                                                                                  tCCACTT‑CGCCCTCGCTCGGCGC‑GTCGACTGCGCTGAGC‑AGGGCGTACACTTGCTGGGCCGCCATCTCATGCCGGGTGTATTC‑a                                                                                                     <  3:147284/84‑1 (MQ=255)
                                                                                                                  tCCACTT‑CGCCCTCGCTCGGCGC‑GTCGACTGCGCTGAGC‑AGGGCGTACACTTGCTGGGCCGCCATCTCATGCCGGGTGTATTC‑AGTGTACTGCGCCCTGTCGAAGCCCGCCTCTGCCAGTCGCCTTTCGAATAACGCCGGGCTGCCaaa                                    >  1:373353/1‑149 (MQ=255)
                                                                                                                  tCCACTT‑CGCCCTCGCTCGGCGC‑GTCGACTGCGCTGAGC‑AGGGCGTACACTTGCTGGGCCGCCATCTCATGCCGGGTGTATTC‑AGTGTACTGCGCCCTGTCGAAGCCCGCCTCTGCCAGTCGCCTTTCGAATAACGCCGGGCTGCCaaa                                    >  7:52960/1‑149 (MQ=255)
                                                                                                                   ccACTT‑CGCCCTCGCTCGGCGC‑GTCGACTGCGCTGAGC‑AGGGCGTACACTTGCTGGGCCGCCATCTCATGCCGGGTGTATTC‑AGTGTACTGCGCCCTGTCGAAGCCCGCCTCTGCCAGTCGCCTTTCGAATAACGCCGGGCTGCCAAAt                                   >  2:216529/1‑149 (MQ=255)
                                                                                                                    cACTT‑CGCCCTCGCTCGGCGC‑GTCGACTGCGCTGAGC‑AGGGCGTACACTTGCTGGGCCGCCATCTCATGCCGGGTGTATTC‑AGTGTACTGCGCCCTGTCGAAGCCCGCCTCTGCCAGTCGCCTTTCGAATAACGCCGGGCTGCCAAATg                                  >  4:196683/1‑149 (MQ=255)
                                                                                                                    cACTT‑CGCCCTCGCTCGGCGC‑GTCGACTGCGCTGAGC‑AGGGCGTACACTTGCTGGGCCGCCATCTCATGCCGGGTGTATTC‑AGTGTACTGCGCACTGTCGAAGCCCGCCTCTGCCAGTCGCCTTTCGAATAACGCCGGTCTGCCAAATg                                  >  6:207848/1‑149 (MQ=255)
                                                                                                                          cGCCCTCGCTCGGCGC‑GTCGACTGCGCTGAGC‑AGGGCGTACACTTGCTGGGCCGCCATCTCATGCCGGGTGTATTC‑AGTGTACTGCGCCCTGTCGAAGCCCGCCTCTGCCAGTCGCCTTTCGAATAACGCCGGGCTGCCAAATGCCGc                              <  4:107302/148‑1 (MQ=255)
                                                                                                                           gCCCTCGCTCGGCGC‑GTCGACTGCGCTGAGC‑AGGGCGTACACTTGCTGGGCCGCCATCTCATGCCGGGTGTATTC‑AGTGTACTGCGCCCTGTCGAAGCCCGCCTCTGCCAGTCGCCTTTCGAATAACGCCGGGCTGCCAAATGCCGCt                             <  2:373353/148‑1 (MQ=255)
                                                                                                                             ccTCGCTCGGCGC‑GTCGACTGCGCTGAGC‑AGGGCGTACACTTGCTGGGCCGCCATCTCATGCCGGGTGTATTC‑AGTGTACTGCGCCCTGTCGAAGCCCGCCTCTGCCATTCGCCTTTCGAATAACGCCGGGCTGCCAAATGCCGCTTc                           >  2:73221/1‑148 (MQ=255)
                                                                                                                              cctgctcaGCGCAGTCGAC‑GCGCCGAGCGAGGGCGTACACTTGCTGGGCCGCCATCTCATGCCGGGTGTATTC‑AGTGTACTGCGCCCTGTCGAAGCCCGCCTCTGCCAGTCGCCTTTCGAATAACGCCGGGCTGCCAAATGCCGCTTc                           >  1:56772/4‑148 (MQ=38)
                                                                                                                                tgctcaGCGCAGTCGAC‑GCGCCGAGC‑AGGGCGTACACTTGCTGGGCCGCCATCTCATGCCGGGTGTATTC‑AGTGTACTGCGCCCTGTCGAAGCCCGCCTCTGCCAGTCGCCTTTCGAATAACGCCGGGCTGCCAAATGCCGCTTCGAtc                       >  6:388222/2‑149 (MQ=255)
                                                                                                                                     ggCGC‑GTCGACTGCGCTGAGC‑AGGGCGTACACTTGCTGGGCCGCCATCTCATGCCGGGTGTATTC‑AGTGTACTGCGCCCTGTCGAAGCCCGCCTCTGCCAGTCGCCTTTCGAATAACGCCGGGCTGCCAAATGCCGCTTCGATCTCGc                   >  5:388402/1‑148 (MQ=255)
                                                                                                                                           gTCGACTGCGCTGAGC‑AGGGCGTACACTTGCTGGGCCGCCATCTCATGCCGGGTGTATTC‑AGTGTACTGCGCCCTGTCGAAGCCCGCCTCTGCCAGTCGCCTTTCGAATAACGCCGGGCTGCCAAATGCCGCTTCGATCTCGCCAACCt             >  7:229245/1‑149 (MQ=255)
                                                                                                                                            tCGACTGCGCTGAGC‑AGGGCGTACACTTGCTGGGCCGCCATCTCATGCCGGGTGTATTC‑AGTGTACTGCGCCCTGTCGAAGCCCGCCTCTGCCAGTCGCCTTTCGAATAACGCCGGGCTGCCAAATGCCGCTTCGATCTCGCCAACCTg            >  7:386978/1‑149 (MQ=255)
                                                                                                                                            tCGACTGCGCTGAGC‑AGGGCGTACACTTGCTGGGCCGCCATCTCATGCCGGGTGTATTC‑AGTGTACTGCGCCCTGTCGAAGCCCGCCTCTGCCAGTCGCCTTTCGAATAACGCCGAGCTGCCAAATGCCGCTTCGATCTCGCCAACCTg            >  1:381245/1‑149 (MQ=255)
                                                                                                                                               aCTGCGCTGAGC‑AGGGCGTACACTTGCTGGGCCGCCATCTCATGCCGGGTGTATTC‑AGTGTACTGCGCCCTGTCGAAGCCCGCCTCTGCCAGTCGCCTTTCGAATAACGCCGGGCTGCCAAATGCCGCTTCGATCTCGCCAACCTgcgc         <  4:103096/149‑1 (MQ=255)
                                                                                                                                                cTGCGCTGAGC‑AGGGCGTACACTTGCTGGGCCGCCATCTCATGCCGGGTGTATTC‑AGTGTACTGCGCCCTGTCGAAGCCCGCCTCTGCCAGTCGCCTTTCGAATAACGCCGGGCTGCCAAATGCCGCTTCGATCTCGCCAACCt             >  2:26244/1‑144 (MQ=255)
                                                                                                                                                cTGCGCTGAGC‑AGGGCGTACACTTGCTGGGCCGCCATCTCATGCCGGGTGTATTC‑AGTGTACTGCGCCCTGTCGAAGCCCGCCTCTGCCAGTCGCCTTTCGAATAACGCCGGGCTGCCAAATGCCGCTTCGATCTCGCCAACCt             <  1:26244/144‑1 (MQ=255)
                                                                                                                                                 tGCGCTGAGC‑AGGGCGTACACTTGCTGGGCCGCCATCTCATGCCGGGTGTATTC‑AGTGTACTGCGCCCTGTCGAAGCCCGCCTCTGCCAGTCGCCTTTCGAATAACGCCGGGCTGCCAAATGCCGCTTCGATCTCGcc                  >  6:55709/1‑138 (MQ=255)
                                                                                                                                                 tGCGCTGAGC‑AGGGCGTACACTTGCTGGGCCGCCATCTCATGCCGGGTGTATTC‑AGTGTACTGCGCCCTGTCGAAGCCCGCCTCTGCCAGTCGCCTTTCGAATAACGCCGGGCTGCCAAATGCCGCTTCGATCTCGcc                  <  5:55709/138‑1 (MQ=255)
                                                                                                                                                   cgcTGAGC‑AGGGCGTACACTTGCTGGGCCGCCATCTCATGCCGGGTGTATTC‑a                                                                                                     >  2:139309/1‑53 (MQ=255)
                                                                                                                                                   cgcTGAGC‑AGGGCGTACACTTGCTGGGCCGCCATCTCATGCCGGGTGTATTC‑a                                                                                                     <  1:139309/53‑1 (MQ=255)
                                                                                                                                                   cgcTGAGC‑AGGGCGTACACTTGCTGGGCCGCCATCTCATGCCGGGTGTATTC‑AGTGTACTGCGCCCTGTCGAAGCCCGCCTCTGCCAGTCGCCTTTCGAATAACGCCGGGCTGCCAAATGCCGCTTCGATCTCGCCAACCTGCGCCGac     <  3:196683/149‑1 (MQ=255)
                                                                                                                                                      tGAGC‑AGGGCGTACACTTGCTGGGCCGCCATCTCATGCCGGGTGTATTC‑AGTGTACTGCGCCCTGTCGAAGCCCGCCTCTGCCAGTCGCCTTTCGAATAACGCCGGGCTGCCAAATGCCGCTTCGATCTCGCCAACCTGCGCCGACACt  >  6:25992/1‑149 (MQ=255)
                                                                                                                                                      |                                                                                                                                                      
TGCCCGATGAGCGTTGCCCGAGCCAAGGCCAGGCCGTGTTCGCGTATGACCAATGGATTATCACTTTGGTTCTGCGCTCCTTGCAGTCTGTCCTGGTTGGCATCAAAAAATGCCTCCACTT‑CGCCCTCGCTCGGCGC‑GTCGACTGCGCTGAGC‑AGGGCGTACACTTGCTGGGCCGCCATCTCATGCCGGGTGTATTC‑AGTGTACTGCGCCCTGTCGAAGCCCGCCTCTGCCAGTCGCCTTTCGAATAACGCCGGGCTGCCAAATGCCGCTTCGATCTCGCCAACCTGCGCCGACACT  >  NC_002947/3948965‑3949261

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