Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,052,096 G→T 6.2% L156M (CTG→ATG)  PP_4452 ← NAD/NADP octopine/nopaline dehydrogenase family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,052,0960GT 6.2% 300.1 / 4.5 94L156M (CTG→ATG) PP_4452NAD/NADP octopine/nopaline dehydrogenase family protein
Reads supporting (aligned to +/- strand):  ref base G (45/43);  new base T (3/3);  total (48/46)
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.77e-01

AATCGGAACGTGGAAAACGGCATTGGTATTATTCAATGTCGACATCAGTGGGTTTTCATCCACTACGAAATCACTACCGAATAAATCACTGCATATACGCTGCGCATCTTCAATTTTTGAAGCTGGCAAGCAGGAGACGCCCAGCGATGGGCGACGCACCATGATATTGACCTTGTCTTCCGAAACCCTTCGTGCGGTAAGTGCAGTTGTGCCAAAACTGGCGACGACAAGATCTACCCCCAGGTTCCTCGCGCGCTCCCGAAGATAAAGCGCCGACAACGAGCCCATTGAAC  >  NC_002947/5051953‑5052245
                                                                                                                                               |                                                                                                                                                     
aaTCGGAACGTGGAAAACGGCATTGGTATTATTCAATGTCGACATCAGTGGGTTTTCATCCACTACGAAATCACTACCGAATAAATCACTGCATATACGCTGCGCATCTTCAATTTTTGAAGCTGGCAAGCAGGAGACGCCCATCGCTgg                                                                                                                                                 <  1:1798816/150‑1 (MQ=255)
aaTCGGAACGTGGAAAACGGCATTGGTATTAATCAATGTAGACATCAGTGGGTTTTCATCCACTACGAAATCACTACCGAATAAATCACTGCATATACGCTGCGCATCATCAATCTTTCAAGCCAGCAAGCAGGCGACGCCCAGcaagag                                                                                                                                                 >  2:1156610/1‑145 (MQ=255)
 aTCGGAACGTGGAAAACGGCATTGGTATTATTCAATGTCGACATCAGTGGGTTTTCATCCACTACGAAATCACTACCGAATAAATCACTGCATATACGCTGCGCATCTTCAATTTTTGAAGCTGGCAAGCAGGAGACGCCCAGCGATggg                                                                                                                                                >  2:1208359/1‑150 (MQ=255)
  tCGGAACGTGGAAAACGGCATTGGTATTATTCAATGTCGACATCAGTGGGTTTTCATCCACTACGAAATCACTACCGAATAAATCACTGCATATACGCTGCGCATCTTCAATTTTTGAAGCTGGCAAGCAGGAGACGCCCAGCGATGGGc                                                                                                                                               >  1:2285357/1‑150 (MQ=255)
  tCGGAACGTGGAAAACGGCATTGGTATTATTCAATGTCGACATCAGTGGGTTTTCATCCACTACGAAATCACTACCGAATAAATCACTGCATATACGCTGCGCATCTTCAATTTTTGAAGCTGGCAAGCAGGAGACGCCCAGCGATGGGc                                                                                                                                               <  1:1456565/150‑1 (MQ=255)
           ggAAAACGGCATTGGTATTATTCAATGTCGACATCAGTGGGTTTTCATCCACTACGAAATCACTACCGAATAAATCACTGCATATACGCTGCGCATCTTCAATTTTTGAAGCTGGCAAGCAGGAGACGCCCAGCGATGGGCGACGCACCa                                                                                                                                      >  1:890017/1‑150 (MQ=255)
                  ggCATTGGTATTATTCAATGTCGACATCAGTGGGTTTTCATCCACTACGAAATCACTACCGAATAAATCACTGCATATACGCTGCGCATCTTCAATTTTTGAAGCTGGCAAGCAGGAGACGCCCAGCGATGGGCGACGCACCATGATAtt                                                                                                                               >  2:826946/1‑150 (MQ=255)
                     aTTGGTATTATTCAATGTCGACATCAGTGGGTTTTCATCCACTACGAAATCACTACCGAATAAATCACTGCATATACGCTGCGCATCTTCAATTTTTGAAGCTGGCAAGCAGGAGACGCCCAGCGATGGGCGACGCACCATGATATTGAc                                                                                                                            <  1:1877299/150‑1 (MQ=255)
                     aTTGGTATTATTCAATGTCGACATCAGTGGGTTTTCATCCACTACGAAATCACTACCGAATAAATCACTGCATATACGCTGCGCATCTTCAATTTTTGAAGCTGGCAAGCAGGAGACGCCCAGCGATGGGCGACGCACCATGATATTGAc                                                                                                                            >  2:1204225/1‑150 (MQ=255)
                        ggTATTATTCAATGTCGACATCAGTGGGTTTTCATCCACTACGAAATCACTACCGAATAAATCACTGCATATACGCTGCGCATCTTCAATTTTTGAAGCTGGCAAGCAGGAGACGCCCAGCGATGGGCGACGCACCATGATATTGACCtt                                                                                                                         <  1:2369726/150‑1 (MQ=255)
                        ggTATTATTCAATGTCGACATCAGTGGGTTTTCATCCACTACGAAATCACTACCGAATAAATCACTGCATATACGCTGCGCATCTTCAATTTTTGAAGCTGGCAAGCAGGAGACGCCCAGCGATGGGCGACGCACCATGATATTGACCtt                                                                                                                         <  2:2285357/150‑1 (MQ=255)
                        ggTATTATTCAATGTCGACATCAGTGGGTTTTCATCCACTACGAAATCACTACCGAATAAATCACTGCATATACGCTGCGCATCTTCAATTTTTGAAGCTGGCAAGCAGGAGACGCCCAGCGATGGGCGACGCACCATGATATTGACCtt                                                                                                                         <  1:1208359/150‑1 (MQ=255)
                          tattatTCAATGTCGACATCAGTGGGTTTTCATCCACTACGAAATCACTACCGAATAAATCACTGCATATACGCTGCGCATCTTCAATTTTTGAAGCTGGCAAGCAGGAGACGCCCAGCGATGGGCGACGCACCATGATATTGACCTTGt                                                                                                                       >  1:1664472/1‑150 (MQ=255)
                               ttCAATGTCGACATCAGTGGGTTTTCATCCACTACGAAATCACTACCGAATAAATCACTGCATATACGCTGCGCATCTTCAATTTTTGAAGCTGGCAAGCAGGAGACGCCCATCGCTGGGCgtctcct                                                                                                                                        >  1:1095185/1‑122 (MQ=255)
                               ttCAATGTCGACATCAGTGGGTTTTCATCCACTACGAAATCACTACCGAATAAATCACTGCATATACGCTGCGCATCTTCAATTTTTGAAGCTGGCAAGCAGGAGACGCCCATCGCTGGGCgtctcct                                                                                                                                        <  2:1095185/128‑7 (MQ=255)
                               ttCAATGTCGACATCAGTGGGTTTTCATCCACTACGAAATCACTACCGAATAAATCACTGCATATACGCTGCGCATCTTCAATTTTTGAAGCTGGCAAGCAGGAGACGCCCAGCGATGGGCGACGCACCATGATATTGACCTTGTCTTcc                                                                                                                  >  2:151970/1‑150 (MQ=255)
                                   aTGTCGACATCAGTGGGTTTTCATCCACTACGAAATCACTACCGAATAAATCACTGCATATACGCTGCGCATCTTCAATTTTTGAAGCTGGCAAGCAGGAGACGCCCAGCGATGGGCGACGCACCATGATATTGACCTTGTCTTCCGaaa                                                                                                              >  2:366036/1‑150 (MQ=255)
                                   aTGTCGACATCAGTGGGTTTTCATCCACTACGAAATCACTACCGAATAAATCACTGCATATACGCTGCGCATCTTCAATTTTTGAAGCTGGCAAGCAGGAGACGCCCAGCGATGGGCGACGCACCATGATATTGACCTTGTCTTCCGaaa                                                                                                              <  1:826946/150‑1 (MQ=255)
                                   aTGTCGACATCAGGGGGTTTTCATCCACTACGAAATCACTACCGAATAAATCACTGCATATACGCTGCGCATCTTCAATTTTTGAAGCTGGCAAGCAGGAGACGCCCAGCGATGGGCGACGCACCATGATATTGACCTTGTCTTCCGaaa                                                                                                              <  2:1333696/150‑1 (MQ=255)
                                      tCGACATCAGTGGGTTTTCATCCACTACGAAATCACTACCGAATAAATCACTGCATATACGCTGCGCATCTTCAATTTTTGAAGCTGGCAAGCAGGAGACGCCCAGCGATGGGCGACGCACCATGATATTGACCTTGTCTTCCGAAAccc                                                                                                           >  2:1921965/1‑150 (MQ=255)
                                         aCATCAGTGGGTTTTCATCCACTACGAAATCACTACCGAATAAATCACTGCATATACGCTGCGCATCTTCACTTTTTGAAGCTGGCAAGCAGGAGACGCCCAGCGATGGGCGACGCACCATGATATTGACCTTGTCTTCCGAAACCCTTc                                                                                                        <  2:1664472/150‑1 (MQ=255)
                                         aCATCAGTGGGTTTTCATCCACTACGAAATCACTACCGAATAAATCACTGCATATACGCTGCGCATCTTCAATTTTTGAAGCTGGCAAGCAGGAGACGCCCAGCGATGGGCGACGCACCATGATATTGACCTTGTCTTCCGAAACCCTTc                                                                                                        >  2:1400506/1‑150 (MQ=255)
                                               gTGGGTTTTCATCCACTACGAAATCACTACCGAATAAATCACTGCATATACGCTGCGCATCTTCAATTTTTGAAGCTGGCAAGCAGGAGACGCCCAGCGATGGGCGACGCACCATGATAtt                                                                                                                               >  2:2349302/1‑121 (MQ=255)
                                               gTGGGTTTTCATCCACTACGAAATCACTACCGAATAAATCACTGCATATACGCTGCGCATCTTCAATTTTTGAAGCTGGCAAGCAGGAGACGCCCAGCGATGGGCGACGCACCATGATAtt                                                                                                                               <  1:2349302/121‑1 (MQ=255)
                                                tGGGTTTTCATCCACTACGAAATCACTACCGAATAAATCACTGCATATACGCTGCGCATCTTCAATTTTTGAAGCTGGCAAGCAGGAGACGCCCAGCGATGGGCGACGCACCATGATATTGACCTTGTCTTCCGAAACCCTTCGTGCGGt                                                                                                 >  2:799647/1‑150 (MQ=255)
                                                    ttttCATCCACTACGAAATCACTACCGAATAAATCACTGCATATACGCTGCGCATCTTCAATTTTTGAAGCTGGCAAGCAGGAGACGCCCAGCGATGGGCGACGCACCATGATATTGACCTTGTCTTCCGAAACCCTTCGTGCGGTAAGt                                                                                             >  2:154371/1‑150 (MQ=255)
                                                     tttCATCCACTACGAAATCACTACCGAATAAATCACTGCATATACGCTGCGCATCTTCAATTTTTGAAGCTGGCAAGCAGGAGACGCCCAGCGATGGGCGACGCACCATGATATTGACCTTGTCTTCCGAAACCCTTCGTGCGGTAAGTg                                                                                            <  1:1204225/150‑1 (MQ=255)
                                                        cATCCACTACGAAATCACTACCGAATAAATCACTGCATATACGCTGCGCATCTTCAATTTTTGAAGCTGGCAAGCAGGAGACGCCCAGCGATGGGCGACGCACCATGATATTGACCTTGTCTTCCGAAACCCTTCGTGCGGTAAGTGCAg                                                                                         >  1:2326250/1‑150 (MQ=255)
                                                            cACTACGAAATCACTACCGAATAAATCACTGCATATACGCTGCGCATCTTCAATTTTTGAAGCTGGCAAGAAGCAGACGCCCAGCGATGGGCGACGCACCATGATATTGACCTTGTCTTCCGAAACCCTTCGTGCGGTAAGTGCAGTtgt                                                                                     <  1:1738355/150‑1 (MQ=255)
                                                             aCTACGAAATCACTACCGAATAAATCACTGCATATACGCTGCGCATCTTCAATTTTTGAAGCTGGCAAGCAGGAGACGCCCAGCGATGGGCGACGCACCATGATATTGACCTTGTCTTCCGAAACCCTTCGTGCGGTAAGTGCAGTtgtg                                                                                    <  1:799647/150‑1 (MQ=255)
                                                             aCTACGAAATCACTACCGAATAAATCACTGCATATACGCTGCGCATCTTCAATTTTTGAAGCTGGCAAGCAGGAGACGCCCAGCGATGGGCGACGCACCATGATATTGACCTTGTCTTCCGAAACCCTTCGTGCGGTAAGTGCAGTtgtg                                                                                    <  1:366036/150‑1 (MQ=255)
                                                               tACGAAATCACTACCGAATAAATCACTGCATATACGCTGCGCATCTTCAATTTTTGAAGCTGGCAAGCAGGAGACGCCCAGCGATGGGCGACGCACCATGATATTGACCTTGTCTTCCGAAACCCTTCGTGCGGTAAGTGCAGTTGTGcc                                                                                  >  1:1303182/1‑150 (MQ=255)
                                                                  gAAATCACTACCGAATAAATCACTGCATATACGCTGCGCATCTTCAATTTTTGAAGCTGGCAAGCAGGAGACGCCCAGCGATGGGCGACGCACCATGATATTGACCTTGTCTTCCGAAACCCTTCGTGCGGTaa                                                                                               >  2:207650/1‑134 (MQ=255)
                                                                  gAAATCACTACCGAATAAATCACTGCATATACGCTGCGCATCTTCAATTTTTGAAGCTGGCAAGCAGGAGACGCCCAGCGATGGGCGACGCACCATGATATTGACCTTGTCTTCCGAAACCCTTCGTGCGGTaa                                                                                               <  1:207650/134‑1 (MQ=255)
                                                                     aTCACTACCGAATAAATCACTGCATATACGCTGCGCATCTTCAATTTTTGAAGCTGGCAAGCAGGAGACGCCCAGCGATGGGCGACGCACCATGATATTGACCTTGTCTTCCGAAACCCTTCGTGCGGTAAGTGCAGTTGTGCCAAAACt                                                                            <  2:1303182/150‑1 (MQ=255)
                                                                        aCTACCGAATAAATCACTGCATATACGCTGCGCATCTTCAATTTTTGAAGCTGGCAAGCAGGAGACGCCCAGCGATGGGCGACGCACCATGATATTGACCTTGTCTTCCGAAACCCTTCGTGCGGTAAGTGCAGTTGTGCCAAAACTGGc                                                                         <  1:1003205/150‑1 (MQ=255)
                                                                          tACCGAATAAATCACTGCATATACGCTGCGCATCTTCAATTTTTGAAGCTGGCAAGCAGGAGACGCCCAGCGATGGGCGACGCACCATGATATTGACCTTGTCTTCCGAAACCCTTCGTGCGGTAAGTGCAGTTGTGCCAAAACTGGcga                                                                       >  2:1912981/1‑150 (MQ=255)
                                                                          tACCGAATAAATCACTGCATATACGCTGCGCATCTTCAATTTTTGAAGCTGGCAAGCAGGAGACGCCCAGCGATGGGCGACGCACCATGATATTGACCTTGTCTTCCGAAACCCTTCGTGCGGTAAGTGCAGTTGTGCCAAAACTGGcga                                                                       <  1:154371/150‑1 (MQ=255)
                                                                          tACCGAATAAATCACTGCATATACGCTGCGCATCTTCAATTTTTGAAGCTGGCAAGCAGGAGACGCCCAGCGATGGGCGACGCACCATGATATTGACCTTGTCTTCCGAAACCCTTCGTGCGGTAAGTGCAGTTGTGCCAAAACTGGcga                                                                       >  1:143790/1‑150 (MQ=255)
                                                                          tACCGAATAAATCACTGCATATACGCTGCGCATCTTCAATTTTTGAAGCTGGCAAGCAGGAGACGCCCAGCGATGGGCGACGCACCATGATATTGACCTTGTCTTCCGAAACCCTTCGTGCGGTAAGTGCAGTTGTGCCAAAACTGGcga                                                                       >  2:927173/1‑150 (MQ=255)
                                                                            ccGAATAAATCACTGCATATACGCTGCGCATCTTCAATTTTTGAAGCTGGCAAGCAGGAGACGCCCAGCGATGGGCGACGCACCATGATATTGACCTTGTCTTCCGAAACCCTTCGTGCGGTAAGTGCAGTTGTGCCAAAACTGGcgacg                                                                     >  1:1505387/1‑150 (MQ=255)
                                                                             cGAATAAATCACTGCATATACGCTGCGCATCTTCAATTTTTGAAGCTGGCAAGCAGGAGACGCCCAGCGATGGGCGACGCACCATGATATTGACCTTGTCTTCCGAAACCCTTCGTGCGGTAAGTGCAGTTGTGCCAAAACTGGcgacga                                                                    <  1:1912981/150‑1 (MQ=255)
                                                                              gAATAAATCACTGCATATACGCTGCGCATCTTCAATTTTTGAAGCTGGCAAGCAGGAGACGCCCAGCGATGGGCGACGCACCATGATATTGACCTTGTCTTc                                                                                                                   >  1:2036641/1‑102 (MQ=255)
                                                                              gAATAAATCACTGCATATACGCTGCGCATCTTCAATTTTTGAAGCTGGCAAGCAGGAGACGCCCAGCGATGGGCGACGCACCATGATATTGACCTTGTCTTc                                                                                                                   <  2:2036641/102‑1 (MQ=255)
                                                                              gAATAAATCACTGCATATACGCTGCGCATCTTCAATTTTTGAAGCTGGCAAGCAGGAGACGCCCAGCGATGGGCGACGCACCATGATATTGACCTTGTCTTCCGAAACCCTTCGTGCGGTAAGTGCAGTTGTGCCAAAACTGGCgacgac                                                                   >  1:153986/1‑150 (MQ=255)
                                                                                aTAAATCACTGCATATACGCTGCGCATCTTCAATTTTTGAAGCTTGCAAGCAGGAGACGCCCAGCGATGGGCGACGCACCATGATATTGACCTTGGCTTCCGAAACCCTTCGTGCGGTAAGTGCAGTTGTGCCAAAACTGGCGACGACaa                                                                 <  1:1033946/150‑1 (MQ=255)
                                                                                    aTCACTGCATATACGCTGCGCATCTTCAATTTTTGAAGCTGGCAAGCAGGAGACGCCCAGCGATGGGCGACGCACCATGATATTGACCTTGTCTTCCGAAACCCTTCGTGCGGTAAGTGCAGTTGTGCCAAAACTGGCGACGACAAGATc                                                             >  2:1913486/1‑150 (MQ=255)
                                                                                    aTCACTGCATATACGCTGCGCATCTTCAATTTTTGAAGCTGGCAAGCAGGAGACGCCCAGCGATGGGCGACGCACCATGATATTGACCTTGTCTTCCGAAACCCTTCGTGCGGTAAGTGCAGTTGTGCCAAAACTGGCGACGACAAGATc                                                             >  1:554885/1‑150 (MQ=255)
                                                                                     tCACTGCATATACGCTGCGCATCTTCAATTTTTGAAGCTGGCAAGCAGGAGACGCCCAGCGATGGGCGACGCACCATGATATTGACCTTGTCTTCCGAAACCCTTCGTGCGGTAAGTGCAGTTGTGCCAAAACTGGCGACGACAAGATCt                                                            >  2:2315977/1‑150 (MQ=255)
                                                                                        cTGCATATACGCTGCGCATCTTCAATTTTTGAAGCTGGCAAGCAGGAGACGCCCAGCGATGGGCGACGCACCATGATATTGACCTTGTCTTCCGAAACCCTTCGTGCGGTAAGTGCAGTtgt                                                                                     >  1:299426/1‑122 (MQ=255)
                                                                                        cTGCATATACGCTGCGCATCTTCAATTTTTGAAGCTGGCAAGCAGGAGACGCCCAGCGATGGGCGACGCACCATGATATTGACCTTGTCTTCCGAAACCCTTCGTGCGGTAAGTGCAGTtgt                                                                                     <  2:299426/122‑1 (MQ=255)
                                                                                        cTGCATATACGCTGCGCATCTTCAATTTTTGAAGCTGGCAAGCAGGAGACGCCCAGCGATGGGCGACGCACCATGATATTGACCTTGTCTTCCGAAACCCTTCGTGCGGTAAGTGCAGTTGTGCCAAAACTGGcgac                                                                      <  1:1146030/137‑1 (MQ=255)
                                                                                        cTGCATATACGCTGCGCATCTTCAATTTTTGAAGCTGGCAAGCAGGAGACGCCCAGCGATGGGCGACGCACCATGATATTGACCTTGTCTTCCGAAACCCTTCGTGCGGTAAGTGCAGTTGTGCCAAAACTGGcgac                                                                      >  2:1146030/1‑137 (MQ=255)
                                                                                          gCATATACGCTGCGCATCTTCAATTTTTGAAGCTGGCAAGCAGGAGACGCCCAGCGATGGGCGACGCACCATGATATTGACCTTGTCTTCCGAAACCCTTCGTGCGGTAAGTGCAGTTGTGCCAAAACTGGCGACGACAAGATCTAcccc                                                       >  1:2263252/1‑150 (MQ=255)
                                                                                           cATATACGCTGCGCATCTTCAATTTTTGAAGCTGGCAAGCAGGAGACGCCCAGCGATGGGCGACGCACCATGATATTGACCTTGTCTTCCGAAACCCTTCGTGCGGTAAGTGCAGTTGTGCCAAAACTGGCGACGACAAGATCTAccccc                                                      <  1:1913486/150‑1 (MQ=255)
                                                                                            atatACGCTGCGCATCTTCAATTTTTGAAGCTGGCAAGCAGGAGACGCCCAGCGATGGGCGACGc                                                                                                                                          >  2:1154963/1‑65 (MQ=255)
                                                                                            atatACGCTGCGCATCTTCAATTTTTGAAGCTGGCAAGCAGGAGACGCCCAGCGATGGGCGACGc                                                                                                                                          <  1:1154963/65‑1 (MQ=255)
                                                                                             tataCGCTGCGCATCTTCAATTTTTGAAGCTGGCAAGCAGGAGACGCCCAGCGATGGGCGACGCACCATGATATTGACCTTGTCTTCCGAAACCCTTCGTGCGGTAAGTGCAGTTGTGCCAAAACTGGCGACGACAAGATCTACCCCCAg                                                    <  2:2263252/150‑1 (MQ=255)
                                                                                                 cGCTGCGCATCTTCAATTTTTGAAGCTGGCAAGCAGGAGACGCCCAGCGATGGGCGACGCACCATGATATTGACCTTGTCTTCCGAAACCCTTCGTGCGGTAAGTGCAGTTGTGCCAAAACTGGCGACGACAAGATCTACCCCCAGGTTc                                                <  1:1921965/150‑1 (MQ=255)
                                                                                                     gcgcATCTTCAATTTTTGAAGCTGGCAAGCAGGAGACGCCCAGCGATGGGCGACGCACCATGATATTg                                                                                                                              <  2:1016116/68‑1 (MQ=255)
                                                                                                     gcgcATCTTCAATTTTTGAAGCTGGCAAGCAGGAGACGCCCAGCGATGGGCGACGCACCATGATATTg                                                                                                                              >  1:1016116/1‑68 (MQ=255)
                                                                                                       gcATCTTCAATTTTTGAAGCTGGCAAGCAGGAGACGCCCAGCGATGGGCGACGCACCATGATATTGACCTTGTCTTCCGAAACCCTTCGTGCGGTAAGTGCAGTTGTGCCAAAACTGGCGACGACAAGATCTACCCCCAGGTTCCTcgcg                                          >  2:2065015/1‑150 (MQ=255)
                                                                                                         aTCTTCAATTTTTGAAGCTGGCAAGCAGGAGACGCCCAGCGATGGGCGACGCACCATGATATTGACCTTGTCt                                                                                                                     >  2:2254361/1‑73 (MQ=255)
                                                                                                         aTCTTCAATTTTTGAAGCTGGCAAGCAGGAGACGCCCAGCGATGGGCGACGCACCATGATATTGACCTTGTCt                                                                                                                     <  1:2254361/73‑1 (MQ=255)
                                                                                                         aTCTTCAATTTTTGAAGCTGGCAAGCAGGAGACGCCCAGCGATGGGCGACGCACCATGATATTGACCTTGTCTTCCGAAACCCTTCGTGCGGTAAGt                                                                                             <  1:1751841/97‑1 (MQ=255)
                                                                                                         aTCTTCAATTTTTGAAGCTGGCAAGCAGGAGACGCCCAGCGATGGGCGACGCACCATGATATTGACCTTGTCTTCCGAAACCCTTCGTGCGGTAAGt                                                                                             >  2:1751841/1‑97 (MQ=255)
                                                                                                            ttCAATTTTTGAAGCTGGCAAGCAGGAGACGCCCAGCGATGGGCGACGCACCATGATATTGACCTTGTCTTCCGAAACCCTTCGTGCGGTAAGTGCAGTTGTGCCAAAACTGGCGACGACAAGATCTACCCCCAGGTTCCTCGCGCGCTc                                     >  2:1592538/1‑150 (MQ=255)
                                                                                                            ttCAATTTTTGAAGCTGGCAAGCAGGAGACGCCCAGCGATGGGCGACGCACCATGATATTGACCTTGTCTTCCGAAACCCTTCGTGCGGTAAGTGCAGTTGTGCCAAAACTGGCGACGACAAGATCTACCCCCAGGTTCCTCGCGCGCTc                                     >  1:2187672/1‑150 (MQ=255)
                                                                                                            ttCAATTTTTGAAGCTGGCAAGCAGGAGACGCCCAGCGATGGGCGACGCACCATGATATTGACCTTGTCTTCCGAAACCCTTCGTGCGGTAAGTGCAGTTGTGCCAAAACTGGCGACGACAAGATCTACCCCCAGGTTCCTCGCGCGCTc                                     >  2:1766647/1‑150 (MQ=255)
                                                                                                             tCAATTTTTGAAGCTGGCAAGCAGGAGACGCCCAGCGATGGGCGACGCACCATGATATTGACCTTGTCTTCCGAAACCCTTCGTGCGGTAAGTGCAGTTGTGCCAAAACTGGCGACGACAAGATCTACCCCCAGGTTCCTCGCGCGCTcc                                    >  1:141547/1‑150 (MQ=255)
                                                                                                             tCAATTTTTGAAGCTGGCAAGCAGGAGACGCCCAGCGATGGGCGACGCACCATGATATTGACCTTGTCTTCCGAAACCCTTCGTGCGGTAAGTGCAGTTGTGCCAAAACTGGCGACGACAAGATCTACCCCCAGGTTCCTCGCGCGCTcc                                    <  1:1766647/150‑1 (MQ=255)
                                                                                                             tCAATTTTTGAAGCTGGCAAGCAGGAGACGCCCAGCGATACGCGAAACAGCATGAAATTGACCTTGTCTTCCGAAACCCTTCGTGCGGTAAGTGCAGTTGTGCCAAAACTGGCGACGACAAGATCTACCCCCAGGTTCCTCGCGCGCTcc                                    >  1:1025585/1‑150 (MQ=255)
                                                                                                                 tttttGAAGCTGGCAAGCAGGAGACGCCCAGCGATGGTCGACGCACCATGATATTGACCTTGTCTTCCGAAACCCTTCGTGCGGTAAGTGCAGTTGTGCCAAAACTGGCGACGACAACATCTACCCCCAGGTTCCTCGCGCGCTCCCGaa                                <  1:822739/150‑1 (MQ=255)
                                                                                                                 tttttGAAGCTGGCAAGCAGGAGACGCCCAGCGATGGGCGACGCACCATGATATTGACCTTGTCTTCCGAAACCCTTCGTGCgg                                                                                                  >  2:959595/1‑84 (MQ=255)
                                                                                                                 tttttGAAGCTGGCAAGCAGGAGACGCCCAGCGATGGGCGACGCACCATGATATTGACCTTGTCTTCCGAAACCCTTCGTGCgg                                                                                                  <  1:959595/84‑1 (MQ=255)
                                                                                                                 tttttGAAGCTGGCAAGCAGGAGACGCCCAGCGATGGGCGACGCACCATGATATTGACCTTGTCTTCCGAAACCCTTCGTGCGGTAAGTGCAGTTGTGCCAAAACTGGCGACGACAAGATCTACCCCCAGGTTCCTCGCGCGCTCCCGaa                                <  1:1400506/150‑1 (MQ=255)
                                                                                                                   tttGAAGCTGGCAAGCAGGAGACGCCCAGCGATGGGCGACGCACCATGATATTGACCTTGTCTTCCGAAACCCTTCGTGCGGTAAGTGCAGTTGTGCCAAAACTGGCGACGACAAGATCTACCCCCAGGTTCCTCGCGCGCTCCCGAAGa                              <  1:2065015/150‑1 (MQ=255)
                                                                                                                      gAAGCTGGCAAGCAGGAGACGCCCAGCGATGGGCGACGCACCATGATATTGACCTTGTCTTCCGaaa                                                                                                              <  1:2089233/67‑1 (MQ=255)
                                                                                                                      gAAGCTGGCAAGCAGGAGACGCCCAGCGATGGGCGACGCACCATGATATTGACCTTGTCTTCCGaaa                                                                                                              >  2:2089233/1‑67 (MQ=255)
                                                                                                                       aaGCTGGCAAGCAGGAGACGCCCAGCGATGGGCGACGCACCATGATATTGACCTTGTCTTCCGAAACCCTTCGTGCGGTAAGTGCAGTTGTGCCAAAACTGGCGACGACAAGATCTACCCCCAGGTTCCTCGCGCGCTCCCGAAGATaaa                          <  2:2187672/150‑1 (MQ=255)
                                                                                                                           tGGCAAGCAGGAGACGCCCAGCGATGGGCGACGCACCATGATATTGACCTTGTCTTCCGAAACCCTTCGTGCGGTAAGTGCAGTTGTGCCAAAACTGGCGACGACAAGATCTACCCCCAGGTTCCTCGCGCGCTCCCGAAGATAAAgcgc                      >  2:1688917/1‑150 (MQ=255)
                                                                                                                              cAAGCAGGAGACGCCCAGCGATGGGCGACGCACCATGATATTGACCTTGTCTTCCGAAACCCTTCGTGCGGTAAGTGCAGTTGTGCCAAAACTGGCGACGACAAGATCTACCCCCAGGTTCCTCGCGCGCTCCCGAAGATAAAGCGCCGa                   <  2:143790/150‑1 (MQ=255)
                                                                                                                               aaGCAGGAGACGCCCAGCGATGGGCGACGCACCATGATATTGACCTTGTCTTCCGAAACCCTTCGTGCGGTAAGTGCAGTTGTGCCAAAACTGGCGACGACAAGATCTACCCCCAGGTTCCTCGCGCGCTCCCGAAGATAAAGCGCCGAc                  >  1:196573/1‑150 (MQ=255)
                                                                                                                                 gCAGGAGACGCCCAGCGATGGGCGACGCACCATGATATTGACCTTGTCTTCCGAAACCCTTCGTGCGGTAAGTGCAGTTGTGCCAAAACTGGCGACGACAAGATCTACCCCCAGGTTCCTcgctc                                         <  2:115620/125‑3 (MQ=255)
                                                                                                                                 gCAGGAGACGCCCAGCGATGGGCGACGCACCATGATATTGACCTTGTCTTCCGAAACCCTTCGTGCGGTAAGTGCAGTTGTGCCAAAACTGGCGACGACAAGATCTACCCCCAGGTTCCTcgctc                                         >  1:115620/1‑123 (MQ=255)
                                                                                                                                  cAGGAGACGCCCAGCGATGGGCGACGCACCATGATATTGACCTTGTCTTCCGAAACCCTTCGTGCGGTAAGTGCAGTTGTGCCAAAACTGGCGACGACAAGATCTACCCCCAGGTTCCTCGCGCGCTCCCGAAGATAAAGCGCCGACAAc               <  2:890017/150‑1 (MQ=255)
                                                                                                                                    tgcgtcGCCCATCGCTGGGCGACGCACCATGATATTGACCTTGTCTTCCGAAACCCTTCGTGCGGTAAGTGCAGTTGTGCCAAAACTGGCGACGACAAGATCTACCCCCAGGTTCCTCGCGCGCTCCCGAAGATAAAGCGCCGACAAc               >  1:293412/6‑148 (MQ=255)
                                                                                                                                    tgcgtcGCCCATCGCTGGGCGACGCACCATGATATTGACCTTGTCTTCCGAAACCCTTCGTGCGGTAAGTGCAGTTGTGCCAAAACTGGCGACGACAAGATCTACCCCCAGGTTCCTCGCGCGCTCCCGAAGATAAAGCGCCGACAAc               <  2:293412/143‑1 (MQ=255)
                                                                                                                                    ggAGACGCCCAGCGATGGGCGACGCACCATGATATTGACCTTGTCTTCCGAAACCCTTCGTGCGGTAAGTGCAGTTGTGCCAAAACTGGCGACGACAAGATCTACCCCCAGGTTCCTCGCGCGCTCCCGAAGATAAAGCGCCGACAACGa             <  1:1795680/150‑1 (MQ=255)
                                                                                                                                       gtcGCCCATCGCTGGGCGACGCACCATGATATTGACCTTGTCTTCCGAAACCCTTCGTGCGGTAAGTGCAGTTGTGCCAAAACTGGCGACGACAAGATCTACCCCCAGGTTCCTCGCGCGCTCCCGAAGATAAAGCGCCGACAACGAGcc          >  2:1512641/3‑150 (MQ=255)
                                                                                                                                       gaCGCCCAGCGATGGGCGACGCACCATGATATTGACCTTGTCTTCCGAAACCCTTCGTGCGGTAAGTGCAGTTGTGCCAAAACTGGCGACGACAAGATCTAc                                                          <  1:33113/102‑1 (MQ=255)
                                                                                                                                       gaCGCCCAGCGATGGGCGACGCACCATGATATTGACCTTGTCTTCCGAAACCCTTCGTGCGGTAAGTGCAGTTGTGCCAAAACTGGCGACGACAAGATCTAc                                                          >  2:33113/1‑102 (MQ=255)
                                                                                                                                           cccAGCGATGGGCGACGCACCATGATATTGACCTTGTCTTCCGAAACCCTTCGTGCGGTAAGTGCAGTTGTGCCAAAACTGGCGACGACAAGATCTACCCCCAGGTTCCTCGCGCGCTCCCGAAGATAAAGCGCCGACAACGAGCCCAtt      <  1:2315977/150‑1 (MQ=255)
                                                                                                                                              aGCGATGGGCGACGCACCATGATATTGACCTTGTCTTCCGAAACCCTTCGTGCGGTAAGTGCAGTTGTGCCAAAACTGGCGACGACAAGATCTACCCCCAGGTTCCTCGCGCGCTCCCGAAGATAAAGCGCCGACAACGAGCCCATTGaa   <  2:153986/150‑1 (MQ=255)
                                                                                                                                               gCGATGGGCGACGCACCATGATATTGACCTTGTCTTCCGAAACCCTTCGTGCGGTAAGGGCAGTTGTGCCAAAACTGGCGACGACAAGATCTACCCCCAGGTTCCTCGCGCGCTCCCGAAGATAAATCGCCGACAACGAGCCCATTGAAc  <  2:554885/150‑1 (MQ=255)
                                                                                                                                               |                                                                                                                                                     
AATCGGAACGTGGAAAACGGCATTGGTATTATTCAATGTCGACATCAGTGGGTTTTCATCCACTACGAAATCACTACCGAATAAATCACTGCATATACGCTGCGCATCTTCAATTTTTGAAGCTGGCAAGCAGGAGACGCCCAGCGATGGGCGACGCACCATGATATTGACCTTGTCTTCCGAAACCCTTCGTGCGGTAAGTGCAGTTGTGCCAAAACTGGCGACGACAAGATCTACCCCCAGGTTCCTCGCGCGCTCCCGAAGATAAAGCGCCGACAACGAGCCCATTGAAC  >  NC_002947/5051953‑5052245

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