Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,447,887 C→G 11.3% intergenic (‑279/‑43) nicC ← / → nicX 6‑hydroxynicotinate 3‑monooxygenase/2,5‑dihydroxypyridine 5,6‑dioxygenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,447,8870CG11.3% 261.7 / 19.0 88intergenic (‑279/‑43)nicC/nicX6‑hydroxynicotinate 3‑monooxygenase/2,5‑dihydroxypyridine 5,6‑dioxygenase
Reads supporting (aligned to +/- strand):  ref base C (44/34);  new base G (4/6);  total (48/40)
Fisher's exact test for biased strand distribution p-value = 5.02e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.74e-01

AATCCGGTTGATTGTTACCGTAAATACTGGTCAATGTGTGGATAGGTAACGTTTGAACAGATCACTTGCGCTCGCTGGCAGGATGCGTTTTCATCTGCACAAGTAGGGTGTACACACTATTATTTTGCACCGAGATGCAGCACCAGGCACTGCCATGGTGCGCCCCTGTCAACGCCACCCAGGAGAACCCCGATGCCGGTGAGCAATGCACAACTG‑‑ACCCAGATGTTCGAGCACGTACTGAAGCTGTCTCGAGTGGACGAAACGCAGAGCGTGGCCGTGCTCAAGAGCCACTACTC  >  NC_002947/4447738‑4448033
                                                                                                                                                      |                                                                                                                                                     
