Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 702,765 A→T 12.8% intergenic (+13/‑119) PP_23SD → / → PP_5SD 23S ribosomal RNA/5S ribosomal RNA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947702,7650AT12.8% 181.8 / 16.2 70intergenic (+13/‑119)PP_23SD/PP_5SD23S ribosomal RNA/5S ribosomal RNA
Reads supporting (aligned to +/- strand):  ref base A (34/27);  new base T (2/7);  total (36/34)
Fisher's exact test for biased strand distribution p-value = 8.02e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.00e-01

CTGGGATCTTGAATCCCCTAAAGGGCCGTCGAAGACTACGACGTTGATAGGTTGGGTGTGTAAGCGCTGTGAGGCGTTGAGCTAACCAATACTAATTGCCCGTGAGGCTTGACCATATAACACCCAAGCAATCTGCTACGCAGATTGCGGTGGTGAAGACGAAAGACCCGAAAATTCGTGCTTTACACAAGAGCTTCGCAAATATCGCATATCCGAATTCGCTGGGCTGTCCATCTGGACATTCTGGCAACAGAATTTCTTGACGACCATAGAGCATTGGAAC  >  NC_002947/702628‑702910
                                                                                                                                         |                                                                                                                                                 
cTGGGATCTTGAATCCCCTAAAGGGCCGTCGAAGACTACGACGTTGATAGGTTGGGTGTGTAAGCGCTGTGAGGCGTTGAGCTAACCAATACTAATTGCCCGTGAGGCTTGACCATATAACACCCAAGCAATCTGCTACGCAGATTGCg                                                                                                                                        <  6:45708/149‑1 (MQ=31)
cTGGGATCTTGAATCCCCTAAAGGGCCGTCGAAGACTACGACGTTGATAGGTTGGGTGTGTAAGCGCTGTGAGGCGTTGAGCTAACCAATACTAATTGCCCGTGAGGCTTGACCATATAACACCCAAGCAATCTGCGTAGCAGATTGct                                                                                                                                        <  2:146256/149‑2 (MQ=2)
 tGGGATCTTGAATCCCCTAAAGGGCCGTCGAAGACTACGACGTTGATAGGTTGGGTGTGTAAGCGCTGTGAGGCGTTGAGCTAACCAATACTAATTGCCCGTGAGGCTTGACCATATAACACCCAAGCAATCTGCTACGCAGATTGCgg                                                                                                                                       >  8:82407/1‑149 (MQ=32)
 tGGGATCTTGAATCCCCTAAAGGGCCGTCGAAGACTACGACGTTGATAGGTTGGGTGTGTAAGCGCTGTGAGGCGTTGAGCTAACCAATACTAATTGCCCGTGAGGCTTGACCATATAACACCCAAGCAATCTGCGTAGCAGATTGctt                                                                                                                                       <  3:4794/149‑3 (MQ=2)
   ggATCTTGAATCCCCTAAAGGGCCGTCGAAGACTACGACGTTGATAGGTTGGGTGTGTAAGCGCTGTGAGGCGTTGAGCTAACCAATACTAATTGCCCGTGAGGCTTGACCATATAACACCCAAGCAATCTGCTACGCAGATTGCggt                                                                                                                                      <  5:329075/148‑1 (MQ=32)
   ggATCTTGAATCCCCTAAAGGGCCGTCGAAGACTACGACGTTGATAGGTTGGGTGTGTAAGCGCTGTGAGGCGTTGAGCTAACCAATACTAATTGCCCGTGAGGCTTGACCATATAACACCCAAGCAATCTGCGTAGCAGATTGcttgg                                                                                                                                     <  8:191710/149‑5 (MQ=2)
   ggATCTTGAATCCCCTAAAGGGCCGTCGAAGACTACGACGTTGATAGGTTGGGTGTGTAAGCGCTGTGAGGCGTTGAGCTAACCAATACTAATTGCCCGTGAGGCTTGACCATATAACACCCAAGCAATCTGCGTAGCAGATTGcttg                                                                                                                                      <  1:119526/148‑4 (MQ=2)
   ggATCTTGAATCCCCTAAAGGGCCGTCGAAGACTACGACGTTGATAGGTTGGGTGTGTAAGCGCTGTGAGGCGTTGAGCTAACCAATACTAATTGCCCGTGAGGCTTGACCATATAACACCCAAGCAATCTGCGTAGCAGATTGcttg                                                                                                                                      <  7:109475/148‑4 (MQ=2)
   ggATCTTGAATCCCCTAAAGGGCCGTCGAAGACTACGACGTTGATAGGTTGGGTGTGTAAGCGCTGTGAGGCGTTGAGCTAACCAATACTAATTGCCCGTGAGGCTTGACCATATAACACCCAAGCAATCTGCGTAGCAGATTGcttg                                                                                                                                      <  5:245710/148‑4 (MQ=2)
   ggATCTTGAATCCCCTAAAGGGCCGTCGAAGACTACGACGTTGATAGGTTGGGTGTGTAAGCGCTGTGAGGCGTTGAGCTAACCAATACTAATTGCCCGTGAGGCTTGACCATATAACACCCAAGCAATCTGCGTAGCAGATTGcttg                                                                                                                                      <  3:149155/148‑4 (MQ=2)
        ttGAATCCCCTAAAGGGCCGTCGAAGACTACGACGTTGATAGGTTGGGTGTGTAAGCGCTGTGAGGCGTTGAGCTAACCAATACTAATTGCCCGTGGGGCTTGACCATATAACACCCAAGCAATCTGCTACGCAGATTGCGgtggtg                                                                                                                                  >  2:213231/1‑147 (MQ=17)
        ttGAATCCCCTAAAGGGCCGTCGAAGACTACGACGTTGATAGGTTGGGTGTGTAAGCGCTGTGAGGCGTTGAGCTAACCAATACTAATTGCCCGTGAGGCTTGACCATATAACACCCAAGCAATCTGCTACGCAGATTGCGgtggtg                                                                                                                                  >  4:74475/1‑147 (MQ=32)
          gAATCCCCTAAAGGGCCGTCGAAGACTACGACGTTGATAGGTTGGGTGTGTAAGCGCTGTGAGGCGTTGAGCTAACCAATACTAATTGCCCGTGAGGCTTGACCATATAACACCCAAGCAATCTGCTACGCAGATTGCGgtggtg                                                                                                                                  <  5:326716/145‑1 (MQ=32)
                ccTAAAGGGCCGTCGAAGACTACGACGTTGATAGGTTGGGTGTGTAAGCGCTGTGAGGCGTTGAGCTAACCAATACTAATTGCCCGTGAGGCTTGACCATATAACACCCAAGCAATCTGCTACGCAGATTGCGGTGGTGAAGACGa                                                                                                                           >  2:304059/1‑146 (MQ=32)
                 cTAAAGGGCCGTCGAAGACTACGACGTTGATAGGTTGGGTGTGTAAGCGCTGTGAGGCGTTGAGCTAACCAATACTAATTGCCCGTGAGGCTTGACCATATAACACCCAAGCAATCTGCTACGCAGATTGCGGTGGTGAAGACGAAAGa                                                                                                                       >  2:330814/1‑149 (MQ=32)
                      gggCCGTCGAAGACTACGACGTTGATAGGTTGGGTGTGTAAGCGCTGTGAGGCGTTGAGCTAACCAATACTAATTGCCCGTGAGGCTTGACCATATAACACCCAAGCAATCTGCTACGCAGATTGCGGTGGTGAAGACGAAAGACCCg                                                                                                                   >  7:30367/1‑148 (MQ=32)
                          cGTCGAAGACTACGACGTTGATAGCTTGGGTGTGTAAGCGCTGTGAGGCGTTGAGCTAACCAATACTAATTGCCCGTGAGGCTTGACCATATAACACCCAAGCAATCTGCTACGCAGATTGCGGTGGTGAAGACGAAAGACCCGAAAAt                                                                                                              <  1:146676/149‑1 (MQ=17)
                               aaGACTACGACGTTGATAGGTTGGGTGTGTAAGCGCTGTGAGGCGTTGAGCTAACCAATACTAATTGCCCGTGAGGCTTGACCATATAACACCCAAGCAATCTGCTACGCAGATTGCGGTGGTGAAGACGAAAGACCCGAAAATTCGTg                                                                                                         <  2:34376/149‑1 (MQ=32)
                                          gTTGATAGGTTGGGTGTGTAAGCGCTGTGAGGCGTTGAGCTAACCAATACTAATTGCCCGTGAGGCTTGACCATATAACACCCAAGCAATCTGCTACGCAGATTGCGGTGGTGAAGACGAAAGAccc                                                                                                                    <  6:62778/127‑1 (MQ=32)
                                          gTTGATAGGTTGGGTGTGTAAGCGCTGTGAGGCGTTGAGCTAACCAATACTAATTGCCCGTGAGGCTTGACCATATAACACCCAAGCAATCTGCTACGCAGATTGCGGTGGTGAAGACGAAAGAccc                                                                                                                    >  5:62778/1‑127 (MQ=32)
                                           ttGATAGGTTGGGTGTGTAAGCGCTGTGAGGCGTTGAGCTAACCAATACTAATTGCCCGTGGGGCTTGACCATATAACACCCAAGCAATCTGCTACGCAGATTGCGGTGGTGAAGACGAAAGACCCGAAAATTCGTGCTTTACACAaga                                                                                             <  1:213231/149‑1 (MQ=17)
                                           ttGATAGGTTGGGTGTGTAAGCGCTGTGAGGCGTTGAGCTAACCAATACTAATTGCCCGTGAGGCTTGACCATATAACACCCAAGCAATCTGCTACGCAGATTGCGGTGGTGAAGACGAAAGACCCGAAAATTCGTGCTTTACACAaga                                                                                             <  7:82407/149‑1 (MQ=32)
                                             gATAGGTTGGGTGTGTAAGCGCTGTGAGGCGTTGAGCTAACCAATACTAATTGCCCGTGAGGCTTGACCATATAACACCCAAGCAATCTGCTACGCAGATTGCGGTGGTGAAGACGAAAGACCCGAAAATTCGTGCATTACACAAGAGc                                                                                           >  3:240224/1‑149 (MQ=17)
                                               tAGGTTGGGTGTGTAAGCGCTGTGAGGCGTTGAGCTAACCAATACTAATTGCCCGTGAGGCTTGACCATATAACACCCAAGCAATCTGCTACGCAGATTGCGGTGGTGAAGACGAAAGACCCGAAAATTCGTGCTTTACACAAGAGCtt                                                                                         >  7:101062/1‑149 (MQ=32)
                                               tAGGTTGGGTGTGTAAGCGCTGTGAGGCGTTGAGCTAACCAATACTAATTGCCCGTGAGGCTTGACCATATAACACCCAAGCAATCTGCTACGCAGATTGCGGTGGTGAAGACGAAAGACCCGAAAATTCGTGCTTTACACAAGAGCtt                                                                                         >  4:174890/1‑149 (MQ=32)
                                                      ggTGTGTAAGCGCTGTGAGGCGTTGAGCTAACCAATACTAATTGCCCGTGAGGCTTGACCATATAACACCCAAGCAATCTGCTACGCAGATTGCGGTGGTGAAGACGAAAGACCCGAAAATTCGTGCTTTACACAAGAGCTTCGCAAat                                                                                  <  3:269998/149‑1 (MQ=32)
                                                            tAAGCGCTGTGAGGCGTTGAGCTAACCAATACTAATTGCCCGTGAGGCTTGACCATATAACACCCAAGCAATCTGCTACGCAGATTGCGGTGGTGAAGACGAAAGACCCGAAAATTCGTGCttt                                                                                                     <  2:123590/124‑1 (MQ=32)
                                                            tAAGCGCTGTGAGGCGTTGAGCTAACCAATACTAATTGCCCGTGAGGCTTGACCATATAACACCCAAGCAATCTGCTACGCAGATTGCGGTGGTGAAGACGAAAGACCCGAAAATTCGTGCttt                                                                                                     >  1:123590/1‑124 (MQ=32)
                                                            tAAGCGCTGTGAGGCGTTGAGCTAACCAATACTAATTGCCCGTGAGGCTTGACCATATAACACCCAAGCAATCTGCTACGCAGATTGCGGTGGTGAAGACGAAAGACCCGAAAATTCGTGCTTTACACAAGAGCTTCGCAAATATCGc                                                                             <  3:74475/148‑1 (MQ=32)
                                                            tAAGCGCTGTGAGGCGTTGAGCTAACCAATACTAATTGCCCGTGAGGCTTGACCATATAACACCCAAGCAATCTGCTACGCAGATTGCGGTGGTGAAGACGAAAGACCCGAAAATTCGTGCTTTACACAAGAGCTTCGCAAATATCGCa                                                                            >  2:31041/1‑149 (MQ=32)
                                                            tAAGCGCTGTGAGGCGTTGAGCTAACCAATACTAATTGCCCGTGAGGCTTGACCATATAACACCCAAGCAATCTGCTACGCAGATTGCGGTGGTGAAGACGAAAGACCCGAAAATTCGTGCTTTACACAAGAGCTTCGCAAATATCGCa                                                                            >  8:245604/1‑149 (MQ=32)
                                                                           gTTGAGCTAACCAATACTAATTGCCCGTGAGGCTTGACCATATAACACCCAAGCAATCTGCTACGCAGATTGCGGTGGTGAAGACGAAAGACCCGAAAATTCGTGCTTTACACAAGAACTTCGCAAATATc                                                                               <  7:18986/131‑1 (MQ=17)
                                                                           gTTGAGCTAACCAATACTAATTGCCCGTGAGGCTTGACCATATAACACCCAAGCAATCTGCTACGCAGATTGCGGTGGTGAAGACGAAAGACCCGAAAATTCGTGCTTTACACAAGAACTTCGCAAATATc                                                                               >  8:18986/1‑131 (MQ=17)
                                                                              gAGCTAACCAATACTAATTGCCCGTGAGGCTTGACCATATAACACCCAAGCAATCTGCTACGCAGATTGCGGTGGTGAAGACGAAAGACCCGAAAATTCGTGCTTTACACAAGAGCTTCGCAAATATc                                                                               >  8:159463/1‑128 (MQ=32)
                                                                              gAGCTAACCAATACTAATTGCCCGTGAGGCTTGACCATATAACACCCAAGCAATCTGCTACGCAGATTGCGGTGGTGAAGACGAAAGACCCGAAAATTCGTGCTTTACACAAGAGCTTCGCAAATATCGCATATCCGAATTCGCTggg                                                           <  1:330814/148‑1 (MQ=32)
                                                                                 cTAACCAATACTAATTGCCCGTGAGGCTTGACCATATAACACCCAAGCAATCTGCTACGCAGATTGCGGTGGTGAAGACGAAAGACCCGAAAATTCGTGCTTTACACAAGAGCTTCGCAAATATCGCATATCCGAATTCGCTgg                                                            <  8:251246/144‑1 (MQ=32)
                                                                                 cTAACCAATACTAATTGCCCGTGAGGCTTGACCATATAACACCCAAGCAATCTGCTACGCAGATTGCGGTGGTGAAGACGAAAGACCCGAAAATTCGTGCTTTACACAAGAGCTTCGCAAATATCGCATATCCGAATTCGCTgg                                                            >  7:251246/1‑144 (MQ=32)
                                                                                    aCCAATACTAATTGCCCGTGAGGCTTGACCATATAACACCCAAGCAATCTGCTACGCAGATTGCGGTGGTGAAGACGAAAGACCCGAAAATTCGTGCTTTACACAAGAGCTTCGCAAATATCGCATATCCGAATTCGCTGGGCTGTCCa                                                    >  6:154245/1‑149 (MQ=32)
                                                                                     ccAATACTAATTGCCCGTGAGGCTTGACAATATAACACCCAAGCAATCTGCTACGCAGATTGCGGTGGTGAAGACGAAAGACCCGAAAATTCGTGCTTTACACAAGAGCTTCGCAAATATCGCATATCCGAATTCGCTGGGCTGTCCAt                                                   <  5:296827/149‑1 (MQ=17)
                                                                                       aaTACTAATTGCCCGTGAGGCTTGACCATATAACACCCAAGCAATCTGCTACGCAGATTGCGGTGGTGAAGACGAAAGACCCGAAAATTCGTGCTTTACACAAGAGCTTCGCAAATATCGCATATCCGAATTCGCTGGGCTGTCCATCt                                                 >  2:317458/1‑149 (MQ=32)
                                                                                            tAATTGCCCGTGAGGCTTGACCATATAACACCCAAGCAATCTGCTACGCAGATTGCGGTGGTGAAGACGAAAGACCCGAAAATTCGTGCTTTACACAAGAGCTTCGCAAATATCGCATATcc                                                                       >  1:330440/1‑122 (MQ=32)
                                                                                            tAATTGCCCGTGAGGCTTGACCATATAACACCCAAGCAATCTGCTACGCAGATTGCGGTGGTGAAGACGAAAGACCCGAAAATTCGTGCTTTACACAAGAGCTTCGCAAATATCGCATATcc                                                                       <  2:330440/122‑1 (MQ=32)
                                                                                            tAATTGCCCGTGAGGCTTGACCATATAACACCCAAGCAATCTGCTACGCAGATTGCGGTGGTGAAGACGAAAGACCCGAAAATTCGTGCTTTACACAAGAGCTTCGCAAATATCGCATATCCGAATTCGCTGGGCTGTCCATCt                                                 <  7:67293/144‑1 (MQ=32)
                                                                                            tAATTGCCCGTGAGGCTTGACCATATAACACCCAAGCAATCTGCTACGCAGATTGCGGTGGTGAAGACGAAAGACCCGAAAATTCGTGCTTTACACAAGAGCTTCGCAAATATCGCATATCCGAATTCGCTGGGCTGTCCATCt                                                 >  8:67293/1‑144 (MQ=32)
                                                                                             aaTTGCCCGTGAGGCTTGACCATATAACACCCAAGCAATCTGCTACGCAGATTGCGGTGGTGAAGACGAAAGACCCGAAAATTCGTGCTTTACACAAGGGCTTCGCAAATATCGCATATCCGAATTCGCt                                                              >  4:277630/1‑130 (MQ=17)
                                                                                             aaTTGCCCGTGAGGCTTGACCATATAACACCCAAGCAATCTGCTACGCAGATTGCGGTGGTGAAGACGAAAGACCCGAAAATTCGTGCTTTACACAAGGGCTTCGCAAATATCGCATATCCGAATTCGCt                                                              <  3:277630/130‑1 (MQ=17)
                                                                                                  cccGTGAGGCTTGACCATATAACACCCAAGCAATCTGCTACGCAGATTGCGGTGGTGAAGACGAAAGACCCGAAAATTCGTGCTTTACACAAGAGCTTCGCAAATATCGCATATCCGAATTCGCTGGGCTGTCCATCTGGACATTCTg                                       <  8:30367/148‑1 (MQ=32)
                                                                                                        aGGCTTGACCATATAACACCCAAGCAATCTGCTACGCAGATTGCGGTGGTGAAGACGAAAGACCCGAAAATTCGTGCTTTACACAAGAGCTTCGCAAATATCGCATATCCGAATTCGCTGGGCTGTCCATCTGGACATTCTGGCaa                                   >  6:187134/1‑146 (MQ=32)
                                                                                                           cTTGACCATATAACACCCAAGCAATCTGCTACGCAGATTGCGGTGGTGAAGACGAAAGACCCGAAAATTCGTGCTTTACACAAGAGCTTCGCAAATATCGCATATCCGAATTCGCTGGGCTGTCCATCTGGACATTCTGGCAACAGAAt                             <  1:317458/149‑1 (MQ=32)
                                                                                                            ttGACCATATAACACCCAAGCAATCTGCTACGCAGATTGCGGTGGTGAAGACGAAAGACCCGAAAATTCGTGCTTTACACAAGAGCTTCGCAAATATCGCATATCCGAATTCGCTGGGCTGTc                                                      >  1:194232/1‑123 (MQ=32)
                                                                                                            ttGACCATATAACACCCAAGCAATCTGCTACGCAGATTGCGGTGGTGAAGACGAAAGACCCGAAAATTCGTGCTTTACACAAGAGCTTCGCAAATATCGCATATCCGAATTCGCTGGGCTGTc                                                      <  2:194232/123‑1 (MQ=32)
                                                                                                              gACCATATAACACCCAAGCAATCTGCTACGCAGATTGCGGTGGTGAAGACGAAAGACCCGAAAATTCGTGCTTTACACAAGAGCTTCGCAAATATCGCATATCCGAATTCGCTGGGCTGTCCATCTGGACATTCTGGCAACAGAATTTc                          >  4:193480/1‑149 (MQ=32)
                                                                                                               aCCATATAACACCCAAGCAATCTGCTACGCAGATTGCGGTGGTGAAGACGAAAGACCCGAAAATTCGTGCTTTACACAAGAGCTTCGCAAATATCGCATATCCGAATTCGc                                                               <  6:299058/111‑1 (MQ=32)
                                                                                                               aCCATATAACACCCAAGCAATCTGCTACGCAGATTGCGGTGGTGAAGACGAAAGACCCGAAAATTCGTGCTTTACACAAGAGCTTCGCAAATATCGCATATCCGAATTCGc                                                               >  5:299058/1‑111 (MQ=32)
                                                                                                               aCCATATAACACCCAAGCAATCTGCTACGCAGATTGCGGTGGTGAAGACGAAAGACCCGAAAATTCGTGCTTTACACAAGAGCTTCGCAAATATCGCATATCCGAATTCGCTGGGCTGTCCATCTGGACATTCTGGCAACAGAATTTc                          <  3:193480/148‑1 (MQ=32)
                                                                                                                ccATATAACACCCAAGCAATCTGCTACGCAGATTGCGGTGGTGAAGACGAAAGACCCGAAAATTCGTGCTTTACACAAGAGCTTCGCAAATATCGCATATCCGAATTCGCTGGGCTGTCCATCTGGACATTCTGGCAACAGAATTTCtt                        <  7:245604/149‑1 (MQ=32)
                                                                                                                 cATATAACACCCAAGCAATCTGCTACGCAGATTGCGGTGGTGAAGACGAAAGACCCGAAAATTCGTGCTTTACACAAGAGCTTCGCAAATATCGCATATCCGAATTCGCTGGGCTGTCCATCTGGACATTCTGGCAACAGAATTTCtt                        >  2:92084/1‑148 (MQ=32)
                                                                                                                 cATATAACACCCAAGCAATCTGCTACGCAGATTGCGGTGGTGAAGACGAAAGACCCGAAAATTCGTGCTTTACACAAGAGCTTCGCAAATATCGCATATCCGAATTCGCTGGGCTGTCCATCTGGACATTCTGGCAACAGAATTTCtt                        >  4:299093/1‑148 (MQ=32)
                                                                                                                           ccAAGCAATCTGCTACGCAGATTGCGGTGGTGAAGACGAAAGACCCGAAAATTCGTGCTTTACACAAGAGCTTCGCAAATATCGCATATCCGAATTCGCTGGGCTGTCCATCTGGACATTCTGGCAACAGAATTTCTTGACGACCa                >  6:50409/1‑146 (MQ=32)
                                                                                                                           ccAAGCAATCTGCTACGCAGATTGCGGTGGTGAAGACGAAAGACCCGAAAATTCGTGCTTTACACAAGAGCTTCGCAAATATCGCATATCCGAATTCGCTGGGCTGTCCATCTGGACATTCTGGCAACAGAATTTCTTGACGACCa                >  5:178482/1‑146 (MQ=32)
                                                                                                                            cAAGCAATCTGCTACGCAGATTGCGGTGGTGAAGACGAAAGACCCGAAAATTCGTGCTTTACACAAGAGCTTCGCAAATATCGCATATCCGAATTCGCTGGGCTGTCCATCTGGACATTCTGGCAACAGAATTTCTTGACGACCATa              >  2:116461/1‑147 (MQ=32)
                                                                                                                             ccgCAATCTGCGTAGCAGATTGCGGTGGTGAAGACGAAAGACCCGAAAATTCGTGCTTTACACAAGAGCTTCGCAAATATCGCATATCCGAATTCGCTGGGCTGTCCATCTGGACATTCTGGCAACAGAATTTCTTGACGACCATagag           >  1:100238/3‑149 (MQ=12)
                                                                                                                               gCAATCTGCGTAGCAGATTGCGGTGGTGAAGACGAAAGACCCGAAAATTCGTGCTTTACACAAGAGCTTCGCAAATATCGCATATCCGAATTCGCTGGGCTGTCTATCTGGACATTCTGGCAACAGAATTTCTTGACGACCATAGAGCa         >  1:74738/1‑149 (MQ=255)
                                                                                                                                cAATCTGCTACGCAGATTGCGGTGGTGAAGACGAAAGACCCGAAAATTCGTGCTTTACACAAGAGCTTCGCAAATATCGCATATCCGAATTCGCTGGGCTGTCCATCTGGACATTCTGGCAACAGAATTTCTTGACGACCATAGAGCAt        <  5:154245/149‑1 (MQ=32)
                                                                                                                                  aTCTGCTACGCAGATTGCGGTGGTGAAGACGAAAGACCCGAAAATTCGTGCTTTACACAAGAGCTTCGCAAATATCGCATATCCGAATTCGCTGGGCTGTCCATCTGGACATTCTGGCAACAGAATTTCTTGACGACCATAGAGCATTg      <  1:92084/149‑1 (MQ=32)
                                                                                                                                   tCTGCTACGCTGATTGCGGTGCTGAAGACGAAAGACCCGAAAATTCGTGCTTTACACAAGAGCTTCGCAAATATCGCATATCCGAATTCGCTGGGCTGTCCATCTGGACATTCTGGCAACAGAATTTCTTGACGACCATAGAGCATTgg     >  5:87298/1‑149 (MQ=17)
                                                                                                                                     tGCTACGCAGATTGCGGTGGTGAAGACGAAAGACCCGAAAATTCGTGCTTTACACAAGAGCTTCGCAAATATCGCATATCCGAATTCGCTGGGCTGTCCATCTGGACATTCTGGCAACAGAATTTCTTGACGACCATAGAGCATTg      <  1:239771/146‑1 (MQ=32)
                                                                                                                                     tGCTACGCAGATTGCGGTGGTGAAGACGAAAGACCCGAAAATTCGTGCTTTACACAAGAGCTTCGCAAATATCGCATATCCGAATTCGCTGGGCTGTCCATCTGGACATTCTGGCAACAGAATTTCTTGACGACCATAGAGCATTg      >  2:239771/1‑146 (MQ=32)
                                                                                                                                     tGCTACGCAGATTGCGGTGGTGAAGACGAAAGACCCGAAAATTAGTGCTTTACACAAGAGCTTCGCAAATATCGAATATCCGAATTCGCTGGGCTGTCCATCTGGAAATTCTGTCAACAGAATTTCATGAACACCATAGAGCATTGGaa   >  4:194130/1‑149 (MQ=11)
                                                                                                                                       cTACGCAGATTGCGGTGGTGAAGACGAAAGACCCGAAAATTCGTGCTTTACACAAGAGCTTCGCAAATATCGCATATCCGAATTCGCTGGGCTGTCCATCTGGACATTCTGGCAACAGAATTTCTTGACGACCATAGAGCATTGGAac  >  1:248819/1‑148 (MQ=31)
                                                                                                                                         |                                                                                                                                                 
CTGGGATCTTGAATCCCCTAAAGGGCCGTCGAAGACTACGACGTTGATAGGTTGGGTGTGTAAGCGCTGTGAGGCGTTGAGCTAACCAATACTAATTGCCCGTGAGGCTTGACCATATAACACCCAAGCAATCTGCTACGCAGATTGCGGTGGTGAAGACGAAAGACCCGAAAATTCGTGCTTTACACAAGAGCTTCGCAAATATCGCATATCCGAATTCGCTGGGCTGTCCATCTGGACATTCTGGCAACAGAATTTCTTGACGACCATAGAGCATTGGAAC  >  NC_002947/702628‑702910

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: