Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,314,685:1 +T 15.2% intergenic (+27/+33) 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,6841.T15.2% 287.3 / 57.3 92intergenic (+26/+34)ompW/yciEouter membrane protein W/putative rubrerythrin/ferritin‑like metal‑binding protein
Reads supporting (aligned to +/- strand):  ref base . (43/35);  new base T (8/6);  total (51/41)
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.69e-02

CGGTGCACAGCAACACGATAGCGTACGCCTCGATCCGTGGGTGTTTATGTTCTCAGCAGGATATCGTTTTTAATTCCGCACAAAAACGACCCCGTAATA‑TACGGGGTCAATAAGGACATGGTATAAAGGGTATTATTTCTTCGCTTCTACGCCATCAGTTTCAGAGCGAATTAAAAATTTCTCAGTTGTTTGCGGAAT  >  NC_000913/1314586‑1314783
                                                                                                   |                                                                                                   
cGGTGCACAGCAACACGATAGCGTACGCCTCGATCCGTGGGTGTTTATGTTCTCAGCAGGATATCGTTTTTAATTCCGCACAAAAACGACCCCGTAata‑t                                                                                                    <  1:1703182/100‑1 (MQ=255)
    ggagaggaaCGCGATAGCGTACGCCGCGAGCCGTGGGTGGGTATGTTCGCAGCAGGATATCGTTTTTGATTCCGCACAAAAACGACCCCGTAATA‑TACgg                                                                                                <  2:2239815/93‑1 (MQ=255)
     cacaGCAACACGATAGCGTACGCCTCGATCCGTGGGTGTTTATGTTCTCAGCAGGATATCGTTTTTAATTCCGCACAAAAACGACCCCGTAATA‑TACggg                                                                                               <  1:784685/100‑1 (MQ=255)
      acaGCAACACGATAGCGTACGCCTCGATCCGTGGGTGTTTATGTTCTCAGCAGGATATCGTTTTTAATTCCGCACAAAAACGACCCCGTAATA‑TACgggg                                                                                              <  2:2054070/100‑1 (MQ=255)
      acaGCAACACGATAGCGTACGCCTCGATCCGTGGGTGTTTATGTTCTCAGCAGGATATCGTTTTTAATTCCGCACAAAAACGACCCCGTAATA‑TACggg                                                                                               <  1:389149/99‑1 (MQ=255)
       caGCAACACGATAGCGTACGCCTCGATCCGTGGGTGTTTATGTTCTCAGCAGGATATCGTTTTTAATTCCGCACAAAAACGACCCCGTAATA‑TACGGGGt                                                                                             >  1:2245037/1‑100 (MQ=255)
       caGCAACACGATAGCGTACGCCTCGATCCGTGGGTGTTTATGTTCTCAGCAGGATATCGTTTTTAATTCCGCACAAAAACGACCCCGTAATA‑TACGGGGt                                                                                             <  1:1128457/100‑1 (MQ=255)
       caGCAACACGATAGCGTACGCCTCGATCCGTGGGTGTTTATGTTCTCAGCAGGATATCGTTTTTAATTCCGCACAAAAACGACCCCGTAATA‑TACGGGGt                                                                                             >  2:523625/1‑100 (MQ=255)
       caGCAACACGATAGCGTACGCCTCGATCCGTGGGTGTTTATGTTCTCAGCAGGATATCGTTTTTAATTCCGCACAAAAACGACCCCGT‑ATATTACGGGGt                                                                                             >  1:303687/1‑100 (MQ=255)
        aGCAACACGATAGCGTACGCCTCGATCCGTGGGTGTTTATGTTCTCAGCAGGATATCGTTTTTAATTCCGCACAAAAACGACCCCGTAATA‑TACGGGGTc                                                                                            <  1:1039255/100‑1 (MQ=255)
         gCAACACGATAGCGTACGCCTCGATCCGTGGGTGTTTATGTTCTCAGCAGGATATCGTTTTTAATTCCGCACAAAAACGACCCCGTAATA‑TACgggg                                                                                              >  2:185143/1‑97 (MQ=255)
         gCAACACGATAGCGTACGCCTCGATCCGTGGGTGTTTATGTTCTCAGCAGGATATCGTTTTTAATTCCGCACAAAAACGACCCCGT‑ATATTACGGGGTcg                                                                                           <  2:303685/100‑2 (MQ=255)
         gCAACACGATAGCGTACGCCTCGATCCGTGGGTGTTTATGTTCTCAGCAGGATATCGTTTTTAATTCCGCACAAAAACGACCCCGT‑ATATTACGGGGTcg                                                                                           <  2:1838060/100‑2 (MQ=255)
           aaCACGATAGCGTACGCCTCGATCCGTGGGTGTTTATGTTCTCAGCAGGATATCGTTTTTAATTCCGCACAAAAACGCACCCGTAATA‑TACGGGGTCAAt                                                                                         >  2:460272/1‑100 (MQ=255)
           aaCACGATAGCGTACGCCTCGATCCGTGGGTGTTTATGTTCTCAGCAGGATATCGTTTTTAATTCCGCACAAAAACGACCCCGT‑ATATTACGGGGTcgtt                                                                                         <  1:772618/100‑4 (MQ=255)
           aaCACGATAGCGTACGCCTCGATCCGTGGGTGTTTATGTTCTCAGCAGGATATCGTTTTTAATTCCGCACAAAAACGACCCCGT‑ATATTACGGGGTcgtt                                                                                         <  2:1495628/100‑4 (MQ=255)
              acGATAGCGTACGCCTCGATCCGTGGGTGTTTATGTTCTCAGCAGGATATCGTTTTTAATTCCGCACAAAAACGACCCCGTAATA‑TACGGGGTCAATAAg                                                                                      <  1:830445/100‑1 (MQ=255)
              acGATAGCGTACGCCTCGATCCGTGGGTGTTTATGTTCTCAGCAGGATATCGTTTTTAATTCCGCACAAAAACGACCCCGTAATA‑TACGGGGTCAATAAg                                                                                      >  2:1130107/1‑100 (MQ=255)
              acGATAGCGTACGCCTCGATCCGTGGGTGTTTATGTTCTCAGCAGGATATCGTTTTTAATTCCGCACAAAAACGACCCCGTAATA‑TACGGGATCAATAAg                                                                                      <  1:472348/100‑1 (MQ=255)
               cGATAGCGTACGCCTCGATCCGTGGGTGTTTATGTTCTCAGCAGGATATCGTTTTTAATTCCGCACAAAAACGACCCCGTAATA‑TACGGGGTCAATAAgg                                                                                     <  2:1874243/100‑1 (MQ=255)
               cGATAGCGTACGCCTCGATCCGTGGGTGTTTATGTTCTCAGCAGGATATCGTTTTTAATTCCGCACAAAAACGACCCCGTAATA‑TACGGGGTCAATAAgg                                                                                     >  1:1339949/1‑100 (MQ=255)
               cGATAGCGTACGCCTCGATCCGTGGGTGTTTATGTTCTCAGCAGGATATCGTTTTTAATTCCGCACAAAAACGACCCCGT‑ATATTACGGGGTcgtttttg                                                                                     <  1:885607/100‑8 (MQ=38)
               cGATAGCGTACGCCTCGATCCGTGGGTGTTTATGTTCTCAGAAGGATATCGTTTTTAATTCAGCACAAAAACGACCCCGTAATA‑TACGGGGGAAATaa                                                                                       >  1:38381/1‑98 (MQ=255)
                gttagcccaCGCCTCGATCCCTGGGTGTTTATGTTCTCAGCAGGATATCGTTTTTAATTCCGCACAAAAACGACCCCGTAATA‑TACGGGGTCAATAAGGa                                                                                    <  2:1613930/92‑1 (MQ=255)
                 aTAGCGTACGCCTCGATCCGTGGGTGTTTATGTTCTCAGCAGGATATCGTTTTTAATTCCGCACAAAAACGACCCCGT‑ATATTACGGGGTcgtttttgtg                                                                                   <  1:2360256/100‑10 (MQ=37)
                    gCGTACGCCTCGATCCGTGGGTGTTTATGTTCTCAGCAGGATATCGTTTTTAATTCCGCACAAAAACGACCCCGTAATA‑TACGGGGTCAATAAGGACATg                                                                                <  2:18921/100‑1 (MQ=255)
                     cGTACGCCTCGATCCGTGGGTGTTTATGTTCTCAGCAGGATATCGTTTTTAATTCCGCACAAAAACGACCCCGTAATA‑TACGGGGTCAATAAGGACATgg                                                                               <  1:44725/100‑1 (MQ=255)
                      gTACGCCTCGATCCGTGGGTGTTTATGTTCTCAGCAGGATATCGTTTTTAATTCCGCACAAAAACGACCCCGTAATA‑TACGGGGTCAATAAGGACATGGt                                                                              >  2:93845/1‑100 (MQ=255)
                       tACGCCTCGATCCGTGGGTGTTTATGTTCTCAGCAGGATATCGTTTTTAATTCCGCACAAAAACGACCCCGTAATA‑TACGGGGTCAATAAGGACATGGta                                                                             <  2:18384/100‑1 (MQ=255)
                          gCCTCGATCCGTGGGTGTTTATGTTCTCAGCAGGATATCGTTTTTAATTCCGAACAAAAACGACCACGTAATG‑TACGGGGTCAATAAGGAAATGGTATaa                                                                          >  2:1943015/1‑100 (MQ=255)
                           ccTCGATCCGTGGGTGTTTATGTTCTCAGCAGGATATCGTTTTTAATTCCGCACAAAAACGACCCCGTAATA‑TACGGGGTCAATAAGGACATGGTATaaa                                                                         >  1:1668882/1‑100 (MQ=255)
                           ccTCGATCCGTGGGTGTTTATGTTCTCAGCAGGATATCGTTTTTAATTCCGCACAAAAACGACCCCGTAATA‑TACGGGGTCAATAAGGACATGGTATaaa                                                                         >  2:1855218/1‑100 (MQ=255)
                           ccTCGATCCGTGGGTGTTTATGTTCTCAGCAGGATATCGTTTTTAATTCCGCACAAAAACGACCCCGTAATA‑TACGGGGTCAATAAGGACATGGTATaaa                                                                         >  1:2312858/1‑100 (MQ=255)
                           ccTCGATCCGTGGGTGTTTATGTTCTCAGCAGGATATCGTTTTTAATTCCGCACAAAAACGACCCCGTAATA‑TACGGGGTCAATAAGGACATGGTATaaa                                                                         >  2:950564/1‑100 (MQ=255)
                            cTCGATCCGTGGGTGTTTATGTTCTCAGCAGGATATCGTTTTAAATACCGCACAAAAACGATCCCGTAATA‑TACGGGGTCAATAAGGACATGGTATAAAg                                                                        >  2:1294103/1‑100 (MQ=255)
                                   cGTGGGTGTTTATGTTCTCAGCAGGGTATCGTTTTTAATTCCGCACAAAAACGACCCCGTAATA‑TACGGGGTCAATAAGGACATGGTATAAAGGGtatta                                                                 >  1:871416/1‑100 (MQ=255)
                                   cGTGGGTGTTTATGTTCTCAGCAGGATATCGTTTTTAATTCCGCACAAAAACGACCCCGTAATA‑TACGGGGTCAATAAGGACATGGTATAAAGGGtatta                                                                 >  1:471991/1‑100 (MQ=255)
                                   cGTGGGTGTTTATGTTCTCAGCAGGATATCGTTTTTAATTCCGCACAAAAACGACCCCGTAATA‑TACGGGGTCAATAAGGACATGGTATAAAGGGtatta                                                                 >  1:505923/1‑100 (MQ=255)
                                   cGTGGGTGTTTATGTTCTCAGCAGGATATCGTTTTTAATTCCGCACAAAAACGACCCCGTAATA‑TACGGGGTCAATAAGGACATGGTACAAAGGGtatta                                                                 >  1:2267788/1‑100 (MQ=255)
                                    gTGGGTGTTTATGTTCTCAGGAGGATATCGTTTTTAATTCCGCACAAAAACGACCCCGTAATA‑TACGGGGTCAATAAGGACATGGTATAAAGGGtat                                                                   <  2:1954691/97‑1 (MQ=255)
                                      gggTGTTTATGTTCTCAGCAGGATATCGTTTTTAATTCCGCACAAAAACGACCCCGTAATA‑TACGGGGTCAATAAGGACATGGTATAAAGGGTATTAttt                                                              >  1:1834483/1‑100 (MQ=255)
                                      gggTGTGTATATTCTCAACAGGATATCGTTTTTAATTCCGCACAAAAACGACCCAGTAATA‑TACGGGGTCAATAAGGACATCGTAGAAAAGCtatt                                                                  >  2:3667/1‑96 (MQ=255)
                                       ggTGTTTATGTTCTCAGCAGGATATCGTTTTTAATTCCGCACAAAAACGACCCCGTAATA‑TACGGGGTCAATAAGGACATGGTATAAAGGGTATTATttc                                                             >  1:840766/1‑100 (MQ=255)
                                          gtTTATGGTCTCAGCAGGATATCGTTTTTAATTCCGCACAAAAACGACCCCGTAATA‑TACGGGGTAAATAAGGACATGGTATAAAGGGTATTATttcttc                                                          <  1:1339926/100‑1 (MQ=255)
                                             tATGTTCTCAGCAGGATATCGTTTTTAATTCCGCACAAAAACGACCCCGTAATA‑TACGGGGTCAATAAGGACATGGTATAAAGGGTATTATTTCTTCGCt                                                       >  1:376717/1‑100 (MQ=255)
                                             tATGTTCTCAGCAGGATATCGTTTTTAATTCCGCACAAAAACGACCCCGTAATA‑TACGGGGTCAATAAGGACATGGTATAAAGGGTATTATTTCTTCGCt                                                       >  1:1853901/1‑100 (MQ=255)
                                              aTGTTCTCAGCAGGATATCGTTTTTAATTCCGCACAAAAACGACCCCGTAATA‑TACGGGGTCAATAAGGACATGGTATAAAGGGTATTATTTCTTCGCtt                                                      <  2:1820282/100‑1 (MQ=255)
                                                  tctcAGCAGGATATCGTTTTTAATTCCGCACAAAAACGACCCCGTAATA‑TACGGGGTCAATAAGGACATGGTATAAAGGGTATTATTTCTTCGCTTCTAc                                                  <  1:1282892/100‑1 (MQ=255)
                                                   ctcAGCAGGATATCGTTTTTAATTCCGCACAAAAACGACCCCGTAATA‑TACGGGGTCAATAAGGACATGGTATAAAGGGTATTATTTCTTCGCTTCTACg                                                 >  1:913641/1‑100 (MQ=255)
                                                     cagcagGATATCGTTTTTAATTCCGCACAAAAACGACCCCGTAATA‑TACGGGGTCAATAAGGACATGGTATAAAGGGTATTATTTCTTCGCTTCTACGc                                                >  1:2033714/1‑99 (MQ=255)
                                                     cagcagGATATCGTTTTTAATTCCGCACAAAAACGACCCCGTAATA‑TACGGGGTCAAAAAGGACATGGTATAAAGGGTATTATTTCTTCGCTTCTAcccc                                               >  2:868853/1‑100 (MQ=255)
                                                      agcagGATATCGTTTTTAATTCCGCACAAAAACGACCCCGTAATA‑TACGGGGTCAATAAGGACATGGTATAAAGGGTATTATTTCTTCGCTTCTACGCCa                                              >  1:1540741/1‑100 (MQ=255)
                                                         agGATATCGTTTTTAATTCCGCACAAAAACGACCCCGTAATA‑TACGGGGTCAATAAGGACATGGTATAAAGGGTATTATTTCTTCGCTTCTACGCCATCa                                           <  1:1107894/100‑1 (MQ=255)
                                                         agGATACCCTTTTTAATGCCGCACACAGACGACCCCGTGACA‑TACGGGGTCCTTCAGGCCATGGTATAACGTGTCATATTTCCTCGCTTCTACGCCAt                                             <  2:1797648/98‑1 (MQ=255)
                                                           gATATCGTATTTAATGCCGGACAAAAACGACCCCGTAATA‑TACGGGGTCAATAAGGGCATGGTATAAAGGGTATTATTTCTTCTCTTCTACGCCATCAGt                                         <  2:2252851/100‑1 (MQ=255)
                                                            atatCGTTTTTAATTCCGCACAAAAACGACCCCGTAATA‑TACGGGGTCAATAAGGACATGGTATAAAGGGTATTATTTCTTCGCTTCTACGCCATCAGtt                                        <  1:720826/100‑1 (MQ=255)
                                                             tatCGTTTTTAATTCCGCACAAAAACGACCCCGTAATA‑TACGGGGTCAATAAGGACATGGTATAAAGGGTATTATTTCTTCGCTTCTACGCCATCAGttt                                       <  1:858618/100‑1 (MQ=255)
                                                             tatCGTTTTTAATTCCGCACAAAAACGACCCCGTAATA‑TACGGGGTCAATAAGGACATGGTATAAAGGGTATTATTTCTTCGCTTCTACGCCATCAGttt                                       <  1:813136/100‑1 (MQ=255)
                                                              atCGTTTTTAATTCCGCACAAAAACGACCCCGTAATA‑TACGGGGTCAATAAGGACAGGGTATAAAGGGTATTATTTCTTCGCTTCTACGCCATCAGTTTc                                      >  1:2216782/1‑100 (MQ=255)
                                                               tCGTTTTTAATTCCGCACAAAAACGACCCCGTAATA‑TACGGGGTCAATAAGGACATGGTATAAAGGGTATTATTTCTTCGCTTCTACGCCATCAGTTTCa                                     <  1:201361/100‑1 (MQ=255)
                                                                cGTTTTTAATTCCGCACAAAAACGACCCCGTAATA‑TACGGGGTCAATAAGGACATGGTATAAAGGGTATTATTTCTTCGCTTCTACGCCATCAGTTTCag                                    >  1:1509314/1‑100 (MQ=255)
                                                                cGTTTTTAATTCCGCACAAAAACGACCCCGTAATA‑TACGGGGTCAATAAGGACATGGTATAAAGGGTATTATTTCTTCGCTTCTACGCCATCAGTTTCag                                    >  1:1067464/1‑100 (MQ=255)
                                                                 gTTTTTAATTCCGCACAAAAACGACCCCGTAATA‑TACGGGGTCAATAAGGACATGGTATAAAGGGTATTATTTCTTCGCTTCTACGCCATCAGTTTCaga                                   >  1:1033648/1‑100 (MQ=255)
                                                                  tttttAATTCCGCACAAAAACGACCCCGTAATA‑TACGGGGTCAATAAGGACATGGTATAAAGGGTATTATTTCTTCGCTTCTACGCCATCAGTTTCagag                                  <  1:657743/100‑1 (MQ=255)
                                                                   ttttAATTCAGAACAAAAACGACCCCGTAATA‑TACGGGGTCAATAAGGACATGGTATAAAGGGTATTATTTCTTCGCTTCTACGCCATCAGTTTCAGAGc                                 <  2:255919/100‑1 (MQ=255)
                                                                    tttAATTCCGCACAAAAACGACCCCGTAATA‑TACGGGGTCAATAAGGACATGGTATAAAGGGTATTATTTCTTCGCTTCTACGCCATCAGTTTCAGAGCg                                <  2:1239332/100‑1 (MQ=255)
                                                                    tttAATTCCGCACAAAAACGACCCCGTAATA‑TACGCGGTCACTAAGGACATGGCATAAAGGGTATTATTTCTTCGCTTCTACGCCATCAGTTTCAGAGCg                                >  2:2120702/1‑100 (MQ=255)
                                                                      taatcacGCACAAAAACGACCCCGTCATA‑TACGGGGTCAATAAGGACATGGTATAAAGGGTATTATTTCTTCGCTTCTACGCCATCAGTTTCAGAGCGa                               <  2:1885442/93‑1 (MQ=255)
                                                                           ccGCACAAAAACGACCCCGTAATA‑TACGGGGTCAATAAGGACATGGTATAAAGGGTATTATTTCTTCGCTTCTACTCCATCAGTTTCAGAGCGAATTaaa                         >  2:1254118/1‑100 (MQ=255)
                                                                           ccGCACAAAAACGACCCCGTAATA‑TACGGGGTCAATAAGGACATGGTATAAAGGGTATTATTTCTTCGCTTCTACGCCATCAGTTTCAGAGCGAATTaaa                         >  2:2148441/1‑100 (MQ=255)
                                                                            cGCACAAAAACGACCCCGTAATA‑TACGGGGTCAATAAGGACATGGTATAAAGGGTATTATTTCTTCGCTTCTACGCCATCAGTTTCAGAGCGAATTaaaa                        >  2:493334/1‑100 (MQ=255)
                                                                            cGCACAAAAACGACCCCGTAATA‑TACGGGGTCAATAAGGACATGGTATAAAGGGTATTATTTCTTCGCTTCTACACCATCAGTTTCAGAGCGAATTaaaa                        >  2:1878578/1‑100 (MQ=255)
                                                                              cacaAAAACGACCCCGTAATA‑TACGGGGTCAATAAGGACATGGTATAAAGGGTATTATTTCTTCGCTTCTACGCCATCAGTTTCAGAGCGAATTAAAAAt                      <  1:1278322/100‑1 (MQ=255)
                                                                               acaAAAACGACCCAGTAATA‑TACGGGGTCAATAAGGACATGGTATAAAGGGTATTATTTCTTCGCTTCTACGCCATCAGTTTCAGAACAAATTAAAAAtt                     <  2:1699285/100‑1 (MQ=255)
                                                                                 ttccgatctcccGTAATA‑TACGGGGTCAATAAGGACATGGTATAAAGGGTATTATTTCTTCGCTTCTACGCCATCAGTTTCAGAGCGAATTAAAAATTtc                   <  1:746763/91‑1 (MQ=255)
                                                                                 aaaaaaGACCGCGTAATA‑TACGGGATCAATAAGGACAGGGTATAAAGTGTAGTATTTCTTCGCGTCTACGCCATCAGTTTCAGAGCGAATTAAAAATTtc                   <  2:2310907/100‑1 (MQ=255)
                                                                                 aaaaaCGACCCCGTAATA‑TACGGGGTCAATAAGGACATGGTATAAAGGGTATTATTTCTTCGCTTCTACGCCATCAGTTTCAGAGCGAATTAAAAATTtc                   <  2:826766/100‑1 (MQ=255)
                                                                                   tattgACCCCGT‑ATATTACGGGGTCAATAAGGACATGGTATAAAGGGTATTATTTCTTCGCTTCTACGCCATCAGTTTCAGAGCGAATTAAAAATTtctc                 >  2:1364238/5‑100 (MQ=255)
                                                                                   nattgACCCCGT‑ATATTACGGGGTCAATAAGGACATGGTATAAAGGGTATTATTTCTTCGCTTCTACGCCATCAGTTTCAGAGCGAATTAAAAATTtctc                 >  2:1747601/5‑100 (MQ=255)
                                                                                   aaaCGACCCCGTAATA‑TACGGGGTCAATAAGGACATGGTATAAAGGGTATTATTTCTTCGCTTCTACGCCATCAGTTTCAGAGCGAATTAAAAATTtctc                 >  1:622159/1‑100 (MQ=255)
                                                                                   aaaCGACCCCGTAATA‑TACGGGGTCAATAAGGACATGGTATAAAGGGTATTATTTCTTCGCTTCTACGCCATCAGTTTCAGAGCGAATTAAAAATTtct                  >  2:842648/1‑99 (MQ=255)
                                                                                    attgACCCCGT‑ATATTACGGGGTCGAT‑AGGACATGGGATAAAGGGTATTATTTTTTCCGTTCTAAGCCATTCATTTCAAGGGGGATTTAAAAttt                    >  2:1691162/4‑95 (MQ=9)
                                                                                    attgACCCCGT‑ATATTACGGGGTCAATAAGGACATGGTATAAAGGGTATTATTTCTTCGCTTCTACGCCATCAGTTTCAGAGCGAATTAAAAATTTCTCa                >  2:2150270/4‑100 (MQ=255)
                                                                                       gACCCCGT‑ATATTACGGGGTCAATAAGGACATGGTATAAAGGGTATTATTTCTTCGCTTCTACGCCATCAGTTTCAGAGCGAATTAAAAATTTCTCAgtt             >  1:1283388/1‑100 (MQ=255)
                                                                                       gACCCCGT‑ATATTACGGGGTCAATAAGGACATGGTATAAAGGGTATTATTTCTTCGCTTCTACGCCATCAGTTTCAGAGCGAATTAAAAATTTCTCAgtt             >  2:187472/1‑100 (MQ=255)
                                                                                        aCCCCGTAATA‑TACGGGGTCAATAAGGACATGGTATAAAGGGTATTATTCCTTCGCTTCTACGCCATCAGTTTCAGAGCGAATTAAAAATTTCTCAgttg            >  2:964392/1‑100 (MQ=255)
                                                                                        aCCCCGT‑ATATTACGGGGTCAATAAGGACATGGTATAAAGGGTATTATTTCTTCGCTTCTACGCCATCAGTTTCAGAGCGAATTAAAAATTTCTCAgttg            >  1:2296278/1‑100 (MQ=255)
                                                                                         ccccGTAATA‑TACGGGGTCAATAAGGACATGGTATAAAGGGTATTATTTCTTCGCTTCTACGCCATCAGTTTCAGAGCGAATTAAAAATTTCTCAgttgt           >  1:1695551/1‑100 (MQ=255)
                                                                                          cccGTAATA‑TACGGGGTCAATAAGGACATGGTATAAAGGGTATTATTTCTTCGCTTCTACGCCATCAGTTTCAGAGCGAATTAAAAATTtcagatcggaa          >  2:746761/1‑91 (MQ=255)
                                                                                           ccGTAATA‑TACGGGGTCAATAAGGACATGGTATAAAGTGTATTATTTCTTCGCTTCTACGCCATCAGTTTCAGAGCGAATTAAAAATTTCTCAGTTGttt         >  2:1269088/1‑100 (MQ=255)
                                                                                              tAATA‑TACGGGGTCAATAAGGACATGGTATAAAGGGTATTATTTCTTCGCTTCTACGCCATCAGTTTCAGAGCGAATTAAAAATTTCTCAGTTGTTTGCg      <  2:997815/100‑1 (MQ=255)
                                                                                               aaTA‑TACGGGGTCAATAAGGACATGGTATAAAGGGTATTATTTCTTCGCTTCTACGCCATCAGTTTCAGAGCGAATTAAAAATTTCTCAGTTGTTTGCgg     <  2:1668869/100‑1 (MQ=255)
                                                                                               aaTA‑TACGGGGTCAATAAGGACATGGTATAAAGGGTATTATTTCTTCGCTTCTACGCCATCAGTTTCAGAGCGAATTAAAAATTTCTCAGTTGTTTGCgg     >  1:2323279/1‑100 (MQ=255)
                                                                                                  a‑taCGGGGTCAATAAGGACATGGTATAAAGGGTATTATTTCTTCGCTTCTACGCCATCAGTTTCAGAGCGAATTAAAAATTTCTCAGTTGTTTGCGGAat  <  1:2230518/100‑1 (MQ=255)
                                                                                                   |                                                                                                   
CGGTGCACAGCAACACGATAGCGTACGCCTCGATCCGTGGGTGTTTATGTTCTCAGCAGGATATCGTTTTTAATTCCGCACAAAAACGACCCCGTAATA‑TACGGGGTCAATAAGGACATGGTATAAAGGGTATTATTTCTTCGCTTCTACGCCATCAGTTTCAGAGCGAATTAAAAATTTCTCAGTTGTTTGCGGAAT  >  NC_000913/1314586‑1314783

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