Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,360,644 G→C 6.5% intergenic (+28/‑532) gltD → / → gltF glutamate synthase, 4Fe‑4S protein, small subunit/periplasmic protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,360,6440GC 6.5% 391.4 / 16.1 123intergenic (+28/‑532)gltD/gltFglutamate synthase, 4Fe‑4S protein, small subunit/periplasmic protein
Reads supporting (aligned to +/- strand):  ref base G (61/53);  new base C (6/2);  total (68/55)
Fisher's exact test for biased strand distribution p-value = 2.92e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.80e-01

Only 100 of 127 total aligned reads displayed.

AGACCAGCAACCCGAAAATCTTTGCTGGCGGCGATATCGTCCGTGGTTCCGATCTGGTGGTGACCGCTATTGCCGAAGGTCGTAAGGCGGCAGACGGTATTATGAACTGGCTGGAAGTTTAAGCGAGGTAACAATGAATTCGCCAGACTGTGTCTGGCGAATCAAATTACTGACGAAGTGCGCGTTGTTTTCTGGCAAATGCGGCATGAGCGTTTGACCTTTCCTATAATAGGGATAATTTTCAATAAATTGTGCTTTGTTCACCAAGTCTGATGAGTGAAGCACATCAAGCAATTTC  >  NC_000913/3360495‑3360792
                                                                                                                                                     |                                                                                                                                                    
