Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,777,401 G→T 5.9% intergenic (+20/+23) 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,4010GT 5.9% 213.7 / 3.4 68intergenic (+20/+23)betA‑II/PP_5065choline dehydrogenase/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (31/32);  new base T (2/2);  total (34/34)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.69e-01

GATCGCCGAGAAGATTTCGGACAAGATCCGTGGTCGCAAGCCACTGCCGCGCAGCACCGCCAAGTACTACGTGGCGGGGGATGCACCAGTGAAAGGCAAGCCGATGCGTGAAGTGAAACAGGGCTGAGGCCTTGATCAGGGGCCGCGCAGCGGCCCCTTTCACGTTTGCTCAGGCGGGCAGCGCAAAAAATGTGCGCATCAACGAATCAT‑‑CCTCGACACCCAGTCGCCCCGCCAGCGCCTGGCGATCGAACAACTGTCCTGTGCGCATTTCACGCACGGCCTTCCCGCCTGGC  >  NC_002947/5777255‑5777547
                                                                                                                                                  |                                                                                                                                                       
gATCGCCGAGAAGATTTCGGACAAGATCCGTGGTCGCAAGCCACTGCCGCGCAGCACCGCCAAGTACTACGTGGCGGGGGATGCACCAGTGAAAGGCAAGCCGATGCGTGAAGTGAAACAGGGCTGAGGCCTTGATCAGGGGCcgcg                                                                                                                                                         >  6:226259/1‑147 (MQ=255)
     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                                                       >  2:164098/1‑149 (MQ=255)
                                                                                           aaaGGCAAGCCGATGCGTGAAGTGAAACAGGGCTGAGGCCTTGATCAGGGGCCGCGCAGCGGCCCCTTTCACGTTTGCTCAGGCGGGCAGCGCAAAAAATGTGCGCATCAACGAATCAT‑‑CCTCGACACCCAGTCGCCCCGCCAGCGCCt                                                       >  1:180231/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       >  2:33248/1‑149 (MQ=255)
                                                                                                                                           gggCCGCGCAGCGGCCCCTTTCACGTTTGCTCAGGCGGGCAGCGCAAAAAATGTGCGCATCAACGAATCAT‑‑CCTCGACACCCAGTCGCCCCGCCAGCGCCTGGCGATCGAACAACTGTCCTGTGCGCATTTCACGCACGGCCTTCCCGc       >  7:37975/1‑149 (MQ=255)
                                                                                                                                           gggCCGCGCAGCGGCCCCTTTCACGTTTGCTCAGGCGGGCAGCGCAAAAAATGTGCGCATCAACGAATCAT‑‑CCTCGACACCCAGTCGCCCCGCCAGCGCCTGGCGATCGAACAACTGTCCTGTGCGCATTTCACGCACGGCCTTCCCGc       >  3:130733/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)
                                                                                                                                                  |                                                                                                                                                       
GATCGCCGAGAAGATTTCGGACAAGATCCGTGGTCGCAAGCCACTGCCGCGCAGCACCGCCAAGTACTACGTGGCGGGGGATGCACCAGTGAAAGGCAAGCCGATGCGTGAAGTGAAACAGGGCTGAGGCCTTGATCAGGGGCCGCGCAGCGGCCCCTTTCACGTTTGCTCAGGCGGGCAGCGCAAAAAATGTGCGCATCAACGAATCAT‑‑CCTCGACACCCAGTCGCCCCGCCAGCGCCTGGCGATCGAACAACTGTCCTGTGCGCATTTCACGCACGGCCTTCCCGCCTGGC  >  NC_002947/5777255‑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: