Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,314,684 A→T 5.6% intergenic (+26/+34) ompW → / ← yciE outer membrane protein W/putative rubrerythrin/ferritin‑like metal‑binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,314,6840AT 5.6% 732.9 / 24.7 210intergenic (+26/+34)ompW/yciEouter membrane protein W/putative rubrerythrin/ferritin‑like metal‑binding protein
Reads supporting (aligned to +/- strand):  ref base A (113/83);  new base T (10/2);  total (124/86)
Fisher's exact test for biased strand distribution p-value = 1.28e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 218 total aligned reads displayed.

CATGTCAGTGTGGTACATGGATATCGATACCACCGCCAATTATAAGCTGGGCGGTGCACAGCAACACGATAGCGTACGCCTCGATCCGTGGGTGTTTATGTTCTCAGCAGGATATCGTTTTTAATTCCGCACAAAAACGACCCCGTAATATACGGGGTCAATAAGGACATGGTATAAAGGGTATTATTTCTTCGCTTCTACGCCATCAGTTTCAGAGCGAATTAAAAATTTCTCAGTTGTTTGCGGAATATTCTGAATCAGCCAGTCGGCCATTTGCTTTTCCTCATTTAAAATCGCTT  >  NC_000913/1314535‑1314833
                                                                                                                                                     |                                                                                                                                                     