aGACCAGCAACCCGAAAATCTTTGCTGGCGGCGATATCGTCCGTGGATCCGAACAGCTGGAGACCGCTATTGCCGAAGGACGAAAGGCGGCAGCCGGTATTATGAACAGGCTGCAAGTTTAAGCGAGGAAACAATGAATACGCCAGACtg                                                                                                                                                      >  1:70145/1‑150 (MQ=255)
      gCAACCCGAAAATCTTTGCTGGCGGCGATATCGTCCGTGGTTCCGATCTGGTGGTGACCGCTATTGCCGAAGGTCGTAAGGCGGCAGACGGTATTATGAACTGGCTGGAAGTTTAAGCGAGGTAACAATGAATTCGCCAGACTGTGTCTg                                                                                                                                                <  2:655466/150‑1 (MQ=255)
        aaCCCGAAAATCTTTGCTGGCGGCGATATCGTCCGTGGTTCCGATCTGGTGGTGCCCGCTATTGCCGAAGGTCGTAAGGCGGCAGACGGTATTATGAACTGGCTGGAAGTTTAAGCGAGGTAACAATGAATTCGCCAGACTGTGTCTTGc                                                                                                                                              >  2:350342/1‑150 (MQ=255)
        aaCCCGAAAATCTTTGCTGGCGGCGATATCGTCCGTGGTTCCGATCTGGTGGTGACCGCTATTGCCGAAGGTCGTAAGGCGGCAGACGGTATTATGAACTGGCTGGAAGTTTAAGCGAGGTAACAATGAATTCGCCAGACTGTGTCTGGc                                                                                                                                              <  2:1742906/150‑1 (MQ=255)
         ccccGAACATCTTAGCGGGCGGCGATATCGTCCGTGGTTCCGATCTGGTGGTGACCGCTATTGCCGAAGGTCGTAAGGCGGCAGACGGTATAATCAACTGGCTGGAAGTTTAAGCGAGGTAACAATGAATTCGCCAGACTGTGTCTGGCg                                                                                                                                             <  2:182631/149‑1 (MQ=255)
         aCCCGAACATCTTTGCTGGCGGCGATATCGTCCGTGGTTCCGATCTGGTGGTGACCGCTATTGCCGAAGGTCGTAAGGCGGCAGACGGTATTATGAACTGGCTGGAAGTTTAAGCGAGGTAACAATGAATTCGCCAGACTGTGTCTGGCg                                                                                                                                             <  2:1593604/150‑1 (MQ=255)
         aCCCGAAAATCTTTGCTGGCGGCGATATCGTCCGTGGTTCCGATCTGGTGGTGACCGCTATTGCCGAAGGTCGTAAGGCGGCAGACGGTATTATGAACTGGCTGGAAGTTTAAGCGAGGTAACAATGAATTCGCCAGACTGTGTCTGGCg                                                                                                                                             <  2:1657206/150‑1 (MQ=255)
         aCCCGAAAATCTTTGCTGGCGGCGATATCGTCCGTGGTTCCGATCTGGTGGTGACCGCTATTGCCGAAGGTCGTAAGGCGGCAGACGGTATTATGAACTGGCTGGAAGTTTAAGCGAGGTAACAATGAATTCGCCAGACACAGTCTGGCg                                                                                                                                             <  1:1289112/150‑1 (MQ=255)
            cGAAAATCTTTGCTGGCGGCGATATCGTCCGTGGTTCCGATCTGGTGGTGACCGCTATTGCCGAAGGTCGTAAGGCGGCAGACGGTATTATGCACGGGCTGGAAGTTTAAGCGAGGTAACAATGAATTCGCCCGACAGTGTCTGGCGaaa                                                                                                                                          >  2:935736/1‑149 (MQ=255)
             gAAAATCTTTGCTGGCGGCGATATCGTCCGTGGTTCCGATCTGGTGGTGACCGCTATTGCCGAAGGTCGTAAGGCGGCAGACGGTATTATGAACTGGCTGGAAGTTTAAGCGAGGTAACAATGAATTCGCCAGACTGTGTCTGGCGAATc                                                                                                                                         <  1:589307/150‑1 (MQ=255)
             gAAAATCTTTGCTGGCGGCGATATCGTCCGTGGTTCCGATCTGGTGGTGACCGCTATTGCCGAAGGTCGTAAGGCGGCAGACGGTATTATGAACTGGCTGGAAGTTTAAGCGAGGTAACAATGAATTCGCCAGACTGTGTCTGGCGAATc                                                                                                                                         <  1:1605333/150‑1 (MQ=255)
               aaaTCTTTGCTGGCGGCGATATCGTCCGTGGTTCCGATCTGGTGGTGACCGCTATTGCCGAAGGTCGTAAGGCGGCAGACGGTATTATGAACTGGCTGGAAGTTTAAGCGAGGTAACAATGAATTCTCCAGACTGTGTCTGGCGAATCaa                                                                                                                                       >  1:886022/1‑150 (MQ=255)
               aaaTCTTTGCTGGCGGCGATATCGTCCGTGGTTCCGATCTGGTGGTGACCGCTATTGCCGAAGGTCGTAAGGCGGCAGACGGTATTATGAACTGGCTGGAAGTTTAAGCGAGGTAACAATGAATTCGCCAGACTGTGTCTGGCGAATCaa                                                                                                                                       <  1:1339880/150‑1 (MQ=255)
                  tCTTTGCTGGCGGCGATATCGTCCGTGGTTCCGATCTGGTGGTGACCGCTATTGCCGAAGGTCGTAAGGCGGCAGACGGTATTATGAACTGGCTGGAAGTTTAAGCGAGGTAACAATGAATTCGCCAGACTGTGTCTGGCGAATCAAAtt                                                                                                                                    <  2:376990/150‑1 (MQ=255)
                   cTTTGCTGGCGGCGATATCGTCCGTGGTTCCGATCTGGTGGTGACCTCTATTGCCGAAGGTCGTAAGGCGGCAGACGGTATTATGAACTGGCTGGAAGTTTAAGCGAGGTAACAATGAATTCGCCAGACTGTGTCTGGCGAATCAAATTa                                                                                                                                   >  1:436131/1‑150 (MQ=255)
                       gCTGGCGGCGATATCGTCCGTGGTTCCGATCTGGTGGTGACCGCTATCGCCGAAGGTCGTAAGGCGGCAGACGGTATTATGAACTGGCTGGAAGTTTAAGCGAGGTAACAATGAATTCGCCAGACTGTGTCTGGCGAATCAAATTACTGa                                                                                                                               <  1:862462/150‑1 (MQ=255)
                        cTGGCGGCGATATCGTCCGTGGTTCCGATCTGGTGGTGCCCGCTATTGCCGAAGGTCGTAAGGCGGCAGACGGTATTATGAACTGGCTGGAAGTTTAAGCGAGGTAACAATGAATTCGCCAGACTGTGTCTGGCGAATCAAATTACTGAc                                                                                                                              >  1:1206999/1‑150 (MQ=255)
                          ggcggcGATATCGTCCGTGGTTCCGATCTGGTGGTGACCGCTATTGCCGAAGGTCGTAAGGCGGCAGACGGTATTATGAACTGGCTGGAAGTTTAAGCGAGGTAACAATGAATTCGCCAGACTGTGTCTGGCGAATCAAATTACTGACGa                                                                                                                            >  2:1247891/1‑150 (MQ=255)
                             ggcgCTATCGTCCGTGGTGCCGATCTGGTGGTGCCCGCTATTGCCGAAGGTCGTAAGGCGGCAGACGGTATTATGAACTGGCTGGAAGTTTAAGCGAGGTAACAATGAATTCGCCAGACGGTGTCTGGCGAATCAAATTACTGCCGAAGt                                                                                                                         >  2:1334239/1‑150 (MQ=255)
                              gcgATATCGTCCGTGGGTCCGATCTGGTGGTGACCGCTATTGCCGAAGGTCGTAAGGCGGCAGACGGTATTATGAACTGGCTGGAAGTTTAAGCGAGGTAACAATGAATTCGCCAGACTGTGTCTGGCGAATCAAATTACTGACGAAGTg                                                                                                                        <  2:142518/150‑1 (MQ=255)
                                         cGTGGTTCCGATCTGGTGGTGACCGCTATTGCCGAAGGTCGTAAGGCGGCAGACGGTATTATGAACTGGCTGGAAGTTTAAGCGAGGTAACAATGAATTCGCCAGACTGTGTCTGGCGAATCAAATTACTGACGAAGTGCGCGTTGtttt                                                                                                             >  1:1655204/1‑150 (MQ=255)
                                          gTGGTTCCGATCTGGTGGTGACCGCTATTGCCGAAGGTCGTAAGGCGGCAGACGGTATTATGAACTGGCTGGAAGTTTAAGCGAGGTAACAATGAATTCGCCAGACTGTGTCTGGCGAATCAAATTACTGACGAAGTGCGCGTTGtttta                                                                                                            >  2:377413/1‑149 (MQ=255)
                                           tGGTTCCGATCTGGTGGTGACCGCTATTGCCGAAGGTCGTAAGGCGGCAGACGGTATTATGAACTGGCTGGAAGTTTAAGCGAGGTAACAATGAATTCGCCAGACTGTGTCTGGCGAATCAAATTACTGACGAAGTGCGCGTTGTTTTCt                                                                                                           <  1:845354/150‑1 (MQ=255)
                                             gTTCCGATCTGGTGGTGACCGCTATTGCCGAAGGTCGTAAGGCGGCAGACGGTATTATGAACTGGCTGGAAGTTTAAGCGAGGTAACAATGAATTCGCCAGACTGTGTCTGGCGAATCAAATTACTGACGAAGTGCGCGTTGTTTTCTgg                                                                                                         >  1:1458344/1‑150 (MQ=255)
                                             gTTCCGATCTGGTGGTGACCGCTATTGCCGAAGGTCGTAAGGCGGCAGACGGTATTATGAACTGGCTGGAAGTTTAAGCGAGGTAACAATGAATTCGCCAGACTGTGTCTGGCGAATCAAATTACTGACGAAGTGCGCGTTGTTTTCTgg                                                                                                         <  1:914208/150‑1 (MQ=255)
                                             gTTCCGATCTGGTGGTCACCGCTATTGCCGAAGGTCGTAAGGCGGCAGACGGTATTATGAACTCGCTGGAAGTTTAAGCGAGGTAACAATGAATTCGCCAGACTGTGTCTGGCGAATCAAATTACTGACGAAGTGCGCGTTGTTTTCTgg                                                                                                         <  2:970834/150‑1 (MQ=255)
                                              ttCCGATCTGGTGGTGACCGCTATTGCCGAAGGTCGTAAGGCGGCAGACGGTATTATGAACTGGCTGGAAGTTTAAGCGAGGTAACAATGAATTCGCCAGACTGTGTCTGGCGAATCAAATTACTGACGAAGTGCGCGTTGTTTTCTTGc                                                                                                        >  1:786106/1‑150 (MQ=255)
                                               tCCGATCTGGGGGTGACCGCTATTGCCGAAGGGCGTAAGGCGGCAGACGGTATTATGAACTGGCTGGAAGTTTAAGCGAGGTAACAATGAATTCGCCAGACTGTGTCTGGCGAATCAAATTACTGCCGAAGTGCGCGTTGTTTTCTGGCa                                                                                                       <  2:1777025/150‑1 (MQ=255)
                                                 cGATCTGGTGGTGACCGCTATTGCCGAAGGTCGTAAGGCGGCAGACGGTATTATGAACTGGCTGGAAGTTTAAGCGAGGTAACAATGAATTCGCCAGACTGTGTCTGGCGAATCAAATTACTGACGAAGTGCGCGTTGTTTTCTGGCaaa                                                                                                     >  1:1145943/1‑150 (MQ=255)
                                                 cGATCTGGTGGTGACCGCTATTGCCGAAGGTCGTAAGGCGGCAGACGGTATTATGAACTGGCTGGAAGTTTAAGCGAGGTAACAATGAATTCGCCAGACTGTGTCTGGCGAATCAAATTACTGACGAAGTGCGCGTTGTTTTCTGGCaaa                                                                                                     >  1:856601/1‑150 (MQ=255)
                                                     cTGGTGGTGACCGCTATTGCCGAAGGTCGTAAGGCGGCAGACGGTATTATGAACTGGCTGGAAGTTTAAGCGAGGTAACAATGAATTCGCCAGACTGTGTCTGGCGAATCAAATTACTGACGAAGTGCGCGTTGTTTTCTGGCAAATGCg                                                                                                 <  1:419027/150‑1 (MQ=255)
                                                          ggtgACCGCTATTGCCGAAGGTCGTAAGGCGGCAGACGGTATTATGAACTGGCTGGAAGTTTAAGCGAGGTAACAATGAATTCGCCAGACTGTGTCTGGCGAATCAAATTACTGACGAAGTGCGCGTTGTTTTCTGGCAAATGCGGCATg                                                                                            <  2:1558540/150‑1 (MQ=255)
                                                          ggtgACCGCTATTGCCGAAGGTCGTAAGGCGGCAGACGGTATTATGAACTGGCTGGAAGTTTAAGCGAGGTAACAATGAATTCGCCAGACTGTGTCTGGCGAATCAAATTACTGACGAAGTGCGCGTTGTTTTCTGGCAAATGCGGCATg                                                                                            >  1:1777969/1‑150 (MQ=255)
                                                           gtgACCGCTATTGCCGAAGGTCGTAAGGCGGCAGACGGTATTATGAACTGGCTGGAAGTTTAAGCGAGGTAACAATGAATTCGCCAGACTGTGTCTGGCGAATCAAATTACTGACGAAGTGCTAGTTGTTTTATGGAAAATGCGGCATGa                                                                                           >  1:1107776/1‑150 (MQ=255)
                                                             gACCGCTATTGCCGAAGGTCGTAAGGCGGCAGACGGTATTATGAACTGGCTGGTAGTTTAAGCGAGGTAACAATGAATTCGCCAGACTGTGTCTGGCGAATCAAATTACGGACGAAGTGCGCGTTGTTTTCTGGCAAATGCGGCATgcgc                                                                                         >  1:155562/1‑150 (MQ=255)
                                                             gACCGCTATTGCCGAAGGTCGTAAGGCGGCAGACGGTATTATGAACTGGCTGGAAGTTTAAGCGAGGTAACAATGAATTCGCCAGACTTTGTCAGGCGAATCAAATTAATGACGAAGTTCGCGTTGTTATATGGCAAAATCGGCATGAGc                                                                                         >  2:1674945/1‑150 (MQ=255)
                                                               ccGCTATTGCCGAAGGTCGTAAGGCGGCAGACGGTATTATGAACTGGCTGGAAGTTTAAGCGAGGTAACAATGAATTCGCCAGACTGTGTCTGGCGAATCAAATTACTGACGAAGTGCGCGTTGTTTTCTGGCAAATGCGGCATGAGCGt                                                                                       >  1:318694/1‑150 (MQ=255)
                                                               ccGCTATTGCCGAAGGTCGTAAGGCGGCAGACGGTATTATGAACTGGCTGGAAGTTTAAGCGAGGTAACAATGAATTCGCCAGACTGTGTCTGGCGAATCAAATTACTGACGAAGTGCGCGTTGTTTTCTGGCAAATGCGGCATGAGCGt                                                                                       >  2:738315/1‑150 (MQ=255)
                                                                 gCTATTGCCGAAGGTCGTAAGGCGGCAGACGGTATTATGAACTGGCTGGAAGTTTAAGCGAGGTAACAATGAATTCGCCAGACTGTGTCTGGCGAATCAAATTACTGACGAAGTGCGCGTTGTTTTCTGGCAAATGCGGCATGAGCGttt                                                                                     <  2:1291877/150‑1 (MQ=255)
                                                                  cTATTGCCTAAGGTCGTAAGGCGGCAGACGGTATTATGTACTGGCGGGAAGTTTAAGCGAGGTAACAATGAATTCGCCAGACTGTGGCTGGCGAATCAAATTACTGACGAAGTGCGCGTTGTTTTCTGGCAAATGCGGCATGTGCGTTTg                                                                                    <  1:341832/150‑1 (MQ=255)
                                                                  cTATTGCCGAAGGTCGTAAGGCGGCAGACGGTATTATGAACTGGCTGGAAGTTTAAGCGAGGTAACAATGAATTCGCCAGACTGGGTCTGGCGAATCAAATTACTGACGAACTGCTGTTTGATTTCTGGCAAATCCGGAATGAGCGtttt                                                                                    >  2:1402687/1‑149 (MQ=255)
                                                                    aTTGCCGAAGGTCGTAAGGCGGCAGACGGTATTATGAACTGGCTGGAAGTTTAAGCGAGGTAACAATGAATTCGCCAGACTGTGTCTGGCGAATCAAATTACTGACGAAGTGCGCGTTGTTTTCTGGCAAATGCGGCATGAGCGTTTGAc                                                                                  >  1:245398/1‑150 (MQ=255)
                                                                     ttGCCGAAGGTCGTAAGGCGGCAGACGGTATTATGAACTGGCTGGAAGTTTAAGCGAGGTAACAATGAATTCGCCAGACTGTGTCTGGCGAATCAAATTACTGACGAAGTGCGCGTTGTTTTCTGGCAAATGCGGCATGAGCGTTTGAcc                                                                                 >  2:1756207/1‑150 (MQ=255)
                                                                      tGCCGAAGGTCGTAAGGCGGCAGACGGTATTATGAACTGGCTGGAAGTTTAAGCGAGGTAACAATGAATTCGCCAGACTGTGTCTGGCGAATCAAATTACTGACGAAGTGCGCGTTGTTTTCTGGCAAATGCGGCATGAGCGTTTGACCt                                                                                >  2:1402299/1‑150 (MQ=255)
                                                                       gCCGAAGGTCGTAAGGCGGCAGACGGTATTATGAACTGGCTGGAAGTTTAAGCGAGGTAACAATGAATTCGCCAGACTGTGTCTGGCGAATCAAATTACTGACGAAGTGCGCGTTGTTTTCTGGCAAATGCGGCATGAGCGTTTGACCtt                                                                               >  1:1004111/1‑150 (MQ=255)
                                                                         cGAAGGTCGTAAGGCGGCAGACGGTATTATGAACTGGCTGGAAGTTTAAGCGAGGTAACAATGAATTCGCCAGACTGTGTCTGGCGAATCAAATTACTGACGAAGTGCGCGTTGTTTTCTGGCAAATGCGGCATGAGCGTTTGACCTTTc                                                                             <  2:1555224/150‑1 (MQ=255)
                                                                         cGAAGGGCGGAATGCGGGAGTCGGTATTTTGCACTGGCTGGAAGTTTAATCGATGTAACAATGAATTGGCCCTTATGTGTCTGGGGTATCAAAATACTGACGTAGTGCGCGGTGGTTTTTGGCAAATGCGGCATGAGCGTTTGACCTTTc                                                                             <  1:897838/150‑1 (MQ=255)
                                                                             ggTCGTAAGGCGGCAGACGGTATTATGAACTGGCTGGAAGTTTAAGCGAGGTAACAATGAATTCGCCAGACTGTGTCTGGCGAATCAAATTACTGACGAAGTGCGCGTTGTTTTCTGGCAAATGCGGCATGAGCGTTTGACCTTTCCtat                                                                         >  1:978392/1‑150 (MQ=255)
                                                                                ggTAAGGCGGCAGACGGTATTATGAACTGGCTGGAAGTTTAAGCGAGGTAACAATGAATTCGCCAGACTGTGTCTGGCGAATCAAATTACTGACGAAGTGCGCGTTGTTTTCTGGCAAATGCGGCATGAGCGTTTGACCTTTCCTataat                                                                      <  1:368496/149‑1 (MQ=255)
                                                                                 gTAAGGAGGCAGACGGTATTATGAACTGTCTGGATGTTTAAGCGAGGTAACAATGAATTCGCCTTTCTGTGTCTGGCGAATCAAATTACTGACGAAGTGCGCGTTGTTTTCTGGCAAATGCGGCATGAGCGTTTGACCTTTCCTataata                                                                     <  2:70145/150‑1 (MQ=255)
                                                                                   aaGGCGGCAGACGGTATTATGAACTGGCTGGAAGTTTAAGCGAGGTAACAATGAATTCGCCAGACTGTGTCTGGCGAATCAAATTACTGACGAAGTGCGCGTTGTTTTCTGGCAAATGCGGCATGAGCGTTTGACCTTTCCTATAATAgg                                                                   <  2:1577559/150‑1 (MQ=255)
                                                                                   aaGGCGGCAGACGGTATTATGAACTGGCTGGAAGTTTAAGCGAGGTAACAATGAATTCGCCAGACTGTGTCTGGCGAATCAAATTACTGACGAAGTGCGCGTTGTTTTCTGGCAAATGCGGCATGAGCGTTTGACCTTTCCTATAATAgg                                                                   >  1:1802534/1‑150 (MQ=255)
                                                                                    aGGCGGCAGACGGTATTATGAACTGGCTGGAAGTTTAAGCGAGGTAACAATGAATTCGCCAGACTGTGTCTGGCGAATCAAATTACTGACGAAGTGCGCGTTGTTTTCTGGCAAATGCGGCATGAGCGTTTGACCTTTCCTATAATAggg                                                                  <  2:1362413/150‑1 (MQ=255)
                                                                                     ggcggcAGACGGTATTATGAACTGGCTGGAAGTTTAAGCGAGGTAACAATGAATTCGCCAGACTGTGTCTGGCGAATCAAATTACTGACGAAGTGCGCGTTGTTTTCTGGCAAATGCGGCATGAGCGTTTGACCTTTCCTATAATAGGGa                                                                 >  2:1442660/1‑150 (MQ=255)
                                                                                     ggcggcAGACGGTATTATGAACTGGCTGGAAGTTTAAGCGAGGTAACAATGAATTCGCCAGACTGTGTCTGGCGAATCAAATTACTGACGAAGTGCGCGTTGTTTTCTGGCAAATGCGGCATGAGCGTTTGACCTTTCCTATAATAGGGa                                                                 <  1:180647/150‑1 (MQ=255)
                                                                                      gcggcAGACGGTATTATGAACTGGCTGGAAGTTTAAGCGAGGTAACAATGAATTCGCCAGACTGTGTCTGGCGAATCAAATTACTGACGAAGTGCGCGTTGTTTTCTGGCAAATGCGGCATGAGCGTTTGACCTTTCCTATAATAGGGAt                                                                >  2:24293/1‑150 (MQ=255)
                                                                                       cggcAGACGGTATTATGAACTGGCTGGAAGTTTAAGCGAGGTAACAATGAATTCGCCAGACTGTGTCTGGCGAATCAAATTACTGACGAAGTGCGCGTTGTTTTCTGGCAAATGCGGCATGAGCGTTTGACCTTTCCTATAATAGGGATa                                                               >  2:1019055/1‑150 (MQ=255)
                                                                                        ggcAGACGGTATTATGAACTGGCTGGAAGTTTAAGCGAGGTAACAATGAATTCGCCAGACTGTGTCTGGCGAATCAAATTACTGACGAAGTGCGCGTTGTTTTCTGGCAAATGCGGCATGAGCGTTTGACCTTTCCTATAATAGGGATaa                                                              >  1:659905/1‑150 (MQ=255)
                                                                                          cAGACGGTATTATGAACTGGCTGGAAGTTTAAGCGAGGTAACAATGAATTCGCCAGACTGTGTCTGGCGAATCAAATTTAAGACGAACTGCGCGTTGTTTTATGGCAAATGCGGCATGAGCGTTTGACCTTTACTATAATAGGGATAAtt                                                            >  1:874360/1‑150 (MQ=255)
                                                                                          cAGACGGTATTATGAACTGGCTGGAAGTTTAAGCGAGGTAACAATGAATTCGCCAGACTGTGTCTGGCGAATCAAATTACTGACGAAGTGCGCGTTGTTTTCTGGCAAATGCGGCATGAGCGTTTGACCTTTCCTATAATAGGGATAAtt                                                            >  1:1873620/1‑150 (MQ=255)
                                                                                            gACGGTATTATGAACTGGCTGGAAGTTTAAGCGAGGTAACAATGAATTCGCCAGACTGTGTCTGGCGAATCAAATTACTGACGAAGTGCGCGTTGTTTTCTGGCAAATGCGGCATGAGCGTTTGACCTTTCCTATAATAGGGATAAtttt                                                          >  1:1589794/1‑150 (MQ=255)
                                                                                            gACGGTATTATGAACTGGCTGGAAGTTTAAGCGAGGTAACAATGAATTCGCCAGACTGTGTCTGGCGAATCAAATTACTGACGAAGTGCGCGATGTTTTCTGGCAAATGCGGCATGAGCGTTTGACCTTTCCTATAATAGGGATAAtttt                                                          >  1:1005484/1‑150 (MQ=255)
                                                                                             aCGTTATTATGAACTGTCTTGAAGTTTAAGAGAGTTAACAATGAATTCGCCAGACTGTGTCTGGCGAATCAAATTACTGACGAAGTGCGCGTTGTTTTCTGGCAAATGCGGCATGAGCGTTTGACCTTTCCTATAATAGGGATAATTTTc                                                         <  1:1313093/150‑1 (MQ=255)
                                                                                             aCGGTATTATGAACTGGCTGGAAGTTTAAGCGAGGTAACAATGAATTCGCCAGACTGTGTCTGGCGAATCAAATTACTGACGAAGTGCGCGTTGTTTTCTGGCAAATGCGGCATGAGCGTTTGACCTTTCCTATAATAGGGATAATTTTc                                                         <  1:738315/150‑1 (MQ=255)
                                                                                               ggTATTATGAACTGGCTGGAAGTTTAAGCGAGGTAACAATGAATTCGCCAGACTGTGTCTGGCGAATCAAATTACTGACGAAGTGCGCGTTGTTTTCTGGCAAATGCGGCATGAGCGTTTGACCTTTCCTATAATAGGGATAATTTTCaa                                                       <  1:377413/150‑1 (MQ=255)
                                                                                                     atGAACTGGCTGGAAGTTTAAGCGAGGTAACAATGAATTCGCCAGACTGTGTCTGGCGAATCAAATTACTGACGAAGTGCGCGTTGTTTTCTGGCAAATGCGGCATGAGCGTTTGACCTTTCCTATAATAGGGATAATTTTCAATAAAtt                                                 >  2:256688/1‑150 (MQ=255)
                                                                                                      tGAACTGGCTGGAAGTTTAAGCGAGGTAACAATGAATTCGCCAGACTGTGTCTGGCGAATCAAATTACTGACGAAGTGCGCGTTGTTTTCTGGCAAATGCGGCATGAGCGTTTGACCTTTCCTATAATAGGGATAATTTTCAATAAATtg                                                <  1:534006/150‑1 (MQ=255)
                                                                                                       tatctggatggAAGTTTAAGCGAGGTCACAATGAATTCGCCAGACTGTGTCTGGCGAATCAAATTACTGCCGAAGTGCGCGGTGGTTTCTGGCAAATGCGGCATGAGCGTTTGACCTTTCCTATAATAGGGATAATTTTCAATAAATtgt                                               <  1:891096/147‑1 (MQ=255)
                                                                                                       gAACTGGCTGGAAGTTTAAGCGAGGTAACAATGAATTCGCCAGACTGTGTCTGGCGAATCAAATTACTGACGAAGTGCGCGTTGTTTTCTGGCAAATGCGGCATGAGCGTTTGACCTTTCCTATAATAGGGATAATTTTCAATAAATtgt                                               >  2:659058/1‑150 (MQ=255)
                                                                                                        aaCTGTCTGGAAGTTTAAGCGAGGTAACAATGAATTCGCCAGACTGGGTCTGGCGAATCAAATTACTGACGAAGTGCGCGTTGTTTTCTGGCAAATGCGGCATGAGCGTTTGACCTTTCCTATAATAGGGATAATTATCAATAAATtgtg                                              <  2:76053/150‑1 (MQ=255)
                                                                                                         aCTGGCTGGAAGTTTAAGCGAGGTAACAATGAATTCGCCAGACTGTGTCTGGCGAATCAAATTACTGACGAAGTGCGCGTTGTTTTCTGGCAAATGCGGCATGAGCGTTTGACCTTTCCTATAATAGGGATAATTTTCAATAAATTGTGc                                             <  1:1410862/150‑1 (MQ=255)
                                                                                                            ggctggAAGTTTAAGCGAGGTAACAATGAATTCGCCAGACTGTGTCTGGCGAATCAAATTACTGACGAAGTGCGCGTTGTTTTCTGGCAAATGCGGCATGAGCGTTTGACCTTTCCTATAATAGGGATAATTTTCAATAAATTGTGCttt                                          >  1:1156195/1‑150 (MQ=255)
                                                                                                             gctggAAGTTTAAGCGAGGTAACAATGAATTCGCCAGACTGTGTCTGGCGAATCAAATTACTGACGAAGTGCGCGTTGTTTTCTGGCAAATGCGGCATGAGCGTTTGACCTTTCCTATAATAGGGATAATTTTCAATAAATTGTGCTTTg                                         <  1:449914/150‑1 (MQ=255)
                                                                                                                 gAAGTTTAAGCGAGGTAACAATGAATTCGCCAGACTGTGTCTGGCGAATCAAATTACTGACGAAGTGCGCGTTGTTTTCTGGCAAATGCGGCATGAGCGTTTGACCTTTCCTATAATAGGGATAATTTTCAATAAATTGTGCTTTGTTCa                                     <  1:1040825/150‑1 (MQ=255)
                                                                                                                  aaGTTTAAGCGAGGTAACAATGAATTCGCCAGACTGTGTCTGGCGAATCAAATTACTGACGAAGTGCGCGTTGTTTTCTGGCAAATGCGGCATGAGCGTTCGACCTTTCCTATAATAGGGATAATTTTCAATAAATTGTGCTTTGTTCAc                                    <  2:758775/150‑1 (MQ=255)
                                                                                                                     tttAAGCGAGGTAACAATGAATTCGCCAGACTGTGTCTGGCGACTCAAATTACTGACGAAGTGCGCGTTGTTTTCTGGCAAATGCGGCATGAGCGTTTGACCTTTCCTATAATAGGGATAATTTTCAATAAATTGTGCTTTGTTCAcctc                                 >  1:312561/1‑148 (MQ=255)
                                                                                                                     tttAAGCGAGGTAACAATGAATTCGCCAGACTGTGTCTGGCGAATCAAATTACTGACGAAGTGCGCGTTGTTTTCTGGCAAATGCGGCATGAGCGTTTGACCTTTACTATAATAGGGATAATTTTCAATAAATTTTGCTTTGTTCACCaa                                 >  1:1393403/1‑150 (MQ=255)
                                                                                                                      ttAAGCGAGGTAACAATGAATTCGCCAGACTGTGTCTGGCGAATCAAATTACTGACGAAGTGCGCGTTGTTTTCTGGCAAATGCGGCATGAGCGTTTGACCTTTCCTATAATAGGGATAATTTTCAATAAATTGTGCTTTGTTCACCAAg                                >  2:657271/1‑150 (MQ=255)
                                                                                                                        aaGCGAGGTAACAAAGAATTCGCCAGACTGTGTCTGGCGAATCAAATTACTGACGAAGTGCGCGTTGTTTTCTGGCAAATGCGGCATGAGCGTTTGACCTTTCCTATAATAGGGATAATTTTCAATAAATTGTGCTTTGTTCACCAAGTc                              <  2:1659759/150‑1 (MQ=255)
                                                                                                                          gggAGGTCACAATGATTTTGCCAGAATGTGTCTGGCGAATCAAATTACTGACGACGTGCGCGTTGTTTTCTGGCTAATGCGTCATGAGCGTTTGACCTTTCCTATAATAGGGATAATTTTCAATAAATTGTGCTTTGTTCACCAAGTCtg                            <  1:935066/148‑1 (MQ=255)
                                                                                                                            gAGGTAACAATGAATTCGCCAGACTGTGTCTGGCGAATCAAATTACTGACGAAGTGCGCGTTGTTTTCTGGCAAATGCGGCATGAGCGTTTGACCTTTCCTATAATAGGGATAATTTTCAATAAATTGTGCTTTGTTCACCAAGTCtgat                          >  1:278165/1‑150 (MQ=255)
                                                                                                                             aGGTAACAATGAATTCGCCAGACTGTGTCTGGCGAATCAAATTACTGACGAAGTGCGCGTTGTTTTCTGGCAAATGCGGCATGAGCGTTTGACCTTTCCTATAATAGGGATAATTTTCAATAAATTGTGCTTTGTTCACCAAGTCtgatg                         >  1:376577/1‑150 (MQ=255)
                                                                                                                             aGGTAACAATGAATTCGCCAGACTGTGTCTGGCGAATCAAATTACTGACGAAGTGCGCGTTGTTTTCTGGCAAATGCGGCATGAGCGTTTGACCTTTCCTATAATAGGGATAATTTTCAATAAATTGTGCTTTGTTCACCAAGTCtgatg                         <  1:68644/150‑1 (MQ=255)
                                                                                                                                tAACAATGAATTCGCCAGACTGTGTCTGGCGAATCAAATTACTGACGAAGTGCGCGTTGTTTTCTGGCAATTGCGGCATGAGCGTTTGACCTTTCCTATAATAGGGATAATTTTCAATAAATTGTGCTTTGTTCACCAAGTCTGATGAGt                      >  2:123620/1‑150 (MQ=255)
                                                                                                                                 aaCAATGAATTCGCCAGACTGTGTCTGGCGAATCAAATTACTGACGAAGTGCGCGTTGTTTTCTGGCAAATGCGGCATGAGCGTTTGACCTTTCCTATAATAGGGATAATTTTCAATAAATTGTGCTTTGTTCACCAAGTCTGATGAGTg                     <  2:96670/150‑1 (MQ=255)
                                                                                                                                 aaCAATGAATTCGCCAGACTGTGTCTGGCGAATCAAATTACTGACGAAGTGCGCGTTGTTTTCTGGCAAATGCGGCATGAGCGTTTGACCTTTCCTATAATAGGGATAATTTTCAATAAATTGTGCTTTGTTCACCAAGTCTGATGAGTg                     <  1:350342/150‑1 (MQ=255)
                                                                                                                                 aaCAATGAATTCGCCAGACTGTGTCTGGCGAATCAAATTACTGACGAAGTGCGCGTTGTTTTCTGGCAAATGCGGCATGAGCGTTTGACCTTTCCTATAATAGGGATAATTTTCAATAAATTGTGCTTTGTTCACCAAGTCTGATGAGTg                     >  2:634470/1‑150 (MQ=255)
                                                                                                                                 aaCAATGAATTCGCCAGACTGTGTCTGGCGAATCAAATTACTGACGAAGTGCGCGTTGTTTTCTGGCAAATGCGGCATGAGCGTTTGACCTTTCCTATAATAGGGATAATTTTCAATAAATTGTGCTTTGTTCACCAAGTCTGATGAGTg                     >  1:503678/1‑150 (MQ=255)
                                                                                                                                    aaTGAATTCGCCAGACTGTGTCTGGCGAATCAAATTACTGACGAAGTGCGCGTTGTTTTCTGGCAAATGCGGCATGAGCGTTTGACCTTTCCTATAATAGGGATAATTTTCAATAAATTGTGCTTTGTTCACCAAGTCTGATGAGTGAAg                  >  1:1001280/1‑150 (MQ=255)
                                                                                                                                       tgaTTCGCCAGACACAGTCTGGCGAATCAAATTACTGACGAAGTGCGCGTTGTTTTCTGGCAAATGCGGCATGAGCGTTTGACCTTTCCTATAATAGGGATAATTTTCATTAAATTGTGCTTTGTTCACCAAGTCTGATGAGTGAAGctc               >  1:217146/3‑148 (MQ=255)
                                                                                                                                       tgaTTCGCCAGACACAGTCTGGCGAATCAAATTACTGACGAAGTGCGCGTTGTTTTCTGGCAAATGCGGCATGAGCGTTTGACCTTTCCTATAATAGGGATAATTTTCAATAAATTGTGCTTTGTTCACCAAGTCTGATGAGTGAAGcac               >  2:1761978/3‑150 (MQ=255)
                                                                                                                                       gAATTCGCCAGACTGTGTCTGGCGAATCAAATTACTGACGAAGTGCGCGTTGTTTTCTGGCAAATGCGGCATGAGCGTTTGACCTTTCCTATAATAGGGATAATTTTCAATAAATTGTGCTTTGTTCACCAAGTCTGATGAGTGAAGcac               <  1:653512/150‑1 (MQ=255)
                                                                                                                                         aTTCGCCAGACACAGTCTGGCGAATCAAATTACTGACGAAGTGCGCGTTGTTTTCTGGCAAATGCGGCATGAGCGTTTGACCTTTCCTATAATAGGGATAATTTTCAATAAATTGTGCTTTGTTCACCAAGTCTGATGAGTGAAGCCCAt             >  2:1226195/1‑150 (MQ=255)
                                                                                                                                          ttCGCCAGACACAGTCTGGCGAATCAAATTACTGACGAAGTGCGCGTTGTTTTCTGGCAAATGCGGCATGAGCGTTTGACCTTTCCTATAATAGGGATAATTTTCAATAAATTGTGCTTTGTTCACCAAGTCTGATGAGTGAAGCACATc            >  2:131429/1‑150 (MQ=255)
                                                                                                                                          ttCGCCAGACACAGTCTGGCGAATCAAATTACTGACGAAGTGCGCGTTGTTTTCTGGCAAATGCGGCATGAGCGTTTGACCTTTCCTATAATAGGGATAATTTTCAATAAATTGTGCTTTGTTCACCAAGTCTGATGAGTGAAGCACATc            >  1:622037/1‑150 (MQ=255)
                                                                                                                                          ttCGCCAGACACAGTCTGGCGAATCAAATTACTGACGAAGTGAGCGTTGTTTTCTGGCAAATGCGGCATGAGCGTTTGACCTTTCCTATAATAGGGATAATTTTCAATAAATTGTGCTTTGTTCACCAAGTCTGATGAGTGAAGCACATc            >  2:1062991/1‑150 (MQ=255)
                                                                                                                                            cGCCAGACTGTGTCTGGCGAATCAAATTACTGACGAAGTGCGCGTTGTTTTCTGGTAAATGCGGCATGAGCGTTTGACCTTTCCTATAATAGGGATAATTTTCAATAAATTGTGCTTTGTTCACCAAGTCTGATGAGTGAAGCACATCaa          >  1:482538/1‑150 (MQ=255)
                                                                                                                                            cGCCAGACTGTGTCTGGCGAATCAAATTACTGACGAAGTGCGCGTTGTTTTCTGGCAAATGCGGCATGAGCGTTTGACCTTTCCTATAATAGGGATAATTTTCAATAAATTGTGCTTTGTTCACCAAGTCTGATGAGTGAAGCACATCaa          <  1:566094/150‑1 (MQ=255)
                                                                                                                                                  aCTGTGTCTGGCGAATCAAATTACTGACGAAGTGCGCGTTGTTTTCTGGCAAATGCGGCATGAGCGTTTGACCTTTCCTATAATAGGGATAATTTTCAATAAATTGTGCTTTGTTCACCAAGTCTGATGAGTGAAGCACATCAAGCAAtt    >  2:672035/1‑150 (MQ=255)
                                                                                                                                                    tgtgTATGGCGAATCAAATTACTGACGAAGTGCGCGTTGTTTTCTGGCAAATGCGGCATGAGCGTTTGACCTTTCCTATAATAGGGATAATTTTCAATAAATTGTGCTTTGTTCACCAAGTCTGATGAGTGAAGCACATCAAGCAATTtc  <  1:1325936/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
AGACCAGCAACCCGAAAATCTTTGCTGGCGGCGATATCGTCCGTGGTTCCGATCTGGTGGTGACCGCTATTGCCGAAGGTCGTAAGGCGGCAGACGGTATTATGAACTGGCTGGAAGTTTAAGCGAGGTAACAATGAATTCGCCAGACTGTGTCTGGCGAATCAAATTACTGACGAAGTGCGCGTTGTTTTCTGGCAAATGCGGCATGAGCGTTTGACCTTTCCTATAATAGGGATAATTTTCAATAAATTGTGCTTTGTTCACCAAGTCTGATGAGTGAAGCACATCAAGCAATTTC  >  NC_000913/3360495‑3360792

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