Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,033,602 T→A 7.3% G73G (GGT→GGA)  PP_0891 → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,033,6020TA 7.3% 207.2 / 7.5 69G73G (GGT→GGA) PP_0891hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (28/35);  new base A (3/2);  total (32/37)
Fisher's exact test for biased strand distribution p-value = 6.53e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.92e-01

CTGGATGGCGAATCTTCACGGCTGTCGTTCATCACCGGCAAGAACGGCGACACCGCCGAAGTGCATCGCTTTCTGGTGTTGCATGGCACGGTCGATCGCAAGGGTATGGCCGGGTTGAGGATCGAGATGGACTCGGTGAGCAGCGGTATTCCGCTGCGTGACG‑AGCAAATGCGCGACAACCTGTTCGAGGTCGAGCGCTTCGCCGAGGCCACAGTGAAGGCACAGATCGACCTGCAGCCGATCAATGACCTGGCCGACGGTGCACAGATCGAATTGCGCCTGCCGCTGACGGT  >  NC_002947/1033456‑1033748
                                                                                                                                                  |                                                                                                                                                   
cTGGATGGCGAATCTTCACGGCTGTCGTTCATCACCGGCAAGAACGGCGACACCGCCGAAGTGCATTGCTTTCTGGTGTTGCATGGCACGGTCGATCGCAAGGGTATGGCCGGGTTGAGGATCGAGATGGACTCGGTGAGCAGCGGGAtt                                                                                                                                                  >  1:343066/1‑150 (MQ=255)
 ctgcggccGAATCTTCACGGCTGTCGTTCATCACCGGCACGAGCGGGGACCCCGCCGAAATGAATCGCTTTCTGGTGTTGCATGGCACGGTCGATCGCAAGGGTATGGCCGGGTTGAGGATCGAGATGGACTCGGTGAGCAGCGGTATTc                                                                                                                                                 <  2:343066/143‑1 (MQ=255)
     tGGCGAATCTTCACGGCTGTCGTTCATCACCGGCAAGAACGGCGACACCGCCGAAGTGCATCGCTTTCTGGTGTTGCATGGCACGGTCGATCGCAAGGGTATGGCCGGGTTGAGGATCGAGATGGACTCGGTGAGCAGCGGTATTCCGCt                                                                                                                                             <  1:222248/150‑1 (MQ=255)
       gCGAATCTTCACGGCTGTCGTTCATCACCGGCAAGAACGGCGACACCGCCGAAGTGCATCGCTTTCTGGTGTTGCATGGCACGGTCGATCGCAAGGGTATGGCCGGGTTGAGGATCGAGATGGACTCGGTGAGCAGCGGTATTCCGCTGc                                                                                                                                           <  2:339361/150‑1 (MQ=255)
         gAATCTTCACGGCTGTCGTTCATCACCGGCAAGAACGGCGACACCGCCGAAGTGCATCGCTTTCTGGTGTTGCATGGCACGGTCGATCGCAAGGGTATGGCCGGGTTGAGGATCGAGATGGACTCGGTGAGCAGCGGTATTCCGCTGCGt                                                                                                                                         <  2:1616955/150‑1 (MQ=255)
         gAATCTTCACGGCTGTCGTTCATCACCGGCAAGAACGGCGACACCGCCGAAGTGCATCGCTTTCTGGTGTTGCATGGCACGGTCGATCGCAAGGGTATGGCCGGGTTGAGGATCGAGATGGACTCGGTGAGCAGCGGAATACCGCTGctc                                                                                                                                         <  1:33836/150‑3 (MQ=255)
            tCTTCACGGCTGTCGTTCATCACCGGCAAGAACGGCGACACCGCCGAAGTGCATCGCTTTCTGGTGTTGCATGGCAGGGTCGATCGCAAGGGTATGGCCGGGTTGAGGATCGAGATGGACTCGGTGAGCAGCGGTATTCCGCTGCGTGAc                                                                                                                                      <  2:388508/150‑1 (MQ=255)
                       gTCGTTCATCACCGGCAAGAACGGCGACACCGCCGAAGTGCATCGCTTTCTGGTGTTGCATGGCACGGTCGATCGCAAGGGTATGGCCGGGTTGAGGATCGAGATGGACTCGGTGAGCAGCGGTATTCCGCTGCGTGACG‑AGCAAAtaga                                                                                                                          >  1:351538/1‑147 (MQ=255)
                          gTTCATCACCGGCAAGAACGGCGACACCGCCGAAGTGCATCGCTTTCTGGTGTTGCATGGCACGGTCGATCGCAAGGGTATGGCCGGGTTGAGGATCGAGATGGACTCGGTGAGCAGCGGTATTCCGCTGCGTGACG‑AGCAAATGCGCGa                                                                                                                       <  1:885238/150‑1 (MQ=255)
                          gTTCATCACCGGCAAAAACGGCGACACCACCGAAGTGCATCGCTTTCTGGTGTTGCATGGCACGGTCGATCGCAAGGGTATGGCCGGGTTGAGGATCGAGATGGACTCGGTGAGCAGCGGTATTCCGCTGCGTGACG‑AGCAAATGCGCGa                                                                                                                       >  2:1106130/1‑150 (MQ=255)
                           ttCATCACCGGCAAGAACGGCGACACCGCCGAAGTGCATCGCTTTCTGGTGTTGCATGGCACGGTCGATCGCAAGGGTATGGCCGGGTTGAGGATCGAGATGGACTCGGTGAGCAGCGGTATTCCGCTGCGTGACG‑AGCAAATGCGCGAc                                                                                                                      >  1:1214118/1‑150 (MQ=255)
                              atcaCCGGCAAGAACGGCGACACCGCCGAAGTGCATCGCTTTCTGGTGTTGCATGGCACGGTCGATCGCAAGGGTATTGCCGGGTTGAGGATCGAGATGGACTCGGTGAGCAGCGGTATTCCGCTGCGTGACG‑AGCAAATGCGCGACAAc                                                                                                                   <  1:1487579/150‑1 (MQ=255)
                               tcaCCGGCAAGAACGGCGACACCGCCGAAGTGCATCGCTTTCTGGTGTTGCATGGCACGGTCGATCGCAAGGGTATGGCCGGGTTGAGGATCGAGATGGACTCGGTGAGCAGCGGTATTCCGCTGCGTGACG‑AGCaaa                                                                                                                              >  2:1436899/1‑138 (MQ=255)
                               tcaCCGGCAAGAACGGCGACACCGCCGAAGTGCATCGCTTTCTGGTGTTGCATGGCACGGTCGATCGCAAGGGTATGGCCGGGTTGAGGATCGAGATGGACTCGGTGAGCAGCGGTATTCCGCTGCGTGACG‑AGCaaa                                                                                                                              <  1:1436899/138‑1 (MQ=255)
                               tcaCCGGCAAGAACGGCGACACCGCCGAAGTGCATCGCTTTCTGGTGTTGCATGGCACGGTCGATCGCAAGGGTATGGCCGGGTTGAGGATCGAGATGGACTCGGTGAGCAGCGGTATTCCGCTGCGTGACG‑AGCAAATGCGCGACAAcc                                                                                                                  >  1:1149391/1‑150 (MQ=255)
                               tcaCCGGCAAGAACGGCGACACCGCCGAAGTGCATCGCTTTCTGGTGTTGCATGGCACGGTCGATCGCAAGGGTATGGCCGGGTTGAGGATCGAGATGGACTCGGTGAGCAGCGGTATTCCGCTGCGCGACT‑AGCAAATGCGCGACAAcc                                                                                                                  >  2:442972/1‑150 (MQ=255)
                                caCCGGCAAGAACGGCGACACCGCCGAAGTGCATCGCTTTCTGGTGTTGCATGGCACGGTCGATCGCAAGGGTATGGCCGGGTTGAGGATCGAGATGGACTCGGTGAGCAGCGGTATTCCGCTGCGTGACG‑AGCAAATGCGCGACAACCt                                                                                                                 >  1:96311/1‑150 (MQ=255)
                                caCCGGCAAAAACGGCGACACCACCGAAGTGCATCGCTTTCTGGTGTTGCATGGCACGGTCGATCGCAAGGGTATGGCCGGGTTGAGGATCGAGATGGACTCGGTGAGCAGCGGTATTCCGCTGCGTGACG‑AGCAAATGCGCGACAACCt                                                                                                                 <  1:1106130/150‑1 (MQ=255)
                                   cGGCAAGAACGGCGACACCGCCGAAGTGCATCGCTTTCTGGTGTTGCATGGCACGGTCGATCGCAAGGGTGGTTCCGGGTTGAGGATCGAGATGGACTCGGTGAGCAGCGGTATTCCGCTGCGTGACG‑AGCAAATACGCGACAACCTGtt                                                                                                              >  1:529134/1‑150 (MQ=255)
                                       aaGAACGGCGACACCGCCGAAGTGCATCGCTTTCTGGTGTTGCATGGCACGGTCGATCGCAAGGGTATGGCCGGGTTGAGGATCGAGATGGACTCGGTGAGCAGCGGTATTCCGCTGCGTGACG‑AGCAAATGCGCGACAACCTGTTCGAg                                                                                                          <  2:690800/150‑1 (MQ=255)
                                        aGAACGGCGACACCGCCGAAGTGCATCGCTTTCTGGTGTTGCATGGCACGGTCGATCGCAAGGGTATGGCCGGGTTGAGGATCGAGATGGACTCGGTGAGCAGCGGTATTCCGCTGCGTGACG‑AGCAAATGCGCGACAACCTGTTCGAgg                                                                                                         <  2:1214118/150‑1 (MQ=255)
                                         gAACGGCGACACCGCCGAAGTGCATCGCTTTCTGGTGTTGCATGGCACGGTCGATCGCAAGGGTATGGCCGGGTTGAGGATCGAGATGGACTCGGTGAGCAGCGGTATTCCGCTGc                                                                                                                                           <  1:367114/116‑1 (MQ=255)
                                         gAACGGCGACACCGCCGAAGTGCATCGCTTTCTGGTGTTGCATGGCACGGTCGATCGCAAGGGTATGGCCGGGTTGAGGATCGAGATGGACTCGGTGAGCAGCGGTATTCCGCTGc                                                                                                                                           >  2:367114/1‑116 (MQ=255)
                                         gAACGGCGACACCGCCGAAGTGCATCGCTTTCTGGTGTTGCATGGCACGGTCGATCGCAAGGGTATGGCCGGGTTGAGGATCGAGATGGACTCGGTGAGCAGCGGTATTCCGCTGCGTGACG‑AGCAAATGCGCGACAAc                                                                                                                   >  1:120087/1‑139 (MQ=255)
                                         gAACGGCGACACCGCCGAAGTGCATCGCTTTCTGGTGTTGCATGGCACGGTCGATCGCAAGGGTATGGCCGGGTTGAGGATCGAGATGGACTCGGTGAGCAGCGGTATTCCGCTGCGTGACG‑AGCAAATGCGCGACAAc                                                                                                                   <  2:120087/139‑1 (MQ=255)
                                           aCGGCGACACCGCCGAAGTGCATCGCTTTCTGGTGTTGCATGGCACGGTCGATCGCAAGGGTATGGCCGGGTTGAGGATCGAGATGGACTCGGTGAGCAGCGGTATTCCGCTGCGTGACG‑AGCAAATGCGCGACAACCTGTTCGAGGTCg                                                                                                      >  1:1048299/1‑150 (MQ=255)
                                                   acCGCCGAAGTGCATCGCTTTCTGGTGTTGCATGGCACGGTCGATCGCAAGGGTATGGCCGGGTTGAGGATCGAGATGGACTCGGTGAGCAGCGGTATTCCGCTGCGTGACG‑AGCAAATGCGCGACAACCTGTTCGAGGTCGAGCGCTTc                                                                                              >  2:797709/1‑150 (MQ=255)
                                                      gccgAAGTGCATCGCTTTCTGGTGTTGCATGGCACGGTCGATCGCAAGGGTATGGCCGGGTTGAGGATCGAGATGGACTCGGTGAGCAGCGGTATTCCGCTGCGTGACG‑AGCAAATGCGCGACAACCTGTTCGAGGTCGAGCGCTTCGcc                                                                                           >  2:342641/1‑150 (MQ=255)
                                                        cgAAGTGCATCGCTTTCTGGTGTTGCATGGCACGGTCGATCGCAAGGGTATGGCCGGGTTGAGGATCGAGATGGACTCGGTGAGCAGCGGTATTCCGCTGCGTGACG‑AGCAAATGCGCGACAACCTGTTCGAGGTCGAGCGCTTCGCCGa                                                                                         <  2:1149391/150‑1 (MQ=255)
                                                               cATCGCTTTCTGGTGTTGCATGGCACGGTCGATCGCAAGGGTATGGCCGGGTTGAGGATCGAGATGGACTCGGTGAGCAGCGGTATTCCGCTGCGt                                                                                                                                         >  2:1373170/1‑96 (MQ=255)
                                                               cATCGCTTTCTGGTGTTGCATGGCACGGTCGATCGCAAGGGTATGGCCGGGTTGAGGATCGAGATGGACTCGGTGAGCAGCGGTATTCCGCTGCGt                                                                                                                                         <  1:1373170/96‑1 (MQ=255)
                                                               cATCGCTTTCTGGTGTTGCATGGCACGGTCGATCGCAAGGGTATGGCCGGGTTGAGGATCGAGATGGACTCGGTGAGCAGCGGTATTCCGCTGCGTGACG‑AGCAAATGCGCGACAACCTGTTCGAGGTCGAGCGCTTCGCCGAGGCcaca                                                                                  <  1:162657/150‑1 (MQ=255)
                                                                 tCGCTTTCTGGTGTTGCATGCCACGGTCGATCGCAAGGGTATGGCCGGGTTGAGGATCGAGATGGACTCGGTGAGCAGCGGTATTCCGCTGCGTGACG‑AGCAAATGCGCGACAACCTGTTCGAGGTCGAGCGCTTCGCCGAGGCCACAGt                                                                                >  2:482271/1‑150 (MQ=255)
                                                                    gtttCGGGTGTTGCATGGCACGGTCGATCGCAAGGGTAGGGCCGGGTTGAGGATCGAGATGGACTCGGTGAGCAGCGGTATTCCGCTGCGTGACG‑AGCAAATGCGCGACAACCTGTTCGAGGTCGAGCGCTTCGCCGAGGCCACAGTGaa                                                                             <  2:546137/149‑1 (MQ=255)
                                                                      ttCTGGTGTTGCATGGCACGGTCGATCGCAAGGGTATGGCCGGGTTGAGGATCGAGATGGACTCGGTGAGCAGCGGTATTCCGCTGCGTGACG‑AGCAAATGCGCGACAACCTGTTCGAGGTCGAGCGCTTCGCCGAGGCCACAGTGAAgg                                                                           <  1:797709/150‑1 (MQ=255)
                                                                       tCTGGTGTTGCATGGCACGGTCGATCGCAAGGGTATGGCCGGGTTGAGGATCGAGATGGACTCGGTGAGCAGCGGTATTCCGCTGCGTGACG‑AGCAAATGCGCGACAACCTGTTCGAGGTCGAGCGCTTCGCCGAGGCCACAGTGAAGGc                                                                          <  2:96311/150‑1 (MQ=255)
                                                                        cTGGTGTTGCATGGCACGGTCGATCGCAAGGGTATGGCCGGGTTGAGGATCGAGATGGACTCGGTGAGCAGCGGTATTCCGCTGCGTGACG‑AGCAAATGCGCGACAACCTGTTCGAGGTCGAGCGCTTCGCCGAGGCCACAGTGAAGGca                                                                         <  2:527984/150‑1 (MQ=255)
                                                                           gtgtTGCATGGCACGGTCGATCGCAAGGGTATGGCCGGGTTGAGGATCGAGATGGACTCGGTGAGCAGCGGTATTCCGCTGCGTGACG‑AGCAAATGCGCGACAACCTGTTCGAGGTCGAGCGCTTCGCCGAGGCCACAGTGAAGGCACAg                                                                      >  1:524724/1‑150 (MQ=255)
                                                                           gtgtTGCATGGCACGGTCGATCGCAAGGGTATGGCCGGGTTGAGGATCGAGATGGACTCGGTGAGCAGCGGTATTCCGCTGCGTGACG‑AGCAAATGCGCGACAACCTGTTCGAGGTCGAGCGCTTCGCCGAGGCCACAGTGAAGGCACAg                                                                      >  2:1411704/1‑150 (MQ=255)
                                                                              ttGCATGGCACGGTCGATCGCAAGGGTATGGCCGGGTTGAGGATCGAGATGGACTCGGTGAGCAGCGGTATTCCGCTGCGTGACG‑AGCAAATGCGCGACAACCTGTTCGAGGTCGAGCGCTTCGCCGAGGCCACAGTGaa                                                                             <  2:1309858/140‑1 (MQ=255)
                                                                              ttGCATGGCACGGTCGATCGCAAGGGTATGGCCGGGTTGAGGATCGAGATGGACTCGGTGAGCAGCGGTATTCCGCTGCGTGACG‑AGCAAATGCGCGACAACCTGTTCGAGGTCGAGCGCTTCGCCGAGGCCACAGTGaa                                                                             >  1:1309858/1‑140 (MQ=255)
                                                                                gCATGGCACGGTCGATCGCAAGGGTATGGCCGGGTTGAGGATCGAGATGGACTCGGTGAGCAGCGGTATTCCGCTGCGTGACG‑AGCAAATGCGCGACAACCTGTTCGAGGTCGAGCGCTTCGCCGAGGCCACAGTGAAGGCACAGATCGa                                                                 <  1:342641/150‑1 (MQ=255)
                                                                                   tGGCACGGTCGATCGCAAGGGTATGGCCGGGTTGAGGATCGAGATGGACTCGGTGAGCAGCGGTATTCCGCTGCGTGACG‑AGCAAATGCGCGACAACCTGTTCGAGGTCGAGCGCTTCGCCGAGGCCACAGTGAAGGCACAGATCGACCt                                                              <  2:623484/150‑1 (MQ=255)
                                                                                       aCGGTCGATCGCGAGGGTATGGCCGGGTTGAGGATCGAGATGGACTCGGTGAGCAGCGGTATTCCGCTGCGTGACG‑AGCAAATGCGCGACAACCTGTTCGAGGTCGAGCGCTTCGCCGAGGCCACAGTGAAGGCACAGATCGACCTGCAg                                                          >  1:536193/1‑150 (MQ=255)
                                                                                       aCGGTCGATCGCAAGGGTATGGCCGGGTTGAGGATCGAGATGGACTCGGTGAGCAGCGGTATTCCGCTGCGTGACG‑AGCAAATGCGCGACAACCTAAAGCGGGTCGATCGCTTCGCCGAGGCCACAGTGAAGGCACAGATCGACCTGCAg                                                          <  1:1350614/150‑1 (MQ=255)
                                                                                        cGGTCGATCGCAAGGGTATGGCCGGGTTGAGGATCGAGATGGACTCGGTGAGCAGCGGTATTCCGCTGCGTGACG‑AGCAAATGCGCGACAACCTGTTCGAGGTCGAGCGCTTCGCCGAGGCCACAGTGAAGGCACAGATCGACCTGCAGc                                                         <  2:529134/150‑1 (MQ=255)
                                                                                           tCGAGCGCAAGGGTATGGCCGGGTTGAGGATCGAGATCGACTCGGTGAGAACCGCTATTCCGCTGCG‑AACGCTGCAAATGCGCGACAACCTGTTCGAGGTCGAGCGCTTCGCCGAGGCCACAGTGAAGGCACAGATCGACCTGCAGCCGa                                                      <  1:442972/150‑1 (MQ=255)
                                                                                             gATCGCAAGGGTATGGCCGGGTTGAGGATCGAGATGGACTCGGTGAGCAGCGGTATTCCGCTGCGTGACG‑AGCAAATGCGCGACAACCTGTTCGAGGTCGAGCGCTTCGCCGAGGCCACAGTGAAGGCACAGATCGACCTGCAGCCGATc                                                    >  2:1405392/1‑150 (MQ=255)
                                                                                             gATCGCAAGGGTATGGCCGGGTTGAGGATCGAGATGGACTCGGTGAGCAGCGGAATACCg                                                                                                                                               <  1:1628208/60‑1 (MQ=255)
                                                                                             gATCGCAAGGGTATGGCCGGGTTGAGGATCGAGATGGACTCGGTGAGCAGCGGAATACCg                                                                                                                                               >  2:1628208/1‑60 (MQ=255)
                                                                                                cGCAAGGGTATGGCCGGGTTGAGGATCGAGATGGACTCGGTGAGCAGCGGTATTCCGCTGCGTGACG‑AGCAAATGCGCGACAACCTGTTCGAGGTCGAGCGCTTCGCCGAGGCCACAGTGAAGGCACAGATCGACCTGCAGCCGATCAAt                                                 <  2:524724/150‑1 (MQ=255)
                                                                                                cGCAAGGGTATGGCCGGGTTGAGGATCGAGATGGACTCGGTGAGCAGCGGTATTCCGCTGCGTGACG‑AGCAAATGCGCGACAACCTGTTCGAGGTCGAGCGCTTCGCCGAGGCCACAGTGAAGGCACAGATCGACCTGCAGCCGATCAAt                                                 <  1:482271/150‑1 (MQ=255)
                                                                                                 gCAAGGGTATGGCCGGGTTGAGGATCGAGATGGACTCGGTGAGCAGCGGTATTCCGCTGCGTGACG‑AGCAAATGCGCGACAACCTGTTCGAGGTCGAGCGCTTCGCCGAGGCCACAGTGAAGGCACAGATCGAc                                                                <  1:295735/134‑1 (MQ=255)
                                                                                                 gCAAGGGTATGGCCGGGTTGAGGATCGAGATGGACTCGGTGAGCAGCGGTATTCCGCTGCGTGACG‑AGCAAATGCGCGACAACCTGTTCGAGGTCGAGCGCTTCGCCGAGGCCACAGTGAAGGCACAGATCGAc                                                                >  2:295735/1‑134 (MQ=255)
                                                                                                       gTATGGCCGGGTTGAGGATCGAGATGGACTCGGTGAGCAGCGGTATTCCGCTGCGTGACG‑AGCAAATGCGCGACAACCTGTTCGAGGTCGAGCGCTTCGCCGAGGCCACAGTGAAGGCACAGATCGACCTGCAGCCGATCAATGACCTgg                                          >  1:419934/1‑150 (MQ=255)
                                                                                                                ggTTGAGGATCGAGATGGACTCGGTGAGCAGCGGTATTCCGCTGCGTGACG‑AGCAAATGCGCGACAACCTGTTCGAGGTCGAGCGCTTCGCCGAGGCCACAGTGAAGGCACAGATCGACCTGCAGCCGATCAATGACCTGGCCGACGGTg                                 <  1:1411704/150‑1 (MQ=255)
                                                                                                                ggTTGAGGATCGAGATGGACTCGGTGAGCAGCGGTATTCCGCTGCGTGACG‑AGCAAATGCGCGACAACCTGTTCGAGGTCGAGCGCTTCGCCGAGGCCACAGTGAAGGCACAGATCGACCTGCAGCCGATCAATGACCTGGCCGACGGTg                                 <  1:1405392/150‑1 (MQ=255)
                                                                                                                 gTTGAGGATCGAGATGGACTCGGTGAGCAGCGGTATTCCGCTGCGTGACG‑AGCAAATGCGCGACAACCTGTTCGAGGTCGAGCGCTTCGCCGAGGCCACAGTGAAGGCACAGATCGACCTGCAGCCGATCAATGACCTGGCCGACGGTGc                                >  2:773145/1‑150 (MQ=255)
                                                                                                                        aTCGAGATGGACTCGGTGAGCAGCGGTATTCCGCTGCGTGACG‑AGCAAATGCGCGACAACCTGTTCGAGGTCGAGCGCTTCGCCGAGGCCACAGTGAAGGCACAGATCGACCTGCAGCCGATCAATGACCTGGCCGACGGTGCACAGATc                         >  1:729890/1‑150 (MQ=255)
                                                                                                                        aTCGAGATGGACTCGGTGAGCAGCGGTATTCCGCTGCGTGACG‑AGCAAATGCGCGACAACCTGTTCGAGGTCGAGCGCTTCGCCGAGGCCACAGTGAAGGCACAGATCGACCTGCAGCCGATCAATGACCTGGCCGACGGTGCACAGATc                         <  2:536193/150‑1 (MQ=255)
                                                                                                                           gagaTGGACTCGGTGAGCAGCGGTATTCCGCTGCGTGACG‑AGCAAATGCGCGACAACCTGTTCGAGGTCGAGCGCTTCGCCGAGGCCACAGTGAAGGCACAGATCGACCTGCAGCCGATCAATGACCTGGCCGACGGTGCACAGATCGaa                      <  2:1048299/150‑1 (MQ=255)
                                                                                                                              aTGGACTCGGTGAGCAGCGGTATTCCGCTGCGTGACG‑AGCAAATGCGCGACAACCTGTTCGAGGTCGAGCGCTTCGCCGAGGCCACAGTGAAGGCACAGATCGACCTGCAGCCGATCAATGACCTGGCCGACGGTGCACAGATCGAATTg                   >  1:1107122/1‑150 (MQ=255)
                                                                                                                              aTGGACTCGGTGAGCAGCGGTATTCCGCTGCGTGACG‑AGCAAATGCGCGACAACCTGTTCGAGGTCGAGCGCTTCGCCGAGGCCACAGTGAAGGCACAGATCGACCTGCAGCCGATCAATGACCTGGCCGACGGTGCACAGATCGAATTg                   >  1:1058381/1‑150 (MQ=255)
                                                                                                                              aTGGACTCGGTGAGCAGCGGTATTCCGCTGCGTAACG‑AGCAAATGCGCGACAACCTGTTCGAGGTCGAGCGCTTCGCCGAGGCCACAGTGAAg                                                                            >  1:1023689/1‑93 (MQ=255)
                                                                                                                              aTGGACTCGGTGAGCAGCGGTATTCCGCTGCGTAACG‑AGCAAATGCGCGACAACCTGTTCGAGGTCGAGCGCTTCGCCGAGGCCACAGTGAAg                                                                            <  2:1023689/93‑1 (MQ=255)
                                                                                                                                      gtcacgcagcGGAATACCGCTGCGTGACG‑AGCAAATGCGCGACAACCTGTTCGAGGTCGAGCGCTTCGCCGAGGCCACAGTGAAGGCACAGATCGACCTGCAGCCGATCAATGACCTGGCCGACGGTGCACAGATCGAATTGCGCCTGcc           >  1:1272003/6‑150 (MQ=255)
                                                                                                                                         acgcagcGGAATACCGCTGCGTGACG‑AGCAAATGCGCGACAACCTGTTCGAGGTCGAGCGCTTCGCCGAGGCCACAGTGAAGGCACAGATCGACCTGCAGCCGATCAATGACCTGGCCGACGGTGCACAGATCGAATTGCGCCTGCCGCt        >  1:481318/3‑150 (MQ=255)
                                                                                                                                             agcGGTATTCCGCTGCGTGACG‑AGCAACTGCAAGCCTTAAACATTGAAGTCTAGCTTTTCGCCGAGGCCCCCGCTTTGGCACAGATCGACCTGCAGCCGATCAATGACCTGGCCGACGGTGCACAGATCGAATTGCGCCTGCCGCTGACg    <  1:773145/150‑1 (MQ=255)
                                                                                                                                               cGGTATTCCGCTGCGTGACG‑AGCAAATGCGCGACAACCTGTTCGAGGTCGAGCGCTTCGCCGAGGCCACAGTGAAGGCACAGATCGACCTGCAGCCGATCAATGACCTGGCCGACGGTGCACAGATCGAATTGCGCCTGCCGCTGACGGt  >  2:779852/1‑150 (MQ=255)
                                                                                                                                                  |                                                                                                                                                   
CTGGATGGCGAATCTTCACGGCTGTCGTTCATCACCGGCAAGAACGGCGACACCGCCGAAGTGCATCGCTTTCTGGTGTTGCATGGCACGGTCGATCGCAAGGGTATGGCCGGGTTGAGGATCGAGATGGACTCGGTGAGCAGCGGTATTCCGCTGCGTGACG‑AGCAAATGCGCGACAACCTGTTCGAGGTCGAGCGCTTCGCCGAGGCCACAGTGAAGGCACAGATCGACCTGCAGCCGATCAATGACCTGGCCGACGGTGCACAGATCGAATTGCGCCTGCCGCTGACGGT  >  NC_002947/1033456‑1033748

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