cATGTCAGTGTGGTACATGGATATCGATACCACCGCCAATTATAAGCTGGGCGGTGCACAGCAACACGATAGCGTACGCCTCGATCCGTGGGTGTTTATGTTCTCAGCAGGATATCGTTTTTAATTCCGCACAAAAACGACCCCGTAata                                                                                                                                                       <  2:164297/150‑1 (MQ=255)
      aGTGTGGTACATGGATATCGATACCACCGCCAATTATAAGCTGGGCGGTGCACAGCAACACGATAGCGTACGCCTCGATCCGTGGGTGTTTATGTTCTCAGCAGGATATCGTTTTTAATTCCGCACAAAAACGACCCCGTAATCTACggg                                                                                                                                                 >  1:1428541/1‑150 (MQ=255)
              aCATGGATATCGATACCACCGCCAATTATAAGCTGGGCGGTGCACAGCAACACGATAGCGTACGCCTCGATCCGTGGGTGTTTATGTTCTCAGCAGGATATCGTTTTTAATTCCGCACAAAAACGACCCCGTAATATACGGGGTCAATaa                                                                                                                                         >  1:848156/1‑150 (MQ=255)
              aCATGGATATCGATACCACCGCCAATTATAAGCTGGGCGGTGCACAGCAACACGATAGCGTACGCCTCGATCCGTGGGTGTTTATGTTCTCAGCAGGATATCGTTTTTAATTCCGCACAAAAACGACCCCGTAATATACGGGGTCAATaa                                                                                                                                         >  1:2576660/1‑150 (MQ=255)
              aCATGGATATCGATACCACCGCCAATTATAAGCTGGGCGGTGCACAGCAACACGATAGCGTACGCCTCGATCCGTGGGTGTTTATGTTCTCAGCAGGATATCGTTTTTAATTCCGCACAAAAACGACCCCGTAATATACGGGGTCAATaa                                                                                                                                         >  1:161439/1‑150 (MQ=255)
              aCATGGATATCGATACCACCGCCAATTATAAGCTGGGCGGTGCACAGCAACACGATAGCGTACGCCTCGATCCGTGGGTGTTTATGTTCTCAGCAGGATATCGTTTTTAATTCCGCACAAAAACGACCCCGTAATATACGGGGTCAATaa                                                                                                                                         >  1:1701094/1‑150 (MQ=255)
               cATGGATAGCGATACCACCGCCAATTTTAAGCTGGGCGGTGCACTGCAACACGATTGCGTACGCCTCTTTGCGTGGGTGTTTATGTTCTCAGCCGGTTATCGTTTTTAATTCCGCACAAACACGACCCCGTAATATACGGGGTCAATAAg                                                                                                                                        <  2:3911799/150‑1 (MQ=255)
                 tGGATATCGATACCACCGCCAATTATAAGCTGGGCGGTGCACAGCAACACGATAGCGTACGCCTCGATCCGTGGGTGTTTATGTTCTCAGCAGTATATCGTTTTTAATTCCGCACAAAAACGACCCCGTAATATACGGGGTCAATAAGGa                                                                                                                                      <  2:2019656/150‑1 (MQ=255)
                    atatCGATACCACCGCCAATTATAAGCTGGGCGGTGCACAGCAACACGATAGCGTACGCCTCGATCCGTGGGTGTTTATGTTCTCAGCAGGATATCGTTTTTAATTCCGCACAAAAACGACCCCGTAATATACGGGGTCAATAAGGACAt                                                                                                                                   >  1:3954702/1‑150 (MQ=255)
                      atCGATACCACCGCCAATTATAAGCTGGGCGGTGCACAGCAACACGATAGCGTACGCCTCGATCCGTGGGTGTTTATGTTCTCAGCAGGATATCGTTTTTAATTCCGCACAAAAACGACCCCGTAATATACGGGGTCAATAAGGACATgg                                                                                                                                 >  2:111236/1‑150 (MQ=255)
                       tCGATACCACCGCCAATTATAAGCTGGGCGGTGCACAGCAACACGATAGCGTACGCCTCGATCCGTGGGTGTTTATGTTCTCAGCAGGATATCGTTTTTAATTCCGCACAAAAACGACCCCGTAATATACGGGGTCAATAAGGACATGGt                                                                                                                                >  1:443160/1‑150 (MQ=255)
                       tCGATACCACCGCCAATTATAAGCTGGGCGGTGCACAGCAACACGATAGCGTACGCCTCGATCCGTGGGTGTTTATGTTCTCAGCAGGATATCGTTTTTAATTCCGCACAAAAACGACCCCGTAATATACGGGGTCAATAAGGACATGGt                                                                                                                                >  1:3292198/1‑150 (MQ=255)
                       tCGATACCACCGCCAATTATAAGCTGGGCGGTGCACAGCAACACGATAGCGTACGCCTCGATCCGTGGGTGTTTATGTTCTCAGCAGGATATCGTTTTTAATTCCGCACAAAAACGACCCCGTAATATACGGGGTCAATAAGGACATGGt                                                                                                                                >  2:2327874/1‑150 (MQ=255)
                        cGATACCACCGCCAATTATAAGCTGGGCGGTGCACAGCAACACGATAGCGTACGCCTCGATCGGTGGGTGTTTATGTTCTCAGCAGGATATCGTTTTTAATTCCGCACAAAAACGACCCCGTAATATACGGGGTCAATAAGGACATGGta                                                                                                                               <  2:4100559/150‑1 (MQ=255)
                        cGATACCACCGCCAATTATAAGCTGGGCGGTGCACAGCAACACGACAGCGTACGCCTCGATCCGTGGGTGTTTATGTTCTCAGCAGGATATCGTTTTTAATTCCGCACAAAAACGACCCCGTAATATACGGGGTCAATAAGGACATGGta                                                                                                                               <  2:461275/150‑1 (MQ=255)
                          aTACCACCGCCAATTATAAGCTGGGCGGTGCACAGCAACACGATAGCGTACGCCTCGATCCGTGGGTGTTTATGTTCTCAGCAGGATATCGTTTTTAATTCCGCACAAAAACGACCCCGTAATATACGGGGTCAATAAGGACATGGtata                                                                                                                             >  1:699101/1‑150 (MQ=255)
                               accGCCAATTATAAGCTGGGCGGTGCACAGCAACACGATAGCGTACGCCTCGATCCGTGGGTGTTTATGTTCTCAGCAGGATATCGTTTTTAATTCCGCACAAAAACGACCCCGTAATATACGGGGTCACTAAGGACATGGTAGAAAggg                                                                                                                        >  2:2463203/1‑150 (MQ=255)
                                 cGCCAATTATAAGCTGGGCGGTGCACAGCAACACGATAGCGTACGCCTCGATCCGTGGGTGTTTATGTTCTCAGCAGGATATCGTTTTTAATTCCGCACAAAAACGACCCCGTAATATACGGGGTCAATAAGGACATGGTATAAAGGGta                                                                                                                      >  1:1908492/1‑150 (MQ=255)
                                  gCCAATTATAAGCTGGGCGGTGCACAGCAACACGATAGCGTACGCCTCGATCCGTGGGTGTTTATGTTCTCAGCAGGATATCGTTTTTAATTCCGCACAAAATCGACCCCGTAATATACGGGGTCAATAAGGACATGGTATAAAGGGtat                                                                                                                     >  1:1680829/1‑150 (MQ=255)
                                  gCCAATTATAAGCTGGGCGGTGCACAGCAACACGATAGCGTACGCCTCGATCCGTGGGTGTTTATGTTCTCAGCAGGATATCGTTTTTAATTCCGCACAAAAACGACCCCGTAATATACGGGGTCAATAAGGACATGGTATAAAGGGtat                                                                                                                     >  1:2186663/1‑150 (MQ=255)
                                  gCCAATTATAAGCTGGGCGGTGCACAGCAACACGATAGCGTACGCCTCGATCCGTGGGTGTTTATGTTCTCAGCAGGATATCGTTTTTAATTCCGCACAAAAACGACCCCGTAATATACGGGGTCAATAAGGACATGGTATAAAGGGtat                                                                                                                     >  1:2922306/1‑150 (MQ=255)
                                  gCCAATTATAAGCTGGGCGGAGCACAGCAACACGATAGCGAACGCCTCGATCCGTGGGTGTTTATGTTCTCAGCAGGATATCGATTTTAATTCCGCACAAAAACGACCCCGTAATATACGGGGTCACTAAGGACATGGTATAAAGGGtat                                                                                                                     >  1:4094900/1‑150 (MQ=255)
                                   ccAATTATAAGCTGGGCGGTGCACAGCAACACGATAGCGTACGCCTCGATCCGTGGGTGTTTATGTTCTCAGCAGGATATCGTTTTTAATTCCGCACAAAAACGACCCCGTAATATACGGGGTCAATAAGGACATGGTATAAAGGGtatt                                                                                                                    <  2:286952/150‑1 (MQ=255)
                                      aTTATAAGCTGGGCGGTGCACAGCAACACGATAGCGTACGCCTCGATCCGTGGGTGTTTATGTTCTCAGCAGGATATCGTTTTTAATTCCGCACAAAAACGACCCCGTAATATACGGGGTCAATAAGGACATGGTATAAAGGGTattatt                                                                                                                 >  1:2724871/1‑150 (MQ=255)
                                            aGCTGGGCGGTGCACAGCAACACGATAGCGTACGCCTCGATCCGTGGGTGTTTCTGTTCTCAGCAGGATATCGTTTTAAATTCCGCACAAAAACGACCCCGTACTATACGGGGTCAATAAGGACATGGCATAAAGGGTATTATTTCTTCg                                                                                                           >  2:3665138/1‑150 (MQ=255)
                                             gCTGGGCGGTGCACAGCAACACGATAGCGTACGCCTCGATCCGTGGGTGTTTATGTTCTCAGCAGGATATCGTTTTTAATTCCGCACAAAAACGACCCCGTAATATACGGGGTCAATAAGGACATGGTATAAAGGGTATTATTTCTTCGc                                                                                                          >  1:240034/1‑150 (MQ=255)
                                                 ggCGGTGCACAGCAACACGATAGCGTACGCCTCGATCCGTGGGTGTTTATGTTCTCAGCAGGATATCGTTTTTAATTCCGCACAAAAACGACCCCGTAATATACGGGGTCAATAAGGACATGGTATAAAGGGTATTATTTCTTCGCTTCt                                                                                                      >  1:1852196/1‑150 (MQ=255)
                                                 ggCGGTGCACAGCAACACGATAGCGTACGCCTCGATCCGTGGGTGTTTATGTTCTCAGCAGGATATCGTTTTTAATTCCGCACAAAAACGACCCCGTAATATACGGGGTCAATAAGGACATGGTATAAAGGGTATTATTTCTTCGCTTCt                                                                                                      >  1:1363891/1‑150 (MQ=255)
                                                  ggggTGCCCAGCAACACTATAGCGTACGCATCGATCCGTGGGTTTTTTTGTTCTCAGCTGTATATCTTTTTTATTTCCGCAAAAACACGACCCCGTAATATACGGGGTCAATAAGGACATGGTATAAAGGGTATTATTTCTTCGCTTCTa                                                                                                     <  1:1201299/148‑1 (MQ=255)
                                                  gCGGTGCACAGAAACACGATAGCGTACGCCTCGATCCGTGGGTGTTTATGTTCTCAGCAGGATATCGTTTTTAATTCCGCACAAAAACGACCCCGTAATATACGGGGTCAATAAGGACATGGTATAAAGGGTATTATTTCTTCGCTTCTa                                                                                                     <  2:3231063/150‑1 (MQ=255)
                                                    ggTGCACAGCAACTCGATAGCGTACGCCTCGATCCGTGGGTGTGTTTGTTCTCAGCAGGATATCGTTTTTAATTCCGCAAAAAAACGACCCCGTAATATACGGGGTCAATAAGGACATGGTATAAAGGGTATTATTTCTTCGCTTCTACg                                                                                                   <  2:3718978/150‑1 (MQ=255)
                                                    ggTGCACAGCAACCGGATAGCGTACGCCTGGATCCGTGGGTGTTTTTGTTCTGTGCTGGATATCGTTTTTAAGTCCGCACAAAAACGACCCCGTAATCTACGGGGTCAATAAGGACATGGTATAAAGGGTATTATTTCTTCGCTTCTACg                                                                                                   <  2:1357651/150‑1 (MQ=255)
                                                        cacaGCAACACGATAGCGTACGCCTCGATCCGTGGGTGTTTATGTTCTCAGCAGGATATCGTTTTTAATTCCGCACAAAAACGACCCCGTAATATACGGGGTCAATAAGGACATGGTATAAAGGGTATTATTTCTTCGCTTCTACGCCAt                                                                                               >  1:3029581/1‑150 (MQ=255)
                                                        cacaGCAACACGATAGCGTACGCCTCGATCCGTGGGTGTTTATGTTCTCAGCAGGATATCGTTTTTAATTCCGCACAAAAACGACCCCGTAATATACGGGGTCAATAAGGACATGGTATAAAGGGTATTATTTCTTCGCTTCTACGCCAt                                                                                               <  2:3764234/150‑1 (MQ=255)
                                                          caGCAACACGATAGCGTACGCCGCGATCCGGGGGTGTTTATGTTCTCAGCCGGATATCGTTTTTAATTCCGCACAAAAACGACCCCGTAATATACGGGGTCAATAAGGACATGGTATAAAGGGTATTATTTCTTCGCTTCTACGCCATCa                                                                                             <  2:1891214/150‑1 (MQ=255)
                                                              aaCACGATAGCGTACGCCTCGATCCGTGGGTGTTTATGTTCTCAGCAGGATATCGTTTTTAATTCCGCACAAAAACGACCCCGTAATATACGGGGTCAATAAGGACATGGTATAAAGGGTATTATTTCTTCGCTTCTACGCCATCAGttt                                                                                         >  1:3539143/1‑150 (MQ=255)
                                                               acacGATAGCGTACGCCTCGATCCGTGGGTGTTTATGTTCTCAGCAGGATATCGTTTTTAATTCCGCACAAAAACGACCCCGTAATATACGGGGTCAATAAGGACATGGTATAAAGGGTATTATTTCTTCGCTTCTACGCCATCAGTTTc                                                                                        >  1:521987/1‑150 (MQ=255)
                                                               acacGATAGCGTACGCCTCGATCCGTGGGTGTTTATGTTCTCAGCAGGATATCGTTTTTAATTCCGCACAAAAACGACCCCGTAATATACGGGGTCAATAAGGACATGGTATAAAGGGTATTATTTCTTCGCTTCTACGCCATCAGTTTc                                                                                        >  1:3755406/1‑150 (MQ=255)
                                                                cacGATAGCGTACGCCTCGATCCGTGGGTGTTTATGTTCTCAGCAGGATATCGTTTTTAATTCCGCACAAACACGACCCCGTAATATACGGGGTCAATAAGGACATGGTATAAAGGGTATTATTTCTTCGCTTCTACGCCATCAGTTTCa                                                                                       <  1:2089488/150‑1 (MQ=255)
                                                                cacGATAGCGTACGCCTCGATCCGTGGGTGTTTATGTTCTCAGCAGGATATCGTTTTTAATTCCGCACAAAAACGACCCCGTAATATACGGGGTCAATAAGGACATGGTATAAAGGGTATTATTTCTTCGCTTCTACGCCATCAGTTTCa                                                                                       <  2:2686623/150‑1 (MQ=255)
                                                                  cGATAGCGTACGCCTCGATCGGTGGGTGTTTATGTTCTCAGCAGGATATCGTTTTTAATTCCGCACAACAACGACCCCGTAATATACGGGGTCAATAAGGACATGGTATAAAGGGTATTATTTCTTCGCTTCTACGCCATCAGTTTCaga                                                                                     <  2:3188476/150‑1 (MQ=255)
                                                                   gATAGCGTACGCCTCGATCCGTGGGTGTTTATGTTCTCAGCAGGATATCGTTTTTAATTCCGCAAAAAAACGACCCCGTAATATACGGGGTCAATAAGGACATGGTATAAAGGGTATTATTTCTTCGCTTCTACGCCATCAGTTTCagag                                                                                    <  2:2000628/150‑1 (MQ=255)
                                                                        cGTACGCCTCGATCCGTGGGTGTTTATGTTCTCAGCAGGATATCGTTTTTAATTCCGCACAAAAACGACCCCGTAATATACGGGGTCAATAAGGACATGGTATAAAGGGTATTATTTCTTCGCTTCTACGCTATCAGTTTCAGAGCGAAt                                                                               <  2:1774487/150‑1 (MQ=255)
                                                                        cGTACGCCTCGATCCGTGGGTGTTTATGTTCTCAGCAGGATATCGTTTTTAATTCCGCACAAAAACGACCCCGTAATATACGGGGTCAATAAGGACATGGTATAAAGGGTATTATTTCTTCGCTTCTACGCCATCAGTTTCAGAGCGAAt                                                                               <  1:3853811/150‑1 (MQ=255)
                                                                            cGCCTCGATCCGTGGGTGTTTATGTTCTCAGCAGGATATCGTTTTTAATTCCGCACAAAAACGACCCCGTAATATACGGGGTCAATAAGGACATGGTATAAAGGGTATTATTTCTTCGCTTCTACGCCATCAGTTTCAGAGCGAATTaaa                                                                           <  2:1852196/150‑1 (MQ=255)
                                                                             gCCTCGATCGGTGGGTGTTTCTGTTCTCCGCAGGATATCGTTTTTAATTCCGCACAAAAACGACCCCGTAATATACGGGGTCAATAAGGACATGGTATAAAGGGTATTATTTCTTCGCTTCTACGCCATCAGTTTCAGAGCGAATTaaaa                                                                          <  1:3902628/150‑1 (MQ=255)
                                                                               cTCGATCCGTGGGTGTTTATGTTCTCAGCAGGATATCGTTTTTAATTCCGCACAAAAACGACCCCGTAATATACGGGGTCAATAAGGACATGGTATAAAGGGTATTATTTCTTCGCTTCTACGCCATCAGTTTCAGAGCGAATTAAAAAt                                                                        >  2:1688465/1‑150 (MQ=255)
                                                                                 ctaTCCGTGGGTGTTTATGTACTCAGCAGGATATCGTTTTTAATTCCGCAAAAAAACGACCCCGTAATATACGGGGTCAATAAGGACATGGTATAAAGGGTATTATTTCTTCGCTTCTACGCCATCAGTTTCAGAGCGAATTAAAAAttt                                                                      <  2:1446828/148‑1 (MQ=255)
                                                                                  gATCCGTGGGTGTTTATGTTCTCAGCAGGATATCGTTTTTAATTCCGCACAAAAACGACCCCGTAATATACGGGGTCAATAAGGACATGGTATAAAGGGTATTATTTCTTCGCTTCTACGCCATCAGTTTCAGAGCGAATTAAAAATTtc                                                                     <  1:2326835/150‑1 (MQ=255)
                                                                                   aTCCGTGGGTGTTTATGTTCTCAGCAGGATATCGTTTTTAATTCCGCACAAAAACGACCCCGTAATATACGGGGTCAATAAGGACATGGTATAAAGGGTATTATTTCTTCGCTTCTACGCCATCAGTTTCAGAGCGAATTAAAAATTtct                                                                    >  1:1791278/1‑150 (MQ=255)
                                                                                       gTGGGTGTTTATGTTCTCAGCAGGATATCGTTTTTAATTCCGCAAAAAAACGACCCCGTAATATAAGGGGTCAATAAGGACATGGTATAAAGGGTATTATTTCTTAGCTTATACGCGATAAGATTCAGAGAGAATTAAAAATTTCTCAgt                                                                >  1:3684243/1‑150 (MQ=255)
                                                                                          ggTGTTTATGTTCTCAGCAGGATATCGTTTTTAATTCCGCACAAAAACGACCCCGTAATATACGGGGTCAATAAGGACATGGTATAAAGGGTATTATTTCTTCGCTTCTACGCCATCAGTTTCAGAGCGAATTAAAAATTTCTCAgttgt                                                             >  2:1094226/1‑150 (MQ=255)
                                                                                          ggTGTTTATGTTCTCAGCAGGATATCGTTTTTAATTCCGCACAAAAACGACCCCGTAATATACGGGGTCAATAAGGACATGGTATAAAGGGTATTATATCTTCGCTTCTACGCCATCAGTTTCAGAGCGAATTAAAAATTTCTCAgttgt                                                             >  2:3927925/1‑150 (MQ=255)
                                                                                              tttATGTTCTCAGCAGGATATCGTTTTTAATTCCGCACAAAAACGACCCCGTAATATACGGGGTCAATAAGGACATGGTATAAAGGGTATTATTTCTTCGCTTCTACGCCATCAGTTTCAGAGCGAATTAAAAATTTCTCAGTTGTTTGc                                                         >  1:4063627/1‑150 (MQ=255)
                                                                                              tttATGTTCTCAGCAGGATATCGTTTTTAATTCCGCACAAAAACGACCCCGTAATATACGGGGTCAATAAGGACATGGTATAAAGGGTATTATTTCTTCGCTTCTACGCCATCAGTTTCAGAGCGAATTAAAAATTTCTCAGTTGTTTGc                                                         >  1:1100157/1‑150 (MQ=255)
                                                                                                taagTTCTCAGCAGGATATCGTTTTTAATTCCGCACAAAAACGACCCCGTAATATACGGGGTCAATAAGGACATGGTATAAAGGGTATTATTTCTACGCTTCTACGCCATCAGTTTCAGAGCGAATAAAAAATTTCTCAGTTGTTTGCgg                                                       >  2:1227103/4‑150 (MQ=255)
                                                                                                tATGTTCTCAGCAGGATATCGTTTTTAATTCCGCACAAAAACGACCCCGTAATATACGGGGTCAATAAGGACATGGTATAAAGGGTATTATTTCTTCGCTTCTACGCCATCAGTTTCAGAGCGAATTAAAAATTTCTCAGTTGTTTGCgg                                                       >  1:3867886/1‑150 (MQ=255)
                                                                                                      ctcAGCAGGATATCGTTTTTAATTCCGCACAAAAACGACCCCGTAATATACGGGGTCAATAAGGACATGGTATAAAGGGTATTATTTCTTCGCTTCTACGCCATCAGTTTCAGAGCGAATTAAAAATTTCTCAGTTGTTTGCGGAATAtt                                                 >  2:2659932/1‑150 (MQ=255)
                                                                                                       tcAGCAGGATATCGTTTTTAATTCCGCACAAAAACGACCCCGTAATATACGGGGTCAATAAGGACATGGTATAAAGGGTATTATTTCTTCGCTTCTACGCCATCAGTTTCAGAGCGAATTAAAAATTTCTCAGTTGTTTGCGGAATATTc                                                <  2:1326934/150‑1 (MQ=255)
                                                                                                       tcAGCAGGATATCGTTTTTAATTCCGCACAAAAACGACCCCGTAATATACGGGGTCAATAAGGACATGGTATAAAGGGTATTATTTCTTCGCTTCTACGCCATCAGTTTCAGAGCGAATTAAAAATTTCTCAGTTGTTTGCGGAATATTc                                                >  2:1740500/1‑150 (MQ=255)
                                                                                                         agcagGATATCGTTTTTAATTCCGCACAAAAACGACCCCGTAATATACGGGGTCAATAAGGACATGGTATAAAGGGTATTATTTCTTCGCTTCTACGCCATCAGTTTCAGAGCGAATTAAAAATTTCTCAGTTGTTTGCGGAATATTCTg                                              >  2:1318725/1‑150 (MQ=255)
                                                                                                            agGATATCGTTTTTAATTCCGCACAAAAACGACCCCGTAATATACGGGGTCAATAAGGACATGGTATAAAGGGTATTATTTCTTCGCTTCTACGCCATCAGTTTCAGAGCGAATTAAAAATTTCTCAGTTGTTTGCGGAATATTCTGAAt                                           <  2:1067553/150‑1 (MQ=255)
                                                                                                                tatCTTTTTTAAGACCGCACAAACACGACCCCGTAAGATACGGGGTCAATAAGGACATGGTATAAAGGGTATTATTTCTTCGCTTCTACGCCATCAGTTTCAGAGCGAATTAAAAATTTCTCAGTTGTTTGCGGAATATTCTGAATCAGc                                       <  2:2593266/150‑1 (MQ=255)
                                                                                                                 atCGTTTTTAATTCCGCACAAAAACGACCCCGTAATATACGGGGTCAATAAGGACATGGTATAAAGGGTATTATTTCTTCGCTTCTACGCCATCAGTTTCAGAGCGAATTAAAAATTTCTCAGTTGTTTGCGGAATATTCTGAATCAGcc                                      >  2:311641/1‑150 (MQ=255)
                                                                                                                 atCGTTTTTAATTCCGCACAAAAACGACCCCGTAATATACGGGGTCAATAAGGACATGGTATAAAGGGTATTATTTCTTCGCTTCTACGCCATCAGTTTCAGAGCGAATTAAAAATTTCTCAGTTGTTGGCGGAATATTCTGAATCAGcc                                      <  2:2411256/150‑1 (MQ=255)
                                                                                                                  tCGTTTTTAATTCCGCACAAAAACGACCCCGTAATATACGGGGTCAATAAGGACATGGTATAAAGGGTATTATTTCTTCGCTTCTACGCCATCAGTTTCAGAGCGAATTAAAAATTTCTCAGTTGTTTGCGGAATATTCTGAATCAGCCa                                     <  2:3292198/150‑1 (MQ=255)
                                                                                                                    gTTTTTAATTCCGCACAAAAACGACCCCGTAATATACGGGGTCAATAAGGACATGGTATAAAGGGTATTTTTTCTTCGCTTCTACGCCTTGAGTTTCAGAGCGAATTTATAATTTGTGTGTTGTTTGCGGAATATTCTGTTTCTGCCAGt                                   >  2:2229504/1‑150 (MQ=255)
                                                                                                                    gTTTTTAATTCCGCACAAAAACGACCCCGTAATATACGGGGTCAATAAGGACATGGTATAAAGGGTATTATTTCTTCGCTTCTACGCCATCAGTTTCAGAGCGAATTAAAAATTTCTCAGTTGTTTGCGGAATATTCTGAATCAGCCAGt                                   >  2:3451685/1‑150 (MQ=255)
                                                                                                                     tttttAATTCCGCACAAAAACGACCCCGTAATATACGGGGTCAATAAGGACATGGTATAAAGGGTATTATTTCTTCGCTTCTACGCCATCAGTTTCAGAGCGAATTAAAAATTTCTCAGTTGTTTGCGGAATATTCTGAATCAGCCAGTc                                  >  2:1028617/1‑150 (MQ=255)
                                                                                                                     tttttAATTCCGCACAAAAACGACCCCGTAATATACGGGGTCAATAAGGACATGGTATAAAGGGTATTATTTCTTCGCTTCTACGCCATCAGTTTCAGAGCGAATTAAAAATTTCTCAGTTGTTTGCGGAATATTCTGAATCAGCCAGTc                                  >  2:592993/1‑150 (MQ=255)
                                                                                                                     tttttAATTCCGCACAAAAACGACCCCGTAATATACGGGGTCAATAAGGACATGGTATAAAGGGTATTATTTCTTCGCTTCTACGCCATCAGTTTCAGAGCGAATTAAAAATTTCTCAGTTGTTTGCGGAATATTCTGAATCAGCCAGTc                                  >  2:1746180/1‑150 (MQ=255)
                                                                                                                             tCCGCACAAACACGACCCCGTAATATACGGGGTCAATAAGGACATGGTATAAAGGGTATTATTTCTTCGCTTCTACGCCATCAGTTTCAGAGCGAATTAAAAATTTCTCAGTTGTTTGCGGAATATTCTGAATCAGCCAGTCGGCCAttt                          <  1:1278172/150‑1 (MQ=255)
                                                                                                                                gCACAAAAACGACCCCGTAATATACGGGGTCAATAAGGACATGGTATAAAGGGTATTATTTCTTCGCTTCTACGCCATCAGTTTCAGAGCGAATTAAAAATTTCTCAGTTGTTTGCGGAATATTCTGAATCAGCCAGTCGGCCATTTGCt                       >  2:418287/1‑150 (MQ=255)
                                                                                                                                 cacaAAAACGACCCCGTAATATACGGGGTCAATAAGGACATGGTATAAAGGGTATTATTTCTTCGCTTCTACGCCATCAGTTTCAGAGCGAATTAAAAATTTCTCAGTTGTTTGCGGAATATTCTGAATCAGCCAGTCGGCCATTTGCtt                      >  2:2128491/1‑150 (MQ=255)
                                                                                                                                 cacaAAAACGACCCCGTAATATACGGGGTCAATAAGGACATGGTATAAAGGGTATTATTTCTTCGCTTCTACGCCATCAGTTTCAGAGCGAATTAAAAATTTCTCAGTTGTTTGCGGAATATTCTGAATCAGCCAGTCGGCCATTTGCtt                      >  2:3254446/1‑150 (MQ=255)
                                                                                                                                  acaAAAACGACCCCGTAATATACGGGGTCAATAAGGACATGGTATAAAGGGTATTATTTCTTCGCTTCTACGCCATCAGTTTCAGAGCGAATTAAAAATTTCTCAGTTGTTTGCGGAATATTCTGAATCAGCCAGTCGGCCATTTGCttt                     >  1:2300236/1‑150 (MQ=255)
                                                                                                                                  acaAAAACGACCCCGTAATATACGGGGTCAATAAGGACATGGTATAAAGGGTATTATTTCTTCGCTTCTACGCCATCAGTTTCAGAGCGAATTAAAAATTTCTCAGTTGTTTGCGGAATATTCTGAATCAGCCAGTCGGCCATTTGCttt                     <  1:3392674/150‑1 (MQ=255)
                                                                                                                                  acaAAAACGACCCCGTAATATACGGGGTCAATAAGGACATGGTATAAAGGGTATTATTTCTTCGCTTCTACGCCATCAGTTTCAGAGCGAATTAAAAATTTCTCAGTTGTTTGCGGAATATTCTGAATCAGCCAGTCGGCCATTTGCttt                     >  2:1838627/1‑150 (MQ=255)
                                                                                                                                  acaAAAACGACCCCGTAATATACGGGGTCAATAAGGACATGGTATAAAGGGTATTATTTCTTCGCTTCTACGCCATCAGTTTCAGAGCGAATTAAAAATTTCTCAGTTGTTTGCGGAATATTCTGAATCAGCCAGTCGGCCATTTGCttt                     <  1:28287/150‑1 (MQ=255)
                                                                                                                                   caAAAACGACCCCGTAATATACGGGGTCAATACGGACATGGGATAATGGGTATTATTTCTTCGCTTCTACGCCATCAGTTTCAGAGCGAATTAAAGATTTCTCAGTTGTTTGCGGAATATTCTGAATCAGCCAGTCGGCCATTTGCtttt                    >  2:2718747/1‑150 (MQ=255)
                                                                                                                                   caAAAACGACCCCGTAATATACGGGGTCAATAAGGACATGGTATAAAGGGTATTATTTCTTCGCTTCTACGCCATCAGTTTCAGAGCGAATTAAAAATTTCTCAGTTGTTTGCGGAATATTCTGAATCAGCCAGTCGGCCATTTGCtttt                    <  1:2622875/150‑1 (MQ=255)
                                                                                                                                    aaaaaCGACCCCGTAATATACGGGGTCAATAAGGACATGGTATAAAGGGTATTATTTCTTCGCTTCTACGCCATCAGTTTCAGAGCGAATTAAAAATTTCTCAGTTGTTTGCGGAATATTCTGAATCAGCCAGTCGGCCATTTGCTTTTc                   <  1:2944671/150‑1 (MQ=255)
                                                                                                                                    aaaaaCGACCCCGTAATATACGGGGTCAATAAGGACATGGTATAAAGGGTATTATTTCTTCGCTTCTACGCCATCAGTTTCAGAGCGAATTAAAAATTTCTCAGTTGTTTGCGGAATATTCTGAATCAGCCAGTCGGCCATTTGCTTTTc                   <  2:1494278/150‑1 (MQ=255)
                                                                                                                                     ttattgACCCCGTATATTACGGGGTCAATAAGGACATGGTATAAAGGGTATTATTTCTTCGCTTCTACGCCATCAGTTTCAGAGCGAATTAAAAATTTCTCAGTTGTTTGCGGAATATTCTGAATCAGCCAGTCGGCCATTTGCTTTTcc                  >  1:2019309/6‑150 (MQ=255)
                                                                                                                                     ttattgACCCCGTATATTACGGGGTCAATAAGGACATGGTATAAAGGGTATTATTTCTTAGCTTCTACTCCATAAGTTTCATAGCGAATTAAAAATTTATCAGTTGTTTGAGGAAAATACTGAATCAGCCAGTCGGCCATTTGCTTTTcc                  >  1:1567107/6‑150 (MQ=255)
                                                                                                                                     ttattgACCCCGTATATTACGGGGTCAATAAGGACATGGTATAAAGGGTATTATTTCATCGCTTCTACGCCATCAGTTTCAGAGCGAATTAAAAATTTCTCAGTTGTTTGCGGAATATTCTGAATCAGCCAGTCGGCCATTTGCTTTTcc                  >  1:2019446/6‑150 (MQ=255)
                                                                                                                                     aaaaCGACCCCGTAATATACGGGGTCAATAAGGACATGGTATAAAGGGTATTATTTCTTCGCTTCTACGCCATCAGTTTCAGAGCGAATTAAAAATTTCTCAGTTGTTTGCGGAATATTCTGAATCAGCCAGTCGGCCATTTGCTTTTcc                  <  2:770063/150‑1 (MQ=255)
                                                                                                                                     aaaaCGACCCCGTAATATACGGGGTCAATAAGGACATGGTATAAAGGGTATTATTTCTTCGCTTCTACGCCATCAGTTTCAGAGCGAATTAAAAATTTCTCAGTTGTTTGCGGAATATTCTGAATCAGCCAGTCGGCCATTTGCTTTTcc                  <  1:3694428/150‑1 (MQ=255)
                                                                                                                                      aaaCGACCCCGTAATATACGGGGTCAATAAGGACATGGTATAAAGGGTATTATTTCTTCGCTTCTACGCCATCAGTTTCAGAGCGAATTAAAAATTTCTCAGTTGTTTGCGGAATATTCTGAATCAGCCAGTCGGCCATTTGCTTTTCCt                 >  1:1003642/1‑150 (MQ=255)
                                                                                                                                        ttgACCCCGTATATTACGGGGTCAATAAGGACATGGTATAAAGGGTATTATTTCTTCGCTTCTACGCCATCAGTTTCAGAGCGAATTAAAAATTTCTCAGTTGTTTGCGGAATATTCTGAATCAGCCAGTCGGCCATTTGCTTTTCCTCa               >  1:4153903/3‑150 (MQ=255)
                                                                                                                                        ttgACCCCGTATATTACGGGGTCAATAAGGACATGGTATAAAGGGTATTATTTCTTCGCTTCTACGCCATCAGTTTCAGAGCGAATTAAAAATTTCTCAGTTGTTTGCGGAATATTCTGAATCAGCCAGTCGGCCATTTGCTTTTCCTCa               >  1:3577801/3‑150 (MQ=255)
                                                                                                                                        ttgACCCCGTATATTACGGGGTCAATAAGGACATGGTATAAAGGGTATTATTTCTTCGCTTCTACGCCATCAGTTTCAGAGCGAATTAAAAATTTCTCAGTTGTTTGCGGAATATTCTGAATCAGCCAGTCGGCCATTTGCTTTTCCTCa               >  1:2299804/3‑150 (MQ=255)
                                                                                                                                        ttgACCCCGTATATTACGGGGTCAATAAGGACATGGTATAAAGGGTATTATTTCTTCGCTTCTACGCCATCAGTTGCAGAGCGAATTAAAAATTTCTCAGTTGTTTGCGGAATATTCTGAATCAGCCAGTCGGCCATTTGCTTTTCCTCa               >  1:1675457/3‑150 (MQ=255)
                                                                                                                                        aCGACCCCGTAATATACGGGGTCAATAAGGACATGGTATAAAGGGTATTATTTCTTCGCTTCTACGCCATCAGTTTCAGAGCGAATTAAAAATTTCTCAGTTGTTTGCGGAATATTCTGAATCAGCCAGTCGGCCATTTGCTTTTCCTCa               <  2:3023445/150‑1 (MQ=255)
                                                                                                                                        aCGACCCCGTAATATACGGGGTCAATAAGGACATGGTATAAAGGGTATTATTTCTTCGCTTCTACGCCATCAGTTTCAGAGCGAATTAAAAATTTCTCAGTTGTTTGCGGAATATTCTGAATCAGCCAGTCGGCCATTTGCTTTTCCTCa               <  2:1791278/150‑1 (MQ=255)
                                                                                                                                              ccGTAATATACGGGGTCAATAAGGACATGGTATAAAGGGTATTATTTCTTCGCTTCTACGCCATCAGTTTCAGAGCGAATTAAAAATTTCTCAGTTGTTTGCGGAATATTCTGAATCAGCCAGTCGGCCATTTGCTTTTCCTCATTTaaa         <  2:1027332/150‑1 (MQ=255)
                                                                                                                                              ccGTAATATACGGGGTCAATAAGGACATGGTATAAAGGGTATTATTTCTTCGCTTCTACGCCATCAGTTTCAGAGCGAATTAAAAATTTCTCAGTTGTTTGCGGAATATTCTGAATCAGCCAGTCGGCCATTTGCTTTTCCTCATTTaaa         >  1:1637254/1‑150 (MQ=255)
                                                                                                                                               cGTAATATACGGGGTCAATAAGGACATGGTATAAAGGGTATTATTTCTTCGCTTCTACGCCATCAGTTTCAGAGCGAATTAAAAATTTCTCAGTTGTTTGCGGAATATTCTGAATCAGCCAGTCGGCCATTTGCTTTTCCTCATTTaaaa        >  2:1939698/1‑150 (MQ=255)
                                                                                                                                                    tataCGGGGTCAATAAGGACATGGTATAAAGGGTATTATTTCTTCGCTTCTACGCCATCTGTTTCAGAGCGAATTAAATATTTCTCAGTTGTTTGCGGAATATTCTGTATCAGACAGTCTGCCATTTGCTTTTCCTCATTTAAAATCGCt   >  1:2808326/1‑150 (MQ=255)
                                                                                                                                                     ataCGGGGTCAATAAGGACATGGTATAAAGGGTATTATTTCTTCGCTTCTACGCCATCAGTTTCAGAGCGAATTAAAAATTTCTCAGTTGTTTGCGGAATATTCTGAATCAGCCAGTCGGCCATTTGCTTTTCCTCATTTAAAATCGCtt  <  1:3060281/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
CATGTCAGTGTGGTACATGGATATCGATACCACCGCCAATTATAAGCTGGGCGGTGCACAGCAACACGATAGCGTACGCCTCGATCCGTGGGTGTTTATGTTCTCAGCAGGATATCGTTTTTAATTCCGCACAAAAACGACCCCGTAATATACGGGGTCAATAAGGACATGGTATAAAGGGTATTATTTCTTCGCTTCTACGCCATCAGTTTCAGAGCGAATTAAAAATTTCTCAGTTGTTTGCGGAATATTCTGAATCAGCCAGTCGGCCATTTGCTTTTCCTCATTTAAAATCGCTT  >  NC_000913/1314535‑1314833

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 11 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: