Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,777,403 A→C 6.5% intergenic (+22/+21) 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,4030AC 6.5% 219.2 / 4.6 72intergenic (+22/+21)betA‑II/PP_5065choline dehydrogenase/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (34/33);  new base C (3/2);  total (37/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 = 9.37e-01

CCGAGAAGATTTCGGACAAGATCCGTGGTCGCAAGCCACTGCCGCGCAGCACCGCCAAGTACTACGTGGCGGGGGATGCACCAGTGAAAGGCAAGCCGATGCGTGAAGTGAAACAGGGCTGAGGCCTTGATCAGGGGCCGCGCAGCGGCCCCTTTCACGTTTGCTCAGGCGGGCAGCGCAAAAAATGTGCGCATCAACGAATCAT‑‑CCTCGACACCCAGTCGCCCCGCCAGCGCCTGGCGATCGAACAACTGTCCTGTGCGCATTTCACGCACGGCCTTCCCGCCTGGCGCCT  >  NC_002947/5777260‑5777551
                                                                                                                                               |                                                                                                                                                         
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)
                                                                                                                                               aGCGGCCCCTTTCACGTTTGCTCAGGCGGGCAGCGCAAAAAATGTGCGCATCAACGAATCAT‑‑CCTCGACACCCAGTCGCCCCGCCAGCGCCTGGCGATCGAACAACTGTCCTGTGCGCATTTCACGCACGGCCTTCCCGCCTGgcgcag  <  7:244429/149‑3 (MQ=255)
                                                                                                                                               |                                                                                                                                                         
CCGAGAAGATTTCGGACAAGATCCGTGGTCGCAAGCCACTGCCGCGCAGCACCGCCAAGTACTACGTGGCGGGGGATGCACCAGTGAAAGGCAAGCCGATGCGTGAAGTGAAACAGGGCTGAGGCCTTGATCAGGGGCCGCGCAGCGGCCCCTTTCACGTTTGCTCAGGCGGGCAGCGCAAAAAATGTGCGCATCAACGAATCAT‑‑CCTCGACACCCAGTCGCCCCGCCAGCGCCTGGCGATCGAACAACTGTCCTGTGCGCATTTCACGCACGGCCTTCCCGCCTGGCGCCT  >  NC_002947/5777260‑5777551

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: