Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,585,063 C→G 5.4% intergenic (+20/+18) yrhB → / ← ggt stable heat shock chaperone/gamma‑glutamyltranspeptidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,585,0630CG 5.4% 310.7 / 7.2 106intergenic (+20/+18)yrhB/ggtstable heat shock chaperone/gamma‑glutamyltranspeptidase
Reads supporting (aligned to +/- strand):  ref base C (50/49);  new base G (3/2);  total (55/51)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.87e-01

Only 100 of 106 total aligned reads displayed.

TTAGCTGGTAACGCACCTTTTATTATTGATAAAGACAGTGGTGAAATTCATTCTCTGGGAACGGCAAAACCGCTGGAAGAATATCTACAGGATTACGAAATAAAAAAGGCTACCTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTTAGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGCCACTTTCTGACCTTTTGCTT  >  NC_000913/3584915‑3585211
                                                                                                                                                    |                                                                                                                                                    
ttAGCTGGTAACGCACATTTTATATTTTAAAAAGACAGTGGTGAAATTCATTCTCTGGGAACGGCAAAACCGCTGGAAGAATATCAACAGGATTACGAAATAAAAAAGGCTACCTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACg                                                                                                                                                     <  2:711979/150‑1 (MQ=255)
   gCTGGTAACGCACCTTTTATTATTGATAAAGACAGTGGTGAAATTCATTCTCTGGGAACGGCAAAACCGCTGGAAGAATATCTACAGGATTACGAAATAAAAAAGGCTACCTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAag                                                                                                                                                  >  1:236646/1‑150 (MQ=255)
   gCTGGTAACGCACCTTTTATTATTGATAAAGACAGTGGTGAAATTCATTCTCTGGGAACGGCAAAACCGCTGGAAGAATATCTACAGGATTACGAAATAAAAAAGGCTACCTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAag                                                                                                                                                  <  2:436072/150‑1 (MQ=255)
       gTAACGCACCTTTTATTATTGATAAAGACAGTGGTGAAATTCATTCTCTGGGAACGGCAAAACCGCTGGAAGAATATCTACAGGATTACGAAATAAAAAAGGCTACCTTCGGCTTGCCCTGACAAAATAGCCCTCTTCGTGGGAAGAggg                                                                                                                                              <  1:1702224/150‑1 (MQ=255)
       gTAACGCACCTTTTATTATTGATAAAGACAGTGGTGAAATTCATTCTCTGGGAACGGCAAAACCGCTGGAAGAATATCTACAGGATTACGAAATAAAAAAGGCTACCTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAggg                                                                                                                                              <  2:1670434/150‑1 (MQ=255)
         aaCGCACCTTTTATTATTTATAAAGACAGTGGTGAAATTCATTCTCTGGGAACGGCAAAACCGCTGGAAGAATATCTACAGGATTACGAAATAAAAAAGGCTACCTTCGGCTTGCCCTGACAAAATAGCCCTCTTCGTGGGAAGAGGGct                                                                                                                                            <  1:1389895/150‑2 (MQ=255)
           cGCACCTTTTATTATTGATAAAGACAGTGGTGAAATTCATTCTCTGGGAACGGCAAAACCGCTGGAAGAATATCTACAGGATTACGAAATAAAAAAGGCTACCTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGc                                                                                                                                          >  2:466817/1‑150 (MQ=255)
           cGCACCTTTTATTATTGATAAAGACAGTGGTGAAATTCATTCTCTGGGAACGGCAAAACCGCTGGAAGAATATCTACAGGATTACGAAATAAAAAAGGCTACCTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGc                                                                                                                                          >  1:1275039/1‑150 (MQ=255)
            gCACCTTTTATTATTGATAAAGACAGTGGTGAAATTCATTCTCTGGGAACGGCAAAACCGCTGGAAGAATATCTACAGGATTACGAAATAAAAAAGGCTACCTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCt                                                                                                                                         <  1:1573803/150‑1 (MQ=255)
            gCACCTTTTATTATTGATAAAGACAGTGGTGAAATTCATTCTCTGGGAACGGCAAAACCGCTGGAAGAATATCTACAGGATTACGAAATAAAAAAGGCTACCTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCAGCt                                                                                                                                         >  1:52448/1‑150 (MQ=255)
               ccTTTTATTATTGATAAAGACAGTGGTGAAATTCATTCTCTGGGAACGGCAAAACCGCTGGAAGAATATCTACAGGATTACGAAATAAAAAAGGCTACCTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAAc                                                                                                                                      >  1:1167794/1‑150 (MQ=255)
               ccTTTTATTATTGATAAAGACAGTGGTGAAATTCATTCTCTGGGAACGGCAAAACCGCTGGAAGAATATCTACAGGATTACGAAATAAAAAAGGCTACCTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAAc                                                                                                                                      <  1:987168/150‑1 (MQ=255)
                 ttttATTATTGATAAAGACAGTGGTGAAATTCATTCTCTGGGAACGGCAAAACCGCTGGAAGAATATCTACAGGATTACGAAATAAAAAAGGCTACCTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCt                                                                                                                                    >  1:1312053/1‑150 (MQ=255)
                  tttATTATTGATAAAGACAGTGGTGAAATTCATTCTCTGGGAACGGCAAAACCGCTGGAAGAATATCTACAGGATTACGAAATAAAAAAGGCTACCTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCtt                                                                                                                                   <  2:1495020/150‑1 (MQ=255)
                    tattatTGATAAAGACAGTGGTGAAATTCATTCTCTGGGAACGGCAAAACCGCTGGAAGAATATCTACAGGATTACGAAATAAAAAAGGCTACCTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTTAg                                                                                                                                 <  1:1696940/150‑1 (MQ=255)
                      ttattGATAAAGACAGTGGTGAAATTCATTCTCTGGGAACGGCAAAACCGCTGGAAGAATATCTACAGGATTACGAAATAAAAAAGGCTACCTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTTAGTa                                                                                                                               >  1:1204599/1‑150 (MQ=255)
                       tattGATAAAGACAGTGGTGAAATTCATTCTCTGGGAACGGCAAAACCGCTGGAAGAATATCTACAGGATTACGAAATAAAAAAGGCTACCTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTTAGTAc                                                                                                                              >  2:1325807/1‑150 (MQ=255)
                        attGATAAAGACAGTGGTGAAATTCATTCTCTGGGAACGGCAAAACCGCTGGAAGAATATCTACAGGATTACGAAATAAAAAAGGCTACCTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTTAGTAcc                                                                                                                             <  2:730744/150‑1 (MQ=255)
                        attGATAAAGACAGTGGTGAAATTCATTCTCTGGGAACGGCAAAACCGCTGGAAGAATATCTACAGGATTACGAAATAAAAAAGGCTACCTTCGGCTTGCACTGACAAAATAGCCCTCTTCCCCCGAAGAGGGCCGCTAGCCTTAGTAcc                                                                                                                             >  2:369112/1‑150 (MQ=255)
                        attGATAAAGACAGTGGTGAAATTCATTCTCTGGGAACGGCAAAACCGCTGGAAGAATATCTAAAGGATTACGAAATAAAAAAGGCTACCTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTTAGTAcc                                                                                                                             >  2:985399/1‑150 (MQ=255)
                          tGATAAAGACAGTGGTGAAATTCATTCTCTGGGAACGGCAAAACCGCTGGAAGAATATCTACAGGATTACGAAATAAAAAAGGCTACCTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTTAGTAcccc                                                                                                                           >  2:756152/1‑150 (MQ=255)
                           gATAAAGACAGTGGTGAAATTCATTCTATGGGACCGGCAAAACCGATGGAAGAATATCTACAGGATTACGAAATAAAAAAGGCTACCTTCGGCTTGCCCTGACAAAATAGCCCACTTCCCACGAAGAGGGCCGCTAACCTTAGTACAccg                                                                                                                          >  1:1428825/1‑150 (MQ=255)
                             tAAAGACAGTGGTGAAATTCATTCTCTGGGAACGGCAAAACCGCTGGAAGAATATCTACAGGATTACGAAATAAAAAAGGCTACCTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTTAGTACCccgcc                                                                                                                        >  2:971423/1‑150 (MQ=255)
                             tAAAGACAGTGGTGAAATTCAGTCTCTGGGAACGGCAAAACCGCTGGAAGCATATCTCCAGGATTACGAAATAAAAAAGGCTACCTTCGTCTTGCCCTGACAAAATAGCCATCTTCCCACGAAGCGGGCCGCTAACCTTAGTCCCccgcc                                                                                                                        >  2:279525/1‑150 (MQ=255)
                              aaaGACAGTGGTGAAATTCATTCTCTGGGAACGGCAAAACCGCTGGAAGAATATCTACAGGATTACGAAATAAAAAAGGCTACCTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTTAGTACCccgccg                                                                                                                       >  2:267733/1‑150 (MQ=255)
                              aaaGACAGTGGTGAAATTCATTCTCTGGGAACGGCAAAACCGCTGGAAGAATATCTACAGGATTACGAAATAAAAAGGCTACCTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTTAGTACCCCGCCGt                                                                                                                      >  2:1638711/1‑150 (MQ=255)
                               aaGACAGTGGTGAAATTCATTCTCTGGGAACGGCAAAACCGCTGGAAGAATATCTACAGGATTACGAAATAAAAAAGGCTACCTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTTAGTACCCCGCCGt                                                                                                                      >  2:1621731/1‑150 (MQ=255)
                                        gtgAAATTCATTCTCTGGGAACGGCAAAACCGCTGGAAGAATATCTACAGGATTACGAAATAAAAAAGGCTACCTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTTAGTACCCCGCCGTTAAatcatc                                                                                                             >  2:803520/1‑150 (MQ=255)
                                         tgAAATTCATTCTCTTGGAACGGCATAACCGCTGGTAGAAGATCTACAGGATTAGGAAATAAAAAATGCTCCCTTCGGCGTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGTCCGCTAAACTTATTACCCCGCAGTTAAATCATcc                                                                                                            <  1:1277571/150‑1 (MQ=255)
                                         tgAAATTCATTCTCTGGGAACGGCAAAACCGCTGGAAGAATATCTACAGGATTACGAAATAAAAAAGGCTACCTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTTAGTACCCCGCCGTTAAATCATcc                                                                                                            <  1:1015547/150‑1 (MQ=255)
                                                     tctGGGAACGGCAAAACCGCTGGAAGAATATCTACAGGATTACGAAATAAAAAAGGCTACCTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTTAGTACCCCGCCGTTAAATCATCCACCGAGCGCggg                                                                                                <  1:999805/150‑1 (MQ=255)
                                                      ctGGGAACGGCAAAACCGCTGGAAGAATATCTACAGGATTACGAAATAAAAAAGGCTACCTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTTAGTACCCCGCCGTTAAATCATCCCCCGAGCGCGGGt                                                                                               >  1:795045/1‑150 (MQ=255)
                                                      ctGGGAACGGCAAAACCGCTGGAAGAATATCTACAGGATTACGAAATAAAAAAGGCTACCTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTTAGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGt                                                                                               <  2:501057/150‑1 (MQ=255)
                                                        gggAACGGCAAAACCGCTGGAAGAATATCTACAGGATTACGAAATAAAAAAGGCTACCTTCGGCTTGCCCTGACAAAATAGCCCTATTCCCACGAAGAGGCCCCCTAACCTTAGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCg                                                                                             >  2:108166/1‑150 (MQ=255)
                                                          gAACGGCAAAACCGCTGGAAGAATATCTACAGGATTACGAAATAAAAAAGGCTACCTTCGGCTTGCCCTGACAAAATAGCCCTCATCCCACGAAGAGGGCCGCTAACCTTAGTACCCCGCCGTTATATCATCCACCGATCGCCGGTCggt                                                                                           >  2:1070767/1‑149 (MQ=255)
                                                              ggCAAAACCGCTGGAAGAATATCTACAGGATTACGAAATAAAAAAGGCTACCTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTTAGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATgcg                                                                                       >  2:1710648/1‑150 (MQ=255)
                                                              ggCAAAACCGCTGGAAGAATAACTACAGGATTACGAAATAAAAACGGCTACCTTCGGCTTGCCCTGCCAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTTAGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATgcg                                                                                       >  2:777297/1‑150 (MQ=255)
                                                               gCAAAACCGCTGGAAGAATATCTACAGGATTACGAAATAAAAAAGGCTACCTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTTAGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATgcgc                                                                                      <  1:633223/150‑1 (MQ=255)
                                                                cAAAACCGCTGGAAGAATATCTACAGGATTACGAAATAAAAAAGGCTACCTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTTAGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGcc                                                                                     <  1:883680/150‑1 (MQ=255)
                                                                   aaCCGCTGGAAGAATATCTACAGGATTACGAAATAAAAAAGGCTACCTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGACGCTAACCTTAGTACCCCGCCGTGAAATCATCCCACGAGCGCGGGTCGCAAGCGCCGTa                                                                                  >  2:296063/1‑150 (MQ=255)
                                                                      cGCTGGAAGAATATCTACAGGATTACGAAATAAAAAAGGCTACCTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTTAGTACCCCGCCGTTAAATCATCCCCCGAGCGCGGGTCGGATGCGCCGTACaa                                                                               >  2:1275763/1‑150 (MQ=255)
                                                                      cGCTGGAAGAATATCTACAGGATTACGAAATAAAAAAGGCTACCTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTTAGTACCACGCCGTTAAATCATCCCCCGCGCGCGGGACGGATGCGCCGTACac                                                                               >  1:895155/1‑149 (MQ=255)
                                                                       gCTGGAAGAATATCTACAGGATTACGAATTAAAAAAGGCTACCTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTTAGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAAc                                                                              <  2:1133290/150‑1 (MQ=255)
                                                                       gCTGGAAGAATATCTACAGGATTACGAAATAAAAAAGGCTACCTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTTAGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAAc                                                                              <  2:152675/150‑1 (MQ=255)
                                                                          ggAAGAATATCTACAGGATTACGAAATAAAAAAGGCTACCTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTTAGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACACCTcc                                                                           >  1:976636/1‑149 (MQ=255)
                                                                           gaagaaTATCTACAGGATTACGAAATAAAAAAGGCTACCTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTTAGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCAc                                                                          <  1:1467396/150‑1 (MQ=255)
                                                                                 tatCTACATGATTACGAAATAAAAAAGGCTACCTTCGGCTTGCGCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTTAGTACCACGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCg                                                                    <  1:108166/150‑1 (MQ=255)
                                                                                   tCTACAGGATTACGAAATAAAAAAGGCTACCTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTTAGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGc                                                                  <  2:1341910/150‑1 (MQ=255)
                                                                                   tCTACAGGATTACGAAATAAAAAAGGCTACCTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTTAGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGc                                                                  <  2:1728357/150‑1 (MQ=255)
                                                                                    caaCAGGATTACGAAATAAAAAAGGCTACCTTCGGCTTGCCCTGACAAAAAAGCCCTCTTCCCACGAAGATGGCCCCTAACCATAGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGACGTACAACTCACCGTCCgggc                                                                 >  2:1546276/3‑148 (MQ=255)
                                                                                    cTACAGGATTACGAAATAAAAAAGGCTACCTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTTAGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGcc                                                                 <  1:1222425/150‑1 (MQ=255)
                                                                                     tACAGGATTACGAAATAAAAAAGGCGACCTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTTAGGACCCCGCCGGTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGccc                                                                <  1:726087/150‑1 (MQ=255)
                                                                                       cAGGATTACGAAATAAAAAAGGCTACCTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTTAGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCaa                                                              >  2:1344465/1‑150 (MQ=255)
                                                                                          gATTACGAAATAAAAAAGGCTACCTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTTAGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCa                                                           <  1:564052/150‑1 (MQ=255)
                                                                                            ttACGAAATAAAAAAGGCTACCTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTTAGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATa                                                         <  1:1579727/150‑1 (MQ=255)
                                                                                             tACGAAATAAAAAAGGCTACCTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTTAGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATaa                                                        <  1:125838/150‑1 (MQ=255)
                                                                                              aCGAAATAAAAAAGGCTACCTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTTAGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAAt                                                       >  2:275604/1‑150 (MQ=255)
                                                                                               cGAAATAAAAAAGGCTACCTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCCCGAAGAGGGCCGCTAACCTTAGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTAGGATGCGCCGTACAACTCACCGTCCGGCCCATCCATAATg                                                      >  1:1565682/1‑150 (MQ=255)
                                                                                               cGAAATAAAAAAGGCTACCTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTTAGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCCGCCCAACCATAATg                                                      >  1:679784/1‑150 (MQ=255)
                                                                                                 aaaTAAAAAAGGCTACCTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCAAGAAGAGGGCCTCTAACCTTAGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGGACAACTCACCGTCCGGCCCAACCCTACTGCt                                                    >  1:135286/1‑150 (MQ=255)
                                                                                                  aaTAAAAAAGGCTACCTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTTAGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCtt                                                   >  2:574892/1‑150 (MQ=255)
                                                                                                  aaTAAAAAAGGCTACCTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTTAGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCtt                                                   <  1:1659415/150‑1 (MQ=255)
                                                                                                    tAAAAAATGCTACCTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTTAGTACCCCGCCGTTAAATCATCCACCGAGCGCTGGTCGGATGCGCCGTACTACTCACCGTCCGGCCCAACCATAATGCTTtg                                                 <  2:22416/150‑1 (MQ=255)
                                                                                                     aaaaaaGGCTACCTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTTAGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTtgt                                                >  2:1059970/1‑150 (MQ=255)
                                                                                                      aaaaaGTCTACCTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTTAGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAAGTCACCGTCCGGCCCAACCATAATGCTTtgtg                                               <  1:1070767/150‑1 (MQ=255)
                                                                                                      aaaaaGGCTACCTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTTAGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTtgtg                                               >  1:1421188/1‑150 (MQ=255)
                                                                                                      aaaaaGGCTACCTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTTAGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTtgtg                                               <  1:1016019/150‑1 (MQ=255)
                                                                                                       aaaaGGCTACCTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCTCGAAGAGGGCCGCTAACCTTAGTACCCCGCCGTTAAATCATCCAGCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTgtgt                                              <  2:1319420/150‑1 (MQ=255)
                                                                                                         aaGGCTACCTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTTAGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTAc                                            <  2:1167794/150‑1 (MQ=255)
                                                                                                             cTACCTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTTAGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGcc                                        <  1:1242333/150‑1 (MQ=255)
                                                                                                             cTACCTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTTAGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGcc                                        <  2:236646/150‑1 (MQ=255)
                                                                                                                ccTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTTAGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGACGTACAACTCACCGTCCGGCCCAACCATAATGCGTTGTGTACTGCCCAt                                     >  1:1539061/1‑150 (MQ=255)
                                                                                                                 cTTCGGCTTGCCCTGACAAAATAGCCCTCATCCCACGAAGAGGGCCGCTAACCTTAGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATc                                    <  1:134969/150‑1 (MQ=255)
                                                                                                                     ggCTTGCCCTGACAAAATATCCCTCTTCCCACGAAGAGGGCCGCTAACCTTAGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCct                                <  1:1094882/150‑1 (MQ=255)
                                                                                                                      gCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTTAGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCctc                               <  1:413637/150‑1 (MQ=255)
                                                                                                                      gCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTTAGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACACACCGTCCGGCCCAACCATAATGCTATGTGTACTGCCCATCGCctc                               >  1:130215/1‑150 (MQ=255)
                                                                                                                         tGCCCTGACAAAATACCCCTCTTCCCACGCAGAGGGCCGCTAACCTTAGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACATCTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCttt                            <  1:296063/150‑1 (MQ=255)
                                                                                                                          gCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTTAGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTCCACCTCCCCGTCCGGCCCGACCATGATGCTTTGGGTCCTGCCCATCGCCTCTTTc                           >  2:1089352/1‑150 (MQ=255)
                                                                                                                          gCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTTAGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTc                           <  1:1059970/150‑1 (MQ=255)
                                                                                                                             cTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTTAGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAgc                        >  1:504778/1‑150 (MQ=255)
                                                                                                                              tGACTAAATAGCCCTCTTCCCACGAAGAGTGCCGCTAACCTTAGTACCGCGCCCTTAACTCATCCACCGAGCGCGGGTTGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAgcg                       <  2:927389/150‑1 (MQ=255)
                                                                                                                                aCAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTTAGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGcc                     <  2:409158/150‑1 (MQ=255)
                                                                                                                                     gcggCCCTCTTCGTGGGAAGAGGGCCGCTAACCTTAGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGCCACttt                >  1:1212802/4‑150 (MQ=255)
                                                                                                                                      cggCCCTCTTCGTGGGAAGAGGGCCGCTAACCTTAGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGCCACTTTc               >  1:1683216/3‑150 (MQ=255)
                                                                                                                                       ggCCCTCTTCGTGTGAAGAGGGCCGCTAACCTTAGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATACTTCTTTGTGTACTGCCACTCGCCTCTATCAGCGCCACTTTCt              >  1:1245952/2‑150 (MQ=255)
                                                                                                                                       aGCCCTCTTCCCACGAAGAGGGCCGCTAACCTTAGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGCCACTTTCt              <  2:1608330/150‑1 (MQ=255)
                                                                                                                                        gCCCTCTTCGTGGGAAGAGGGCCGCTAACCTTAGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGGACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGCCACTTTCTg             >  1:779298/1‑150 (MQ=255)
                                                                                                                                        gCCCTCTTCCCACGAAGAGGGCCGCTAACCTTAGTACCCCGCCGTTACATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGCCACTTTCTg             <  1:985399/150‑1 (MQ=255)
                                                                                                                                        gCCCTCTTCCCACGAAGAGGGCCGCTAACCTTAGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGCCACTTTCTg             >  2:1481330/1‑150 (MQ=255)
                                                                                                                                        gCCCTCTTCCCACGAAGAGGGCCGCTAACCTTAGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGCCACTTTCTg             <  2:785255/150‑1 (MQ=255)
                                                                                                                                            tctTCCCACGAAGAGGGCCGCTAACCTTAGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGCCACTTTCTGACCt         <  2:52448/150‑1 (MQ=255)
                                                                                                                                                cctaCGAAGAGGGCCGCTAACCTTAGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGCCACTTTCTGACCTTTTg     >  1:644058/4‑150 (MQ=255)
                                                                                                                                                cccACGAAGAGGGCCGCTAACCTTAGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGCCACTTTCTGACCTTTTg     <  1:1621731/150‑1 (MQ=255)
                                                                                                                                                cccACGAAGAGGGCCGCTAACCTTAGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGCCACTTTCTGACCTTTTg     >  1:454965/1‑150 (MQ=255)
                                                                                                                                                cccACGAAGAGGGCCGCTAACCTTAGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGCCACTTTCTGACCTTTTg     >  1:583952/1‑150 (MQ=255)
                                                                                                                                                 ccACGAAGAGGGCCGCTAACCTTAGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGCCACTTTCTGACCTTTTgg    >  1:536862/1‑149 (MQ=255)
                                                                                                                                                 ccACGAAGAGGGCCGCTAACCTTAGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACGGCCCATCGCCTCGTTCAGCGCCACGTTCTGACCTTTTGc    >  1:989683/1‑150 (MQ=255)
                                                                                                                                                 ccACGAAGAGGGCCGCTAACCTTAGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCAAACCATAATGCTTTGTGTACTGCACATCGCCTCTTTCAGCGCCACTTTCGGACCTTTTGc    >  2:1552741/1‑150 (MQ=255)
                                                                                                                                                  cACGAAGAGGGCCGCTAACCTTAGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGCCACTTTCTGACCTTTTGCt   <  2:1421188/150‑1 (MQ=255)
                                                                                                                                                   aCGAAGAGGGCCGCTAACCTTAGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGCCACTTTCTGACCTTTTGCtt  <  1:1638711/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
TTAGCTGGTAACGCACCTTTTATTATTGATAAAGACAGTGGTGAAATTCATTCTCTGGGAACGGCAAAACCGCTGGAAGAATATCTACAGGATTACGAAATAAAAAAGGCTACCTTCGGCTTGCCCTGACAAAATAGCCCTCTTCCCACGAAGAGGGCCGCTAACCTTAGTACCCCGCCGTTAAATCATCCACCGAGCGCGGGTCGGATGCGCCGTACAACTCACCGTCCGGCCCAACCATAATGCTTTGTGTACTGCCCATCGCCTCTTTCAGCGCCACTTTCTGACCTTTTGCTT  >  NC_000913/3584915‑3585211

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: