Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,725,088:1 +C 5.1% intergenic (+27/+41) PP_5661 → / ← PP_4183 putative lipoprotein/membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,725,0881.C 5.1% 334.6 / 10.7 78intergenic (+27/+41)PP_5661/PP_4183putative lipoprotein/membrane protein
Reads supporting (aligned to +/- strand):  ref base . (36/38);  new base C (2/2);  total (38/40)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.32e-01

GAAGCCGCCGAGGAAAACGCCGAAGCTGCCAAAGATCAGGCCGAAGCCGAGCAGAAGGCTGCCGAAGAAGCCGCTCCGGCTACTCCGGCTCCTACTACACCGGCAGAACCTGCCAAGCAGTAATTGGCTAAAACAAAAAAACCCGACGATGTCGGGTTTTTTTGTGCCTGCGAACTTTCGATTAATCGTTTCAGTTGGAAACAACCGAGCTTTCCTTGGCCGGCTCCGACAGTACACTTTCGGTCGGTGTGAACACCATCACATCGAGTACATCCGAATGGAACTCTCGACGGTA  >  NC_002947/4724939‑4725233
                                                                                                                                                      |                                                                                                                                                 
gAAGCCGCCGAGGAAAACGCCGAAGCTGCCAAAGATCAGGCCGAAGCCGAGCAGAAGGCTGCCGAAGAAGCCGCTCCGGCTACTCCGGCTCCTACTACACCGGCAGAACCTGCCAAGCAGTAATTGGCTAAAACAAAAAAACCCGACGAt                                                                                                                                                    >  2:1598030/1‑150 (MQ=255)
gAAGCCGCCGAGGAAAACGCCGAAGCTGCCAAAGATCAGGCCGAAGCCGAGCAGAAGGCTGCCGAAGAAGCCGCTCCGGCTACTCCGGCTCCTACTAAAACGGCAGAACCTGCCAAGCAGTAATTGGCTAAAACAAAAAAACCCGACGAt                                                                                                                                                    <  1:1806442/150‑1 (MQ=255)
    cggccgAGGAAAACGCCGAAGATGCCAAAGATCAGGCCAAAGCCGAGCAGAAGGCTGCCGAAGAAGCCGCTCCGGCTACGACGGCTCCTACTACACCGGCAGAACCTGCCAAGCAGTAATTGGCTAAAACAAAAAAACCCGACGATGTCg                                                                                                                                               <  2:160247/148‑1 (MQ=255)
     cgccgAGGAAAACGCCGAAGCTGCCAAAGATCAGGCCGAAGCCGAGCAGAAGGCTGCCGAAGAAGCCGCTCCGGCTACTCCGGCTCCTACTACACCGGCAGAACCTGCCAAGCAGTAATTGGCTAAAACAAAAAAACCCGACGATGTCgg                                                                                                                                              <  1:467338/150‑1 (MQ=255)
         gAGGAAAACGCCGAAGCTGCCAAAGATCAGGCCGAAGCCGAGCAGAAGGCTGCCGAAGAAGCCGCTCCGGCTACTCCGGCTCCTACTACACCGGCAGAACCTGCCAAGCAGTAATTGGCTAAAACAAAAAAACCCGACGATGTCGGGttt                                                                                                                                          >  1:1337160/1‑150 (MQ=255)
            gAAAACGCCGAAGCTGCCAAAGATCAGGCCGAAGCCGAGCAGAAGGCTGCCGAAGAAGCCGCTCCGGCTACTCCGGCTCCTACTACACCGGCAGAACCTGCCAAGCAGTAATTGGCTAAAACAAAAAAACCCGACGATGTCGGGtttttt                                                                                                                                       >  1:1912937/1‑150 (MQ=255)
            gAAAACGCCGAAGCTGCCAAAGATCAGGCCGAAGCCGAGCAGAAGGCTGCCGAAGAAGCCGCTCCGGCTACTCCGGCTCCTACTACACCGGCAGAACCTGCCAAGCAGTAATTGGCTAAAACAAAAAAACCCGACATCGTCGGGtttttt                                                                                                                                       <  1:1396500/150‑1 (MQ=255)
             aaaaCGCCGAAGCTGCCAAAGATCAGGCCGAAGCCGAGCAGAAGGCTGCCGAAGAAGCCGCTCCGGCTACTCCGGCTCCTACTACACCGGCAGAACCTGCCAAGCAGTAATTGGCTAAAACAAAAAAACCCGACGATGTCGGGttttttt                                                                                                                                      <  1:450199/150‑1 (MQ=255)
             aaaaCGCCGAAGCTGCCAAAGATCAGGCCGAAGCCGAGCAGAAGGCTGCCGAAGAAGCCGCTCCGGCTACTCCGGCTCCTACTACACCGGCAGAACCTGCCAAGCAGTAATTGGCTAAAACAAAAAAACCCGACGATGTCGGGttttttt                                                                                                                                      <  2:1912937/150‑1 (MQ=255)
             aaaaCGCCGAAGCTGCCAAAGATCAGGCCGAAGCCGAGCAGAAGGCTGCCGAAGAAGCCGCTCCGGCTACTCCGGCTCCTACTACACCGGCAGAACCTGCCAAGCAGTAATTGGCTAAAACAAAAAAACCCGACGATGTCGGGttttttt                                                                                                                                      <  1:1100279/150‑1 (MQ=255)
              aaaCGCCGAAGCTGCCAAAGATCAGGCCGAAGCCGAGCAGAAGGCTGCCGAAGAAGCCGCTCCGGCTACTCCGGCTCCTACTACACCGGCAGAACCTGCCAAGCAGTAATTGGCTAAAACAAAAAAACCCGACGATGTCGGGTTTTTTtg                                                                                                                                     <  2:231344/150‑1 (MQ=255)
                aCGCCGAAGCTGCCAAAGATCAGGCCGAAGCCGAGCAGAAGGCTGCCGAAGAAGCCGCTCCGGCTACTCCGGCTCCTACTACACCGGCAGAACCTGCCAAGCAGTAATTGGCTAAAACAAAAAAACCCGACGATGTCGGGTTTTTTtgtg                                                                                                                                   >  2:1499632/1‑150 (MQ=255)
                         cTGCCAAAGATCAGGCCGAAGCCGAGCAGAAGGCTGCCGAAGAAGCCGCTCCGGCTACTCCGGCTCCTACTACACCGGCAGAACCTGCCAAGCAGTAATTGGCTAAAACAAAAAAACCCGACGATGTCGGGTTTTTTTGTGCCTGCGAAc                                                                                                                          <  2:1312348/150‑1 (MQ=255)
                             cAAAGATCAGGCCGAAGCCGAGCAGAAGGCTGCCGAAGAAGCCGCTCCGGCTACTCCGGCTCCTACTACACCGGCAGAACCTGCCAAGCAGTAATTGGCTAAAACAAAAAAACCCGACGATGTCGGGTTTTTTTGTGCCTGCGAACTTTc                                                                                                                      >  1:1515505/1‑150 (MQ=255)
                              aaaGATCAGGCCGAAGCCGAGCAGAAGGCTGCCGAAGAAGCCGCTCCGGCTACTCCGGCTCCTACTACACCGGCAGAACCTGCCAAGCAGTAATTGGCTAAAACAAAAAAACCCGACGATGTCGGGttttttt                                                                                                                                      <  2:1099390/133‑1 (MQ=255)
                              aaaGATCAGGCCGAAGCCGAGCAGAAGGCTGCCGAAGAAGCCGCTCCGGCTACTCCGGCTCCTACTACACCGGCAGAACCTGCCAAGCAGTAATTGGCTAAAACAAAAAAACCCGACGATGTCGGGttttttt                                                                                                                                      >  1:1099390/1‑133 (MQ=255)
                                  aTCAGGCCGAAGCCGAGCAGAAGGCTGCCGAAGAAGCCGCTCCGGCTACTCCGGCTCCTACTACACCGGCAGAACCTGCCAAGCAGTAATTGGCTAAAACAAAAAAACCCGACGATGTCGGGTTTTTTTGTGCCTGCGAACTTTCGATTa                                                                                                                 >  2:1528969/1‑150 (MQ=255)
                                    cAGGCCGAAGCCGAGCAGAAGGCTGCCGAAGAAGCCGCCCCGGCTACTCCGGCTCCTACTACACCGGCAGAACCTGCCAAGCAGTAATTGGCTAAAACAAAAAAACCCGACGATGTCGGGTTTTTTTGTGCCTGCGAACTTTCGATTAAt                                                                                                               <  2:1323440/150‑1 (MQ=255)
                                       gCCGAAGCCGAGCAGAAGGCTGCCGAAGGAGCCGTTCCGGCTACTCCGGCTCCTACTACACCGGCAGAACCTGCCAAGCAGTAATTGGCTAAAACAAAAAAACCCGACGATGTCGGGTTTTTTTGTGCCTGCGAACTTTCGATTAATCGt                                                                                                            <  2:413541/150‑1 (MQ=255)
                                         cGAAGCCGAGCAGAAGGCTGCCGAAGAAGCCGCTCCGGCTACTCCGGCTCCTACTACACCGGCAGAACCTGCCAAGCAGTAATTGGCTAAAACAAAAAAACCCGACGATGTCGGGTTTTTTTGTGCCTGCGAACTTTCGATTAATCGttt                                                                                                          <  1:1330982/150‑1 (MQ=255)
                                           aaGCCGAGCAGAAGGCTGCCGAAGAAGCCGCTCCGGCTACTCCGGCTCCTACTACACCGGCAGAACCTGCCAAGCAGTAATTGGCTAAAACAAAAAAACCCGACGATGTCGGGTTTTTTTGTGCCTGCGAACTTTCGATTAATCGTTTCa                                                                                                        >  1:1062929/1‑150 (MQ=255)
                                             gCCGAGCAGAAGGCTGCCGAAGAAGCCGCTCCGGCTACTCCGGCTCCTACTACACCGGCAGAACCTGCCAAGCAGTAATTGGCTAAAACAAAAAAACCCGACGATGTCGGGTTTTTTTGTGCCTGCGAACTTTCGATTAATCGTTTCAGt                                                                                                      <  1:1285231/150‑1 (MQ=255)
                                                gAGCAGAAGGCTGCCGAAGAAGCCGCTCCGGCTACTCCGGCTCCTACTACACCGGCAGAACCTGCCAAGCAGTAATTGGCTAAAACAAAAAAACCCGACGATGTCGGGTTTTTTTGTGCCTGCGAACTTTCGATTAATCGTTTCAGTTgg                                                                                                   >  2:1806422/1‑150 (MQ=255)
                                                gAGCAGAAGGCTGCCGAAAAAGCCGCTCCGGCTACTCCGGCTCCTACTACACCGGCAGAACCTGCCAAGCAATAATTGGCTAAAACAAAAAAACCCGACGATGTCGGGTTTTTTGTGCCTGCGAACTTTCGa                                                                                                                    >  1:1197427/1‑132 (MQ=255)
                                                gAGCAGAAGGCTGCCGAAAAAGCCGCTCCGGCTACTCCGGCTCCTACTACACCGGCAGAACCTGCCAAGCAATAATTGGCTAAAACAAAAAAACCCGACGATGTCGGGTTTTTTGTGCCTGCGAACTTTCGa                                                                                                                    <  2:1197427/132‑1 (MQ=255)
                                                      aaGGCTGCCGAAGAAGCCGCTCCGGCTACTCCGGCTCCTACTACACCGGCAGAACCTGCCAAGCAGTAATTGGCTAAAACAAAAAAACCCGACGATGTCGGGTTTTTTTGTGCCTGCGAACTTTCGATTAATCGTTTCAGTTGGAaacaa                                                                                             <  2:1496417/150‑1 (MQ=255)
                                                        ggCTGCCGAAGAAGCCGCTCCGGCTACTCCGGCTCCTACTACACCGGCAGAACCTGCCAAGCAGTAATTGGCTAAAACAAAAAAACCCGACGATGTCGGGTTTTTTTGTGCCTGCGAACTTTCGATTAATCGTTTCAGTTGGAAACAAcc                                                                                           <  2:1062929/150‑1 (MQ=255)
                                                        ggCTGCCGAAGAAGCCGCTCCGGCTACTCCGGCCCCCATCACACCGGCAGAACCTGCCAAGCAGTAATTGGCTAAAACAAAAAAACCCGACGATGTCGGGTTTTTTTGTGCCTGCGAACTTTCGATTAATCGTTTCAGTTGGAAACAAcc                                                                                           <  1:101575/150‑1 (MQ=255)
                                                           tGCCGAAGAAGCCGCTCCGGCTACTCCGGCTCCTACTACACCGGCAGAACCTGCCAAGCAGTAATTGGCTAAAACAAAAAAACCCGACGATGTCGGGTTTTTTTGTGCCTGCGAACTTTCGATTAATCGTTTCAGTTGGAAACAACCGa                                                                                         >  2:557749/1‑149 (MQ=255)
                                                           tGCCGAAGAAGCCGCTCCGGCTACTCCGGCTCCTACTACACCGGCAGAACCTGCCAAGCAGTAATTGGCTAAAACAAAAAAACCCGACGATGTCGGGTTTTTTTGTGCCTGCGAACTTTCGATTAATCGTTTCAGTTGGAAACAACCGa                                                                                         <  1:557749/149‑1 (MQ=255)
                                                           tGCCGAAGAAGCCGCTCCGGCTACTCCGGCTCCTACTACACCGGCAGAACCTGCCAAGCAGTAATTGGCTAAAACAAAAAAACCCGACGATGTCGGGTTTTTTTGTGCCTGCGAACTTTCGATTAATCGTTTCAGTTGGAAACAACCGAg                                                                                        <  2:1337160/150‑1 (MQ=255)
                                                           tGCCGAAGAAGCCGCTCCGGCTACTCCGGCTCCTACTACACCGGCAGAACCTGCCAAGCAGTAATTGGCTAAAACAAAAAAACCCGACGATGTCGGGTTTTTTTGTGCCTGCGAACTTTCGATTAATCGTTTCAGTTGGAAACAACCGAg                                                                                        >  1:699555/1‑150 (MQ=255)
                                                            gCCGAAGAAGCCGCTCCGGCTACTCCGGCTCCTACTACACCGGCAGAACCTGCCAAGCAGTAATTGGCTAAAAAAAAAAAACCCGACGATGTCGGGAATCTTCGGTCGACCGGATTGGCTGTGAATCGTTTCAGTTGGAAACAACCGAGc                                                                                       >  1:1399100/1‑150 (MQ=255)
                                                               gaagaaGCCGCTCCGGCTACTCCGGCTCCTACTACACCGGCAGAACCTGCCAAGCAGTAATTGGCTAAAACAAAAAAACCCGACGATGTCGGGTTTTTTTGTGCCTGCGAACTTTCGATTAATCGTTTCAGTTGGAAACAACCGAGCttt                                                                                    <  1:1621189/150‑1 (MQ=255)
                                                               gaagaaGCCGCTCCGGCTACTCCGGCTCCTACTACACCGGCAGAACCTGCCAAGCAGTAATTGGCTAAAACAAAAAAACCCGACATCGTCGGGTTTTTTtgttt                                                                                                                                  >  2:972488/1‑102 (MQ=38)
                                                               gaagaaGCCGCTCCGGCTACTCCGGCTCCTACTACACCGGCAGAACCTGCCAAGCAGTAATTGGCTAAAACAAAAAAACCCGACATCGTCGGGTTTTTTtgttt                                                                                                                                  <  1:972488/104‑3 (MQ=38)
                                                                   aagCCGCTCCGGCTACTCCGGCTCCTACTACACCGGCAGAACCTGCCAAGCAGTAATTGGCTAAAACAAAAAAACCCGACGATGTCGGGTTTTTTTGTGCCTGCGAACTTTCGATTAATCGTTTCAGTTGGAAACAACCGAGCTTTCCtt                                                                                <  1:1170389/150‑1 (MQ=255)
                                                                    agCCGCTCCGGCTACTCCGGCTCCTACTAAACCGGCAGAACCTGCCAAGCAGTAATTGGCTAAAACAAAAAAACCCGACGATGTCGGGTTTTTTTGTGCCTGCGAACTTTCGATTAATCGTTTCAGTTGGAAACAACCGAGCTTTCCTTg                                                                               >  2:583009/1‑150 (MQ=255)
                                                                     gccgccCCGGCTACTCCGGCTCCTACTACACCGGCAGAACCTGCCAAGCAGTAATTGGCTAAAACAAAAAAACCCGACGATGTCg                                                                                                                                               >  1:1230472/1‑85 (MQ=255)
                                                                     gccgccCCGGCTACTCCGGCTCCTACTACACCGGCAGAACCTGCCAAGCAGTAATTGGCTAAAACAAAAAAACCCGACGATGTCg                                                                                                                                               <  2:1230472/85‑1 (MQ=255)
                                                                       cGCTCCGGCTACTCCGGCTCCTACTACACCGGCAGAACCTGCCAAGCAGTAATTGGCTAAAACAAAAAAACCCGACGATGTCGGGTTTTTTTGTGCCTGCGAACTTTCGATTAATCGTTTCAGtt                                                                                                     <  2:1559828/125‑1 (MQ=255)
                                                                       cGCTCCGGCTACTCCGGCTCCTACTACACCGGCAGAACCTGCCAAGCAGTAATTGGCTAAAACAAAAAAACCCGACGATGTCGGGTTTTTTTGTGCCTGCGAACTTTCGATTAATCGTTTCAGtt                                                                                                     >  1:1559828/1‑125 (MQ=255)
                                                                       cGCTCCGGCTACTCCGGCTCCTACTACACCGGCAGAACCTGCCAAGCAGTAATTGGCTAAAACAAAAAAACCCGACGATGTCGGGTTTTTTTGTGCCTGCGAACTTTCGATTAATCGTTTCAGTTGGAAACAACCGAGCTTTCCTTGGcc                                                                            <  2:1515505/150‑1 (MQ=255)
                                                                         cTCCGGCTACTCCGGCTCCTACTACACCGGCAGAACCTGCCAAGCAGTAATTGGCTAAAACAAAAAAACCCGACGATGTCGGGTTTTTTTGTGCCTGCGAACTTTCGATTAATCGTTTCAGTTGGAAACAACAGAGCTTTCCTTGGCCgg                                                                          >  2:1109241/1‑150 (MQ=255)
                                                                               cTACTCCGGCTCCTACTACACCGGCAGAACCTGCCAAGCAGTAATTGGCTAAAACAAAAAAACCCGACGATGTCGGGTTTTTTTGTGCCTGCGAACTTTCGATTAATCGTTTCAGTTGGAAACAACCGAGCTTTCCTTGGCCGGCTCCGa                                                                    <  2:1777826/150‑1 (MQ=255)
                                                                                 gcTCCGGGTGCTACTACACCGGCAGAACCTGCGAAGCAGTAATTGGCTAAAACAAAAAAACCCGACGATGTCGGGTTTTTTTGTGCCTGCGAACTTTCGATTAATCGTTTCAGTTGGAAACAACCGAGCTTTCCTTGGCCGGCTCCGACa                                                                  <  2:1535031/149‑1 (MQ=255)
                                                                                 aCTCCGGCTCCTACTACACCGGCAGAACCTGCCAAGCAGTAATTGGCTAAAACAAAAAAACCCGACGATGTCGGGTTTTTTTGTGCCTGCGAACTTTCGATTAATCGTTTCAGTTGGAAACAACCGAGCTTTCCTTGTCCGGCTCCGACa                                                                  <  2:699555/150‑1 (MQ=255)
                                                                                   tCCGGCTCCTACTACACCGGCAGAACCTGCCAAGCAGTAATTGGCTAAAACAAAAAAACCCGACGATGTCGGGTTTTTTTGTGCCTGCGAACTTTCGATTAATCGTTTCAGTTGGAAACAACCGAGCTTTCCTTGGCCGGCTCCGACAGt                                                                <  1:1598030/150‑1 (MQ=255)
                                                                                     cGGCTCCTACTACACCGGCAGAACCTGCCAAGCAGTAATTGGCTAAAACAAAAAAACCCGACGATGTCGGGTTTTTTTGTGCCTGCGAACTTTCGATTAATCGTTTCAGTTGGAAACAACCGAGCTTTCCTTGGCCGGCTCCGACAGTac                                                              >  2:762668/1‑150 (MQ=255)
                                                                                     cGGCTCCTACTACACCGGCAGAACCTGCCAAGCAGTAATTGGCTAAAACAAAAAAACCCGACGATGTCGGGTTTTTTTGTGCCTGCGAACTTTCGATTAATCGTTTCAGTTGGAAACAACCGAGCTTTCCTTGGCCGGCTCCGACAGTac                                                              >  1:1430242/1‑150 (MQ=255)
                                                                                           ctactaCACCGGCAGAACCTGCCAAGCAGTAATTGGCTAAAACAAAAAAACCCGACGATGTCGGGTTTTTTTGTGCCTGCGAACTTTCGATTAATCGTTTCAGTTGGAAACAACCGAGCTTTCCTTGGCCGGCTCCGACAGTACACTTtc                                                        <  1:762668/150‑1 (MQ=255)
                                                                                              ctacACCGGCAGAACCTGCCAAGCAGTAATTGGCTAAAACAAAAAAACCCGACGATGTCGGGTTTTTTTGTGCCTGCGAACTTTCGATTAATCGTTTCAGTTGGAAACAACCGAGCTTTCCTTGGCCGGCTCCGACAGTACACTTtcggt                                                     >  1:370797/1‑150 (MQ=255)
                                                                                               tacACCGCCGGAACCTCCCCCGGAGTACTTGGCTAAGACAAAAAGAGCCTGCGGGGTCTGGATGTTTTTTGCCTGCGTACGTCCGATTAATCGTTTCAGTTGTAAACAACCGAGCTTTACTTGGCAGGCCCCGGCACTACTCTTtcgctc                                                    >  2:1801681/1‑150 (MQ=255)
                                                                                                acacCGGCAGAACCTGCCAAGCAGTAATTGGCTAAAACAAAAAAACCCGACGATGTCGGGTTTTTTTGTGCCTGCGAACTTTCGATTAATCGTTTCAGTTGGAAACAACCGAGCTTTCCTTGGCCGGCTCCGACAGTACACTTtcggtcg                                                   >  1:1550365/1‑150 (MQ=255)
                                                                                                  acCGGCAGAACCTGCCAAGCAGTAATTGGCTAAAACAAAAAAACCCGACGATGTCGGGTTTTTTTGTGCCTGCGAACTTTCGATTAATCGTTTCAGTTGGAAACAACCGAGCTTTCCTTGGCCGGCTCCGACAGTACACTTTcggtcggt                                                 >  1:1589005/1‑150 (MQ=255)
                                                                                                  acCGGCAGAACCTGCCAAGCAGTAATTGGCTAAAACAAAAAAACCCGACGATGTCGGGTTTTTTTGTGCCTGCGAACTTTCGATTAATCGTTTCAGTTGGAAACAACCGAGCTTTCCTTGGCCGGCTCCGACAGTACACTTTcggtcggt                                                 >  2:1730535/1‑150 (MQ=255)
                                                                                                   ccGGCAGAACCTGCCAAGCAGTAATTGGCTAAAACAAAAAAACCCGACGATGTCGGGTTTTTTTGTGCCTGCGAACTTTCGATTAATCGTTTCAGTTGGAAACAACCGAGCTTTCCTTGGCCGGCTCCGACAGTACACTTTCGGTCGgtg                                                >  1:250071/1‑150 (MQ=255)
                                                                                                    cGGCAGAACCTGCCAAGCAGTAATTGGCTAAAACAAAAAAACCCGACGATGTCGGGTTTTTTTGTGCCTGCGAACTTTCGATTAATCGTTTCAGTTGGAAACAACCGAGCTTTCCTTGGCCGGCTCCGACAGTACACTTTCGGTCGgtgt                                               >  2:1409066/1‑150 (MQ=255)
                                                                                                     ggCAGAACCTGCCAAGCAGTAATTGGCTAAAACAAAAAAACCCGACGATGTCGGGTTTTTTTGTGCCTGCGAACTTTCAATTAATCGTTTCAGTTGGAAACAACCGAGCATTCCTTGGCTGGCTCCGACAGTACACTTTCGGTCGGtgtg                                              <  1:1499632/150‑1 (MQ=255)
                                                                                                      gCAGAACCTGCCAAGCAGTAATTGGCTAAAACAAAAAAACCCGACGATGTCGGGTTTTTCTGTGCCTGCGAACTTTCGATTAATCGTTTCAGTTGGAAACAACCGAGCTTTCCTTGGCCGGCTCCGACAGTACACTTTcggtcggt                                                 <  2:1102917/146‑1 (MQ=255)
                                                                                                      gCAGAACCTGCCAAGCAGTAATTGGCTAAAACAAAAAAACCCGACGATGTCGGGTTTTTCTGTGCCTGCGAACTTTCGATTAATCGTTTCAGTTGGAAACAACCGAGCTTTCCTTGGCCGGCTCCGACAGTACACTTTcggtcggt                                                 >  1:1102917/1‑146 (MQ=255)
                                                                                                          aaCCTGCCAAGCAGTAATTGGCTAAAACAAAAAAACCCGACGATGTCGGGTTTTTTTGTGCCTGCGAACTTTCGATTAATCGTTTCAGTTGGAAACAACCGAGCTTTCCTTGGCCGGCTCCGACAGTACACTTTCGGTCGGTGTGAacac                                         >  2:1907068/1‑150 (MQ=255)
                                                                                                            ccTGCCAAGCAGTAATTGGCTAAAACAAAAAAACCCGACGATGTCGGGTTTTTTTGTGCCTGCGAACTTTCGATTAATCGTTTCAGTTGGAAACAACCGAGCTTTCCTTGGCCGGCTCCGACAGTACACTTTCGGTCGATGTGatggtgt                                       <  1:1109241/150‑7 (MQ=255)
                                                                                                             cTGCCAAGCAGTAATTGGCTAAAACAAAAAAACCCGACGATGTCGGGTTTTTTTGTGCCTGCAAACTATCGATTAATCGTTTCAGTAGGAAACAACCGTGCGTGTCGTGGTCGGCTCCGACAGTACACTTTCGGTCGGTGTGAACACCAt                                      >  2:707647/1‑150 (MQ=255)
                                                                                                                 cAAGCAGTAATTGGCTAAAACAAAAAAACCCGACGATGTCGGGTTTTTTTGTGCCTGCGAACTTTCGATTAATCGTTTCAGTTGGAAACAACCGAGCTTTCCTTGGCCGGCTCCGACAGTACACTTTCGGTCGGTGTGAACACCATcaca                                  >  1:1254100/1‑150 (MQ=255)
                                                                                                                  aaGCAGTAATTGGCTAAAACAAAAAAACCCGACGATGTCGGGTTTTTTTGTGCCTGCGAACTTTCGATTAATCGTTTCAGTTGGAAACAACCGAGCTTTCCTTGGCCGGCTCCGACAGTACACTTTCGGTCGGTGTGAACACCATCACAt                                 >  2:1623802/1‑150 (MQ=255)
                                                                                                                     cAGTAATTGGCTAAAACAAAAAAACCCGACGATGTCGGGTTTTTTTGTGCCTGCGAACTTTCGATTAATCGTTTCAGTTGGAAACAACCGAGCTTTCCTTGGCCGGCTCCGACAGTACACTTTCGGTCGGTGTGAACACCATCACATCGa                              <  1:1222288/150‑1 (MQ=255)
                                                                                                                      aGTAATTGGCTAAAACAAAAAAACCCGACGATGTCGGGTTTTTTTGTGCCTGCGAACTTTCGATTAATCGTTTCAGTTGGAAACAACCGAGCTTTCCTTGGCCGGCTCCGACAGTACACTTTCGGTCGGTGTGAACACCATCACATCGAg                             <  2:1589005/150‑1 (MQ=255)
                                                                                                                           ttGGCTAAAACAAAAAAACCCGACGATGTCGGGTTTTTTTGTGCCTGCGAACTTTCGATTAATCGTTTCAGTTGGAAACAACCGAGCTTTCCTCGGCCGGCTCCGACAGTACACTTTCGGTCGGTGTGAACACCATCACATCGAGTACAt                        <  2:1254100/150‑1 (MQ=255)
                                                                                                                              gCTAAAACAAAAAAACCCGACGATGTCGGGTTTTTTTGTGCCTGCGAACTTTCGATTAATCGTTTCAGTTGGAAACAACCGAGCTTTCCTTGGCCGGCTCCGACAGTACACTTTCGGTCGGTGTGAACACCATCACATCGAGTACATCCg                     >  2:1649777/1‑150 (MQ=255)
                                                                                                                                tAAAACAAAAAAACCCGACGATGTCGGGTTTTTTTGTGCCTGCGAACTTTCGATTAATCGTTTCAGTTGGAAACAACCGAGCTTTCCTTGGCCGGCTCCGACAGTACACTTTCGGTCGGTGTGAACACCATCACATCGAGTACATCCGaa                   <  1:216247/150‑1 (MQ=255)
                                                                                                                                       aaaaaaCCCGACGATGTCGGGTTTTTTTGTGCCTGCGAACTTTCGATTAATCGTTTCAGTTGGAAACAACCGAGCTTTCCTTGGCCGGCTCCGACAGTACACTTTCGGTCGGTGTGAACACCATCACATCGAGTACATCCGAATGGAAct            >  1:753673/1‑150 (MQ=255)
                                                                                                                                        aaaaaCCCGACGATGTCGGGTTTTTTTGTGCCTGCGAACTTTCGATTAATCGTTTCAGTTGGAAACAACCGAGCTTTCCTTGGCCGGCTCCGACAGTACACTTTCGGTCGGTGTGAACACCATCACATCGAGTACATCCGAATGGAActc           >  2:1630260/1‑150 (MQ=255)
                                                                                                                                        aaaaaCCCGACATCGTCGGGTTTTTTTGTGCCTGCGAACTTTCGATTAATCGTTTCAGTTGGAAACAACCGAGCTTTCCTTGGCCGGCTCCGACAGTACACTTTCGGTCGGTGTGAACACCATCACATCGa                              <  1:127357/131‑1 (MQ=255)
                                                                                                                                        aaaaaCCCGACATCGTCGGGTTTTTTTGTGCCTGCGAACTTTCGATTAATCGTTTCAGTTGGAAACAACCGAGCTTTCCTTGGCCGGCTCCGACAGTACACTTTCGGTCGGTGTGAACACCATCACATCGa                              >  2:127357/1‑131 (MQ=255)
                                                                                                                                         aaaaCCCGACGATGTCGGGTTTTTTTGTGCCTGCGAACTTTCGGTTAATCGTTTCAGTTGGAAACAACCGAGCTTTCCTTGGCCGGCTCCGACAGTACACTTTCGGTCGGTGTGAACACCATCACATCGAGTACATCCGAATGGAActct          >  1:1559945/1‑150 (MQ=255)
                                                                                                                                         aaaaCCCGACATCGTCGGGTTTTTTTGTGCCTGCGAACTTTCGATTAATCGTTTCAGTTGGAAACAACCGAGCTTTCCTTGGCCGGCTCCGACAGTACACTTTCGGTCGGTGTGAACACCATCACATCGAGTACATCCGAATGGAActct          >  2:69535/1‑150 (MQ=255)
                                                                                                                                          aaaCCCGACGATGTCGGGTTTTTTTGTGCCTGCGAACTTTCGATTAATCGTTTCAGTTGGAAACAACCGAGCTTTCCTTGGCCGGCTCCGACAGTACACTTTCGGTCGGTGTGAACACCATCACATCGAGTACATCCGAATGGAACtctc         <  2:1314313/150‑1 (MQ=255)
                                                                                                                                              ccGACGATGTCGGGTTTTTTTTTTCCTGCGAACTTTCGATTAATCGTTTCAGTTGGAAACAACCGAGCTTTCCTTGGCCGGCTCCGACAGTACACTTTCGGTCGGTGTGAACACCATCACATCGAGTACATCCGAATGGAACTCTCGACg     >  2:1504041/1‑150 (MQ=255)
                                                                                                                                                 acgaTGTCGGGTTTTTTTGTGCCTGCGAACTTTCGATTAATCTTTTCATGTGGAAACAACCGAGCTTTCCTTGGCCGGCTCCGACAGTACACTTTCGGTCGGTGTGAACACCATCACATCGAGTACATCCGAATGGAACTCTCGACGGTa  <  1:916509/150‑1 (MQ=255)
                                                                                                                                                      |                                                                                                                                                 
GAAGCCGCCGAGGAAAACGCCGAAGCTGCCAAAGATCAGGCCGAAGCCGAGCAGAAGGCTGCCGAAGAAGCCGCTCCGGCTACTCCGGCTCCTACTACACCGGCAGAACCTGCCAAGCAGTAATTGGCTAAAACAAAAAAACCCGACGATGTCGGGTTTTTTTGTGCCTGCGAACTTTCGATTAATCGTTTCAGTTGGAAACAACCGAGCTTTCCTTGGCCGGCTCCGACAGTACACTTTCGGTCGGTGTGAACACCATCACATCGAGTACATCCGAATGGAACTCTCGACGGTA  >  NC_002947/4724939‑4725233

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: