Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,023,295 G→T 5.1% P660P (CCG→CCT xdhD → putative hypoxanthine oxidase, molybdopterin‑binding/Fe‑S binding

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,023,2950GT 5.1% 271.4 / 3.1 77P660P (CCG→CCTxdhDputative hypoxanthine oxidase, molybdopterin‑binding/Fe‑S binding
Reads supporting (aligned to +/- strand):  ref base G (40/33);  new base T (2/2);  total (42/35)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.94e-01

CCTCCGTTCCTTCCGCCGCCAGCTGTGCACTGGAAGAGATCCTGCGTCAGGGTCGCGAGATGATCCAATGGTCGTCACCGAAACCGCAAAATGGTGACTGGCACATCGGTCGCGGCGTCGCCATTATCATGCAGAAATCGGGTATCCCGGATATCGATCAGGCTAACTGCATGATCAAACTGGAATCGGACGGCACCTTTATCGTTCATTCTGGCGGTGCGGATATTGGTACTGGTCTGGATACTGTAGTGACGAAACTGGCAGCAGAAGTGCTGCACTGCCCACCGCAGGACGTGC  >  NC_000913/3023147‑3023443
                                                                                                                                                    |                                                                                                                                                    
ccTCCGTTCCTTCCGCCGCCAGCTGTGCACTGGAAGAGATCCTGCGTCAGGGTCGCGAGATGATCCAATGGTCGTCACCGAAACCGCAAAATGGTGACTGGCACATCGGTCGCGGCGTCGCCATTATCATGCAGAAATCGGGTATCCCg                                                                                                                                                      <  2:462059/149‑1 (MQ=255)
  tCCGTTCCTTCCGCCGCCAGCTGTGCACTGGAAGAGATCCTGCGTCAGGGTCGCGAGATGATCCAATGGTCGTCACCGAAACCGCAAAATGGTGACTGGCACATCGGTCGCGGCGTCGCCATTATCATGCAGAAATCGGGTATCCCGGa                                                                                                                                                    <  4:64281/149‑1 (MQ=255)
  tCCGTTCCTTCCGCCGCCAGCTGTGCACTGGAAGAGATCCTGCGTCAGGGTCGCGAGATGATCCAATGGTCGTCACCGAAACCGCAAAATGGTGACTGGCACATCGGTCGCGGCGTCGCCATTATCATGCAGAAATCGGGTATCCCGGa                                                                                                                                                    <  1:266347/149‑1 (MQ=255)
  tCCGTTCCTTCCGCCGCCAGCTGTGCACTGGAAGAGATCCTGCGTCAGGGTCGCGAGATGATCCAATGGTCGTCACCGAAACCGCAAAATGGTGACTGGCACATCGGTCGCGGCGTCGCCATTATCATGCAGAAATCGGGTATCCCGGa                                                                                                                                                    >  2:342191/1‑149 (MQ=255)
  tCCGTTCCTTCCGCCGCCAGCTGTGCACTGGAAGAGATCCTGCGTCAGGGTCGCGAGATGATCCAATGGTCGTCACCGAAACCGCAAAATGGTGACTGGCACATCGGTCGCGGCGTCGCCATTATCATGCAGAAATCGGGTATCCCGGa                                                                                                                                                    <  1:395495/149‑1 (MQ=255)
   ccGTTCCTTCCGCCGCCAGCTGTGCACTGGAAGAGATCCTGCGTCAGGGTCGCGAGATGATCCAATGGTCGTCACCGAAACCGCAAAATGGTGACTGGCACATCGGTCGCGGCGTCGCCATTATCATGCAGAAATCGGGTATCCCGGat                                                                                                                                                   >  3:163948/1‑149 (MQ=255)
         cttccGCCGCCAGCTGTGCACTGGAAGAGATCCTGCGTCAGGGTCGCGAGATGATCCAATGGTCGTCACCGAAACCGCAAAATGGTGACTGGCACATCGGTCGCGGCGTCGCCATTATCATGCAGAAATCGGGTATCCCGGATATCGAt                                                                                                                                             >  3:74443/1‑149 (MQ=255)
          ttccGCCGCCAGCTGTGCACTGGAAGAGATCCTGCGTCAGGGTCGCGAGATGATCCAATGGTCGTCACCGAAACCGCAAAATGGTGACTGGCACATCGGTCGCGGCGTCGCCATTATCATGCAGAAATCGGGTATCCCGGATATCGATc                                                                                                                                            >  3:225432/1‑149 (MQ=255)
           tccGCCGCCAGCTGTGCACTGGAAGAGATCCTGCGTCAGGGTCGCGAGATGATCCAATGGTCGTCACCGAAACCGCAAAATGGTGACTGGCACATCGGTCGCGGCGTCGCCATTATCATGCAGAAATCGGGTATCCCGGATATCGATCa                                                                                                                                           <  2:131493/149‑1 (MQ=255)
                 gccAGCTGTGCACTGGAAGAGATCCTGCGTCAGGGTCGCGAGATGATCCAATGGTCGTCACCGAAACCGCAAAATGGTGACTGGCACATCGGTCGCGGCGTCGCCATTATCATGCAGAAATCGGGTATCCCGGATATCGATCAGGCTaa                                                                                                                                     >  3:142979/1‑149 (MQ=255)
                 gccAGCTGTGCACTGGAAGAGATCCTGCGTCAGGGTCGCGAGATGATCCAATGGTCGTCACCGAAACCGCAAAATGGTGACTGGCACATCGGTCGCGGCGTCGCCATTATCATGCAGAAATCGGGTATCCCGGATATCGATCAGGCTaa                                                                                                                                     >  2:99641/1‑149 (MQ=255)
                     gCTGTGCACTGGAAGAGATCCTGCGTCAGGGTCGCGAGATGATCCAATGGTCGTCACCGAAACCGCAAAATGGTGACTGGCACATCGGTCGCGGCGTCGCCATTATCATGCAGAAATCGGGTATCCCGGATATCGATCAGGCTAACTGc                                                                                                                                 <  4:67428/149‑1 (MQ=255)
                       tgtgCACTGGAAGAGATCCTGCGTCAGGGTCGCGAGATGATCCAATGGTCGTCACCGAAACCGCAAAATGGTGACTGGCACATCGGTCGCGGCGTCGCCATTATCATGCAGAAATCGGGTATCCCGGATATCGATCAGGCTAACTGCAt                                                                                                                               <  2:21809/149‑1 (MQ=255)
                          gCACTGGAAGAGATCCTGCGTCAGGGTCGCGAGATGATCCAATGGTCGTCACCGAAACCGCAAAATGGTGACTGGCACATCGGTCGCGGCGTCGCCATTATCATGCAGAAATCGGGTATCCCGGATATCGATCAGGCTAACTGCATGAt                                                                                                                            <  1:143763/149‑1 (MQ=255)
                            aCTGGAAGAGATCCTGCGTCAGGGTCGCGAGATGATCCAATGGTCGTCACCGAAACCGCAAAATGGTGACTGGCACATCGGTCGCGGCGTCGCCATTATCATGCAGAAATCGGGTATCCCTGTTATCGATCAGGCTAACTGCATGATCa                                                                                                                          >  2:346157/1‑149 (MQ=255)
                            aCTGGAAGAGATCCTGCGTCAGGGTCGCGAGATGATCCAATGGTCGTCACCGAAACCGCAAAATGGTGACTGGCACATCGGTCGCGGCGTCGCCATTATCATGCAGAAATCGGGTATCCCTGATATCGATCAGGCTAACTGCATGATCa                                                                                                                          >  2:170352/1‑149 (MQ=255)
                              tGGAAGAGATCCTGCGTCAGGGTCGCGAGATGATCCAATGGTCGTCACCGAAACCGCAAAATGGTGACTGGCACATCGGTCGCGGCGTCGCCATTATCATGCAGAAATCGGGTATCCCGGATATCGATCAGGCTAACTGCATGATCaaa                                                                                                                        <  3:52726/149‑1 (MQ=255)
                                gAAGAGATCCTGCGTCAGGGTCGCGAGATGATCCAATGGTCGGCACCGAAACCGCGAAATGGTGACTGGCACATCGGTCGCGGCGTCGCCATTATCATGCAGAAATCGGGTATCCCGGATATCGATCAGGCTAACTGCATGATCAAACt                                                                                                                      <  3:531185/149‑1 (MQ=255)
                                  agagATCCTGCGTCAGGGTCGCGAGATGATCCAATGGTCGTCACCGAAACCGCAAAATGGTGACTGGCACATCGGTCGCGGCGTCGCCATTATCATGCAGAAATCGGGTATCCCGGATATCGATCAGGCTAACTGCATGATCAAACTgg                                                                                                                    >  1:110366/1‑149 (MQ=255)
                                  agagATCCTGCGTCAGGGTCGCGAGATGATCCAATGGTCGTCACCGAAACCGCAAAATGGTGACTGGCACATCGGTCGCGGCGTCGCCATTATCATGCAGAAATCGGGTATCCCGGATATCGATCAGGCTAACTGCATGATCAAACTgg                                                                                                                    >  3:402588/1‑149 (MQ=255)
                                      aTCCTGCGTCAGGGTCGCGAGATGATCCAATGGTCGTCACCGAAACCGCAAAATGGTGACTGGCACATCGGTCGCGGCGTCGCCATTATCATGCAGAAATCGGGTATCCCGGATATCGATCAGGCTAACTGCATGATCAAACTGGAATc                                                                                                                >  3:111559/1‑149 (MQ=255)
                                      aTCCTGCGTCAGGGTCGCGAGATGATCCAATGGTCGTCACCGAAACCGCAAAATGGTGACTGGCACATCGGTCGCGGCGTCGCCATTATCATGCAGAAATCGGGTATCCCGGATATCGATCAGGCTAACTGCATGATCAAACTGGAATc                                                                                                                >  2:269007/1‑149 (MQ=255)
                                        ccTGCGTCAGGGTCGCGAGATGATCCAATGGTCGTCACCGAAACCGCAAAATGGTGACTGGCACATCGGTCGCGGCGTCGCCATTATCATGCAGAAATCGGGTATCCCGGATATCGATCAGGCTAACTGCATGATCAAACTGGAATCgg                                                                                                              >  4:348659/1‑149 (MQ=255)
                                         cTGCGTCAGGGTCGCGAGATGATCCAATGGTCGTCACCGAAACCGCAAAATGGTGACTGGCACATCGGTCGCGGCGTCGCCATTATCATGCAGAAATCGGGTATCCCGGATATCGATCAGGCTAACTGCATGATCAAACTGGAATCGGa                                                                                                             >  2:73820/1‑149 (MQ=255)
                                         cTGCGTCAGGGTCGCGAGATGATCCAATGGTCGTCACCGAAACCGCAAAATGGTGACTGGCACATCGGTCGCGGCGTCGCCATTATCATGCAGAAATCGGGTATCCCGGATATCGATCAGGCTAACTGCATGATCAAACTGGAATCGGa                                                                                                             >  2:92533/1‑149 (MQ=255)
                                         cTGCGTCAGGGTCGCGAGATGATCCAATGGTCGTCACCGAAACCGCAAAATGGTGACTGGCACATCGGTCGCGGCGTCGCCATTATCATGCAGAAATCGGGTATCCCGGATATCGATCAGGCTAACTGCATGATCAAACTGGAATCGGa                                                                                                             <  1:372778/149‑1 (MQ=255)
                                            cGTCAGGGTCGCGAGATGATCCAATGGTCGTCACCGAAACCGCAAAATGGTGACTGGCACATCGGTCGCGGCGTCGCCATTATCATGCAGAAATCGGGTATCCCTGATATCGATCAGGCTAACTGCATGATCAAACTGGAATCGGACgg                                                                                                          <  1:346157/149‑1 (MQ=255)
                                            cGTCAGGGTCGCGAGATGATCCAATGGTCGTCACCGAAACCGCAAAATGGTGACTGGCACATCGGTCGCGGCGTCGCCATTATCATGCAGAAATCGGGTATCCCTGATATCGATCAGGCTAACTGCATGATCAAACTGGAATCGGACgg                                                                                                          <  1:170352/149‑1 (MQ=255)
                                                    tCGCGAGATGATCCAATGGTCGTCACCGAAACCGCAAAATGGTGACTGGCACATCGGTCGCGGCGTCGCCATTATCATGCAGAAATCGGGTATCCCGGATATCGATCAGGCTAACTGCATGATCAAACTGGAATCGGACGGCACCTTTa                                                                                                  <  4:20346/149‑1 (MQ=255)
                                                       cgAGATGATCCCATGGTCGTCACCGAAACCGCAAAATGGTGACTGGCACATCGGTCGCGGCGTCGCCATTATCATGCAGAAATCGGGTATCCCGGATATCGATCAGGCTAACTGCATGATCAAACTGGAATCGGACGGCACCTTTATCg                                                                                               >  1:434409/1‑149 (MQ=255)
                                                        gagaTGATCCAATGGTCGTCACCGAAACCGCAAAATGGTGACTGGCACATCGGTCGCGGCGTCGCCATTATCATGCAGAAATCGGGTATCCCGGATATCGATCAGGCTAACTGCATGTTCAAACTGGAATCGGACGGCACCTTTATCGt                                                                                              >  4:498933/1‑149 (MQ=255)
                                                           atgatCCAATGGGCGTCACCGAAACCGCCAAATGGTGACTGGCACATCGGTCGCGGCGTCGCCATTATCACGCAGAAATCGGGTAGCCCGGATATCGATCAGGCTAACTGCATGATCAAACTGGAATCGGACGGAACCTTTATCGTTCa                                                                                           <  4:240594/149‑1 (MQ=255)
                                                              atCCAATGGTCGTCACCGAAACCGCAAAATGGTGACTGGCACATCGGTCGCGGCGTCGCCATTATCATGCAGAAATCGGGTATCCCGGATATCGATCAGGCTAACTGCATGATCAAACTGGAATCGGACGGCACCTTTATCGTTCATTc                                                                                        >  4:269594/1‑149 (MQ=255)
                                                               tCCAATGGTCGTCACCGAAACCGCAAAATGGTGACTGGCACATCGGTCGCGGCGTCGCCATTATCATGCAGAAATCGGGTATCCCGGATATCGATCAGGCTAACTGCATGATCAAACTGGAATCGGACGGCACCTTTATCGTTCATTCt                                                                                       >  2:236970/1‑149 (MQ=255)
                                                               tCCAATGGTCGTCACCGAAACCGCAAAATGGTGACTGGCACATCGGTCGCGGCGTCGCCATTATCATGCAGAAATCGGGTATCCCGGATATCGATCAGGCTAACTGCATGATCAAACTGGAATCGGACGGCACCTTTATCGTTCATTCt                                                                                       >  1:329555/1‑149 (MQ=255)
                                                                ccAATGGTCGTCACCGAAACCGCAAAATGGTGACTGGCACATCGGTCGCGGCGTCGCCATTATCATGCAGAAATCGGGTATCCCGGATATCGATCAGGCTAACTGCATGATCAAACTGGAATCGGACGGCACCTTTATCGTTCATTCTg                                                                                      <  1:412674/149‑1 (MQ=255)
                                                                     ggTCGTCACCGAAACCGCAAAATGGTGACTGGCACATCGGTCGCGGCGTCGCCATTATCATGCAGAAATCGGGTATCCCGGATATCGATCAGGCTAACTGCATGATCAAACTGGAATCGGACGGCACCTTTATCGTTCATTCTGGCGGt                                                                                 <  2:41132/149‑1 (MQ=255)
                                                                       tcgtcACCGAAACCGCAAAATGGTGACTGGCACATCGGTCGCGGCGTCGCCATTATCATGCAGAAATCGGGTATCCCGGATATCGATCAGGCTAACTGCATGATCAAACTGGAATCGGACGGCACCTTTATCGTTCATTCTGGCGGTGc                                                                               >  1:103458/1‑149 (MQ=255)
                                                                          tcACCGAAACCGCAAAATGGTGACTGGCACATCGGTCGCGGCGTCGCCATTATCATGCAGAAATCGGGTATCCCGGATATCGATCAGGCTAACTGCATGATCAAACTGGAATCGGACGGCACCTTTATCGTTCATTCTGGCGGTGCGGa                                                                            >  1:551937/1‑149 (MQ=255)
                                                                            aCCGAAACCGCAAAATGGTGACTGGCACATCGGTCGCGGCGTCGCCATTATCATGCAGAAATCGGGTATCCCGGATATCGATCAGGCTAACTGCATGATCAAACTGGAATCGGACGGCACCTTTATCGTTCATTCTGGCGGTGCGGata                                                                          >  3:201945/1‑149 (MQ=255)
                                                                            aCCGAAACCGCAAAATGGTGACTGGCACATCGGTCGCGGCGTCGCCATTATCATGCAGAAATCGGGTATCCCGGATATCGATCAGGCTAACTGCATGATCAAACTGGAATCGGACGGCACCTTTATCGTTCATTCTGGCGGTGCGGata                                                                          >  4:60324/1‑149 (MQ=255)
                                                                              cGAAACCGCAAAATGGTGACTGGCACATCGGTCGCGGCGTCGCCATTATCATGCAGAAATCGGGTATCCCGGATATCGATCAGGCTAACTGCATGATCAAACTGGAATCGGACGGCACCTTTATCGTTCATTCTGGCGGTGCGGATAtt                                                                        >  3:396616/1‑149 (MQ=255)
                                                                               gAAACCGCAAAATGGTGACTGGCACATCGGTCGCGGCGTCGCCATTATCATGCAGAAATCGGGTATCCCGGATATCGATCAGGCTAACTGCATGATCAAACTGGAATCGGACGGCACCTTTATCGTTCATTCTGGCGGTGCGGATATTg                                                                       <  1:321480/149‑1 (MQ=255)
                                                                               gAAACCGCAAAATGGTGACTGGCACATCGGTCGCGGCGTCGCCATTATCATGCAGAAATCGGGTATCCCGGATATCGATCAGGCTAACTGCATGATCAAACTGGAATCGGACGGCACCTTTATCGTTCATTCTGGCGGTGCGGATATTg                                                                       >  4:142757/1‑149 (MQ=255)
                                                                                aaaCCGCAAAATGGTGACTGGCACATCGGTCGCGGCGTCGCCATTATCATGCAGAAATCGGGTATCCCGGATATCGATCAGGCTAACTGCATGATCAAACTGGAATCGGACGGCACCTTTATCGTTCATTCTGGCGGTGCGGATATTgg                                                                      <  4:402588/149‑1 (MQ=255)
                                                                                aaaCCGCAAAATGGTGACTGGCACATCGGTCGCGGCGTCGCCATTATCATGCAGAAATCGGGTATCCCGGATATCGATCAGGCTAACTGCATGATCAAACTGGAATCGGACGGCACCTTTATCGTTCATTCTGGCGGTGCGGATATTgg                                                                      >  1:283113/1‑149 (MQ=255)
                                                                                aaaCCGCAAAATGGTGACTGGCACATCGGTCGCGGCGTCGCCATTATCATGCAGAAATCGGGTATCCCGGATATCGATCAGGCTAACTGCATGATCAAACTGGAATCGGACGGCACCTTTATCGTTCATTCTGGCGGTGCGGATATTgg                                                                      <  2:110366/149‑1 (MQ=255)
                                                                                        aaaTGGTGACTGGCACATCGGTCGCGGCGTCGCCATTATCATGCAGAAATCGGGTATCCCGGATATCGATCAGGCTAACTGCATGATCAAACTGGAATCGGACGGCACCTTTATCGTTCATTCTGGCGGTGCGGATATTGGTACTGGTc                                                              >  4:438858/1‑149 (MQ=255)
                                                                                        aaaTGGTGACTGGCACATCGGTCGCGGCGTCGCCATTATCATGCAGAAATCGGGTATCCCGGATATCGATCAGGCTAACTGCATGATCAAACTGGAATCGGACGGCACCTTTATCGTTCATTCTGGCGGTGCGGATATTGGTACTGGTc                                                              >  3:22667/1‑149 (MQ=255)
                                                                                               gACTGGCACATCGGTCGCGGCGTCGCCATTATCATGCAGAAATCGGGTATCCCGGATATCGATCAGGCTAACTGCATGATCAAACTGGAATCGGACGGCACCTTTATCGTTCATTCTGGCGGTGCGGATATTGGTACTGGTCTGGATAc                                                       <  1:73820/149‑1 (MQ=255)
                                                                                               gACTGGCACATCGGTCGCGGCGTCGCCATTATCATGCAGAAATCGGGTATCCCGGATATCGATCAGGCTAACTGCATGATCAAACTGGAATCGGACGGCACCTTTATCGTTCATTCTGGCGGTGCGGATATTGGTACTGGTCTGGATAc                                                       <  1:92533/149‑1 (MQ=255)
                                                                                                aCTGGCACATCGGTCGCGGCGTCGCCATTATCATGCAGAAATCGGGTATCCCGGATATCGATCAGGCTAACTGCATGATCAAACTGGAATCGGACGGCACCTTTATCGTTCATTCTGGCGGTGCGGATATTGGTACTGGTCTGGATACt                                                      <  3:264091/149‑1 (MQ=255)
                                                                                                  tGGCACATCGGTCGCGGCGTCGCCATTATCATGCAGAAATCGGGTATCCCGGATATCGATCAGGCTAACTGCATGATCAAACTGGAATCGGACGGCACCTTTATCGTTCATTCTGGCGGTGCGGATATTGGTACTGGTCTGGATACTGt                                                    >  1:9888/1‑149 (MQ=255)
                                                                                                   ggCACATCGGTCGCGGCGTCGCCATTATCATGCAGAAATCGGGTATCCCGGATATCGATCAGGCTAACTGCATGATCAAACTGGAATCGGACGGCACCTTTATCGTTCATTCTGGCGGTGCGGATATTGGTACTGGTCTGGATACTGTa                                                   <  2:103458/149‑1 (MQ=255)
                                                                                                   ggCACATCGGTCGCGGCGTCGCCATTATCATGCAGAAATCGGGTATCCCGGATATCGATCAGGCTAACTGCATGATCAAACTGGAATCGGACGGCACCTTTATCGTTCATTCTGGCGGTGCGGATATTGGTACTGGTCTGGATACTGTa                                                   <  2:503262/149‑1 (MQ=255)
                                                                                                   ggCACATCGGTCGCGGCGTCGCCATTATCATGCAGAAATCGGGTATCCCGGATATCGATCAGGCTAACTGCATGATCAAACTGGAATCGGACGGCACCTTTATCGTTCATTCTGGCGGTGCGGATATTGGTACTGGTCTGGATACTGTa                                                   <  2:35982/149‑1 (MQ=255)
                                                                                                    gCACATCGGTCGCGGCGTCGCCATTATCATGCAGAAATCGGGTATCCCGGATATCGATCAGGCTAACTGCATGATCAAACTGGAATCGGACGGCACCTTTATCGTTCATTCTGGCGGTGCGGATATTGGTACTGGTCTGGATACTGTAg                                                  >  3:201453/1‑149 (MQ=255)
                                                                                                      acaTCGGTCGCGGCGTCGCCATTATCATGCAGAAATCGGGTATCCCGGATATCGATCAGGCTAACTGCATGATCAAACTGGAATCGGACGGCACCTTTATCGTTCATTCTGGCGGTGCGGATATTGGTACTGGTCTGGATACTGTAGTg                                                >  3:220703/1‑149 (MQ=255)
                                                                                                       caTCGGTCGCGGCGTCGCCATTATCATGCAGAAATCGGGTATCCCGGATATCGATCAGGCTAACTGCATGATCAAACTGGAATCGGACGGCACCTTTATCGTTCATTCTGGCGGTGCGGATATTGGTACTGGTCTGGATACTGTAGTGa                                               <  1:18218/149‑1 (MQ=255)
                                                                                                       caTCGGTCGCGGCGTCGCCATTATCATGCAGAAATCGGGTATCCCGGATATCGATCAGGCTAACTGCATGATCAAACTGGAATCGGACGGCACCTTTATCGTTCATTCTGGCGGTGCGGATATTGGTACTGGTCTGGATACTGTAGTGa                                               <  1:67809/149‑1 (MQ=255)
                                                                                                       caTCGGTCGCGGCGTCGCCATTATCATGCAGAAATCGGGTATCCCGGATATCCATCAGGCTAACTGCATGATCAAACTGGAATCGGACGGCACCTTAATCGTTCATTCTGGCGGTGCGGATATTGGTACAGGTCTGGATACCGTAGTGa                                               <  4:480425/149‑1 (MQ=255)
                                                                                                           ggTCGCGGCGTCGCCATTATCATGCAGAAATCGGGTATCCCGGATATCGATCAGGCTAACTGCATGATCAAACTGGAATCGGACGGCACCTTTATCGTTCATTCTGGCGGTGCGGATATTGGTACTGGTCTGGATACTGTAGTGACGaa                                           <  1:236970/149‑1 (MQ=255)
                                                                                                                  gcgTCGCCATTATCATGCAGAAATCGGGTATCCCGGATATCGATCAGGCTAACTGCATGATCAAACTGGAATCGGACGGCACCTTTATCGTTCATTCTGGCGGTGCGGATATTGGTACTGGTCTGGATACTGTAGTGACGAAACTGgca                                    <  3:60324/149‑1 (MQ=255)
                                                                                                                     tCGCCATTATCATGCAGAAATCGGGTATCCCGGATATCGATCAGGCTAACTGCATGATCAAACTGGAATCGGACGGCACCTTTATCGTTCATTCTGGCGGTGCGGATATTGGTACTGGTCTGGATACTGTAGTGACGAAAc                                         <  3:33411/141‑1 (MQ=255)
                                                                                                                     tCGCCATTATCATGCAGAAATCGGGTATCCCGGATATCGATCAGGCTAACTGCATGATCAAACTGGAATCGGACGGCACCTTTATCGTTCATTCTGGCGGTGCGGATATTGGTACTGGTCTGGATACTGTAGTGACGAAAc                                         >  4:33411/1‑141 (MQ=255)
                                                                                                                       gCCATTATCATGCAGAAATCGGGTATCCCGGATATCGATCAGGCTAACTGCATGATCAAACTGGAATCGGACGGCACCTTTATCGTTCATTCTGGCGGTGCGGATATTGGTACTGGTCTGGATACTGTAGTGACGAAACTGGCAGCAGa                               >  2:195014/1‑149 (MQ=255)
                                                                                                                       gCCATTATCATGCAGAAATCGGGTATCCCGGATATCGATCAGGCTAACTGCATGATCAAACTGGAATCGGACGGCACCTTTATCGTTCATTCTGGCGGTGCGGATATTGGTACTGGTCTGGATACTGTAGTGACGAAACTGGCAGCAGa                               >  3:18274/1‑149 (MQ=255)
                                                                                                                       gCCATTATCATGCAGAAATCGGGTATCCCGGATATCGATCAGGCTAACTGCATGATCAAACTGGAATCGGACGGCACCTTTATCGTTCATTCTGGCGGTGCGGATATTGGTACTGGTCTGGATACTGTAGTGACGAAACTGGCAGCAGa                               >  4:346006/1‑149 (MQ=255)
                                                                                                                        ccATTATCATGCAGAAATCGGGTATCCCGGATATCGATCAGGCTAACTGCATGATCAAACTGGAATCGGACGGCACCTTTATCGTTCATTCTGGCGGTGCGGATATTGGTACTGGTCTGGATACTGTAGTGACGAAACTGGCAGCAGaa                              <  1:430413/149‑1 (MQ=255)
                                                                                                                          aTTATCATGCAGAAATCGGGTATCCCGGATATCGATCAGGCTAACTGCATGATCAAACTGGAATCGGACGGCACCTTTATCGTTCATTCTGGCGGTGCGGATATTGGTACTGGTCTGGATACTGTAGTGACGAAACTGGCAGCAGAAGt                            >  3:202794/1‑149 (MQ=255)
                                                                                                                             aTCATGCAGAAATCGGGTAACCCGGATATCGATCAGGCTAACTGCATGATCAAACTGGAATCGGACGGCACCTTTATCGTTCATTCTGGCGGTGCGGATATTGGTACTGGTCTGGATACTGTAGTGACGAAACTGGCAGCAGAAGtgct                         <  4:448459/149‑1 (MQ=255)
                                                                                                                                aTGCAGAAATCGGGTATCCCGGATATCGATCAGGCTAACTGCATGATCAAACTGGAATCGGACGGCACCTTTATCGTTCATTCTGGCGGTGCGGATATTGGTACTGGTCTGGATACTGTAGTGACGAAACTGGCAGCAGAAGTGCTGCa                      >  3:301972/1‑149 (MQ=255)
                                                                                                                                aTGCAGAAATCGGGTATCCCGGATATCGATCAGGCTAACTGCATGATCAAACTGGAATCGGACGGCACCTTTATCGTTCATTCTGGCGGTGCGGATATTGGTACTGGTCTGGATACTGTAGTGACGAAACTGGCAGCAGAAGTGCTGCa                      >  4:524309/1‑149 (MQ=255)
                                                                                                                                   cAGAAATCGGGTATCCCGGATATCGATCAGGCTAACTGCATGATCAAACTGGAATCGGACGGCACCTTTATCGTTCATTCTGGCGGTGCGGATATTGGTACTGGTCTGGATACTGTAGTGACGAAACTGGCAGCAGAAGTGCTGCACTg                   <  4:163948/149‑1 (MQ=255)
                                                                                                                                    aGAAATCGGGTATCCCGGATATCGATCAGGCTAACTGCATGATCAAACTGGAATCGGACGGCACCTTTATCGTTCATTCTGGCGGTGCGGATATTGGTACTGGTCTGGATACTGTAGTGACGAAACTGGCAGCAGAAGTGCTGCACTGc                  >  3:357972/1‑149 (MQ=255)
                                                                                                                                         tCGGGTATCCCGGATATCGATCAGGCTAACTGCATGATCAAACTGGAATCGGACGGCACCTTTATCGTTCATTCTGGCGGTGCGGATATTGGTACTGGTCTGGATACTGTAGTGACGAAACTGGCAGCAGAAGTGCTGCACTGCCCAcc             >  4:140986/1‑149 (MQ=255)
                                                                                                                                                 cccGGATATCGATCAGGCTAACTGCATGATCAAACTGGAATCGGACGGCACCTTTATCGTTCATTCTGGCGGTGCGGATATTGGTACTGGTCTGGATACTGTAGTGACGAAACTGGCAGCAGAAGTGCTGCACTGCCCACCGCAGGACg     >  4:433779/1‑149 (MQ=255)
                                                                                                                                                  ccGGATATCGATCAGGCTAACTGCATGATCAAACTGGAATCGGACGGCACCTTTATCGTTCATTCTGGCGGTGCGGATATTGGTACTGGTCTGGATACTGTAGTGACGAAACTGGCAGCAGAAGTGCTGCACTGCCCACCGCAGGACGt    <  4:120282/149‑1 (MQ=255)
                                                                                                                                                    ggATATCGATCAGGCTAACTGCATGATCAAACTGGAATCGGACGTCACCTTTATCGTTCATTCTGGCGGTGCGGATATTGGTACTGGTCTGGATACTGTAGTGACGAAACTGGCAGCAGAAGTGCTGCACTGCCCACCGCCGGACGTGc  >  2:61578/1‑149 (MQ=255)
                                                                                                                                                    ggATATCGATCAGGCTAACTGCATGATCAAACTGGAATCGGACGGCACCTTTATCGTTCATTCTGGCGGTGCGGATATTGGTACTGGTCTGGATACTGTAGTGACGAAACTGGCAGCAGAAGTGCTGCACTGCCCACCGCAGGACGTGc  <  3:348659/149‑1 (MQ=255)
                                                                                                                                                    ggATATCGATCAGGCTAACTGCATGATCAAACTGGAATCGGACGGCACCTTTATCGTTCATTCTGGCGGTGCGGATATTGGTACTGGTCTGGATACTGTAGTGACGAAACTGGCAGCAGAAGTGCTGCACTGCCCACCGCAGGACGTGc  >  3:551965/1‑149 (MQ=255)
                                                                                                                                                    ggATATCGATCAGGCTAACTGCATGATCAAACTGGAATCGGACGGCACCTTTATCGTTCATGCTGGCGGTGCGGATATTGGTACTGGTCTGGATACTGTGGTGACGAAACTGGCAGCAGAAGTGCTCCACTGCCCACCGCAGGACGTGc  <  4:265231/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
CCTCCGTTCCTTCCGCCGCCAGCTGTGCACTGGAAGAGATCCTGCGTCAGGGTCGCGAGATGATCCAATGGTCGTCACCGAAACCGCAAAATGGTGACTGGCACATCGGTCGCGGCGTCGCCATTATCATGCAGAAATCGGGTATCCCGGATATCGATCAGGCTAACTGCATGATCAAACTGGAATCGGACGGCACCTTTATCGTTCATTCTGGCGGTGCGGATATTGGTACTGGTCTGGATACTGTAGTGACGAAACTGGCAGCAGAAGTGCTGCACTGCCCACCGCAGGACGTGC  >  NC_000913/3023147‑3023443

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