Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,777,402 C→G 5.6% intergenic (+21/+22) betA‑II → / ← PP_5065 choline dehydrogenase/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,777,4020CG 5.6% 233.5 / 2.8 71intergenic (+21/+22)betA‑II/PP_5065choline dehydrogenase/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (34/33);  new base G (2/2);  total (36/35)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.10e-01

CCGAGAAGATTTCGGACAAGATCCGTGGTCGCAAGCCACTGCCGCGCAGCACCGCCAAGTACTACGTGGCGGGGGATGCACCAGTGAAAGGCAAGCCGATGCGTGAAGTGAAACAGGGCTGAGGCCTTGATCAGGGGCCGCGCAGCGGCCCCTTTCACGTTTGCTCAGGCGGGCAGCGCAAAAAATGTGCGCATCAACGAATCAT‑‑CCTCGACACCCAGTCGCCCCGCCAGCGCCTGGCGATCGAACAACTGTCCTGTGCGCATTTCACGCACGGCCTTCCCGCCTGGC  >  NC_002947/5777260‑5777547
                                                                                                                                              |                                                                                                                                                      
ccGAGAAGATTTCGGACAAGATCCGTGGTCGCAAGCCACTGCCGCGCAGCACCGCCAAGTACTACGTGGCGGGGGATGCACCAGTGAAAGGCAAGCCGATGCGTGAAGTGAAACAGGGCTGAGGCCTTGATCAGGGGCCGCGCAGCGGc                                                                                                                                                  >  2:80298/1‑149 (MQ=255)
 cGAGAAGATTTCGGACAAGATCCGTGGTCGCAAGCCACTGCCGCGCAGCACCGCCAAGTACTACGTGGCGGGGGATGCACCAGTGAAAGGCAAGCCGATGCGTGAAGTGAACCAGGGCTGAGGCCTTGCTCAGGGGCCGCGCCGCGGcc                                                                                                                                                 >  8:124425/1‑149 (MQ=255)
  gagaAGTTTTCGGACAAGATGCGTGGTCGCAAGCCACTGCCGCGCAGCACCGCCAAGTACTACGTGGCGGGGGATGCACCAGTGAAAGGCAAGCCGATGCGTGAAGTGAAACAGGGCTGAGGCCTTGATCAGGGGCCGCGCAGCGGccc                                                                                                                                                <  1:7852/149‑1 (MQ=255)
  gagaAGATTTCGGACAAGATCCGTGGTCGCAAGCCACTGCCGCGCAGCACCGCCAAGTACTACGTGGCGGGGGATGCACCAGTGAAAGGCAAGCCGATGCGTGAAGTGAAACAGGGCTGAGGCCTTGATCAGGGGCCGCTGCGCGGccc                                                                                                                                                <  4:175928/149‑1 (MQ=255)
  gagaAGATTTCGGACAAGATCCGTGGTCGCAAGCCACTGCCGCGCAGCACCGCCAAGTACTACGTGGCGGGGGATGCACCAGTGAAAGGCAAGCCGATGCGTGAAGTGAAACAGGGCTGAGGCCTTGATCAGGGGCCGCGCAGCGGccc                                                                                                                                                >  7:153278/1‑149 (MQ=255)
      agaTTTCGGACAAGATCCGTGGTCGCAAGCCACTGCCGCGCAGCACCGCCAAGTACTACGTGGCGGGGGATGCACCAGTGAAAGGCAAGCCGATGCGTGAAGTGAAACAGGGCTGAGGCCTTGATCAGGGGCCGCTGCGCGGCCCCtga                                                                                                                                            <  1:119788/149‑3 (MQ=255)
          ttCGGACAAGATCCGTGGTCGCAAGCCACTGCCGCGCAGCACCGCCAAGTACTACGTGGCGGGGGATGCACCAGTGAAAGGCAAGCCGATGCGTGAAGTGAAACAGGGCTGAGGCCTTGATCAGGGGCCGCGCAGCGGCCCCTTTCACg                                                                                                                                        >  1:199943/1‑149 (MQ=255)
             ggACAAGCTCCGTGGTCGCAAGCCACTGCCGCGCAGCACCGCCAAGTACTACGTGGCGGGGGATGCACCAGTGAAAGGCAAGCCGATGCGTGAAGTGAAACAGGGCTGAGGCCTTGATCAGGGGCCGCGCAGCGGCCCCTTTCACGtt                                                                                                                                      <  7:249083/148‑1 (MQ=255)
              gACAAGATCCTTGTTGTCAAGCCTCTGCCGCGCATCACCGCCAAGTACTACGTTGCGGGGGATTCATCAGTGAAAGGCAAGCCGATGCGTTAAGTGAAACAGGGCTGAGGCCTTGATCAGGGGCCGCGCAGCGGCCCCTTTCACGttt                                                                                                                                     <  2:343935/148‑1 (MQ=255)
               aCAAGATCCGTGGTCGCAAGCCACTGCCGCGCAGCACCGCCAAGTACTACGTGGCGGGGGATGCACCAGTGAAAGGCAAGCCGATGCGTGAAGTGAAACAGGGCTGAGGCCTTGATCAGGGGCCGCGCAGCGGCCCCTTTCACGTTTGc                                                                                                                                   <  2:126021/149‑1 (MQ=255)
                cAAGATCCGTGGTCGCAAGCCACTGCCGCGCAGCACCGCCAAGTACTACGTGGCGGGGGATGCACCAGTGAAAGGCAAGCCGATGCGTGAAGTGAAACAGGGCTGAGGCCTTGATCAGGGGCCGCGCAGCGGCCCCTTTCACGTTTGCt                                                                                                                                  <  3:203597/149‑1 (MQ=255)
                 aaGATCCGTGGTCGCAAGCCACTGCCGCGCAGCACCGCCAAGTACTACGTGGCGGGGGATGCACCAGTGAAAGGCAAGCCGATGCGTGAAGTGAAACAGGGCTGAGGCCTTGATCAGGGGCCGCGCAGCGGCCCCTTTCACGTTTGCt                                                                                                                                  >  8:101139/1‑148 (MQ=255)
                 aaGATCCGTGGTCGCAAGCCACTGCCGCGCAGCACCGCCAAGTACTACGTGGCGGGGGATGCACCAGTGAAAGGCAAGCCGATGCGTGAAGTGAAACAGGGCTGAGGCCTTGATCAGGGGCCGCGCAGCGGCCCCTTTCACGTTTGCt                                                                                                                                  >  2:107776/1‑148 (MQ=255)
                  aGATCCGTGGTCGCAAGCCACTGCCGCGCAGCACCGCCAAGTACTACGTGGCGGGGGATGCACCAGTGAAAGGCAAGCCGATGCGTGAAGTGAAACAGGGCTGAGGCCTTGATCAGGGGCCGCGCAGCGGCCCCTTTCACGTTTGCTCa                                                                                                                                <  1:222493/149‑1 (MQ=255)
                       cGTGGTCGCAAGCCACTGCCGCGCAGCACCGCCAAGTACTACGTGGCGGGGGATGCACCAGTGAAAGGCAAGCCGATGCGTGAAGTGAAACAGGGCTGAGGCCTTGATCAGGGGCCGCGCAGCGGCCCCTTTCACGTTTGCTCAGGCgg                                                                                                                          <  1:107776/149‑1 (MQ=255)
                        gTGGTCGCAAGCCACTGCCGCGCAGCACCGCCAAGTACTACGTGGCGGGGGATGCACCAGTGAAAGGCAAGCCGATGCGTGAAGTGAAACAGGGCTGAGGCCTTGATCAGGGGCCGCGCAGCGGCCCCTTTCACGTTTGCTCAGGCggg                                                                                                                        >  4:119813/1‑149 (MQ=255)
                          ggTCGCAAGCCACTGCCGCGCCGCCCCGCCAAGTACTACGGGGCGGGGGATGCACCCGGGAAAGGCAAGCCGATGCGTGAAGTGCAACAGGGCTGAGGCCTTGATCCGGGGCCGCGCAGCGGCCCCTTTCACGTTTGCTCAGGCGGGCa                                                                                                                      <  7:124425/149‑1 (MQ=255)
                                  gCCACTGCCGCGCAGCACCGCCAAGTACTACGTGGCGGGGGATGCACCAGTGAAAGGCAAGCCGATGCGTGAAGTGAAACAGGGCTGAGGCCTTGATCAGGGGCCGCGCAGCGGCCCCTTTCACGTTTGCTCAGGCGGGCAGCGCaaa                                                                                                               >  5:352006/1‑148 (MQ=255)
                                     aCTGCCGCGCAGCACCGCCAAGTACTACGTGGCGGGGGATGCACCAGTGAAAGGCAAGCCGATGCGTGAAGTGAAACAGGGCTGAGGCCTTGATCAGGGGCCGCGCAGCGGCCCCTTTCACGTTTGCTCAGGCGGGCAGCGCAAAAAAt                                                                                                           <  3:77947/149‑1 (MQ=255)
                                      cTGCCGCGCAGCACCGCCAAGTACTACGTGGCGGGGGATGCACCAGTGAAAGGCAAGCCGATGCGTGAAGTGAAACAGGGCTGAGGCCTTGATCAGGGGCCGCGCAGCGGCCCCTTTCACGTTTGCTCAGGCGGGCAGCGCAAAAAAtg                                                                                                          <  3:148364/149‑1 (MQ=255)
                                         ccGCGCAGCACCGCCAAGTACTACGTGGCGGGGGATGCACCAGTGAAAGGCAAGCCGATGCGTGAAGTGAAACAGGGCTGAGTCCTTGATCAGGGGCCGCGCAGCGGCCCCTTTCACGTTTGCTCAGGCGGGCAGCGCAAAAAATGTgc                                                                                                       <  5:270803/149‑1 (MQ=255)
                                           gcgcAGCACCGCCAAGTACTACGTGGCGGGGGATGCACCAGTGAAAGGCAAGCCGATGCGTGAAGTGAAACAGGGCTGAGGCCTTGATCAGGGGCCGCGCAGCGGCCCCTTTCACGTTTGCTCAGGCGGGCAGCGCAAAAAATGTgcgc                                                                                                     <  2:310083/149‑1 (MQ=255)
                                             gcagcaCCGCCAAGTACTACGTGGCGGGGGATGCACCAGTGAAAGGCAAGCCGATGCGTGAAGTGAAACAGGGCTGAGGCCTTGATCAGGGGCCGCGCAGCGGCCCCTTTCACGTTTGCTCAGGCGGGCAGCGCAAAAAATGTGCGCAt                                                                                                   <  3:355817/149‑1 (MQ=255)
                                               agcaCCGCCAAGTACTACGTGGCGGGGGATGCACCAGTGAAAGGCAAGCCGATGCGTGAAGTGAAACAGGGCTGAGGCCTTGATCAGGGGCCGCGCAGCGGCCCCTTTCACGTTTGCTCAGGCGGGCAGCGCAAAAAATGTGCGCATc                                                                                                  >  1:89623/1‑148 (MQ=255)
                                                 caCCGCCAAGTACTACGTGGCGGGGGATGCACCAGTGAAAGGCAAGCCGATGCGTGAAGTGAAACAGGGCTGAGGCCTTGATCAGGGGCCGCGCAGCGGCCCCTTTCACGTTTGCTCAGGCGGGCAGCGCAAAAAATGTGCGCATCAAc                                                                                               >  3:92627/1‑149 (MQ=255)
                                                   ccGCCACCTAGTACTTGGCGGGGGATGCACCAGTGAAAGGCAAGCCGATGCGTGAAGTGAAACAGGGCTGAGGCCTTGATCAGAGGCCGCGCAGCGGCCCCTTTCACGTTTGCTCAGGCGGGCAGCGCAAAAAATGTGCGCATCAACg                                                                                              >  7:251782/1‑148 (MQ=255)
                                                    cGCCAAGTACTACGTGGCGGGGGATGCACCAGTGAAAGGCAAGCCGATGCGTGAAGTGAAACAGGGCTGAGGCCTTGATCAGGGGCCGCGCAGCGGCCCCTTTCACGATTGCTCAGGCGGGCAGCGCAAAAAATGTGCGCATCAACGa                                                                                             >  7:156268/1‑148 (MQ=255)
                                                     gCCAAGTACTACGTGGCGGGGGATGCACCAGTGAAAGGCAAGCCGATGCGTGAAGTGAAACAGGGCTGAGGCCTTGATCAGGGGCCGCGCAGCGGCCCCTTTCACGTTTGCTCAGGCGGGCAGCGCAAAAAATGTGCGCATCAACGAat                                                                                           <  2:181332/149‑1 (MQ=255)
                                                     gCCAAGTACTACGTGGCGGGGGATGCACCAGTGAAAGGCAAGCCGATGCGTGAAGTGAAACAGGGCTGAGGCCTTGATCAGGGGCCGCGCAGCGGCCCCTTTCACGTTTGCTCAGGCGGGCAGCGCAAAAAATGTGCGCATCAACGAat                                                                                           <  3:119813/149‑1 (MQ=255)
                                                       cAAGTACTACGTGGCGGGGGATGCACCAGTGAAAGGCAAGCCGATGCGTGAAGTGAAACAGGGCTGAGGCCTTGATCAGGGGCCGCGCAGCGGCCCCTTTCACGTTTGCTCAGGCGGGCAGCGCAAAAAATGTGCGCATCAACGAatca                                                                                         <  7:316565/149‑1 (MQ=255)
                                                         aGTACTACGTGGCGGGGGATGCACCAGTGAAAGGCAAGCCGATGCGTGAAGTGAAACAGGGCTGAGGCCTTGATCAGGGGCCGCGCAGCGGCCCCTTTCACGTTTGCTCAGGCGGGCAGCGCAAAAAATGTGCGCATCAACGAatcat‑‑c                                                                                     <  2:199943/149‑1 (MQ=255)
                                                         aGTACTACGTGGCGGGGGATGCACCAGTGAAAGGCAAGCCGATGCGTGAAGTGAAACAGGGCTGAGGCCTTGATCAGGGGCCGCGCAGCGGCCCCTTTCACGTTTGCTCAGGCGGGCAGCGCAAAAAATGTGCGCATCAACGAatcat‑‑c                                                                                     >  2:112029/1‑149 (MQ=255)
                                                           tactacGTGGCGGGGGATGCACCAGTGAAAGGCAAGCCGATGCGTGAAGTGAAACAGGGCTGAGGCCTTGATCAGGGGCCGCGCAGCGGCCCCTTTCACGTTTGCTCAGGCGGGCAGCGCAAAAAATGTGCGCATCAACGAATCAT‑‑CCt                                                                                   >  1:66482/1‑149 (MQ=255)
                                                             ctacGTGGCGGGGGATGCACCAGTGAAAGGCAAGCCGATGCGTGAAGTGAAACAGGGCTGAGGCCTTGATCAGGGGCCGCGCAGCGGCCCCTTTCACGTTTGCTCAGGCGGGCAGCGCAAAAAATGTGCGCATCAACGAATCAT‑‑CCTc                                                                                  <  5:60743/148‑1 (MQ=255)
                                                                cGTGGCGGGGGATGCACCAGTGAAAGGCAAGCCGATGCGTGAAGTGAAACAGGGCTGAGGCCTTGATCAGGGGCCGCGCAGCGGCCCCTTTCACGTTTGCTCAGGCGGGCAGCGCAAAAAATGTGCGCATCAACGAATCAT‑‑CCTCGaca                                                                             >  6:271268/1‑149 (MQ=255)
                                                                cGTGGCGGGGGATGCACCAGTGAAAGGCAAGCCGATGCGTGAAGTGAAACAGGGCTGAGGCCTTGATCAGGGGCCGCGCAGCGGCCCCTTTCACGTTTGCTCAGGCGGGCAGCGCAAAAAATGTGCGCATCAACGAATCAT‑‑CCTCGac                                                                              <  8:261453/148‑1 (MQ=255)
                                                                             gCACCAGTGAAAGGCAAGCCGATGCGTGAAGTGAAACAGGGCTGAGGCCTTGATCAGGGGCCGCGCAGCGGCCCCTTACGGGGGGGTTCAGGGGAGAGCAGCGCAAAAAATGTGAGCATAAACGAATAAT‑‑CCTCGACACCCAGTCGccc                                                                  >  1:126650/1‑149 (MQ=255)
                                                                              cACCAGTGAAAGGCAAGCCGATGCGTGAAGTGAAACAGGGCTGAGGCCTTGATCAGGGGCCGCGCAGCGGCCCCTTTCACGTTTGCTCAGGCGGGCAGCGCAAAAAATGTGCGCATCAACGAATCAT‑‑CCTCGACACCCAGTCGCCCCGc                                                               >  6:183652/1‑149 (MQ=255)
                                                                                 cAGTGAAAGGCAAGCCGATGCGTGAAGTGAAACAGGGCTGAGGCCTTGATCAGGGGCCGCGCAGCGGCCCCTTTCACGTTTGCTCAGGCGGGCAGCGCAAAAAATGTGCGCATCTACGAATCAT‑‑CCTCGACACCCAGTCGCCCCGCCAg                                                            >  8:132295/1‑149 (MQ=255)
                                                                                    tGAAAGGCAAGCCGATGCGTGAAGTGAAACAGGGCTGAGGCCTTGATCAGGGGCCGCGCAGCGGCCCCTTTCACGTTTGCTCAGGCGGGCAGCGCAAAAAATGTGCGCATCAACGAATCAT‑‑CCTCGACACCCAGTCGCCCCGCCAgcgc                                                         >  2:77177/1‑149 (MQ=255)
                                                                                    tGAAAGGCAAGCCGATGCGTGAAGTGAAACAGGGCTGAGGCCTTGATCAGGGGCCGCGCAGCGGCCCCTTTCACGTTTGCTCAGGCGGGCAGCGCAAAAAAAAAGCGCATAAACGAATAAT‑‑CCTCGACACCCAGTCGCCCCGCCAgcgc                                                         >  1:130467/1‑149 (MQ=255)
                                                                                      aaaGGCAAGCCGATGCGTGAAGTGAAACAGGGCTGAGGCCTTGATCAGGGGCCGCGCAGCGGCCCCTTTCACGTTTGCTCAGGCGGGCAGCGCAAAAAATGTGCGCATCAACGAATCAT‑‑CCTCGACACCCAGTCGCCCCGCCAGCGCCt                                                       >  1:180231/1‑149 (MQ=255)
                                                                                      aaaGGCAAGCCGATGCGTGAAGTGAAACAGGGCTGAGGCCTTGATCAGGGGCCGCGCAGCGGCCCCTTTCACGTTTGCTCAGGCGGGCAGCGCAAAAAATGTGCGCATCAACGAATCAT‑‑CCTCGACACCCAGTCGCCCCGCCAGCGCCt                                                       >  2:164098/1‑149 (MQ=255)
                                                                                           cAAGCCGATGCGTGAAGTGAAACAGGGCTGAGGCCTTGATCAGGGGCCGCGCAGCGGCCCCTTTCACGTTTGCTCAGGCGGGCAGCGCAAAAAATGTGCGCATCAACGAATCAT‑‑CCTCGACACCCAGTCGCCCCGCCAGCGCCTGGCGa                                                  <  2:89623/149‑1 (MQ=255)
                                                                                            aaGCCGATGCGTGAAGTGAAACAGGGCTGAGGCCTTGATCAGGGGCCGCGCAGCGGCCCCTTTCACGTTTGCTCAGGCGGGCAGCGCAAAAAATGTGCGCATCAACGAATCAT‑‑CCTCGACACCCAGTCGCCCCGCCAGGCGCTGGCggg                                                 <  5:183652/149‑3 (MQ=255)
                                                                                            aaGCCGATGCGTGAAGTGAAACAGGGCTGAGGCCTTGATCAGGGGCCGCGCAGCGGCCCCTTTCACGTTTGCTCAGGCGGGCAGCGCAAAAAATGTGCGCATCAACGAATCAT‑‑CCTCGACACCCAGTCGCCCCGCCAGCGCCTGGCGAt                                                 <  2:66482/149‑1 (MQ=255)
                                                                                             aGCCGATGCGTGAAGTGAAACAGGGCTGAGGCCTTGATCAGGGGCCGCGCAGCGGCCCCTTTCACGTTTGCTCAGGCGGGCAGCGCAAAAAATGTGCGCATCAACGAATCAT‑‑CCTCGACACCCAGTCGCCCCGCCAGCGCCTGGCGATc                                                <  8:89488/149‑1 (MQ=255)
                                                                                                     cGTGAAGTGAAACAGGGCTGAGGCCTTGATCAGGGGCCGCGCAGCGGCCCCTTTCACGTTTGCTCAGGCGGGCAGCGCAAAAAATGTGCGCATCAACGAATCAT‑‑CCTCGACACCCAGTCGCCCCGCCAGCGCCTGGCGATCGaacaac                                         <  1:77177/148‑1 (MQ=255)
                                                                                                     cGTGAAGTGAAACAGGGCTGAGGCCTTGATCAGGGGCCGCGCAGCGGCCCCTTTCACGTTTGCTCAGGCGGGCAGCGCAAAAAATGTGCGCATCAACGAATCAT‑‑CCTCGACACCCAGTCGCCCCGCCAGCGCCTGGCGATCGAACAACt                                        <  8:153278/149‑1 (MQ=255)
                                                                                                       tgaagtgaaACAGGGCTGAGGCCTTGATCAGGGGCCGCGCAGCGGCCCCTTTCACGTTTGCTCAGGCGGGCAGCGCAAAAAATGTGCGCATCAACGAATCAT‑‑CCTCGACACCCAGTCGCCCCGCCAGCGCCTGGCGATCGAACAACTGt                                      <  6:352006/149‑1 (MQ=255)
                                                                                                         aagtgaaACAGGGCTGAGGCCTTGATCAGGGGCCGCGCAGCGGCCCCTTTCACGTTTGCTCAGGCGGGCAGCGCAAAAAATGTGCGCATCAACGAATCAT‑‑CCTCGACACCCAGTCGCCCCGCCAGCGCCTGGCGATCGAACAACTGTc                                     <  5:155414/148‑1 (MQ=255)
                                                                                                           gtgaaACAGGGCTGAGGCCTTGATCAGGGGCCGCGCAGCGGCCCCTTTCACGTTTGCTCAGGCGGGCAGCGCAAAAAATGTGCGCATCAACGAATCAT‑‑CCTCGACCACCCAGTCGCCCCGCCAGCGCCTGGCGATCGAACAACTGTCCt                                   <  1:80298/149‑1 (MQ=255)
                                                                                                             gaaACAGGGCTGAGGCCTTGATCAGGGGCCGCACAGCGGCCCCTTTCACGTTTGCTCAGGCGGGCAGCGCAAAAAATGTGCGCATCAACGAATCAT‑‑CCTCGACACCCAGTCGCCCCGCCAGCGCCTGGCGATCGAACAACTGTCCtgtg                                >  5:24454/1‑149 (MQ=255)
                                                                                                              aaaCAGGGCTGAGGCCTTGATCAGGGGCCGCGCAGCGGCCCCTTTCACGTTTGCTCAGGCGGGCAGCGCAAAAAATGTGCGCATCAACGAATCAT‑‑CCTCGACACCCAGTCGCCCCGCCAGCGCCTGGCGATCGAACAACTGTCCTGTgc                               >  7:39848/1‑149 (MQ=255)
                                                                                                               aaCAGGGCTGAGGCCTTGATCAGGGGCCGCGCAGCGGCCCCTTTCACGTTTGCTCAGGCGGGCAGCGCAAAAAATGTGCGCATCAACGAATCAT‑‑CCTCGACACCCAGTCGCCCCGCCAGCGCCTGGCGATCGAACAACTGTCCTGTgcg                              <  6:167301/149‑1 (MQ=255)
                                                                                                                aCAGGGCTGAGGCCTTGATCAGGGGCCGCGCAGCGGCCCCTTTCACGTTTGCTCAGGCGGGCAGCGCAAAAAATGTGCGCATCAACGAATCAT‑‑CCTCGACACCCAGTCGCCCCGCCAGCGCCTGGCGATCGAACAACTGTCCTGTgcg                              >  3:2627/1‑148 (MQ=255)
                                                                                                                    ggCTGAGGCCTTGATCAGGGGCCGCGCAGCGGCCCCTTTCACGTTTGCTCAGGCGGGCAGCGCAAAAAATGTGCGCATCAACGAATCAT‑‑CCTCGACACCCAGTCGCCCCGCCAGCGCCTGGCGATCGAACAACTGTCCTGTGCGCAttt                         <  4:92627/149‑1 (MQ=255)
                                                                                                                         aGGCCTTGATCAGGGGCCGCGCAGCGGCCCCTTTCACGTTTGCTCAGGCGGGCAGCGCAAAAAATGTGCGCATCAACGAATCATCCCCTCGACACCCAGTCGCCCCGCCAGCGCCTGGCGATCGAACAACTGTCCTGTGCGCATTTca                       <  1:112029/148‑1 (MQ=255)
                                                                                                                             cTTGATCAGGGGCCGCGCAGCGGCCCCTTTCACGTTTGCTCAGGCGGGCAGCGCAAAAAATGTGCGCATCAACGAATCAT‑‑CCTCGACACCCAGTCGCCCCGCCAGCGCCTGGCGATCGAACAACTGTCCTGTGCGCATTTCACGCACgg                >  5:21236/1‑149 (MQ=255)
                                                                                                                             cTTGATCAGGGGCCGCGCAGCGGCCCCTTTCACGTTTGCTCAGGCGGGCAGCGCAAAAAATGTGCGCATCAACGAATCAT‑‑CCTCGACACCCAGTCGCCCCGCCAGCGCCTGGCGATCGAACAACTGTCCTGTGCGCATTTCACGCACgg                >  1:118857/1‑149 (MQ=255)
                                                                                                                              ttGATCAGGGGCCGCGCAGCGGCCCCTTTCACGTTTGCTCAGGCGGGCAGCGCAAAAAATGTGCGCATCAACGAATCAT‑‑CCTCGACACCCAGTCGCCCCGCCAGCGCCTGGCGATCGAACAACTGTCCTGTGCGCATTTCACGCACGGc               >  6:354970/1‑149 (MQ=255)
                                                                                                                                gATCAGGGGCCGCGCAGCGTCCCCTTTCACGTTTGCTCAGGCGGGCAGCGCAAAAAATGTGCGCATCAACGAATCAT‑‑CCTCGACTCCCAGTCGCCCCGCCAGCGCCTGGCGATCGAACAACTGTCCTGTGCGCATTTCACGCACGGcc              <  8:156268/148‑1 (MQ=255)
                                                                                                                                   aaGGGGCCGCTGCGCGGCCCCTTTCACGTTTGCTCAGGCGGGCAGCGCAAAAAATGTGCGCATCCACGACTCAT‑‑CCTCGACACCCAGTCGCCCCGCCAGCGCCTGGCGATCGAACAACTGTCCTGTGCGCATTTCACGCACGGCCTTcc          >  7:362457/2‑149 (MQ=255)
                                                                                                                                     ggggCCGCTGCGCGGCCCCTTTCACGTTTGCTCAGGCGGGCAGCGCAAAAAATGTGCGCATCAACGAATCAT‑‑CCTCGACACCCAGTCGCCCCGCCAGCGCCTGGCGATCGAACAACTGTCCTGTGCGCATTTCACGCACGGCCTTCCCg        >  1:346604/1‑149 (MQ=255)
                                                                                                                                      gggCCGCGCAGCGGCCCCTTTCACGTTTGCTCAGGCGGGCAGCGCAAAAAATGTGCGCATCAACGAATCAT‑‑CCTCGACACCCAGTCGCCCCGCCAGCGCCTGGCGATCGAACAACTGTCCTGTGCGCATTTCACGCACGGCCTTCCCGc       >  3:130733/1‑149 (MQ=255)
                                                                                                                                      gggCCGCGCAGCGGCCCCTTTCACGTTTGCTCAGGCGGGCAGCGCAAAAAATGTGCGCATCAACGAATCAT‑‑CCTCGACACCCAGTCGCCCCGCCAGCGCCTGGCGATCGAACAACTGTCCTGTGCGCATTTCACGCACGGCCTTCCCGc       >  2:33248/1‑149 (MQ=255)
                                                                                                                                      gggCCGCGCAGCGGCCCCTTTCACGTTTGCTCAGGCGGGCAGCGCAAAAAATGTGCGCATCAACGAATCAT‑‑CCTCGACACCCAGTCGCCCCGCCAGCGCCTGGCGATCGAACAACTGTCCTGTGCGCATTTCACGCACGGCCTTCCCGc       >  7:37975/1‑149 (MQ=255)
                                                                                                                                      gggCCGCGCAGCGGCCCCTTTCACGTTTGCTCAGGCGGGCAGCGCAAAAAATGTGCGCATCAACGAATCAT‑‑CCTCGACACCCAGTAGCCCCGCCAGCGCCGGGCGATCGAACAACAGACCTGTGCGCATTTCCCGCACGGCCTTCCCGc       >  8:244429/1‑149 (MQ=255)
                                                                                                                                          cgcgCAGCGGCCCCTTTCACGTTTGCTCAGGCGGGCAGCGCAAAAAATGTGCGCATCAACGAATCAT‑‑CCTCGACACCCAGTCGCCCCGCCAGCGCCTGGCGATCGAACAACTGTCCTGTGCGCATTTCACGCACGGCCTTCCCGCCTgg   >  7:348809/1‑149 (MQ=255)
                                                                                                                                          cgcgCAGCGGCCCCTTTCACGTTTGCTCAGGCGGGCAGCACAAAAAATGTGCGCATCAATGAATCAT‑‑CCTCGACACCCAGTCGCCCCGCCAGCGCCTGGCGATCGAACAACTGTCCTGTGCGCAttt                         >  6:379284/1‑127 (MQ=255)
                                                                                                                                          cgcgCAGCGGCCCCTTTCACGTTTGCTCAGGCGGGCAGCACAAAAAATGTGCGCATCAATGAATCAT‑‑CCTCGACACCCAGTCGCCCCGCCAGCGCCTGGCGATCGAACAACTGTCCTGTGCGCAttt                         <  5:379284/127‑1 (MQ=255)
                                                                                                                                           gcgcAGCGGCCCCTTTCACGTTTGCTCAGGCGGGCAGCGCAAAAAATGTGCGCATCAACGAATCAT‑‑CCTCGCCACCCAGTCGCCCCGCCAGCGCCTGGCGATCGAACAACTGTCCTGTGCGCATTTCACGCACGGCCTTCCCGCCTGgc  >  5:324031/1‑149 (MQ=255)
                                                                                                                                              |                                                                                                                                                      
CCGAGAAGATTTCGGACAAGATCCGTGGTCGCAAGCCACTGCCGCGCAGCACCGCCAAGTACTACGTGGCGGGGGATGCACCAGTGAAAGGCAAGCCGATGCGTGAAGTGAAACAGGGCTGAGGCCTTGATCAGGGGCCGCGCAGCGGCCCCTTTCACGTTTGCTCAGGCGGGCAGCGCAAAAAATGTGCGCATCAACGAATCAT‑‑CCTCGACACCCAGTCGCCCCGCCAGCGCCTGGCGATCGAACAACTGTCCTGTGCGCATTTCACGCACGGCCTTCCCGCCTGGC  >  NC_002947/5777260‑5777547

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