Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,360,644 G→C 7.8% 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 7.8% 432.2 / 20.3 138intergenic (+28/‑532)gltD/gltFglutamate synthase, 4Fe‑4S protein, small subunit/periplasmic protein
Reads supporting (aligned to +/- strand):  ref base G (65/62);  new base C (6/5);  total (71/67)
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.76e-01

Only 100 of 142 total aligned reads displayed.

AGACCAGCAACCCGAAAATCTTTGCTGGCGGCGATATCGTCCGTGGTTCCGATCTGGTGGTGACCGCTATTGCCGAAGGTCGTAAGGCGGCAGACGGTATTATGAACTGGCTGGAAGTTTAAGCGAGGTAACAATGAATTCGCCAGACTGTGTCTGGCGAATCAAATTACTGACGAAGTGCGCGTTGTTTTCTGGCAAATGCGGCATGAGCGTTTGACCTTTCCTATAATAGGGATAATTTTCAATAAATTGTGCTTTGTTCACCAAGTCTGATGAGTGAAGCACATCAAGCAATTTC  >  NC_000913/3360495‑3360792
                                                                                                                                                     |                                                                                                                                                    
aGACCAGCAACCCGAAAATCTTTGCTGGCGGCGATATCGTCCGGGGATCCGATCTCGTGGTGACCGCTATTGCCGAAGGTCGTAAGGCGGCAGACGGTATTATGAACTGGCTGGAAGTTTAAGCGAGGTAACAATGAATTCGCCAGACtg                                                                                                                                                      <  2:1314234/150‑1 (MQ=255)
   ccAGCAACCCGAAAATCTTTGCTGGCGGCGATATCGTCCGTGGTTCCGATCTGGTGGTGACCGCTATTGCCGAAGGTCGTAAGGCGGCAGCCGGTATTATGAACTGGCTGGAAGTTTAAGCGAGGTAACAATGAATTCGCCAGACtgtgg                                                                                                                                                   >  2:977010/1‑149 (MQ=255)
   ccAGCAACCCGAAAATCTTTGCTGGCGGCGATATCGTCCGTGGTTCCGATCTGGTGGTGACCGCTATTGCCGAAGGTCGTAAGGCGGCAGACGGTATTATGATCTGGCTGGAAGTTTAAGCGAGGTAACAATGAATTCGCCAGACTgtgt                                                                                                                                                   >  1:714382/1‑150 (MQ=255)
    cAGCAACCCGAAAATCTTTGCTGGCGGCGATATCGTCCGTGGTTCCGATCTGGTGGTGACCGCTATTGCCGAAGGTCGTAAGGCGGCAGACGGTATTATGAACTGGCTGGAAGTTTAAGCGAGGTAACAATGAATTCGCCAGACTGTGTc                                                                                                                                                  <  1:1708608/150‑1 (MQ=255)
     aGCAACCCGAAAATCTTTGCTGGCGGCGATATCGTCCGTGGTTCCGATCTGGTGGTGACCGCTATTGCCGAAGGTCGTAAGGCGGCAGACGGTATTATGAACTGGCTGGAAGTTTAAGCGAGGTAACAATGAATTCGCCAGACTGTGTCt                                                                                                                                                 >  2:1356797/1‑150 (MQ=255)
         aCCCGAAAATCTTTGCTGGCGGCGATATCGTCCGTGGTTCCGATCTGGTGGTGACCGCTATTGCCGAAGGTCGTAAGGCGGCAGACGGTATTATGAACTGGCTGGAAGTTTAAGCGAGGTAACAATGAATTCGCCAGACACAGTCTGGCg                                                                                                                                             <  2:1136804/150‑1 (MQ=255)
          cccGAAAATCATTGCTGGCGGCGATATCGTCCGTGGTTCCCAACTGGTGGTGACCGCTATTCCCGAAGGTCGTAACGCGGCACACGCTATAAAGAACTGGCTGGAAGTATAAGCGAGGTAACAATCAATTCGCCAGACAGTGTCTGGCGa                                                                                                                                            >  1:1422187/1‑150 (MQ=255)
            cgtaaaTCTTTGCTGGCGGCGATATCGGCCGCGGTTCCGATCTGGTGGGGACCGCTATTGCCGAAGGCCGTAAGGCGGCAGACTGTATTATGAACTGGCTGGAAGTTTAAGCGAGGTAACAATGAATTCGCCAGACACAGGCTGGCGAAt                                                                                                                                          <  1:1254324/147‑1 (MQ=255)
             gAAAATCTTTGCTGGCGGCGATATCGTCCGTGGTTCCGATCTGGTGGTGACCGCTATTGCCGAAGGTCGTAAGGCGGCAGACGGTATTCTGAACTGGCTGGAAGTTTAAGCGAGGTAACAATGAATTCGCCAGACTGTGTCTGGCGAATc                                                                                                                                         <  1:741790/150‑1 (MQ=255)
             gAAAATCTTTGCTGGCGGCGATATCGTCCGTGGTTCCGATCTGGTGGTGACCGCTATTGCCGAAGGGCGTAAGGCGGCAGACGGTATTATGAACTGGCTGGAAGTTTAAGCGAGGTAACAATGAATTCGCCAGACTGTGTCTGGCGAATc                                                                                                                                         <  1:984261/150‑1 (MQ=255)
             gAAAATCTTTGCTGGCGGCGATATCGTCCGTGGTTCCGAGCTGGGGGTGACCGCTATTGCCGAAGGTCGTAAGGCGGCAGACGGTATTATGAACTGGCTGGAAGTTTAAGCGAGGTAACAATGAATTCGCCAGACTGTGTCTGGCGAATc                                                                                                                                         <  2:1534659/150‑1 (MQ=255)
             caaatTCTATGCTGGCGGCGATATCGTCCGTGGTTCCGATCTGGTGGTGCCCGCTATTGCCGACGGTCGTAAGCGGGCAGACGGTAATATGAACTGGCTGGCAGTTTAAGCGAGGTAACAATGAATTCGCCAGACACAGTCTGGCGAAtt                                                                                                                                         <  2:1819015/149‑2 (MQ=255)
                  tCTTTGCTGGCGGCGATATCGTCCGTGGTTCCGATCTGGTGGTGACCGCTATTGCCGAAGGGCGTAAGGCGGCAGACGGTAATATGAACTGGCTGGAAGTGTAAGCGAGGTAACAATGAATTCGCCAGACACAGCCTGGCGAAttcattg                                                                                                                                    <  1:276888/150‑7 (MQ=255)
                   cTTTGCTGGCGGCGATATCGTCCGTGGTTCCGATCTGGTGGTGACCGCTATTGCCGAAGGTCGTAAGGCGGCAGACGGTATTATGAACTGGCTGGAAGTTTAAGCGAGGTAACAATGAATTCGCCAGACTGTGTCTGGCGAATCAAATTa                                                                                                                                   >  1:901978/1‑150 (MQ=255)
                     ttGCTGGCGGCGATATCGTCCGTGGTTCCGATCTGGTGGTGACCGCTATTGCCGAAGGTCGTAAGGCGGCAGACGGTATTATGAACTGGCTGGAAGTTTAAGCGAGGTAACAATGAATTCGCCAGACTGTGTCTGGCGAATCAAATTACt                                                                                                                                 <  2:1369547/150‑1 (MQ=255)
                      tGCTGGCGGAGATATCGTCCGTGGTTCCGATCTGGTGGTGACCGCTATTGCCGAAGGTCGTAAGGCGGCAGACGGTATTATGAACTGGCTGGAAGTTTAAGCGAGGTAACAATGAATTCGCCAGACTGTGTCTGGCGAATCAAATTACTg                                                                                                                                <  2:1428086/150‑1 (MQ=255)
                        cTGGCGGCGATATCGTCCGTGGTTCCGATCTGGTGGTGACCGCTATTGCCGAAGGTCGTAAGGCGGCAGACGGTATTATGAACTGGCTGGAAGTTTAAGCGAGGTAACAATGAATTCGCCAGACTGTGTCTGGCGAATCAAATTACTGAc                                                                                                                              <  1:758125/150‑1 (MQ=255)
                        cTGGCGGCGATATCGTCCGTGGTTCCGATCTGGTGGTGACCGCTATTGCCGAAGGTCGTAAGGCGGCAGACGGTATTATGAACTGGCTGGAAGTTTAAGCGAGGTAACAATGAATTCGCCAGACTGTGTCTGGCGAATCAAATTACTGAc                                                                                                                              <  1:27539/150‑1 (MQ=255)
                        cTGGCGGCGATATCGTCCGTGGTTCCGATCTGGTGGTGACCGCTATTGCCGAAGGTCGTAAGGCGGCAGACGGTATTATGAACTGGCTGGAAGTTTAAGCGAGGTAACAATGAATTCGCCAGACTGTGTCTGGCGAATCAAATTACTGAc                                                                                                                              <  1:1074671/150‑1 (MQ=255)
                        cTGGCGGCGATATCGTCCGTGGTTCCGATCTGGTGGTGACCGCTATTGCCGAAGGTCGTAAGGCGGCAGACGGTATTATGAACTGGCTGGAAGTTTAAGCGAGGTAACAATGAATTCGCCAGACTGTGTCTGGCGAATCAAATTACTGAc                                                                                                                              >  1:1649026/1‑150 (MQ=255)
                        cTGGCGGCGATATCGGCCGTGGTTCCGCGCTGGTGGTGACCGCTATTGCCGAAGGTCGTAAGGCGGCAGACTGTATTATGAACTGTCTTGAATTTAAAGCGTGGTAACAATGAATTCGCGAGACTGTGTCTGGCGAATCATATTACTGAc                                                                                                                              <  1:291210/150‑1 (MQ=255)
                          ggcggcGATATCGTCCGTGGTTCCGATCTGGTGGTGACCGCTATTGCCGAAGGTCGTAAGGCGGCAGACGGTATTATGAACTGGCTGGAAGTTTAAGCGAGGTAACAATGAATTCGCCAGACTGTGTCTGGCGAATCAAATTACTGACGa                                                                                                                            >  1:940783/1‑150 (MQ=255)
                          ggcggcGATATCGTCCGTGGTTCCGATCTGGTGGTGACCGCTATTGCCGAAGGTCGTAAGGCGGCAGACGGTATTATGAACTGGCTGGAAGTTTAAGCGAGGTAACAATGAATTCGCCAGACTGTGTCTGGCGAATCAAATTACTGACGa                                                                                                                            >  2:110067/1‑150 (MQ=255)
                              gcgATATCGTCCGTGGTTCCGATCTGGTGGTGACCGCTATTGCCGAAGGTCGTAAGGCGGCAGACGGTATTATGAACTGGCTGGAAGTTTAAGCGAGGTAACAATGAATTCGCCAGACTGTGTCTGGCGAATCAAATTACTGACGAAGTg                                                                                                                        <  1:642988/150‑1 (MQ=255)
                              gcgATATCGTCCGTGGTTCCGATCTGGTGGTGACCGCTATTGCCGAAGGTCGTAAGGCGGCAGACGGTATTATGAACTGGCTGGAAGTTTAAGCGAGGTAACAATGAATTCGCCAGACTGTGTCTGGCGAATCAAATTACTGACGAAGTg                                                                                                                        <  1:210030/150‑1 (MQ=255)
                                gttatCGTCCGTGGTTCCGATCTGGTGGTGACCGCTATTGCCGAAGGTCGTAAGGCGGCAGACGGTATTATGAACTGGCTGGAAGTTTAAGCGAGGTAACAATGAATTCGCCAGACTGTGTCTGGCGAATCAAATTACTGACGAAGTgcg                                                                                                                      <  1:880119/148‑1 (MQ=255)
                                      gTCCGTGGTTCCGATCTGGTGGTGACCGCTATTGCCGAAGGTCGTAAGGCGGCAGACGGTATTATGAACTGGCTGGAAGTTTAAGCGAGGTAACAATGAATTCGCCAGACTGTGTCTGGCGAATCAAATTACTGACGAAGTGCGCgttgt                                                                                                                <  1:1311631/150‑1 (MQ=255)
                                      gTCCGTGGTTCCGATCTGGTGGTGACCGCTATTGCCGAAGGTCGTAAGGCGGCAGACGGTATTATGAACTGGCTGGAAGTTTAAGCGAGGTAACAATGAATTCGCCAGACTGTGTCTGGCGAATCAAATTACTGACGAAGTGCGCgttgt                                                                                                                <  1:1635974/150‑1 (MQ=255)
                                         cGTGGTTCCGATCTGGTGGTGACCGCTATTGCCGAAGGTCGTAAGGCGGCAGACGGTATTATGAACTGGCTGGAAGTTTAAGCGAGGTAACAATGAATTCGCCAGACTGTGTCTGGCGAATCAAATTACTGACGAAGTGCGCGTTGtttt                                                                                                             <  1:1200633/150‑1 (MQ=255)
                                         cGTGGTTCCGATCTGGTGGTGACCGCTATTGCCGAAGGTCGTAAGGCGGCAGACGGTATTATGAACTGGCTGGAAGTTTAAGCGAGGTAACAATGAATTCGCCAGACTGTGTCTGGCGAATCAAATTACTGACGAAGTGCGCGTTGtttt                                                                                                             >  1:156114/1‑150 (MQ=255)
                                           tGGTTCCGATCTGGTGGTGACCGCTATTGCCGAAGGTCGTAAGGCGGCAGACGGTATTATGAACTGGCTGGAAGTTTAAGCGAGGTAACAATGAATTCGCCAGACTGTGTCTGGCGAATCAAATTACTGACGAAGTGCGCGTTGTTTTCt                                                                                                           <  1:1155678/150‑1 (MQ=255)
                                           tGGTTCCGATCTGGTGGTGACCGCTATTGCCGAAGGTCGTAAGGCGGAAGACGGTATTATGAACTGGCTGGAAGTTTAAGCGAGGTAACATTGAATTCGCCAGACTGTGTCTGGCGAATCAAATTACTGACGAAGTGCGCGTTGTTTTCt                                                                                                           <  2:1271741/150‑1 (MQ=255)
                                               tCCGATCTGGTGGTGACCGCTATTGCCGAAGGTCGTAAGGCGGCAGACGGTATTATGAACTGGCTGGAAGTTTAAGCGAGGTAACAATGAATTCGCCAGACTGTGTCTGGCGAATCAAATTACTGACGAAGTGCGCGTTGTTTTCTGGCa                                                                                                       >  1:968330/1‑150 (MQ=255)
                                                 cGATCTGGTGGTGACCGCTATTGCCGAAGGTCGTAAGGCGGCAGACGGTATTATGAACTGGCTGGAAGTTTAAGCGAGGTAACAATGAATTCGCCAGACTGTGTCTGGCGAATCAAATTACTGACGAAGTGCGCGTTGTTTTCTGGCaaa                                                                                                     >  2:180466/1‑150 (MQ=255)
                                                 cGATCTGGTGGTGACCGCTATTGCCGAAGGTCGTAAGGCGGCAGACGGTATTATGAACTGGCTGGAAGTTTAAGCGAGGTAACAATGAATTCGCCAGACTGTGTCTGGCGAATCAAATTACTGACGAAGTGCGCGTTGTTTTCTGGCaaa                                                                                                     >  1:1386688/1‑150 (MQ=255)
                                                      tggtggTGACCGCTATTGCCGAAGGTCGTAAGGCGGCAGACGGTATTATGAACTGGCTGGAAGTTTAAGCGAGGTAACAATGAATTCGCCAGACTGTGTCTGGCGAATCAAATTACTGACGAAGTGCGCGTTGTTTTCTGGCAAATGCgg                                                                                                >  2:615857/1‑150 (MQ=255)
                                                      tggtggTGACCGCTATTGCCGAAGGTCGTAAGGCGGCAGACGGTATTATGAACTGGCTGGAAGTTTAAGCGAGGTAACAATGAATTCGCCAGACTGTGTCTGGCGAATCAAATTACTGACGAAGTGCGCGTTGTTTTCTGGCAAATGCgg                                                                                                >  2:904865/1‑150 (MQ=255)
                                                       ggtggtGACCGCTATTGCCGAAGGTCGTAAGGCGGCAGACGGTATTATGAACTGGCTGGAAGTTTAAGCGAGGTAACAATGAATTCGCCAGACTGTGTCTGGCGAATCAAATTACTGACGAAGTGCGCGTTGTTTTCTGGCAAATGCGGc                                                                                               <  1:1022932/150‑1 (MQ=255)
                                                       ggtggtGACCGCTATTGCCGAAGGTCGTAAGGCGGCAGACGGTATTATGAACTGGCTGGAAGTTTAAGCGAGGTAACAATGAATTCGCCAGACTGTGTCTGGCGAATCAAATTACTGACGAAGTGCGCGTTGTTTTCTGGCAAATGCGGc                                                                                               >  1:1609843/1‑150 (MQ=255)
                                                        gtggtgACCGCTATTGCCTAAGGTCGTAAGGCGGCAGACGGTATTATGAACTGGCTGGAAGTTTAAGCGAGGTAACAATGAATTCGCCAGACTGTGTCTGGCGAATCAAATTACTGACGAAGTGCGCGTTGTTTTCTGGCAAATGCGGCa                                                                                              <  2:727897/150‑1 (MQ=255)
                                                          ggtgACCGCTATTGCCGAAGGTCGTAAGGCGGCAGACGGTATTATGAACTGGCTGGAAGTTTAAGCGAGGTAACAATGAATTCGCCAGACTGTGTCTGGCGAATCAAATTACTGACGAAGTGCGCGTTGTTTTCTGGCAAATGCGGCATg                                                                                            <  2:1601233/150‑1 (MQ=255)
                                                                 gCTATTGCCGAAGGTCGTAAGGCGGCAGACGGTATTCTGAACTGGCTGGAAGTTTAAGCGAGGTAACAATGAATTCGCCAGACTGTGTCTGGCGAATCAAATTACTGACGAAGTGCGCGTTGTTTTCTGGCAAATGCGGCATGAGCGttt                                                                                     <  2:915396/150‑1 (MQ=255)
                                                                 gCTATTGCCGAAGGTCGTAAGGCGGCAGACGGTATTATGAACTGGCTGGAAGTTTAAGCGAGGTAACAATGAATTCGCCAGACTGTGTCTGGCGAATCAAATTACTGACGAAGTGCGCGTTGTTTTCTGGCAAATGCGGCATGAGCGttt                                                                                     <  1:581422/150‑1 (MQ=255)
                                                                 gCTATTGCCGAAGGTCGTAAGGCGGCAGACGGTATTATGAACTGGCTGGAAGTTTAAGCGAGGTAACAATGAATTCGCCAGACTGTGTCTGGCGAATCAAATTACTGACGAAGTGCGCGTTGTTTTCTGGCAAATGCGGCATGAGCGttt                                                                                     >  2:885/1‑150 (MQ=255)
                                                                  cTATTGCCGAAGGTCGTAAGGCGGCAGACGGTATTATGAACTGGCTGGAAGTTTAAGCGAGGTAACAATGAATTCGCCAGACTGTGTCTGGCGAATCAAATTACTGACGAAGTGCGCGTTGTTTTCTGGCAAATGCGGCATGAGCGTTTg                                                                                    <  1:839486/150‑1 (MQ=255)
                                                                  cTATTGCCGAAGGTCGTAAGGCGGCAGACGGTATTATGAACTGGCTGGAAGTTTAAGCGAGGTAACAATGAATTCGCCAGACTGTGTCTGGCGAATCAAATTACTGACGAAGTGCGCGTTGTTTTCTGGCAAATGCGGCATGAGCGTTTg                                                                                    >  2:479553/1‑150 (MQ=255)
                                                                       gCCGAAGGTCGTAAGGCGGCAGACGGTATTATGAACTGGCTGGAAGTTTAAGCGAGGTAACAATGAATTCGCCAGACTGTGTCTGGCGAATCAAATTACTGACGAAGTGCGCGTTGTTTTCTGGCAAATGCGGCATGAGCGTTTGACCtt                                                                               >  2:1632269/1‑150 (MQ=255)
                                                                       gCCGAAGGTCGTAAGGCGGCAGACGGTATTATGAACTGGCTGGAAGTTTAAGCGAGGTAACAATGAATTCGCCAGACTGTGTCTGGCGAATCAAATTACTGACGAAGTGCGCGTTGTTTTCTGGCAAATGCGGCATGAGCGTTTGACCtt                                                                               >  1:339121/1‑150 (MQ=255)
                                                                        ccGAAGGTCGTAAGGCGGCAGACGGTATTATGAACTGGCTGGAAGTTTAAGCGAGGTAACAATGAATTCGCCAGACTGTGTCTGGCGAATCAAATTACTGACGAAGTGCGCGTTGTTTTCTGGCAAATGCGGCATGAGCGTTTGACCttt                                                                              <  1:1810502/150‑1 (MQ=255)
                                                                         cGAAGGTCGTAAGGCGGCAGACGGTATTATGAACTGGCTGGAAGTTTAAGCGAGGTAACAATGAATTCGCCAGACAGTGTCTGGCGAATCAAATTACTGACGAAGTGCGCGTTGTTTTCTGGCAAATGCGGCATGCGCGTTTGACCTTTc                                                                             >  2:344784/1‑150 (MQ=255)
                                                                            aGGTCGTAAGGCGGCAGACGGTATTATGAACTGGCTGGAAGTTTAAGCGAGGTAACAATGAATTCGCCAGACTGTGTCTGGCGAATCAAATTACTGACGAAGTGCGCGTTGTTTTCTGGCAAATGCGGCATGAGCGTTTGACCTTTCCta                                                                          <  2:1559625/150‑1 (MQ=255)
                                                                             ggTCGTAAGGGGGCAGACGGTATTATGAACTGGCTGGAAGTTTAAGCGAGGTAACAATGAATTCGGCAGACTGTGTCTGGCGAATCTAATTACTGAATAAGTGCTCGTTGTTTTCTTGCAAATGCGGCATGAGCGTTTGACCTTTCCtat                                                                         <  1:480712/150‑1 (MQ=255)
                                                                              gTCGTAAGGCGGCAGACGGTATTATGAACTGGCTGGAAGTTTAAGCGAGGTAACAATGAATTCGCCAGACTGTGTCTGGCGAATCAAATTACTGACGAAGTGCGGGTTGTTTTCTGGCAAATGCGGCATGAGCGTTTGCCCTTTCCtata                                                                        >  2:754166/1‑150 (MQ=255)
                                                                                  tAAGGCGGCAGACGGTATTATGAACTGGCTGGAAGTTTAAGCGAGGTAACAATGAATTCGCCAGACTGTGTCTGGCGAATCAAATTACTGACGAAGTGCGCGTTGTTTTCTGGCAAATGCGGCATGAGCGTTTGACCTTTCCTATAATAg                                                                    <  2:730370/150‑1 (MQ=255)
                                                                                  tAAGGCGGCAGACGGTATTATGAACTGGCTGGAAGTTTAAGCGAGGTAACAATGAATTCGCCAGACTGTGTCTGGCGAATCAAATTACTGACGAAGTGCGCGTTGTTTTCTGGCAAATGCGGCATGAGCGTTTGACCTTTCCTATAATAg                                                                    >  2:422135/1‑150 (MQ=255)
                                                                                   aaGGCGGCAGACGGTATTATGAACTGGCTGGAAGTTTAAGCGAGGTAACAATGAATTCGCCAGACTGTGTCTGGCGAATCAAATTACTGACGAAGTGCGCGTTGTTTTCTGGCAAATGCGGCATGAGCGTTTGACCTTTCCTATAATAgg                                                                   >  2:1540204/1‑150 (MQ=255)
                                                                                   aaGGCGGCAGACGGTATTATGAACTGGCTGGAAGTTTAAGCGAGGTAACAATGAATTCGCCAGACTGTGTCTGGCGAATCAAATTACTGACGAAGTGCGCGTTGTTTTCTGGCAAATGCGGCATGAGCGTTTGACCTTTCCTATAATAgg                                                                   <  1:593070/150‑1 (MQ=255)
                                                                                     ggcggcAGACGGTATTATGAACTGGCTGGAAGTTTAAGCGAGGTAACAATGAATTCGCCAGACTGTGTCTGGCGAATCAAATTACTGACGAAGTGCGCGTTGTTTTCTGTCAAATGCTGAATGAGCGTTTGACCTTTCCTATAATAGGGa                                                                 >  2:281651/1‑150 (MQ=255)
                                                                                     ggcggcAGACGGTATTATGAACTGGCTGGAAGTTTAAGCGAGGTAACAATGAATTCGCCAGACTGTGTCTGGCGAATCAAATTACTGACGAAGTGCGCGTTGTTTTCTGGCAAATGCGGCATGAGCGGTTGACCTTTCCTATAATAGGGa                                                                 <  2:1312532/150‑1 (MQ=255)
                                                                                       cggcAGACGGTATTATGAACTGGCTGGAAGTTTAAGCGAGGTAACAATGAATTCGCCAGACTGTGTCTGGCGAACCAAATTACTGACGAAGTGCGCGTTGTTTTCTGGCAAATGCGGCATGAGCGTTTGACCTTTCCTATAATAGGGATa                                                               >  2:1237807/1‑150 (MQ=255)
                                                                                             aCGGTATTATGAACTGGCTGGAAGTTTAAGCGAGGTAACAATGAATTCGCCAGACTGTGTCTGGCGAATCAAATTACTGACGAAGTGCGCGTTGTTTTCTGGCAAATGCGGCATGAGCGTTTGACCTTTCCTATAATAGGGATAATTTTc                                                         <  1:41569/150‑1 (MQ=255)
                                                                                             aCGGTATTATGAACTGGCTGGAAGTTTAAGCGAGGTAACAATGAAATCGCCAGACTGTGTCTGGCGAATCAAATTACTGACGAAGTGCGCGTTGTTTTCTGGCAAATGCGGCATGAGCGTTTGACCTTTCCTATAATAGGGATAATTTTc                                                         <  1:946663/150‑1 (MQ=255)
                                                                                              cGGTATTATGAACTGGCTGGAAGTTTAAGCGAGGTAACAATGAATTCGCCAGACTGTGTCTGGCGAATCAAATTACTGACGAAGTGCGCGTTGTTTTCTGGCAAATGCGGCATGAGCGTTTGACCTTTCCTATAATAGGGATAATTTTCa                                                        <  2:1723500/150‑1 (MQ=255)
                                                                                              cGGTATTATGAACTGGCTGGAAGTTTAAGCGAGGTAACAATGAATTCGCCAGACTGTGTCTGGCGAATCAAATTACTGACGAAGTGCGCGTTGTTTTCTGGCAAATGCGGCATGAGCGTTTGACCTTTCCTATAATAGGGATAATTTTCa                                                        >  1:1721017/1‑150 (MQ=255)
                                                                                               ggTATTATGAACTGGCTGGAAGTTTAAGCGAGGTAACAATGAATTCGCCAGACTGTGTCTGGCGAATCAAATTACTGACGAAGTGCGCGTTGTTTTCTGGCAAATGCGGCATGAGCGTTTGACCTTTCCTATAATAGGGATAATTTTCaa                                                       >  2:544952/1‑150 (MQ=255)
                                                                                                      tGAACTGGCTGGAAGTTTAAGCGAGGTAACAATGAATTCGCCAGACTGTGTCTTGCGAATCAAATTACTGACGAAGTGCGCGTTGTTTTCTGGCAAATGCGGCATGAGCGTTTGACCTTTCCTATAATAGGGATAATTTTCAATAAATtg                                                <  1:1574146/150‑1 (MQ=255)
                                                                                                        aaCTGGCTGGAAGTTTAAGCGAGGTAACAATGAATTCGCCAGACTGTGTCTGGCGAATCAAATTACTGACGAAGTGCGCGTTGTTTTCTGGCAAATGCGGCATGAGCGTTTGACCTTTCCTATAATAGGGATAATTTTCAATAAATtgtg                                              <  1:885/150‑1 (MQ=255)
                                                                                                            ggctggAAGTTTAAGCGAGGTAACAATGAATTCGCCAGACTGTGTCTGGCGAATCAAATTACTGACGAAGTGCGCGTTGTTTTCTGGCAAATGCGGCATGAGCGTTTGACCTTTCCTATAATAGGGATAATTTTCAATAAATTGTGCttt                                          >  2:1219134/1‑150 (MQ=255)
                                                                                                             gctggAAGTTTTAGGGATGTTACAATGAATTCGCCAGTCTGTTTCTGGCGAATCAATTTACTGACGAAGTGCGCTTTGTTTTCTGGCTAATGCGGCATGAGCGTTTGACCTTTCCTATATTAGGGATAATTTTCAATAAATTGTGCTTTg                                         <  1:1804583/150‑1 (MQ=255)
                                                                                                              ctggAAGTTTAAGCGAGGTAACAATGAATTCGCCAGACTGTGTCTGGCGAATCAAATTACTGACGAAGTGCGCGTTGTTTTCTGGCAAATGCGGCATGAGCGTTTGACCTTTCCTATAATAGGGATAATTTTCAATAAATTGTGCTTTGt                                        >  1:1676102/1‑150 (MQ=255)
                                                                                                                ggAAGTTTAAGCGAGGTAACAATGAATTCGCCAGACTGTGTCTGGCGAATCAAATTACTGACGAAGTGCGCGTTGTTTTCTGGCAAATGCGGCATGAGCGTTTGACCTTTCCTATAATAGGGATAATTTTCAATAAATTGTGCTTTGtta                                      >  1:1484215/1‑149 (MQ=255)
                                                                                                                ggAAGTTTAAGCGAGGTAACAATGAATTCGCCAGACTGTGTCTGGCGAATCAAATTACTGACGAAGTGCGCGTTGTTTTCTGGCAAATGCGGCATGAGCGTTTGACCTTTCCTATAATAGGGATAATTTTCAATAAATTGTGCTTTGTTc                                      <  1:615857/150‑1 (MQ=255)
                                                                                                                 gAAGTTTAAGCGAGGTAACAATGAATTCGCCAGACTGTGTCTGGCGAATCAAATTACTGACGAAGTGCGCGTTGTTTTCTGGCAAATGCGGCATGAGCGTTTGACCTTTCCTATAATAGGGATAATTTTCAATAAATTGTGCTTTGTTCa                                     >  2:1569242/1‑150 (MQ=255)
                                                                                                                  aaGTTTAAGCGAGGTAACAATGAATTCGCCAGACTGTGTCTGGCGAATCAAATTACTGACGAAGTGCGCGTTGTTTTCTGGCAAATGCGGCATGAGCGTTTGACCTTTCCTATAATAGGGATAATTTTCAATAAATTGTGCTTTGTTCAc                                    >  2:549645/1‑150 (MQ=255)
                                                                                                                   aGTTTAAGCGAGGTAACAATGAATTCGCCAGACTGTGTCTGGCGAATCAAATTACTGACGAAGTGCGCGTTGTTTTCTGGCAAATGCGGCATGAGCGTTTGACCTTTCCTATAATAGGGATAATTTTCAATAAATTGTGCTTTGTTCAcc                                   <  1:614464/150‑1 (MQ=255)
                                                                                                                   aGTTTAAGCGAGGTAACAATGAATTCGCCAGACTGTGTCTGGCGAATCAAATTACTGACGAAGTGCGCGTTGTTTTCTGGCAAATGCGGCATCACCGTTTGCCCTTTCCAATAAAAGGTATAATTTTCAATAAATTGTGATTTGTTCAca                                   >  2:276004/1‑149 (MQ=255)
                                                                                                                   aGTTTAAGCGAGGTAAAAATGAATTCGCCAGACTGTGTCTGGCGAATCAAATTACTGACGAAGTGCGCGTTGTTTTCTGGCAAATGCGGCATGAGCGTTTGACCTTTCCTATAATAGGGATAATTTTCAATAAATTGTGCTTTGTTCAcc                                   <  1:23247/150‑1 (MQ=255)
                                                                                                                    gTTTAAGCGAGGTAACAATGAATTCGCCAGACTGTGTCTGGCGAATCAAATTACTGACGAAGTGCGCGTTGTTTTCTGGCAAATGCGGCATGAGCGTTTGACCTTTCCTATAATAGGGATAATTTTCAATAAATTGTGCTTTGTTCACCa                                  >  1:516569/1‑150 (MQ=255)
                                                                                                                         aGCGAGGTAACAATGAATTCGCCAGACTGTGTCTGGCGAATCAAATTACTGACGAAGTGCGCGTTGTTTTCTGGCAAATGCGGCATGAGCGTTTGACCTTTCCTATAATAGGGATAATTTTCAATAAATTGTGCTTTGTTCACCAAGTCt                             <  1:904865/150‑1 (MQ=255)
                                                                                                                            gAGGTAACAATGAATTCGCCAGACTGTGTCTGGCGAATCAAATTACTGACGAAGTGCGCGTTGTTTTCTGGCAAATGCGGCATGAGCGTTTGACCTTTCCTATAATAGGGATAATTTTCAATAAATTGTGCTTTGTTCACCAAGTCtgat                          <  2:1752545/150‑1 (MQ=255)
                                                                                                                            gAGGTAACAATGAATTCGCCAGACTGTGTCTGGCGAATCAAATTAATGACGAAGTGCGCGTTGTTTTCTGGCAAATGCGGCATGAGCGTTTTACCTTTACTATAAAAGGGATAATTTTAAATAAATTGTGCTTTTTTCACCAAGTCtgat                          >  2:1388382/1‑150 (MQ=255)
                                                                                                                                  aCAATGAATTCGCCAGACTGTGTCTGGCGAATCAAATTACTGACGAAGTGCGCGTTGTTTTCTGGCAAATGCGGCATGAGCGTTTGACCTTTCCTATAATAGGGATAATTTTCAATAAATTGTGCTTTGTTCACCAAGTCTGATGAGTGa                    >  1:775479/1‑150 (MQ=255)
                                                                                                                                   cAATGAATTCGCCAGACTGTGTCTGGCGAATCAAATTACTGACGAAGTGCGCGTTGTTTTCTGGCAAATGCGGCATGAGCGTTTGACCTTTCCTATAATTGGGATAATTTTCAATAAATTGTGCTTTGTTCACCAAGTCTGATGAGTGaa                   >  1:1149599/1‑150 (MQ=255)
                                                                                                                                    aaTGAATTCGCCAGACTGTGTCTGGCGAATCAAATTTCTGACGAAGTGCGCGTTGTTTTCTGGCAAATGCAGCATGAGCGTTTGACCTTTCCTATAATAGGGATAATTTTCAATAAATTGTGCTTTGTTCACCAAGTCTGATGAGTGAAg                  >  1:1683577/1‑150 (MQ=255)
                                                                                                                                    aaTGAATTCGCCAGACTGTGTCTGGCGAATCAAATTACTGACGAAGTGCGCGTTGTTTTCTGGCAAATGCGGCATGAGCGTTTGACCTTTCCTATAATAGGGATAATTTTCAATAAATTGTGCTTTGTTCACCAAGTCTGATGAGTGAAg                  >  2:814/1‑150 (MQ=255)
                                                                                                                                    aaTGAATTCGCCAGACTGTGTCTGGCGAATCAAATTACTGACGAAGTGCGCGTTGTTTTCTGGCAAATGCGGCATGAGCGTTTGACCTTTCCTATAATAGGGATAATTTTCAATAAATTGTGCTTTGTTCACCAAGTCTGATGAGTGAAg                  >  1:1366358/1‑150 (MQ=255)
                                                                                                                                      tGAATTCGCCAGACTGTGTCTGGCGAATCAAATTACTGACGAAGTGCGCGTTGTTTTCTGGCAAATGCGGCTTGAGCGTTTGACCTTTCCTATAATAGGGATAATTTTCAATAAATTGTGCTTTGTTCACCACGTCTGATGAGTGAAGca                >  1:985131/1‑150 (MQ=255)
                                                                                                                                        aaTTCGCCAGACTGTGTCTGGCGAATCAAATTACTGACGAAGTGCGCGTTGTTTTCTGGCAAATGCGGCATGAGCGTTTGACCTTTCCTATAATAGGGATAATTTTCAATAAATTGTGCTTTGTTCACCAAGTCTGATGAGTGAAGcaca              <  2:1626668/150‑1 (MQ=255)
                                                                                                                                        aaTTCGCCAGACTGTGTCTGGCGAATCAAATTACTGACGAAGTGCGCGTTGTTTTCTGGCAAATGCGGCATGAGCGTTTGACCTTTCCTATAATAGGGATAATTTTCAATAAATTGTGCTTTGTTCACCAAGTCTGATGAGTGAAGcaca              <  1:110067/150‑1 (MQ=255)
                                                                                                                                        aaTTCGCCAGACTGTGTCTGGCGAATCAAATTACTGACGAAGTGCGCGTTGTTTTCTGGCAAATGCGGCATGAGCGTTTGACCTTTCCTATAATAGGGATAATTTTCAATAAATTGTGCTTTGTTCACCAAGTCTGATGAGTGAAGcaca              >  1:1798293/1‑150 (MQ=255)
                                                                                                                                         aTTCGCCAGACACAGTCTGGCGAATCAAATTACTGACGAAGTGCGCGTTGTTTTCTGGCAAATGCGGCATGAGCGTTTGACCTTTCCTATAATAGGGATAATTTTCAATAAATTGTGCTTTGTTCACCAAGTCTGATGAGTGAAGCACAt             >  1:1301172/1‑150 (MQ=255)
                                                                                                                                          ttCGCCAGACACAGTCTGGCGAATCAAATTACTGACGAAGTGCGCGTTGTTTTCTGGCAAATGCGGCATGAGCGTTTGACCTTTCCTATAATAGGGATAATTTTCAATAAATTGTGCTTTGTTCACCAAGTCTGATGAGTGAAGCACAta            >  2:580189/1‑149 (MQ=255)
                                                                                                                                          ttCGCCAGACACAGTCTGGCGAATCAAATTACTGACGAAGTGCGCGTTGTTTTCTGGCAAATGCGGCATGAGCGTTTGACCTTTCCTATAATAGGGATAATTTTCAATAAATTGTGCTTTGTTCACCAAGTCTGATGAGTGAAGCACATc            >  1:721287/1‑150 (MQ=255)
                                                                                                                                          ttCGCCAGACACAGTCTGGCGAATCAAATTACTGACGAAGTGCGCGTTGTTTTCTGGCAAATGCGGCATGAGCGTTTGACCTTTCCTATAATAGGGATAATTTTCAATAAATTGTGCTTTGTTCACCAAGTCTGATGAGTGAAGCACATc            >  1:575796/1‑150 (MQ=255)
                                                                                                                                          ttCGACAGACTGTGTCTGGCGAATCAAATTACTGACGAAGTGCGCGTTGTTTTCTGGCAAATGCGGCATGAGCGTTTGACCTTTCCTATAATAGGGATAATTTTCAATAAATTGTGCTTTGTTCACCAAGTCTGATGAGTGAAGCACATc            <  1:647939/150‑1 (MQ=255)
                                                                                                                                           tCGCCAGACTGTGTCTGGCGAATCAAATTACTGACGAAGTGCGCGTTGTTTTCTGGCAAATGCGGCATGAGCGTTTGACCTTTCCTATAATAGGGATAATTTTCAATAAATTGTGCTTTGTTCACCAAGTCTGATGAGTGAAGCACATCa           <  2:771619/150‑1 (MQ=255)
                                                                                                                                            cGCCAGACTGTGTCTGGCGAATCAAATTACTGACGAAGTGCGCGTTGTTTTCTGGCAAATGCGGCATGAGCGTTTGACCTTTCCTATAATAGGGATAATTTTCAATAAATTGTGCTTTGTTCACCAAGTCTGATGAGTGAAGCACATCaa          >  1:1226431/1‑150 (MQ=255)
                                                                                                                                            cGCCAGACACAGTCTGGCGAATCAAATTACTGACGAAGTGCGCGTTGTTTTCTGGCAAATGCGGCATGAGCGTTTGACCTTTCCTATATTAGGGATAATTTTCAATAAATTGTGCTTTGTTCACCAAGTCTGATGAGTGAAGCACATCaa          >  2:1462503/1‑150 (MQ=255)
                                                                                                                                             gCCAGACTGTGTCTGGCGAATCAAATTACTGACGAAGTGCGCGTTGTTTTCTGGCAAATGCGGCATGAGCGTTTGACCTTTCCTATAATAGGGATAATTTTCAATAAATTGTGCTTTGTTCACCAAGTCTGATGAGTGAAGCACATCAAg         >  2:271498/1‑150 (MQ=255)
                                                                                                                                              ccAGACTGTGTCTGGCGAATCAAATTACTGACGAAGTGCGCGTTGTTTTCTGGCAAATGCGGCATGAGCGTTTGACCTTTCCTATAATAGGGATAATTTTCAATAAATTGTGCTTTGTTCACCAAGTCTGATGAGTGAAGCACATCAAGc        <  1:1473171/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: