Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,348,666:1 +C 6.1% intergenic (‑31/+443) PP_3820 ← / ← galU group II intron‑encoding maturase/UTP‑glucose‑1‑phosphate uridylyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,348,6661.C 6.1% 186.1 / 8.7 65intergenic (‑31/+443)PP_3820/galUgroup II intron‑encoding maturase/UTP‑glucose‑1‑phosphate uridylyltransferase
Reads supporting (aligned to +/- strand):  ref base . (32/29);  new base C (2/2);  total (34/31)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.82e-01

TCTGCATTACGGTTACCAACGAGACGGCGACGCCTACTGGCGGCATAACCTGTTCGGCCCTTGGTGGCGCGGTTAAC‑‑CGCCACTTACTACGGCTTCGGCTGACTTCTGCACGCTCATCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTCGCGCGACCTTCCTGCTTATGCCTGTCGGATCTACGTCACAGCGTTCC  >  NC_002947/4348590‑4348814
                                                                             |                                                                                                                                                     
tCTGCATTACGGTTACCAACGAGACGGCGACGCCTACTGGCGGCATAACCTGTTCGGCCCTTGGTGGCGCGGTTAAC‑‑CGCCACTTACTACGGCTTCGGCTGACTTCTGCACGCTCGTCCCGTCGCCTTTCGACGCGCGGTAGCACAGt                                                                               <  2:88771/148‑1 (MQ=14)
 cagCATTACGGTTACCAACGAGACGGCGACGCCTACTGGCGGCATAACCTGTTCGGCCCTTGGTGGCGCGGTTAAC‑‑CGCCACTTACTACGGCTTCGGCTGACTTCTGCACGCTCGTCCCGTCGCCTTTCGACGCGCGGTAGCACAGTg                                                                              <  8:296230/146‑1 (MQ=17)
 cTGCATTACGGTTACCAACGAGACGGCGACGCCTACTGGCGGCATAACCTGTTCGGCCCTTGGTGGCGCGGTTAAC‑‑CGCCACTTACTACGTCTTCGGCTGACTTCTGCACGCTCGTCCCGTCGCCTTTCGACGCGCGGTAGCACAGTg                                                                              <  1:290202/148‑1 (MQ=17)
 cTGCATTACGGTTACCAACGAGACGGCGACGCCTACTGGCGGCATAACCTGTTCGGCCCTTGGTGGCGCGGTTAAC‑‑CGCCACTTACTACGGCTTCGGCTGACTTCTGCACGCTCATCCCGTCGCCTTTCGACGCGCGGTAGCACAGTg                                                                              <  5:87320/148‑1 (MQ=32)
    cATTACGGTTACCAACGAGACGGCGACGCCTACTGGCGGCATAACCTGTTCGGCCCTTGGTGGCGCGGTTAAC‑‑CGCCACTTACTACGGCTTCGGCTGACTTCTGCACGCTCATCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCaa                                                                          <  4:36011/149‑1 (MQ=32)
     aTTACGGTTACCAACGAGACGGCGACGCCTACTGGCGGCATAACCTGTTCGGCCCTTGGTGGCGCGGTTAAC‑‑CGCCACTTACTACGGCTTCGGCTGACTTCTGCACGCTCATCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAg                                                                         >  3:151186/1‑149 (MQ=32)
      ttACGGTTACCAACGAGACGGCGACGCCTACTGGCGGCATAACCTGTTCGGCCCTTGGTGGCGCGGTTAAC‑‑CGCCACTTACTACGGCTTCGGCTGACTTCTGCACGCTCATCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGc                                                                        <  4:231501/149‑1 (MQ=33)
       tACGGTTACCAACGAGACGGCGACGCCTACTGGCGGCATAACCTGTTCGGCCCTTGGTGGCGCGGTTAAC‑‑CGCAACTTCCTACGGCTTCGGCTGACTTCTGCACGCTCGTCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCg                                                                       >  5:265027/1‑149 (MQ=255)
         cGGTTACCAACGAGACGGCGACGCCTACTGGCGGCATAACCTGTTCGGCCCTTGGTGGCGCGGTTAAC‑‑CGCCACTTACTACGGCTTCGGCTGACTTCTGCACGCTCATCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGt                                                                      <  5:124250/148‑1 (MQ=34)
         cGGTTACCAACGAGACGGCGACGCCTACTGGCGGCATAACCTGTTCGGCCCTTGGTGGCGCGGTTAAC‑‑CGCCACTTACTACGGCTTCGGCTGACTTCTGCACGCTCATCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTa                                                                     <  8:93447/149‑1 (MQ=34)
           gTTACCAACGAGACGGCGACGCCTACTGGCGGCATAACCTGTTCGGCCCTTGGTGGCGCGGTTAAC‑‑CGCCACTTACTACGGCTTCGGCTGACTTCTGCACGCTCGTCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACa                                                                   <  4:265449/149‑1 (MQ=21)
           gTTACCAACGAGACGGCGACGCCTACTGGCGGCATAACCTGTTCGGCCCTTGGTGGCGCGGTTAAC‑‑CGCCACTTACTACGGCTTCGGCTGACTTCTGCACGCTCATCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTa                                                                     >  5:115088/1‑147 (MQ=34)
                cAACGAGACGGCGACGCCTACTGGCGGCATACCCTGTTCGGCCCTTGGTGGCGCGGTTAAC‑‑CGCCACGTACTCCGGCGTCGGCGGACTTCTGCACGCTCATCCCGTCGCCGTTCGACGCGCGGGAGCACCGTGGCAAGCGTCCCGAtct                                                              >  6:348911/1‑149 (MQ=11)
                cAACGAGACGGCGACGCCTACTGGCGGCATAACCTGTTCGGCCCTTGGTGGCGCGGTTAAC‑‑CGCCACTTACTACGGCTTCGGCTGACTTCTGCACGCTCGTCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGAtct                                                              >  8:332748/1‑149 (MQ=21)
                  aCGAGACGGCGACGCCTACTGGCGGCATAACCTGTTCGGCCCTTGGTGGCGCGGTTAAC‑‑CGCCACTTACTACGGCTTCGGCTGACTTCTGCACGCTCATCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTcc                                                            >  4:322812/1‑149 (MQ=34)
                  aCGAGACGGCGACGCCTACTGGCGGCATACCTGTTCGGCCCTTGGTGGCGCGGTTAAC‑‑CGCCACTTACTACGGCTTCGGCTGACTTCTGCACGCTCGTCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTccc                                                           <  8:344879/149‑1 (MQ=19)
                     agaCGGCGACGCCTACTGGCGGCATAACCTGTTCGGCCCTTGGTGGCGCGGTTAAC‑‑CGCCACTTACTACGGCTTCGGCTGACTTCTGCACGCTCATCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAg                                                         <  4:337681/149‑1 (MQ=34)
                      gaCGGCGACGCCTACTGGCGGCATAACCTGTTCGGCCCTTGGTGGCGCGGTTAAC‑‑CGCCACTTACTACGGCTTCGGCTGACTTCTGCACGCTCATCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAgg                                                        >  2:342588/1‑149 (MQ=34)
                           cGACGCCTACTGGCGGCATAACCTGTTCGGCCCTTGGTGGCGCGGTTAAC‑‑CGCCACTTACTACGGCTTCGGCTGACTTCTGCACGCTCGTCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAAt                                                   >  3:162143/1‑149 (MQ=21)
                            gACGCCTACTGGCGGCATAACCTGTTCGGCCCTTGGTGGCGCGGTTAAC‑‑CGCCACTTACTACGGCTTCGGCTGACTTCTGCACGCTCGTCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAAtt                                                  >  3:270527/1‑149 (MQ=21)
                                 cTACTGGCGGCATAACCTGTTCGGCCCTTGGTGGCGCGGTTAAC‑‑CGCCACTTACTACGGCTTCGGCTGACTTCTGCACGCTCATCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTcgcgc                                             >  1:328284/1‑149 (MQ=34)
                                 cTACTGGCGGCATAACCTGTTCGGCCCTTGGTGGCGCGGTTAAC‑‑CGCCACTTACTACGGCTTAGGCTGACTTCTGCACGCTCGTCCCGTCGACTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTcgcgc                                             >  2:167764/1‑149 (MQ=255)
                                  tACTGGCGGTATAACCTGTTCGGCCCTTGGTGGCGCGGTTAAC‑‑CGCCACTTACTACGGCTTCGGCTGACTTCTGCACGCTCGTCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTcgcgc                                             >  5:129988/1‑148 (MQ=21)
                                    cTGGCGGCATAACCTGTTCGGCCCTTGGTGGCGCGGTTAAC‑‑CGCCACTTACTACGGCTTCGGCTGACTTCTGCACGCTCATCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTCGCGCGAc                                          >  3:135326/1‑149 (MQ=34)
                                     tGGCGGCATAACCTGTTCGGCCCTTGGTGGCGCGGTTAAC‑‑CGCCACTTACTACGGCTTCGGCTGACTTCTGCACGCTCATCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTCGCGCGAcc                                         >  6:216074/1‑149 (MQ=34)
                                         ggcATAACCTGTTCGGCCCTTGGTGGCGCGGTTAAC‑‑CGCCACTTACTACGGCTTCGGCTGACTTCTGCACGCTCGTCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTCGCGCGACCTTcc                                     <  3:308015/149‑1 (MQ=21)
                                         ggcATAACCTGTTCGGCCCTTGGTGGCGCGGTTAAC‑‑CGCCACTTACTACGGCTTCGGCTGACTTCTGCACGCTCATCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTCGCGCGACCTTcc                                     <  4:135326/149‑1 (MQ=34)
                                          gcATAACCTGTTCGGCCCTTGGTGGCGCGGTTAAC‑‑CGCCACTTACTACGGCTTCGGCTGACTTCTGCACGCTCATCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTCGCGCGACCTTCCt                                    >  7:79024/1‑149 (MQ=34)
                                            aTAACCTGTTCGGCCCTTGGTGGCGCGGTTAAC‑‑CGCCACTTACTACGGCTTCGGCTGACTTCTGCACGCTCATCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTCGCGCGACCTTCCTGc                                  <  1:342588/149‑1 (MQ=34)
                                            aTAACCTGTTCGGCCCTTGGTGGCGCGGTTAAC‑‑CGCCACTTACTACGGCTTCGGCTGACTTCTGCACGCTCATCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTCGCGCGACCTTCCTGc                                  >  8:39579/1‑149 (MQ=34)
                                             tAACCTGTTCGGCCCTTGGTGGCGCGGTTAAC‑‑CGCCACTTACTACGGCTTCGGCTGACTTCTGCACGCTCATCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAAtt                                                  <  6:308004/132‑1 (MQ=34)
                                             tAACCTGTTCGGCCCTTGGTGGCGCGGTTAAC‑‑CGCCACTTACTACGGCTTCGGCTGACTTCTGCACGCTCATCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAAtt                                                  >  5:308004/1‑132 (MQ=34)
                                              aaCCTGTTCGGCCCTTGGTGGCGCGGTTAAC‑‑CGCCACTTACTACGGCTTCGGCTGACTTCTGCACGCTCATCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTCGCGCGACCTTCCTGCt                                 >  8:222618/1‑148 (MQ=34)
                                              aaCCTGTTCGGCCCTTGGTGGCGCGGTTAAC‑‑CGCCACTTACTACGGCGTCGGCTGACTTCTGCACGCTCATCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTCATTCGCGCGACCTTCCTGCt                                 >  8:216921/1‑148 (MQ=21)
                                               aCCTGTTCGGCCCTTGGTGGCGCGGTTAAC‑‑CGCCACTTACTACGGCTTCGGCTGACTTCTGCACGCTCATCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTCGCGCGACCTTCCTGCTTa                               >  7:326883/1‑149 (MQ=34)
                                               aCCTGTTCGGCCCTTGGTGGCGCGGTTAAC‑‑CGCCACTTACTACGGCTTCGGCTGACTTCTGCACGCTCATCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTCGCGCGACCTTCCTGCTTa                               >  2:187877/1‑149 (MQ=34)
                                                ccTGTTCGGCCCTTGGTGGCGCGGTTAAC‑‑CGCTACTTACTACGGCTTCGGCTGACTTCTGCACGCTCGTCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTCGCGCGACCTTCCTGCTTAt                              <  6:265027/149‑1 (MQ=21)
                                                ccTGTTCGGCCCTTGGTGGCGCGGTTAAC‑‑CGCCACTTACTACGGCTTCGGCTGACTTCTGCACGCTCGTCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTCGCGCGACCTTCCTGCTTAt                              >  1:252800/1‑149 (MQ=21)
                                                 cTGTTCGGCCCTTGGTGGCGCGGTTAAC‑‑CGCCACTTACTACGGCTTCGGCTGACTTCTGCACGCTCGTCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCa                                                          >  3:142186/1‑120 (MQ=25)
                                                 cTGTTCGGCCCTTGGTGGCGCGGTTAAC‑‑CGCCACTTACTACGGCTTCGGCTGACTTCTGCACGCTCGTCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCa                                                          <  4:142186/120‑1 (MQ=25)
                                                 cTGTTCGGCCCTTGGTGGCGCGGTTAAC‑‑CGCCACTTACTACGGCTTCGGCTGACTTCTGCACGCTCATCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTCGCGCGACCTTCCTGCTTATg                             <  5:216074/149‑1 (MQ=34)
                                                   gTTCGGCCCTTGGTGGCGCGGTTAAC‑‑CGCCACTTACTACGGCTTCGGCTGACTTCTGCACGCTCGTCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTCGCGCGACCTTCCTGCTTATGcc                           >  8:11059/1‑149 (MQ=21)
                                                    ttCGGCCCTTGGTGGCGCGGTTAAC‑‑CGCCACTTACTACGGCTTCGGCTGACTTCTGCACGCTCGTCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTcgcg                                              <  7:349397/129‑1 (MQ=25)
                                                    ttCGGCCCTTGGTGGCGCGGTTAAC‑‑CGCCACTTACTACGGCTTCGGCTGACTTCTGCACGCTCGTCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTcgcg                                              >  8:349397/1‑129 (MQ=25)
                                                    ttCGGCCCTTGGTGGCGCGGTTAAC‑‑CGCCACTTACTACGGCTTCGGCTGACTTCTGCACGCTCGTCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTCGCGCGACCTTCCTGCTTATGc                            <  4:270527/147‑1 (MQ=21)
                                                       ggCCCTTGGGGGCGCGGTTAAC‑‑CGCCACTGACTACCGCTTCGGCTGACATCTGCACGCTCATCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTCGCGCGACCTTCCTGCTTATGCCTGTc                       <  4:109085/149‑1 (MQ=255)
                                                        gCCCTTGGTGGCGCGGTTAAC‑‑CGCCACTTACTACGGCTTCGGCTGACTTCTGCACGCTCGTCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTCGCGCGACCTTCCTGCTTATGCCTGTCg                      >  7:143780/1‑149 (MQ=21)
                                                         cccTTGGTGGCGCGGTTAAC‑‑CGCCACTTACTACGGCTTCGGCTGACTTCTGCACGCTCATCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTCGCGCGACCTTCCTGCTTATGCCTGTCgg                     <  8:79024/149‑1 (MQ=34)
                                                          ccTTGGTGGCGCGGTTAAC‑‑CGCCACTTACTACGGCTTCGGCTGACTTCTGCACGCTCATCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTCGCGCGACCTTCCTGCTTATGCCTGTCGGa                    >  8:99370/1‑149 (MQ=34)
                                                               gtggCGCGGTTAAC‑‑CGCCACTTACTACGGCTTCGGCTGACTTCTGCACGCTCATCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTCGCGCGACCTTCCTGCTTATGCCTGTCGGATCTAc               <  5:154646/149‑1 (MQ=34)
                                                               gtggCGCGGTTAAC‑‑CGCCACTTACTACGGCTTCGGCTGACTTCTGCACGCTCATCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTCGCGCGACCTTCCTGCTTATGCCTGTCGGATCTAc               <  6:115088/149‑1 (MQ=34)
                                                                 gtggcgGTTAACCGCGCCACTTACTACGGCTTCGGCTGACTTCTGCACGCTCATCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTCGCGCGACCTTCCTGCTTATGCCTGTCGGAt                   <  3:74991/142‑1 (MQ=255)
                                                                 gtggcgGTTAACCGCGCCACTTACTACGGCTTCGGCTGACTTCTGCACGCTCATCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTCGCGCGACCTTCCTGCTTATGCCTGTCGGAt                   >  4:74991/4‑145 (MQ=255)
                                                                 gtggcgGTTAACCGCGCCACTTACTACGGCTTCGGCTGACTTCTGCACGCTCATCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTCGCGCGACCTTCCTGCTTATGCCTGTCGGATCTAc               >  2:254847/4‑149 (MQ=255)
                                                                  gcgcGGTTAAC‑‑CGCCACTTACTACGGCTTCGGCTGACTTCTGCACGCTCGTCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTCGCGCGACCTTCCTGCTTATGCCTGTCGGATCTACGt             <  3:303141/148‑1 (MQ=21)
                                                                  gcgcGGTTAC‑‑CGCCACTTACTACGGCTTCGGCTGACTTCTGCACGCTCATCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTCGCGCGACCTTCCTGCTTATGCCTGTCGGATCTACGTca           >  8:155309/1‑149 (MQ=21)
                                                                    gcgGTTAACCGCGCCACTTACTACGGCTTCGGCTGACTTCTGCACGCTCATCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTCGCGCGACCTTCCTGCTTATGCCTGTCGGATCTACGTc            <  1:254847/149‑1 (MQ=255)
                                                                    gcgGTTAAC‑‑CGCCACTTACTACGGCTTCGGCTGACTTCTGCACGCTCATCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTCGCGCGACCTTc                                      <  2:292317/121‑1 (MQ=35)
                                                                    gcgGTTAAC‑‑CGCCACTTACTACGGCTTCGGCTGACTTCTGCACGCTCATCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTCGCGCGACCTTc                                      >  1:292317/1‑121 (MQ=35)
                                                                    gcgGTTAAC‑‑CGCCACTTACTACGGCTTCGGCTGACTTCTGCACGCTCATCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTCGCGCGACCTTCCTGCTTATGCCTGTCGGATCTACGTca           >  4:45579/1‑148 (MQ=34)
                                                                    gcgGTTAAC‑‑CGCCACTTACTACGGCTTCGGCTGACTTCTGCACGCTCATCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTCGCGCGACCTTCCTGCTTATGCCTGTCGGATCTACGTca           >  7:34914/1‑148 (MQ=34)
                                                                     cgGTTAAC‑‑CGCCACTTACTACGGCTTCGGCTGACTTCTGCACGCTCGTCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTCGCGCGACCTTCCTGCTTATGCCTGTCGGATCTACGTcaca         <  4:174485/149‑1 (MQ=21)
                                                                     cgGTTAAC‑‑CGCCACTTACTACGGCTTCGGCTGACTTCTGCACGCTCGTCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTCGCGCGACCTTCCTGCTTATGCCTGTCGGATCTACGTcaca         <  4:162143/149‑1 (MQ=21)
                                                                     cgGTTAAC‑‑CGCCACTTACTACGGCTTCGGCTGACTTCTGCACGCTCATCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTCGCGCGACCTTCCTGCTTATGCCTGTCGGATCTACGTcaca         <  2:196676/149‑1 (MQ=34)
                                                                          aaC‑‑CGCCACTTACTACGGCTTCGGCTGACTTCTGCACGCTCGTCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTCGCGCGACCTTCCTGCTTATGCCTGTCGGATCTACGTCACAGCGtt    >  3:279105/1‑149 (MQ=21)
                                                                            c‑‑cGCCACTTACTACGGCTTCGGCTGACTTCTGCACGCTCGTCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAGGCGTACAGATCTCCCAGGGTAATTCGCGCGACCTTCCTGCTTATGCCTGTCGGATCTACGTCACAGCGTTcc  <  7:93420/149‑1 (MQ=21)
                                                                            c‑‑cGCCACTTACTACGGCTTCGGCTGACTTCTGCACGCTCATCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTCGCGCGACCTTCCTGCTTATGCCTGTCGGATCTACGTCACAGCGTTcc  >  3:190239/1‑149 (MQ=34)
                                                                             |                                                                                                                                                     
TCTGCATTACGGTTACCAACGAGACGGCGACGCCTACTGGCGGCATAACCTGTTCGGCCCTTGGTGGCGCGGTTAAC‑‑CGCCACTTACTACGGCTTCGGCTGACTTCTGCACGCTCATCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTCGCGCGACCTTCCTGCTTATGCCTGTCGGATCTACGTCACAGCGTTCC  >  NC_002947/4348590‑4348814

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