aaTCCGGTTGATTGTTACCGTAAATACTGGTCAATGTGTGGATAGGTAACGTTTGAACAGATCACTTGCGCTCGCTGGCAGGATGCGTTTTCATCTGCACAAGTAGGGTGTACACACTATTATTTTGCACCGAGATGCAGCACCAGGCAc                                                                                                                                                       <  1:752519/150‑1 (MQ=255)
      gTTGATTGTTACCGTAAATACTGGTCAATGTGTGGATAGGTAACGTTTGAACAGATCACTTGCGCTCGCTGGCAGGATGCGTTTTCATCTGCACAAGTAGGGTGTACACACTATTATTTTGCACCGAGATGCAGCACCAGGCACTGCCAt                                                                                                                                                 >  2:716115/1‑150 (MQ=255)
       ttGATTGTTACCGTAAATACTGGTCAATGTGTGGATAGGTAACGTTTGAACAGATCACTTGCGCTCGCTGGCAGGATGCGTTTTCATCTGCACAAGTAGGGTGTACACACTATTATTTTGCACCGAGATGCAGCACCAGGCACTGCCATg                                                                                                                                                <  2:128243/150‑1 (MQ=255)
        tGATTGTTACCGTAAATACTGGTCAATGTGTGGATAGGTAACGTTTGAACAGATCACTTGCGCTCGCTGGCAGGATGCGTTTTCATCTGCACAAGTAGGGTGTACACACTATTATTTTGCACCGAGATGCAGCACCATGGCAGTGcctgg                                                                                                                                                >  2:1063656/1‑147 (MQ=255)
        tGATTGTTACCGTAAATACTGGTCAATGTGTGGATAGGTAACGTTTGAACAGATCACTTGCGCTCGCTGGCAGGATGCGTTTCATCTGCACAAGTAGGGTGTACACACTATTATTTTGCACCGAGATGCAGCACCATGGCAGTGcctggt                                                                                                                                               >  1:714087/1‑146 (MQ=255)
         gATTGTTACCGTAAATACTGGTCAATGTGTGGATAGGTAACGTTTGAATAGATCACTTGCGCTCGCTGGCAGGATGCGTTTTCATCTGCACAAGTAGGGTGTACACACTATTATTTTGCACCGAGATGCAGCACCAGGCACTGCCATGGt                                                                                                                                              <  2:2008464/150‑1 (MQ=255)
         gATTGTTACCGTAAATACTGGTCAATGTGTGGATAGGTAACGTTTGAACAGATCACTTGCGCTCGCTGGCAGGATGCGTTTTCATCTGCACAAGTAGGGTGTACACACTATTATTTTGCACCGAGATGCAGCACCAGGCACTGCCATGGt                                                                                                                                              <  2:1428962/150‑1 (MQ=255)
          aTTGTTACCGTAAATACTGGTCAATGTGTGGATAGGTAACGTTTGAACAGATCACTTGCGCTCGCTGGCAGGATGCGTTTCATCTGCACAAGTAGGGTGTACACACTATTATTTTGCACCGAGATGCAGCACCATGGCAGTGCCTGGTgc                                                                                                                                            <  2:714087/150‑1 (MQ=255)
           ttGTTACCGTAAATACTGGTCAATGTGTGGATAGGTAACGTTTGAACAGATCACTTGCGCTCGCTGGCAGGATGCGTTTTCATCTGCACAAGTAGGGTGTACACACTATTATTTTGCACCGAGATGCAGCACCAGGCACTGCCATGGTgc                                                                                                                                            >  1:1920413/1‑150 (MQ=255)
            tGTTACCGTAAATACTGGTCAATGTGTGGATAGGTAACGTTTGAACAGATCACTTGCGCTCGCTGGCAGGATGCGTTTTCATCTGCACAAGTAGGGTGTACACACTATTATTTTGCACCGAGATGCAGCACCATGGCAGTGCCTGGTgct                                                                                                                                           <  1:1063656/150‑2 (MQ=255)
            tGTTACCGTAAATACTGGTCAATGTGTGGATAGGTAACGTTTGAACAGATCACTTGCGCTCGCTGGCAGGATGCGTTTTCATCTGCACAAGTAGGGTGTACACACTATTATTTTGCACCGAGATGCAGCACCATGGCAGTGCCTGGTgct                                                                                                                                           <  1:1165664/150‑2 (MQ=255)
                 ccGTAAATACTGGTCAATGTGTGGATAGGTAACGTTTGAACAGATCACTTGCGCTCGCTGGCAGGATGCGTTCTACAGAGCACAAGTAAGTTTTACACACTATTATTTTGCACCGAGATGCAGCACCATGGCAGTGCCTGGTgctgcatc                                                                                                                                      <  1:2099839/150‑7 (MQ=255)
                   gTAAATACTGGTCAATGTGTGGATAGGTAACGTTTGAACAGATCACTTGCGCTCGCTGGCAGGATGCGTTTTCATCTGCACAAGTAGGGTGTACACACTATTATTTTGCACCGAGATGCAGCACCAGGCACTGCCATGGTGCGCCCCTGt                                                                                                                                    >  2:1913478/1‑150 (MQ=255)
                   gTAAATACTGGTCAATGTGTGGATAGGTAACGTTTGAACAGATCACTTGCGCTCGCTGGCAGGATGCGTTTTCATCTGCACAAGTAGGGTGTACACACTATTATTTTGCACCGAGATGCAGCACCAGGCACTGCCATGGTGCGCCCCTGt                                                                                                                                    >  2:1303191/1‑150 (MQ=255)
                       tccgatctCAATGTGTGGATAGGTAACGTTTGAACAGATCACTTGCGCTCGCTGGCAGGATGCGTTTTCATCTGCACAAGTAGGGTGTACACACTATTATTTTGCACCGAGATGCAGCACCAGGCACTGCCATGGTGCGCCCCTGTCAAc                                                                                                                                <  1:571444/143‑1 (MQ=255)
                               cAATGTGTGGATAGGTAACGTTTGAACAGATCACTTGCGCTCGCTGGCAGGATGCGTTTTCATCTGCACAAGTAGGGTGTACACACTATTATTTTGCACCGAGATGCAGCACCAGGCACTGCCATGGTgcgc                                                                                                                                          <  2:345062/132‑1 (MQ=255)
                               cAATGTGTGGATAGGTAACGTTTGAACAGATCACTTGCGCTCGCTGGCAGGATGCGTTTTCATCTGCACAAGTAGGGTGTACACACTATTATTTTGCACCGAGATGCAGCACCAGGCACTGCCATGGTgcgc                                                                                                                                          >  1:345062/1‑132 (MQ=255)
                                aaTGTGTGGATAGGTAACGTTTGAACAGATCACTTGCGCTCGCTGGCAGGATGCGTTTTCATCTGCACAAGTAGGGTGTACACACTATTATTTTGCACCGAGATGCAGCACCAGGCACTGCCATGGTGCGCCCCTGTCAACGCCACCCAg                                                                                                                       <  1:1183311/150‑1 (MQ=255)
                                              tAACGTTTGAACAGATCACTTGCGCTCGCTGGCAGGATGCGTTTTCATCTGCACAAGTAGGGTGTACACACTATTATTTTGCACCGAGATGCAGCACCAGGCACTGCCATGGTGCGCCCCTGTCAACGCCACCCAGGAGAACCCCGATGc                                                                                                         >  2:1671177/1‑150 (MQ=255)
                                                aCGTTTGAACAGATCACTTGCGCTCGCTGGCAGGATGCGTTTTCATCTGCACAAGTAGGGTGTACACACTATTATTTTGCACCGAGATGCAGCACCAGGCACTGCCATGGTGCGCCCCTGTCAACGCCACCCAGGAGAACCCCGATGCCg                                                                                                       <  2:2223658/150‑1 (MQ=255)
                                                   tttGAACAGATCACTTGCGCTCGCTGGCAGGATGCGTTTTCATCTGCACAAGTAGGGTGTACACACTATTATTTTGCACCGAGATGCAGCACCAGGCACTGCCATGGTGCGCCCCTGTCAACGCCACCCAGGAGAACCCCGATGCCGGTg                                                                                                   >  1:1794783/1‑150 (MQ=255)
                                                        aCAGATCACTTGCGCTCGCTGGCAGGATGCGTTTTCATCTGCACAAGTAGGGTGTACACACTATTATTTTGCACCGAGATGCAGCACCATGGCAGTGCCTGGTgc                                                                                                                                            >  1:851139/1‑105 (MQ=255)
                                                        aCAGATCACTTGCGCTCGCTGGCAGGATGCGTTTTCATCTGCACAAGTAGGGTGTACACACTATTATTTTGCACCGAGATGCAGCACCATGGCAGTGCCTGGTgc                                                                                                                                            <  2:851139/105‑1 (MQ=255)
                                                         cAGATCACTTGCGCTCGCTGGCAGGATGCGTTTTCATCTGCACAAGTAGGGTGTACACACTATTATTTTGCACCGAGATGCAGCACCAGGCACTGCCATGGTGCGCCCCTGTCAACGCCACCCAGGAGAACCCCGATGCCGGTGAGCAAt                                                                                             <  1:1303191/150‑1 (MQ=255)
                                                          aGATCACTTGCGCTCGCTGGCAGGATGCGTTTTCATCTGCACAAGTAGGGTGTACACACTATTATTTTGCACCGAGATGCAGCACCAGGCACTGCCATGGTGCGCCCCTGTCAACGCCACCCAGGAGAACCCCGATGCCGGTGAGCAATg                                                                                            <  2:1920413/150‑1 (MQ=255)
                                                               aCTTGCGCTCGCTGGCAGGATGCGTTTTCATCTGCACAAGTAGGGTGTACACACTATTATTTTGCACCGAGATGCAGCACCAGGCACTGCCATGGTGCGCCCCTGTCAACGCCACCCAGGAGAACCCCGATGCCGGTGAGCAATGCACaa                                                                                       <  1:920558/150‑1 (MQ=255)
                                                                cTTGCGCTCGCTGGCAGGATGCGTTTTCATCTGCACAAGTAGGGTGTACACACTATTATTTTGCACCGAGATGCAGCACCAGGCACTGCCATGGTGCGCCCCTGTCAACGCCACCCAGGAGAACCCCGATGCCGGTGa                                                                                                  <  2:1618527/138‑1 (MQ=255)
                                                                cTTGCGCTCGCTGGCAGGATGCGTTTTCATCTGCACAAGTAGGGTGTACACACTATTATTTTGCACCGAGATGCAGCACCAGGCACTGCCATGGTGCGCCCCTGTCAACGCCACCCAGGAGAACCCCGATGCCGGTGa                                                                                                  >  1:1618527/1‑138 (MQ=255)
                                                                 ttGCGCTCGCTGGCAGGATGCGTTTTCATCTGCACAAGTAGGGTGTACACACTATTATTTTGCACCGAGATGCAGCACCAGGCACTGCCATGGTGCGCCCCTGTCAACGCCACCCAGGAGAACCCCGATGCCGGTGAGCAATGCACAACt                                                                                     >  1:1941814/1‑150 (MQ=255)
                                                                 ttGCGCTCGCTGGCAGGATGCGTTTTCATCTGCACAAGTAGGGTGTACACACTATTATTTTGCACCGAGATGCAGCACCAGGCACTGCCATGGTGCGCCCCTGTCAACGCCACCCAGGAGAACCCCGATGCCGGTGAGCAATGCACAACt                                                                                     >  1:504190/1‑150 (MQ=255)
                                                                 ttGCGCTCGCTGGCAGGATGCGTTTTCATCTGCACAAGTAGGGTGTACACACTATTATTTTGCACCGAGATGCAGCACCAGGCACTGCCATGGTGCGCCCCTGTCAACGCCACCCAGGAGAACCCCGATGCCGGTGAGCAATGCACAACt                                                                                     >  1:2155773/1‑150 (MQ=255)
                                                                       tcgctGTCAGGATGCGTTTTCATCTGCACAAGTAGGGTGTACACACTATTATTTTGCACCGATATGCAGCACCAGGCACTGCCATGGTGCGCCCCTGTCAACGCCACCCAGGAGAACACCGATGCCGGTGAGCAATGCACAACTG‑‑Accgt                                                                             >  2:1189167/1‑148 (MQ=255)
                                                                         gctGGCAGGATGCGTTTTCATCTGCACAAGTAGGGTGTACACACTATTATTTTGCACCGAGATGCAGCACCATGGCAGTGCCTGGTgctgc                                                                                                                                         >  2:27996/1‑88 (MQ=38)
                                                                         gctGGCAGGATGCGTTTTCATCTGCACAAGTAGGGTGTACACACTATTATTTTGCACCGAGATGCAGCACCATGGCAGTGCCTGGTgctgc                                                                                                                                         <  1:27996/91‑4 (MQ=38)
                                                                            ggCAGGATGCGTTTTCATCTGCACAAGTAGGGTGTACACACTATTATTTTGCACCGAGATGCAGCACCAGGCACTGCCATGGTGCGCCCCTGTCAACGCCACCCAGGAGAACCCCGATGCCGGTGAGCAATGCACAACTG‑‑ACCCAGATGt                                                                        >  1:1099254/1‑150 (MQ=255)
                                                                              cAGGATGCGTTTTCATCTGCACAAGTAGGGTGTACACACTATTATTTTGCACCGAGATGCAGCACCAGGCACTGCCATGGTGCGCCCCTGTCAACGCCACCCAGGAGAACCCCGATGCCGGTGAGCAATGCACAACTG‑‑ACCCAGATGTTc                                                                      <  1:2014197/150‑1 (MQ=255)
                                                                                 gATGCGTTTTCATCTGCACAAGTAGGGTGTACACACTATTATTTTGCACCGAGATGCAGCACCAGGCACTGCCATGGTGCGCCCCTGTCAACGCCACCCAGGAGAACCCCGATGCCGGTGAGCAATGCACAACTG‑‑Accc                                                                              >  1:2146170/1‑139 (MQ=255)
                                                                                 gATGCGTTTTCATCTGCACAAGTAGGGTGTACACACTATTATTTTGCACCGAGATGCAGCACCAGGCACTGCCATGGTGCGCCCCTGTCAACGCCACCCAGGAGAACCCCGATGCCGGTGAGCAATGCACAACTG‑‑Accc                                                                              <  2:2146170/139‑1 (MQ=255)
                                                                                  aTGCGTTTTCATCTGCACAAGTAGGGTGTACACACTATTATTTTGCACCGAGATGCAGCACCAGGCACTGCCATGGTGCGCCCCTGTCAACGCCACCCAGGAGAACCCCGAt                                                                                                           <  2:1678215/112‑1 (MQ=255)
                                                                                  aTGCGTTTTCATCTGCACAAGTAGGGTGTACACACTATTATTTTGCACCGAGATGCAGCACCAGGCACTGCCATGGTGCGCCCCTGTCAACGCCACCCAGGAGAACCCCGAt                                                                                                           >  1:1678215/1‑112 (MQ=255)
                                                                                    gCGTTTTCATCTGCACAAGTAGGGTGTACACACTATTATTTTGCACCGAGATGCAGCACCAGGCACTGCCATGGTGCGCCCCTGTCAACGCCACCCAGGAGAACCCCGATGCCGGTGAGCAATGCACAACTG‑‑ACCCAGATGTTCGAGCAc                                                                >  1:1785660/1‑150 (MQ=255)
                                                                                     cGTTTTCATCTGCACAAGTAGGGTGTACACACTATTATTTTGCACCGAGATGCAGCACCAGGCACTGCCATGGTGCGCCCCTGTCAACGCCACCCAGGAGAACCCCGATGCCGGTGAGCAATGCACAACTG‑‑ACCCAGATGTTCGAGCACg                                                               <  1:1189167/150‑1 (MQ=255)
                                                                                            aTCTGCACAAGTAGGGTGTACACACTATTATTTTGCACCGAGATGCAGCACCAGGCACTGCCATGGTGCGCCCCTGTCAACGCCACCCAGGAGAACCCCGATGCCGGTGAGCAATGCACAACTG‑‑ACCCAGATGTTCGAGCACGTACTGaa                                                        <  2:1785660/150‑1 (MQ=255)
                                                                                            aTCTGCACAAGTAGGGTGTACACACTATTATTTTGCACCGAGATGCAGCACCAGGCACTGCCATGGTGCGCCCCTGTCAACGCCACCCAGGAGAACCCCGATGCCGGTGAGCAATGCACAACTG‑‑ACCCAGATGTTCGAGCACGTACTGaa                                                        >  2:730038/1‑150 (MQ=255)
                                                                                            aTCTGCACAAGTAGGGTGTACACACTATTATTTTGCACCGAGATGCAGCACCAGGCACTGCCATGGTGCGCCCCTGTCAACGCCACCCAGGAGAACCCCGATGCCGGTGAGCAATGCACAACTG‑‑ACCCAGATGTTCGAGCACGTACTGaa                                                        <  1:716115/150‑1 (MQ=255)
                                                                                             tCTGCACAAGTAGGGTGTACACACTATTATTTTGCACCGAGATGCAGCACCAGGCACTGCCATGGt                                                                                                                                              <  1:774611/66‑1 (MQ=255)
                                                                                             tCTGCACAAGTAGGGTGTACACACTATTATTTTGCACCGAGATGCAGCACCAGGCACTGCCATGGt                                                                                                                                              >  2:774611/1‑66 (MQ=255)
                                                                                             tCTGCACAAGTAGGGTGTACACACTATTATTTTGCACCGAGATGCAGCACCAGGCACTGCCATGGTGCGCCCCTGTCAACGCCACCCAGGAGAACCCCGATGCCGGTGAGCAATGCACAACTG‑‑ACCCAGATGTTCGAGCACGTACTGAAg                                                       <  2:2155773/150‑1 (MQ=255)
                                                                                             tCTGCACAAGTAGGGTGTACACACTATTATTTTGCACCGAGATGCAGCACCAGGCACTGCCATGGTGCGCCCCTGTCAACGCCACCCAGGAGAACCCCGATGCCGGTGAGCAATGCACAACTG‑‑ACCCAGATGTTCGAGCACGTACTGAAg                                                       >  1:2155658/1‑150 (MQ=255)
                                                                                             ccttcacaAGTAGGGTGTACACACTATTATTTTGCACCGAGATGCAGCACCAGGCACTGCCATGGTGCGCCCCTGTCAACGCCACCCAGGAGAACCCCGATGCCGGTGAGCAATGCAAAACTG‑‑ACCCAGATGTTCGAGCACGTACTGAAg                                                       >  1:1891849/5‑150 (MQ=255)
                                                                                               tGCACAAGTAGGGGGTACACACTGTTATTTTGCACCGAGAGGCAGCACCAGGCACTGCCATGGTGCGACCCTGTCAACGCAACCCAGGAGAACCCCGATGCCGGTGAGCAATGCACAACTG‑‑ACCCAGATGTTCGAGCCCGTACTCAAGcc                                                     >  2:486892/1‑149 (MQ=255)
                                                                                                gCACAAGTAGGGTGTACACACTATTATTTTGCACCGAGATGCAGCACCAGGCACTGCCATGGTGCGCCCCTGTCAACGCCACCCAGGAGAACCCCGATGCCGGTGAGCAATGCACAACTG‑‑ACCCAGATGTTc                                                                      >  2:1216989/1‑132 (MQ=255)
                                                                                                gCACAAGTAGGGTGTACACACTATTATTTTGCACCGAGATGCAGCACCAGGCACTGCCATGGTGCGCCCCTGTCAACGCCACCCAGGAGAACCCCGATGCCGGTGAGCAATGCACAACTG‑‑ACCCAGATGTTc                                                                      <  1:1216989/132‑1 (MQ=255)
                                                                                                    aaGTAGGGTGTACACACTATTATTTTGCACCGAGATGCAGCACCAGGCACTGCAATGGTTCGCCCCTGTCAACGCCACACAAGCGAACCCCGATGCCGCTGAACAATACACAACTG‑‑ACCCAGAAGTTCGAGCACGTACTGAAGCTGtctc                                                >  2:1027196/1‑150 (MQ=255)
                                                                                                      gTAGGGTGTACACACTATTATTTTGCACCGAGATGCAGCACCAGGCACTGCCATGGTGCGCCCCTGTCAACGCCACCCAGGAGAACCCCGATGCCGGTGAGCAATGCACAACTG‑‑ACCCAGATGTTCGAGCACGTACTGAAGCTGTCTCGa                                              >  2:2065233/1‑150 (MQ=255)
                                                                                                       tAGGGTGTACACACTATTATTTTGCACCGAGATGCAGCACCAGGCACTGCCATGGTGCGCCCCTGTCAACGCCACCCAGGAGAACCCCGATGCCGGTGAGCAATGCACAACTG‑‑ACCCAGATGTTCGAGCACGTACTGAAGCTGTCTCGAg                                             <  2:1099254/150‑1 (MQ=255)
                                                                                                       tAGGGTGTACACACTATTATTTTGCACCGAGATGCAGCACCAGGCACTGCCATGGTGCGCCCCTGTCAACGCCACCCAGGAGAACCCCGATGCCGGTGAGCAATGCACAACTG‑‑ACCCAGATGTTCGAGCACGTACTGAAGCTGTCTCGAg                                             <  1:730038/150‑1 (MQ=255)
                                                                                                        aGGGTGTACACACTATTATTTTGCACCGAGATGCAGCACCAGGCACTGCCATGGTGCGCCCCTGTCAACGCCACCCAGGAGAACCCCGATGCCGGTGAGCAATGCACAACTG‑‑ACCCAGATGTTCGAGCACGTACTGAAGCTGTCTCGAGt                                            <  2:1891849/150‑1 (MQ=255)
                                                                                                         gggTGTACACACTATTATTTTGCACCGAGATGCAGCACCAGGCACTGCCATGGTGCGCCCCTGTCAACGCCACCCAGGAGAACCCCGATGCCGGTGAGCAATGCACAACTG‑‑ACCCAGATGTTCGAGCACGTACTGAAGCTGTCTCGAGTg                                           >  2:291677/1‑150 (MQ=255)
                                                                                                          ggTGTACACACTATTATTTTGCACCGAGATGCAGCACCAGGCACTGCCATGGTGCGCCCCTGTCAACGCCACCCAGGAGAACCCCGATGCCGGTGAGCAATGCACAACTG‑‑ACCCAGATGTTCGAGCACGTACTGAAGCTGTCTCGAGTgg                                          >  2:1062012/1‑150 (MQ=255)
                                                                                                           gtgtACACACTATTATTTTGCACCGAGATGCAGCACCAGGCACTGCCATGGTGCGCCCCTGTCAACGCCACCCAGGAGAACCCCGATGCCGGTGAGCAATGCACAACTG‑‑ACCCAGATGTTCGAGCACGTACTGAAGCTGTCTCGAGTGGa                                         <  1:2065233/150‑1 (MQ=255)
                                                                                                           gtgtACACACTATTATTTTGCACCGAGATGCAGCACCAGGCACTGCCATGGTGCGCCCCTGTCAACGCCACCCAGGAGAACCCCGATGCCGGTGAGCAATGCACAACTG‑‑ACCCAGATGTTCGAGCACGTACTGAAGCTGTCTCGAGTGGa                                         >  2:581682/1‑150 (MQ=255)
                                                                                                                   acTATTATTTTGCACCGAGATGCAGCACCAGGCACTGCCATGGTTCGCCCCTGTCAACGCCAACCAGGAGAACCCCGCTGACCGTGTGCAAGGCACAACGG‑‑AGACAGATCTTCGAGCACGTACTGAAGCTGTCTCGCGTGGACGAACCGc                                 >  2:642024/1‑150 (MQ=255)
                                                                                                                     tattatTTTGCACCGAGATGCAGCACCAGGCACTGCCATGGTGCTCCCCTGTCAACGCACACCAGGAGAAAACCAAGATCGGAAGAGCA‑‑‑CACGTCTGAAATCCAGATGTTCGAGCACGTACTGAAGCTGTCTAGAATAGACGAAACGCag                               >  1:551063/1‑150 (MQ=255)
                                                                                                                     tattatTTTGCACCGAGATGCAGCACCAGGCACTGCCATGGTGCGCCCCTGTCAACGCCACCCAGGAGAACCCCGATGCCGGTGAGCAACGCACAACTG‑‑ACCCAGATGTTCGAGCACGTACTGAAGCTGTCTCGAGTGGACGAAACGCag                               >  2:1277542/1‑150 (MQ=255)
                                                                                                                      attattTTGCACCGAGATGCAGCACCAGGCACTGCCATGGTGCGCCCCTGTCAACGCCACCCAGGAGAACCCCGATGCCGGTGAGCAATGCACAACTG‑‑ACCCAGATGTTCGAGCACGTACTGAAGCTGTCTCGAGTGGACGAAACGCaga                              >  2:1108640/1‑150 (MQ=255)
                                                                                                                       ttattTTGCACCGAGATGCAGCACCAGGCACTGCCATGGTGCGCCCCTGTCAACGCCACCCAGGAGAACCCCGATGCCGGTGAGCAATGca                                                                                          >  2:558510/1‑91 (MQ=255)
                                                                                                                       ttattTTGCACCGAGATGCAGCACCAGGCACTGCCATGGTGCGCCCCGGTCAACGCCACCCAGGAGAACCCCGATGCCGGTGAGCAATGca                                                                                          <  1:558510/91‑1 (MQ=255)
                                                                                                                           tttGCACCGAGATGCAGCACCAGGCACTGCCATGGTGCGCCCCTGTCAACGCCACCCAGGAGAACCCCGATGCCGGTGAGCAATGCACAACTG‑‑ACCCAGATGTTCGAGCACGTACTGAAGCTGTCTCGAGTGGACGAAACGCAGAGCGTg                         <  2:2034483/150‑1 (MQ=255)
                                                                                                                             tGCACCGAGATGCAGCACCAGGCACTGCCATGGTGCGCCCCTGTCAACGCCACCCAGGAGAACCCCGATGCCGGTGAGCAATGCACAACTG‑‑ACCCAGATGTTCGAGCACGTACTGAAGCTGTCTCGAGTGGACGAAACGCAGAGCGTGGc                       <  2:504190/150‑1 (MQ=255)
                                                                                                                               cACCGAGATGCAGCACCAGGCACTGCCATGGTGCGCCCCTGTCAACGCCACCCAGGAGAACCCCGATGCCGGTGAGCAATGCACAACTG‑‑ACCCAGATGTTCGAGCACGTACTGAAGCTGTCTCGAGTGGACGAAACGCAGAGCGTGGccc                     >  2:1521556/1‑149 (MQ=255)
                                                                                                                                aCCGAGATGCAGCACCAGGCACTGCCATGGTGCGCCCCTGTCAACGCCACCCAGGAGAACCCCGATGCCGGTGAGCAATGCACAACTG‑‑ACCCAGATGTTCGAGCACGTACTGAAGCTGTCTCGAGTGGACGAAACGCAGAGCGTGGCCGt                    <  1:581682/150‑1 (MQ=255)
                                                                                                                                  ttataTGCAGCACCAGGCACTGCCATGGTGCGCCCCTGTCAACGCCACCCAGGAGAACCCCGATGCCGGTGAGCAATGCACAACTG‑‑ACCCAGATGTTCGAGCACGTACTGAAGCTGTCTCGAGTGGACGAAACGCAGAGCGTGGCCGTGc                  <  2:551063/146‑1 (MQ=255)
                                                                                                                                  cggggtGCAGCACTACGCGCTGGCATGGTGCGCCCCTGTCAACGCAACCCAGGAGAACCCCGAGGCCGGGGAGCAATGCACAACTG‑‑ACCCAGATGTTCGAGCACGTACTGAAGCTGTCTCGAGTGGACGAAACGCAGAGCGTGGCCGTGc                  >  2:1296928/6‑150 (MQ=255)
                                                                                                                                        gcagcaCCAGGCACTGCCATGGTGCGCCCCTGTCAACGCCACCCAGGAGAACCCCGATGCCGGTGAGCAATGCACAACTG‑‑ACCCAGATGTTCGAGCACGTACTGAAGCTGTCTCGa                                              <  1:1105775/116‑1 (MQ=255)
                                                                                                                                        gcagcaCCAGGCACTGCCATGGTGCGCCCCTGTCAACGCCACCCAGGAGAACCCCGATGCCGGTGAGCAATGCACAACTG‑‑ACCCAGATGTTCGAGCACGTACTGAAGCTGTCTCGa                                              >  2:1105775/1‑116 (MQ=255)
                                                                                                                                         cagcaCCAGGCACTGCCATGGTGCGCCCCTGTCAACGCCACCCAGGAGAACCCCGATGCCGGTGAGCAATGCACAACTG‑‑ACCCAGATGTTCGAGCACGTACTGAAGCTGTCTCGAGTGGACGAAACGCAGAGCGTGGCCGTGCTCAagag           <  1:291677/150‑1 (MQ=255)
                                                                                                                                          agcaCCAGGCACTGCCATGGTGCGCCCCTGTCAACGCCACCCAGGAGAACCCCGATGCCGGTGAGCAATGCACAACTG‑‑ACCCAGATGTTCGAGCACGTACTGAAGCTGTCTCGAGTGGACGAAACGCAGAGCGTGGCCGTGCTCAAGAGc          >  1:1703710/1‑150 (MQ=255)
                                                                                                                                           gcaCCAGGCACTGCCATGGTGCGCCCCTGTCAACGCCACCCAGGAGAACCCCGATGCCGGTGAGCAATGCACAACTG‑‑ACCCAGATGTTCGAGCACGTACTGAAGCTGTCTCGAGTGGACGAAACGCAGAGCGTGGCCGTGCTCAAGAGcc         >  1:574244/1‑150 (MQ=255)
                                                                                                                                            caCCAGGCACTGCCATGGTGCGCCCCTGTCAACGCCACCCAGGAGAACCCCGATGCCGGTGAGCAATGCACAACTG‑‑ACCCAGATGTTCGAGCACGTACTGAAGCTGTCTCGAGTGGACGAAACGCAGAGAGTGGCCGTGCTCAAGAGCCa        >  1:1929660/1‑150 (MQ=255)
                                                                                                                                             cccAGGCACTGCCATGGTGCGCCCCTGTCACCGCCACCCAGGAGAACCCCGATGCCGGTGAGCAATGCACAACTG‑‑ACCCAGATGTTCGAGCACGTACTGAAGCTGTCTCGAGTGGACGAAACGCAGAGCGTGGCCGTGCTCAAGAGCCac       <  2:1638057/149‑1 (MQ=255)
                                                                                                                                             aCCAGGCACTGCCATGGTGCGCCCCTGTCAACGCCACCCAGGAGAACCCCGATGCCGGTGAGCAATGCACAACTG‑‑ACCCAGATGTTCGAGCACGTACTGAAGCTGTCTCGAGTGGACGAAACGCAGAGCGTTGCCGTGCTCAAGAGCCac       <  1:1062012/150‑1 (MQ=255)
                                                                                                                                             aCCAGGCACTGCCATGGTGCGCCCCTGTCAACGCCACCCAGGAGAACCCCGATGCCGGTGAGCAATGCACAACTG‑‑ACCCAGATGTTCGAGCACGTACTGAAGCTGTCTCGAGTGGACGAAACGCAGAGCGTGGCCGTGCTCAAGAGCCac       >  1:1638057/1‑150 (MQ=255)
                                                                                                                                               cAGGCACTGCCATGGTGCGCCCCTGTCAACGCCACCCAGGAGAACCCCGATGCCGGTGAGCAATGCACAACTG‑‑ACCCAGATGTTCGAGCACGTACTGAAGCTGt                                                   >  2:335025/1‑104 (MQ=255)
                                                                                                                                               cAGGCACTGCCATGGTGCGCCCCTGTCAACGCCACCCAGGAGAACCCCGATGCCGGTGAGCAATGCACAACTG‑‑ACCCAGATGTTCGAGCACGTACTGAAGCTGt                                                   <  1:335025/104‑1 (MQ=255)
                                                                                                                                                aGGCACTGCCATGGTGCGCCCCTGTCAACGCCACCCAGGAGAACCCCGATGCCGGTGAGCAATGCACAACTG‑‑ACCCAGATGTTCGAGCACGTACTGAAGCTGTCTCGAGTGGACGAAACGCAGAGCGTGGCCGTGCTCAAGAGCCactac    >  2:2128038/1‑150 (MQ=255)
                                                                                                                                                 agacACTGCCATGGTGCCCCCCTGTCAAAGCCACCCAGAAGAACTTCGATTAAACAGAATTAAGCACAACTG‑‑ACCCAGATGTTCGAGCACGTACTAACGCTGTCTCGAATTGACAAACCTCCACTCGTGCCCGTGCTAAAAAGCCaatac    >  2:329235/4‑150 (MQ=255)
                                                                                                                                                  ggCACTGCCATGGTGCGCCCCTGTCAACGCCACCCAGGAGAACCCCGATGCCGGTGAGCAATGCACAACTG‑‑ACCCGGATGTTCGAGCACGTACTGAAGCTGTCTCGAGTGGACGAAACGCAGAGCGTGGCCGTGCTCAAGAGCCactact   >  2:1027048/1‑150 (MQ=255)
                                                                                                                                                   gCACTGCCATGGTGCGCCCCTGTCACCGCTGCCCAGGAGAACCCCGATGCCGGTGAGCAATGCACAACTG‑‑ACCCAGCTGTTCGAGCACGTACTGAAGCTGTCTCGAGTGGACGAAACGCAGAGCGTGGCCGTGCTCAAGAGCCACTACTc  >  1:2245749/1‑150 (MQ=255)
                                                                                                                                                      |                                                                                                                                                     
AATCCGGTTGATTGTTACCGTAAATACTGGTCAATGTGTGGATAGGTAACGTTTGAACAGATCACTTGCGCTCGCTGGCAGGATGCGTTTTCATCTGCACAAGTAGGGTGTACACACTATTATTTTGCACCGAGATGCAGCACCAGGCACTGCCATGGTGCGCCCCTGTCAACGCCACCCAGGAGAACCCCGATGCCGGTGAGCAATGCACAACTG‑‑ACCCAGATGTTCGAGCACGTACTGAAGCTGTCTCGAGTGGACGAAACGCAGAGCGTGGCCGTGCTCAAGAGCCACTACTC  >  NC_002947/4447738‑4448033

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