Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,535,980 C→G 7.9% L333V (CTG→GTG) ‡ pncB → nicotinate phosphoribosyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,535,9800CG 7.9% 279.0 / 11.3 89L333V (CTG→GTG) ‡pncBnicotinate phosphoribosyltransferase
Reads supporting (aligned to +/- strand):  ref base C (41/41);  new base G (4/3);  total (45/44)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.53e-01

TGCGCCATGATTCGGGCGAGCCGGTAGCCTGGGCGGAGAAGGCCATTGCCCATTACCAGAAACTGGGTATCGACCCGATGACCAAGACGCTGGTGTTTTCCGACGGCCTCAACCTGACCCGCTCGCTGGAAATTTTCCGCGCCCTGCGCGGACGCATCAACGTCAGCTTTGGCATCGGCACCAACCTGACCTGCGACATCCCAGGTGTGGCCCCCATGAACATCGTGCTTAAAATGACCGACTGCAACGGCGCGCCGGTGGCCAAGATCTCTGACGAAGCCGCCAA  >  NC_002947/5535837‑5536122
                                                                                                                                               |                                                                                                                                              
tGCGCCATGATTCGGGCGAGCCGGTAGCCTGGGCGGAGAAGGCCATTGCCCATTACCAGAAACTGGGTATCGACCCGATGACCAAGACGCTGGTGTTTTCCGACGGCCTCAACCTGACCCGCTCGCTGGAAATTTTCCGCGCCCTgcgc                                                                                                                                           >  2:223645/1‑149 (MQ=255)
 gcgcCATGATTCGGGCGAGCCGGTAGCCTGGGCGGAGAAGGCCATTGCCCATTACCAGAAACTGGGTATCGACCCGATGACCAAGACGCTGGTGTTTTCCGACGGCCTCAACCTGACCCGCTCGCTGGAAATTTTCCGCGCCCTGcgcg                                                                                                                                          >  2:141315/1‑149 (MQ=255)
   gcCATGATTCGGGCGAGCCGGTAGCCTGGGCGGAGAAGGCCATTGCCCATTACCAGAAACTGGGTATCGACCCGATGACCAAGACGCTGGTGTTTTCCGACGGCCTCAACCTGACCCGCTCGCTGGAAATTTTCCGCGCCCTGCGCGGa                                                                                                                                        <  8:277661/149‑1 (MQ=255)
   gcCATGATTCGGGCGAGCCGGTAGCCTGGGCGGAGAAGGCCATTGCCCATTACCAGAAACTGGGTATCGACCCGATGACCAAGACGCTGGTGTTTTCCGACGGCCTCAACCTGACCCGCTCGCTGGAAATTTTCCGCGCCCTGCGCGGa                                                                                                                                        <  4:131563/149‑1 (MQ=255)
   gcCATGATTCGGGCGAGCCGGTAGCCTGGGCGGAGAAGGCCATTGCCCATTACCAGAAACTGGGTATCGACCCGATGACCAAGACGCTGGTGTTTTCCGACGGCCTCAACCTGACCCGCTCGCTGGAAATTTTCCGCGCAGGGCGCGGa                                                                                                                                        <  7:405775/149‑1 (MQ=255)
    cctactttCGGGCGTGCCGGTAGCATCGTCGGTGAAGCCCATTTCCCATTACCAGAAACTGGGTATCTACCCGATGTCCAAGACGCTGGTGTTTTCCGACGGCCTCAACCTGACCCGCTCGCTGGAAATTTTCCGCGCCCTGCGCGGAc                                                                                                                                       <  2:390378/143‑1 (MQ=255)
     cATGATTCGGGCGAGCCGGTAGCCTGGGCGGAGAAGGCCATTGCCCATTACCAGAAACTGGGTATCGACCCGATGACCAAGACGCTGGTGTTTTCCGACGGCCTCAACCTGACCCGCTCGCTGGAAATTTTCCGCGCAGGGCGCGGaaa                                                                                                                                      <  6:412260/149‑3 (MQ=255)
      aTGATTCGGGCGAGCCGGTAGCCTGGGCGGAGAAGGCCATTGCCCATTACCAGAAACTGGGTATCGACCCGATGACCAAGACGCTGGTGTTTTCCGACGGCCTCAACCTGACCCGCTCGCTGGAAATTTTCCGCGCCCTGCGCGGACGc                                                                                                                                     <  1:185910/149‑1 (MQ=255)
         aTTCGGGCGAGCCGGTAGCCTGGGCGGAGAAGGCCATTGCCCATTACCAGAAACTGGGTATCGACCCGATGACCAAGACGCTGGTGTTTTCCGACGGCCTCAACCTGACCCGCTCGCTGGAAATTTTCCGCGCCCTGCGCGGACGCATc                                                                                                                                  <  7:409127/149‑1 (MQ=255)
             gggCGAGCCGGTAGCCTGGGCGGAGAAGGCCATTGCCCATTACCAGAAACTGGGTATCGACCCGATGACCAAGACGCTGGTGTTTTCCGACGGCCTCAACCTGACCCGCTCGCTGGAAATTTTCCGCGCCCTGCGCGGACGCATCAAc                                                                                                                               <  5:336459/148‑1 (MQ=255)
               gCGAGCCGGTAGCCTGGGCGGAGAAGGCCATTGCCCATTACCAGAAACTGTGTATCGACCCGATGACCAAGACGCTGGTGTTTTCCGACGGCCTCAACCTGACCCGCTCGCTGGAAATTTTCCGCGCCCTGCGCGGACGCATCAACGt                                                                                                                             >  1:326382/1‑148 (MQ=255)
               gCGAGCCGGCAGCCTGGGCGGAGAAGGCCATTGCCCATTACCAGAAACTGGGTATCGACCCGATGACCAAGACGCTGGTGATTTCCGACGGCCTCAACCTGACCCGCTCGCTGGAAATTTTCCGCGCCCTGCGCGGACGCATCAACGt                                                                                                                             >  1:186766/1‑148 (MQ=255)
                 gAGCCGGTAGCCTGGGCGGAGAAGGCCATTGCCCATTACCAGAAACTGGGTATCGACCCGATGACCAAGACGCTGGTGTTTTCCGACGGCCTCAACCTGACCCGCTCGCTGGAAATTTTCCGCGCCCTGCGCGGACGCATCAACGTCAg                                                                                                                          <  3:13773/149‑1 (MQ=255)
                     cGGTAGCCTGGGCGGAGAAGGCCATTGCCCATTACCAGAAACTGGGTATCGACCCGATGACCAAGACGCTGGTGTTTTCCGACGGCCTCAACCTGACCCGCTCGCTGGAAATTTTCCGCGCCCTGCGCGGACGCATCAACGTCAGCttt                                                                                                                      <  5:142094/149‑1 (MQ=255)
                      ggTAGCCTGGGCGGAGAAGGCCATTGCCCATTACCAGAAACTGGGTATCGACCCGATGACCAAGACGCTGGTGTTTTCCGACGGCCTCAACCTGACCCGCTCGCTGGAAATTTTCCGCGCCCTGCGCGGACGCATCAACGTCAGCTTTg                                                                                                                     <  3:399724/149‑1 (MQ=255)
                       gTAGCCTGGGCGGAGAAGGCCATTGCCCATTACCAGAAACTGGGTATCGACCCGATGACCAAGACGCTGGTGTTTTCCGACGGCCTCAACCTGACCCGCTCGCTGGAAATTTTCCGCGCCCTGCGCGGACGCATCAACGTCAGCTTTgg                                                                                                                    >  6:404280/1‑149 (MQ=255)
                            cTGGGCGGAGAAGGCCATTGCCCATTACCAGAAACTGGGTATCGACCCGATGACCAAGACGCTGGTGTTTTCCGACGGCCTCAACCTGACCCGCTCGCTGGAAATTTTCCGCGCCCTGCGCGGACGCATCAACGTCAGCTTTGGCATCg                                                                                                               >  5:394012/1‑149 (MQ=255)
                            cTGGGCGGAGAAGGCCATTGCCCATTACCAGAAACTGGGTATCGACCCGATGACCAAGACGCTGGTGTTTTCCGACGGCCTCAACCTGACCCGCTCGCTGGAAATTTTCCGCGCCCTGCGCGGACGCATCAACGTCAGCTTTGGCATCg                                                                                                               >  7:406029/1‑149 (MQ=255)
                              gggCGGAGAAGGCCATTGCCCATTACCAGAAACTGGGTATCGACCCGATGACCAAGACGCTGGTGTTTTCCGACGGCCTCAACCTGACCCGCTCGCTGGAAATTTTCCGCGCCCTGCGCGGACGCATCAACGTCAGCTTTGGCATCgg                                                                                                              >  5:212527/1‑148 (MQ=255)
                               ggCGGAGAAGGCCATTGCCCATTACCAGAAACTGGGTATCGACCCGATGACCAAGACGCTGGTGTTTTCCGACGGCCTCAACCTGACCCGCTCGCTGGAAATTTTCCGCGCCCTGCGCGGACGCATCAACGTCAGCTTTGGCATCGGCa                                                                                                            <  3:112408/149‑1 (MQ=255)
                                    agaAGGCCATTGCCCATTACCAGAAACTGGGTATCGACCCGATGACCAAGACGCTGGTGTTTTCCGACGGCCTCAACCTGACCCGCTCGCTGGAAATTTTCCGCGCCCTGCGCGGACGCATCAACGTCAGCTTTGGCATCGGCACCAac                                                                                                       <  2:23252/149‑1 (MQ=255)
                                        ggCCATTGCCCATTACCAGAAACTGGGTATCGACCCGATGACCAAGACGCTGGTGTTTTCCGACGGCCTCAACCTGACCCGCTCGCTGGAAATTTTCCGCGCCCTGCGCGGACGCATCAACGTCAGCTTTGGCATCGGCACCAacctga                                                                                                   <  1:223645/149‑1 (MQ=255)
                                         gCCATTGCCCATTACCAGAAACTGGGTATCGACCCGATGACCAAGACGCTGGTGTTTTCCGACGGCCTCAACCTGACCCGCTCGCTGGAAATTTTCCGCGCCCTGCGCGGACGCATCAACGTCAGCTTTGGCATCGGCACCAacctgac                                                                                                  <  2:277045/149‑1 (MQ=255)
                                         gCCATTGCCCATTACCAGAAACTGCGTATCGACCCGATGACCAAGACGCTGGTGTTTTCCGACGGCCTCCACCTGACCCGCTCGCTTGAAATATTACGCGCCCTGCGCGGACGCATCAACGTCACATATGGCATCGGCACCAacctgac                                                                                                  >  2:19294/1‑149 (MQ=255)
                                             ttGCCCATTACCAGAAACTGTGTATCGACCCGATGACCAAGACGCTGGTGTTTTCCGACGGCCTCAACCTGACCCGCTCGCTGGAAATTTTCCGCGCCCTGCGCGGACGCATCAACGTCAGCTTTGGCATCGGCACCAACCTGACCTGc                                                                                              <  5:134309/149‑1 (MQ=255)
                                                cccATTACCAGAAACTGGGTATCGACCCGATGACCAAGACGCTGGTGTTTTCCGACGGCCTCAACCTGACCCGCTCGCTGGAAATTTTCCGCGCCCTGCGCGGACGCATCAACGTCAGCTTTGGCATCGGCACCAACCTGACCTGCGa                                                                                            <  6:212527/148‑1 (MQ=255)
                                                cccATTACCAGAAACTGGGTATCGACCCGATGACCAAGACGCTGGTGTTTTCCGACGGCCTCAACCTGACCCGCTCGCTGGAAATTTTCCGCGCCCTGCGCGGACGCATCAACGTCAGCTTTGGCATCGGCACCAACCTGACCTGCGa                                                                                            <  7:136687/148‑1 (MQ=255)
                                                cccATTACCAGAAACTGGGTATCGACCCGATGACCAAGACGCTGGTGTTTTCCGACGGCCTCAACCTGACCCGCTCGCTGGAAATTTTCCGCGCCCTGCGCGGACGCATCAACGTCAGCTTTGGCATCGGCACCAACCTGACCTGCGa                                                                                            <  2:33609/148‑1 (MQ=255)
                                                cccATTACCAGAAACTGGGTATCGACCCGATGACCAAGACGCTGGTGTTTTCCGACGGCCTCAACCTGACCCGCTCGCTGGAAATTTTCCGCGCCCTGCGCGGACGCATCAACGTCAGCTTTGGCATCGGCACCAACCTGACCTGCGa                                                                                            <  7:46176/148‑1 (MQ=255)
                                                  cATTACCAGAAACTGGGTATCGACCCGATGACCAAGACGCTGGTGTTTTCCGACGGCCTCAACCTGACCCGCTCGCTGGAAATTTTCCGCGCCCTGCGCGGACGCATCAACGTCAGCTTTGGCATCGGCACCAACCTGACCTGCGACAt                                                                                         >  6:149982/1‑149 (MQ=255)
                                                       ccAGAAACTGGGTATCGACCCGATGACCAAGACGCTGGTGTTTTCCGACGGCCTCAACCTGACCCGCTCGCTGGAAATTTTCCGCGCCCTGCGCGGACGCATCAACGTCAGCTTTGGCATCGGCACCAACCTGACCTGCGACATCCCAg                                                                                    >  1:170797/1‑149 (MQ=255)
                                                       ccAGAAACTGGGTATCGACCCGATGACCAAGACGCTGGTGTTTTCCGACGGCCTCAACCTGACCCGCTCGCTGGAAATTTTCCGCGCCCTGCGCGGACGCATCAACGTCAGCTTTGGCATCGGCACCAACCTGACCTGCGACATCCCAg                                                                                    <  2:326382/149‑1 (MQ=255)
                                                        cAGAAACTGGGTATCGACCCGATGACCAAGACTCTGGTGTTTTCAGACGGCCTCAACCTGACCAGCTCGCTGGAAATTTTCCGCGCCCTGCGCGGACGCATCAACGTCAGCTTTTGCATCGGCACCAAACTGACATGCGACATCCCAgg                                                                                   >  5:439953/1‑149 (MQ=255)
                                                        cAGAAACTGGGTATCGACCCGATGACCAAGACGCTGGTGTTTTCCGACGGCCTCAACCTGACCCGCTCGCTGGAAATTTTCCGCGCCCTGCGCGGACGCATCAACGTCAGCttt                                                                                                                      >  4:163824/1‑114 (MQ=255)
                                                        cAGAAACTGGGTATCGACCCGATGACCAAGACGCTGGTGTTTTCCGACGGCCTCAACCTGACCCGCTCGCTGGAAATTTTCCGCGCCCTGCGCGGACGCATCAACGTCAGCttt                                                                                                                      <  3:163824/114‑1 (MQ=255)
                                                        cAGAAACTGGGTATCGACCCGATGACCAAGACGCTGGTGTTTTCCGACGGCCTCAACCTGACCCGCTCGCTGGAAATTTTCCGCGCCCTGCGCGGACGCATCAACGTCAGCTTTGGCATCGGCACCAACCTGACCTGCGACATCCCAgg                                                                                   >  1:295272/1‑149 (MQ=255)
                                                        cAGAAACTGGGTATCGACCCGATGACCAAGACGCTGGTGTTTTCCGACGGCCTCAACCTGACCCGCTCGCTGGAAATTTTCCGCGCCCTGCGCGGACGCATCAACGTCAGCTTTGGCATCGGCACCAACCTGACCTGCGACATCCCAgg                                                                                   >  7:63103/1‑149 (MQ=255)
                                                         aGAAACTGGGTATCGACCCGATGACCAAGACGCTGGTGTTTTCCGACGGCCTCAACCTGACCCGCTCGCTGGAAATTTTCCGCGCCCTGCGCGGACGCATCAACGTCAGCTTTGGCATCGGCACCAACCTGACCTGCGACATCCCAGgt                                                                                  >  3:253742/1‑149 (MQ=255)
                                                          gAAACTGGGTATCGACCCGATGACCAAGACGCTGGTGTTTTCCGACGGCCTCAACCTGACCCGCTCGCTGGAAATTTTCCGCGCCCTGCGCGGACGCATCAACGTCAGCTTTGGCATCGGCACCAACCTGACCTGCGACATCCCAGgtg                                                                                 >  5:271057/1‑149 (MQ=255)
                                                          gAAACTGGGTATCGACCCGATGACCAAGACGCTGGTGTTTTCCGACGGCCTCAACATGACCCGCTCGCTGGAAATTTTCCGCGCCCTGCGCGGACGCATCAACGTCAGCTTTGGCATCGGCACCAACCTGACCTGCGACATCCCAGgtg                                                                                 >  1:277896/1‑149 (MQ=255)
                                                               tGGGTATCGACCCGATGACCAAGACGCTGGTGTTTTCCGACGGCCTCAACCTGACCCGCTCGCTGGAAATTTTCCGCGCCCTGCGCGGACGCATCAACGTCAGCTTTGGCATCGGGAGGAACCTGACCTGGGACATc                                                                                        >  2:173991/1‑137 (MQ=255)
                                                               tGGGTATCGACCCGATGACCAAGACGCTGGTGTTTTCCGACGGCCTCAACCTGACCCGCTCGCTGGAAATTTTCCGCGCCCTGCGCGGACGCATCAACGTCAGCTTTGGCATCGGCACCAACCTGACCTGCGACATc                                                                                        <  1:173991/137‑1 (MQ=255)
                                                                     tCGACCCGATGACCAAGACGCTGGTGTTTTCCGACGGCCTCAACCTGACCCGCTCGCTGGAAATTTTCCGCGCCCTGCGCGGACGCATCAACGTCAGCTTTGGCATCGGCACCAACCTGACCTGCGACATCCCAGGTGTGGCCCCCAt                                                                       >  6:145795/1‑148 (MQ=255)
                                                                     tCGACCCGATGACCAAGACGCTGGTGTTTTCCGACGGCCTCAACCTGACCCGCTCGCTGGAAATTTTCCGCGCCCTGCGCGGACGCATCAACGTCAGCTTTGGCATCGGCACCAACCTGACCTGCGACATCCCAGGTGTGGCCCCCATg                                                                      <  5:404280/149‑1 (MQ=255)
                                                                     tCGACCCGATGACCAAGACGCTGGTGATTTCCGACGGCCTCAACCTGACCCGCTCGCTGGAAATTTTCCGCGCCCTGCGCGGACGCATCAACGTCAGCTTTGGCATCGGCACCAACCTGACCTGCGACATCCCAGGTGTGGCCCCCAt                                                                       <  2:186766/148‑1 (MQ=255)
                                                                          ccGATGACCAAGACGCTGGTGTTTTCCGACGGCCTCAACCTGACCCGCTCGCTGGAAATTTTCCGCGCCCTGCGCGGACGCATCAACGTCAGCTTTGGCATCGGCACCAACCTGACCTGCGACATCCCAGGTGTGGCCCCCATGAACAt                                                                 >  6:353859/1‑149 (MQ=255)
                                                                           cGATGACCAAGTCGCTGGTGTTTTCCGACGGCCTCAACCTGTCCCGCTCGCTGGAAATTTTCCGCGCCCTGCGCGGACGCATCAACGTCAGCTTTGGCATCGGCACCAACCTGACCTGCGACATCCCAGGTGTGGCCCCCATGAACATc                                                                <  2:277896/149‑1 (MQ=255)
                                                                           cGATGACCAAGACGCTGGTGTTTTCCGACGGCCTCAACCTGACCCGCTCGCTGGAAATTTTCCGCGCCCTGCGCGGACGCATCAACGTCAGCTTTGGCATCGGCACCAACCTGACCTGCGACATCCCAGGTGTGGCCCCCATGAACATc                                                                <  2:295272/149‑1 (MQ=255)
                                                                            gATGACCAAGACGCTGGTGTTTTCCGACGGCCTCAACCTGACCCGCTCGCTGGAAATTTTCCGCGCCCTGCGCGGACGCATCAACGTCAGCTTTGGCATCGGCACCAACCTGACCTGCGACATCCCAGGTGTGGCCCCCATGAACATCg                                                               >  5:380675/1‑149 (MQ=255)
                                                                             aTGACCAAGACGCTGGTGTTTTCCGACGGCCTCAACCTGACCCGCTCGCTGGAAATTTTCCGCGCCCTGCGCGGACGCATCAACGTCAGCTTTGGCATCGGCACCAACCTGACCTGCGACATCCCAGGTGTGGCCCCCATGAACATCGt                                                              >  7:336503/1‑149 (MQ=255)
                                                                                  cAAGACGCTGGTGTTTTCCGACGGCCTCAACCTGACCCGCTCGCTGGAAATTTTCCGCGCCCTGCGCGGACGCATCAACGTCAGCTTTGGCATCGGCACCAACCTGACCTGCGACATCCCAGGTGTGGCCCCCATGAACATCGt                                                              >  2:8086/1‑144 (MQ=255)
                                                                                  cAAGACGCTGGTGTTTTCCGACGGCCTCAACCTGACCCGCTCGCTGGAAATTTTCCGCGCCCTGCGCGGACGCATCAACGTCAGCTTTGGCATCGGCACCAACCTGACCTGCGACATCCCAGGTGTGGCCCCCATGAACATCGt                                                              <  1:8086/144‑1 (MQ=255)
                                                                                   aaGACGCTGGTGTTTTCCGACGGCCTCAACCTGACCCGCTCGCTGGAAATTTTCCGCGCCCTGCGCGGACGCATCAACGTCAGCTTTGGCATCGGCACCAACCTGACCTGCGACATCCCAGGTGTGGCCCCCATGAACATCGTGCTTaa                                                        <  5:149982/149‑1 (MQ=255)
                                                                                   aaGACGCTGGTGTTTTCCGACGGCCTCAACCTGACCCGCTCGCTGGAAATTTTCCGCGCCCTGCGCGGACGCATCAACGTCAGCTTTGGCATCGGCACCAACCTGACCTGCGACATCCCAGGTGTGGCCCCCATGAACATCGTGCTTa                                                         >  8:77223/1‑148 (MQ=255)
                                                                                   aaGACGCTGGTGTTTTCCGACGGCCTCAACCTGACCCGCTCGCTGGAAATTTTCCGCGCCCAGCGCGGACGCATCAACGTCAGCTTTGGCATCGGCACCAACCTGACCTGCGACATCCCAGGTGTGGCCCCCATGAACATCGTGc                                                            >  5:194985/1‑145 (MQ=255)
                                                                                   aaGACGCTGGTGTTTTCCGACGGCCTCAACCTGACCCGCTCGCTGGAAATTTTCCGCGCCCAGCGCGGACGCATCAACGTCAGCTTTGGCATCGGCACCAACCTGACCTGCGACATCCCAGGTGTGGCCCCCATGAACATCGTGc                                                            <  6:194985/145‑1 (MQ=255)
                                                                                            gtgtATTCCGACGGCCTCAACCTGCCCCGCTCGCTGGAAATTTTCCGCGCCCTGCGCGGACGCATCAACGTCAGCGTTGGCATCGGCACCAACCTGACCTGCGACATCCCAGGTGTGGCCCCCATGAACATCGTCCTTAAAATGACCGa                                               <  6:333338/149‑1 (MQ=255)
                                                                                               ttttCCGACGGCCTCAACCTGACCCGCTCGCTGGAAATTTTCCGCGCCCTGCGCGGACGCATCAACGTCAGCTTTGGCATCGGCACCAACCTGACCTGCGACATCCCAGGTGTGGCCCCCATGAACATCGTGCtt                                                          >  3:163233/1‑135 (MQ=255)
                                                                                               ttttCCGACGGCCTCAACCTGACCCGCTCGCTGGAAATTTTCCGCGCCCTGCGCGGACGCATCAACGTCAGCTTTGGCATCGGCACCAACCTGACCTGCGACATCCCAGGTGTGGCCCCCATGAACATCGTGCtt                                                          <  4:163233/135‑1 (MQ=255)
                                                                                                  tCCGACGGCCTCAACCTGACCCGCTCGCTGGAAATTTTCCGCGCCCTGCGCGGACGCATCAACGTCAGCTTTGGCATCGTCACCAACCTGACCTGCGACATCCCAGGTGTGGCCCCCATGAACATCGTGCTTAAAATGACCGACTGCaa                                         >  6:215582/1‑149 (MQ=255)
                                                                                                    cGACGGCCTCAACCTGACCCGCTCGCTGGAAATTTTCCGCGCCCTGCGCGGACGCATCAACGTCAGCTTTGGCATCGGCACCAACCTGACCTGCGACATCCCAGGTGTGGCCCCCATGAACATCGTGCTTAAAATGACCGACTGCAACg                                       >  1:437432/1‑149 (MQ=255)
                                                                                                     gACGGCCTCAACCTGACCCGCTCGCTGGAAATTTTCCGCGCCCTGCGCGGACGCATCAACGTCAGCTTTGGCATCGGCACCAACCTGACCTGCGACATCCCAGGTGTGGCCCCCATGAACATCGTGCTTAAAATGa                                                   <  5:148040/136‑1 (MQ=255)
                                                                                                     gACGGCCTCAACCTGACCCGCTCGCTGGAAATTTTCCGCGCCCTGCGCGGACGCATCAACGTCAGCTTTGGCATCGGCACCAACCTGACCTGCGACATCCCAGGTGTGGCCCCCATGAACATCGTGCTTAAAATGa                                                   >  6:148040/1‑136 (MQ=255)
                                                                                                      aCGGCCTCAACCTGACCCGCTCGCTGGAAATTTTCCGCGCCCTGCGCGGACGCATCAACGTCAGCTTTGGCATCGGCACCAACCTGACCTGCGACATCCCAGGTGTGGCCCCCa                                                                        <  8:164937/114‑1 (MQ=255)
                                                                                                      aCGGCCTCAACCTGACCCGCTCGCTGGAAATTTTCCGCGCCCTGCGCGGACGCATCAACGTCAGCTTTGGCATCGGCACCAACCTGACCTGCGACATCCCAGGTGTGGCCCCCa                                                                        >  7:164937/1‑114 (MQ=255)
                                                                                                      aCGGCCTCAACCTGACCCGCTCGCTGGAAATTTTCCGCGCCCTGCGCGGACGCATCAACGTCAGCTTTGGCATCGGCACCAACCTGACCTGCGACATCCCAGGTGTGGCCCCCATGAACATCGTGCTTAAAATGACCGACTGCAACGgc                                     >  7:34438/1‑149 (MQ=255)
                                                                                                          ccTCAACCTGACCCGCTCGCTGGAAATTTTCCGCGCCCTGCGCGGACGCATCAACGTCAGCTTTGGCATCGGCACCAACCTGACCTGCGACATCCCAGGTGTGGCCCCCATGAACATCGTGCTTAAAATGACCGACTGCAACGgcgcgc                                 >  1:82737/1‑149 (MQ=255)
                                                                                                           cTCAACCTGACCCGCTCGCTGGAAATTTTCCGCGCCCTGCGCGGACGCATCAACGTCAGCTTTGGCATCGGCACCAACCTGACCTGCGACATCCCAGGTGTGGCCCCCATGAACATCGTGCTTAAAATGACCGACTGCAACGGCGCGcc                                >  8:323835/1‑149 (MQ=255)
                                                                                                              aaCCTGACCCGCTCGCTGGAAATTTTCCGCGCCCTGCGCGGACGCATCAACGTCAGCTTTGGCATCGGCACCAACCTGACCTGCGACATCCCAGGTGTGGCCCCCATGAACATCGTGCTTAAAATTACCGACTGCAACGGCGCGCCGGt                             >  6:71301/1‑149 (MQ=255)
                                                                                                                ccTGACCCGCTCGCTGGAAATTTTCCGCGCCCTGCGCGGACGCATCAACGTCAGCTTTGGCATCGGCACCAACCTGACCTGCGACATCCCAGGTGTGGCCCCCATGAACATCGTGCTTAAAATGACCGACTGCAACGGCGCGCCGGTgg                           >  1:405187/1‑149 (MQ=255)
                                                                                                                      ccGCTCGCTGGAAATTTTCCGCGCCCTGCGCGGACGCATCAACGTCAGCTTTGGCATCGGCACCAACCTGACCTGCGACATCCCAGGTGTGGCCCCCATGAACATCGTGCTTAAAATGACCGACTGCAACGGCGCGCCGGTGGCCAAGa                     <  7:77223/149‑1 (MQ=255)
                                                                                                                        gctcgctGGAAATTTTCCGCGCCCTGCGCGGACGCATCAACGTCAGCTTTGGCATCGGCACCAACCTGACCTGCGACATCCCAGGTGTGGCCCCCATGAACATCGTGCTTAAAATGACCGACTGCAACGGCGCGCCGGTGGCCAAGAtc                   <  5:353859/149‑1 (MQ=255)
                                                                                                                        gctcgctGGAAATTTTCCGCGCCCTGCGCGGACGCATCAACGTCAGCTTTGGCATCGGCACCAACCTGACCTGCGACATCCCAGGTGTGGCCCCCATGAACATCGTGCTTAAAATGACCGACTGCAACGGCGCGCCGGTGGCCAAGAt                    >  6:124434/1‑148 (MQ=255)
                                                                                                                           cgctGGAAATTTTCCGCGCCCTGCGCGGACGCATCAACGTCAGCTTTGGCATCGGCACCAACCTGACCTGCGACATCCCAGGTGTGGCCCCCATGAACATCGTGCTTAAAATCACCGACTGCAACGGCGCGCCGGTGGCCAAGATctct                <  5:71301/149‑1 (MQ=255)
                                                                                                                             ctGGAAATTTTCCGCGCCCTGCGCGGACGCATCAACGTCAGCTTTGGCATCGGCACCAACCTGACCTGCGACATCCCAGGTGTGGCCCCCATGAACATCGTGCTTAAAATGACCGACTGCAACGGCGCGCCGGTGGCCAAGATCTCTGa              >  1:221337/1‑149 (MQ=255)
                                                                                                                             ctGGAAATTTTCCGCGCCCTGCGCGGACGCATCAACGTCAGCTTTGGCATCGGCACCAACCTGACCTGCGACATCCCAGGTGTGGCCCCCATGAACATCGTGCTTAAAATGACCGACTGCAACGGCGCGCCGGTGGCCAAGATCTCTGa              <  6:394012/149‑1 (MQ=255)
                                                                                                                              tGGAAATTTTCCGCGCCCTGCGCGGACGCATCAACGTCAGCTTTGGCATCGGCACCAACCTGACCTGCGACATCCCAGGTGTGGCCCCCATGAACATCGTGCTTAAAATGACCGACTGCAACGGCGCGCCGGTGGCCAAGATctct                <  3:317335/146‑1 (MQ=255)
                                                                                                                              tGGAAATTTTCCGCGCCCTGCGCGGACGCATCAACGTCAGCTTTGGCATCGGCACCAACCTGACCTGCGACATCCCAGGTGTGGCCCCCATGAACATCGTGCTTAAAATGACCGACTGCAACGGCGCGCCGGTGGCCAAGATctct                >  4:317335/1‑146 (MQ=255)
                                                                                                                              tGGAAATTTTCCGCGCCCTGCGCGGACGCATCAACGTCAGCTTTGGCATCGGCACCAACCTGACCTGCGACATCCCAGGTGTGGCCCCCATGAACATCGTGCTTAAAATGACCGACTGCAACGGCGCGCCGGTGGCCAAGATCTCTGAc             >  2:399596/1‑149 (MQ=255)
                                                                                                                                gAAATTTTCCCCGCCCTGCGCGGACGCATCAACGTCAGCTTTGGCATCGGCACCAACCTGACCTGCGACATCCCAGGTGTGGCCCCCATGAACATCGTGCTTAAAATGACCGACTGCAACGGCGCGCCGGTGGCCAAGATCTCTGAc             <  8:406029/147‑1 (MQ=255)
                                                                                                                                  aaTTTTCCGCGCCCTGCGCGGACGCATCAACGTCAGCTTTGGCATCGGCACCAACCTGACCTGCGACATCCCAGGTGTGGCCCCCATGAACATCGTGCTTAAGATGACCGACTGCAACGGCGCGCCGGTGGCCGAGATCTCTGACGAAg         >  8:78924/1‑149 (MQ=255)
                                                                                                                                      gtCCGCGCAGGGCGCGGACGCATCAACGTCAGCTTTGGCATCGGCACCAACCTGCCCTGCGACATCCCAGGTGTGGCCCCCATGAACATCGTGCTTAAAATGACCGACTGCAACGGCGCGCCGGTGGCCAAGATCTCTGACGAAgccgc     >  1:122479/2‑149 (MQ=255)
                                                                                                                                      gtCCGCGCAGGGCGCGGACGCATCAACGTCAGCTTTGGCATCGGCACCAACCTGACCTGCGACATCCCAGGTGTGGCCCCCATGAACATCGTGCTTAAAATGACCGACTGCAACGGCGCGCCGGTGGCCAAGATCTCTGACGAAgccgc     >  6:160537/2‑149 (MQ=255)
                                                                                                                                      gtCCGCGCAGGGCGCGGACGCATCAACGCCAGCTTTGGCATCGGCACCAACCTGACCTGCGACATCCCAGGTGTGGCCCCCATGAACATCGTGCTTAAAATGACCGAGTGCAACGGCGCGCCGGt                             >  5:376240/2‑125 (MQ=255)
                                                                                                                                      gtCCGCGCAGGGCGCGGACGCATCAACGCCAGCTTTGGCATCGGCACCAACCTGACCTGCGACATCCCAGGTGTGGCCCCCATGAACATCGTGCTTAAAATGACCGACTGCAACGGCGCGCCGGt                             <  6:376240/124‑1 (MQ=255)
                                                                                                                                       tCCTCGCCCTGCGCGGACGCATCAACGTCAGCTTTGGCATCGTCACCAACCTGTCCTGCGACATCCCAGGTGTGGCCCCCATGAACATCGTGCTTAAAATGACCGACTTCAACGGCGCGCCGGTGGCCAAGATCTCTGACGAAgccgc     <  6:439953/148‑1 (MQ=255)
                                                                                                                                       tCCGCGCAGGGCGCGGACGCATCAACGTCAGCTTTGGCATCGGCACCAACCTGACCTGCGACATCCCAGGTGTGGCCCCCATGAACATCGTGCTTAAAATGACCGACTGCAACGGCGCGCCGGTGGCCAAGATCTCTGACGAAgccgc     >  5:281260/1‑148 (MQ=255)
                                                                                                                                         cgcgCCCTGCGCGGACGCATCAACGTCAGCTTTGGCATCGTCACCAACCTGACCTGCGACATCCCAGGTGTGGCCCCCATGAACATCGTGCTTAAAATGACCGACTGCAACGGCGCGCCGGTGGCCAAGATCTCTGACGAAGCCGCCaa  <  5:215582/149‑1 (MQ=255)
                                                                                                                                         cgcgCCCTGCGCGGACGCATCAACGTCAGCTTTGGCATCGGCACCAACCTGACCTGCGACATCCCAGGTGTGGCCCCCATGAACATCGTGCTTAAAATGACCGACTGCAACGGCGCGCCGGTGGCCAAGATCTCTGACGAAGCCGCCaa  <  6:324139/149‑1 (MQ=255)
                                                                                                                                               |                                                                                                                                              
TGCGCCATGATTCGGGCGAGCCGGTAGCCTGGGCGGAGAAGGCCATTGCCCATTACCAGAAACTGGGTATCGACCCGATGACCAAGACGCTGGTGTTTTCCGACGGCCTCAACCTGACCCGCTCGCTGGAAATTTTCCGCGCCCTGCGCGGACGCATCAACGTCAGCTTTGGCATCGGCACCAACCTGACCTGCGACATCCCAGGTGTGGCCCCCATGAACATCGTGCTTAAAATGACCGACTGCAACGGCGCGCCGGTGGCCAAGATCTCTGACGAAGCCGCCAA  >  NC_002947/5535837‑5536122

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: ‑