Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,840,547 T→C intergenic (+39/+2) yicS → / ← nepI putative periplasmic protein/putative transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,840,5470TC91.9% 305.8 / 12.7 99intergenic (+39/+2)yicS/nepIputative periplasmic protein/putative transporter
Reads supporting (aligned to +/- strand):  ref base T (3/5);  new base C (43/48);  total (46/53)
Fisher's exact test for biased strand distribution p-value = 7.21e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.72e-01

TAATAATCAGAGCAATTATTGGGAAGCGTTGGGTAAGGTAGAATGCCCGGATATGTAAGAACAGCAACAACTACTGACCAACGCTTCTTTACCTGATATCAGGATTTCTTCATTTTCACCTTTGCAGTAACCAACAATGCAGTCAGCAACATCAATGTGCCGGACAACATCAACGGCGAAGTCAGACCAATATT  >  NC_000913/3840451‑3840644
                                                                                                |                                                                                                 
taataaTCAGAGCAATTATTGGGAAGCGTTGGGTAAGGTATAATGCCCGGATATGTAAGAACAGCAACAACTACTGACCAACGCTTCTTTACCTGacat                                                                                                 >  1:52430/1‑99 (MQ=255)
taataaTCAGAGCAATTATTGGGAAGCGTTGGGTAAGGTAGAATGCCCGGATATGTAAGAACAGCAACAACTACTGACCAACGCTTCTTTACCTGACATc                                                                                                >  1:407434/1‑100 (MQ=255)
 aataatCAGAGCAATTATTGGGAAGCGTTGGGTAAGGTAGAATGCCCGGATATGTAAGAACAGCAACAACTACTGACCAACGCTTCTTTACCTGACATCa                                                                                               >  1:113020/1‑100 (MQ=255)
    aatCAGAGCAATTATTGGGAAGCGTTGGGTAAGGTAGAATGCCCGGATATGTAAGAACAGCAACAACTACTGACCAACGCTTCTTTACCTGACATCAGGa                                                                                            >  1:407098/1‑100 (MQ=255)
    aatCAGAGCAATTATTGGGAAGCGTTGGGTAAGGTAGAATGCCCGGATATGTAAGAACAGCAACAACTACTGACCAACGCTTCTTTACCTGACATCAGGa                                                                                            >  1:299677/1‑100 (MQ=255)
       cAGAGCAATTATTGGGAAGCGTTGGGTAAGGTAGAATGCCCGGATATGTAAGAACAGCAACAACTACTGACCAACGCTTCTTTACCTGACATCAGGAttt                                                                                         <  2:1668487/100‑1 (MQ=255)
       cAGAGCAATTATTGGGAAGCGTTGGGTAAGGTAGAATACCCGGATATGTAAGAACAGCAACAACTACTGACCAACGCTTCTTTACCTGACATCAGGAttt                                                                                         >  1:1040890/1‑100 (MQ=255)
        agagCAATTATTGGGAAGCGTTGGGTAAGGTAGAATGCCCGGATATGTAAGAACAGCAACAACTACTGACCAACGCTTCTTTACCTGACATCAGGATttc                                                                                        >  2:1559402/1‑100 (MQ=255)
          agCAATTATTGGGAAGCGTTGGGTAAGGTAGAATGCCCGGATATGTAAGAACAGCAACAACTACTGACCAACGCTTCTTTACCTGACATCAGGATttctt                                                                                      <  2:775416/100‑1 (MQ=255)
           gCAATTATTGGGAAGCGTTGGGTAAGGTAGAATGCCCGGATATGTAAGAACAGCAACAACTACTGACCAACGCTTCTTTACCTGACATCAGGATttcttc                                                                                     >  1:1156234/1‑100 (MQ=255)
            cAATTATTGGGAAGCGTTGGGTAAGGTAGAATGCCCGGATATGTAAGAACAGCAACAACTACTGACCAACGCTTCTTTACCTGACATCAGGATTTCTTCa                                                                                    <  2:1740581/100‑1 (MQ=255)
            cAATTATTGGGAAGCGTTGGGTAAGGTAGAATGCCCGGATATGTAAGAACAGCAACAACTACTGACCAACGCTTCTTTACCTGACATCAGGATTTCTTCa                                                                                    <  1:242189/100‑1 (MQ=255)
            cAATTATTGGGAAGCGTTGGGTAAGGTAGAATGCCCGGATATGTAAGAACAGCAACAACTACTGACCAACGCTTCTTTACCTGACATCAGGATTTCTTCa                                                                                    >  1:1794527/1‑100 (MQ=255)
             aaTTATTGGGAAGCGTTGGGTAAGGTAGAATGCCCGGATATGTAAGAACAGCAACAACTACTGACCAACGCTTCTTTACCTGACATCAGGATTTCTTCAt                                                                                   >  2:837300/1‑100 (MQ=255)
              attattGGGAAGCGTTGGGTAAGGTAGAATGCCCGGATATGTAAGAACAGCAACAACTACTGACCAACGCTTCTTTACCTGACATCAGGATTTCTTCAtt                                                                                  >  2:643800/1‑100 (MQ=255)
               ttattGGGAAGCGTTGGGTAAGGTAGAATGCCCGGATATGTAAGAACAGCAACAACTACTGACCAACGCTTCTTTACCTGACATCAGGATTTCTTCAttt                                                                                 <  2:748851/100‑1 (MQ=255)
               ttattGGGAAGCGTTGGGTAAGGTAGAATGCCCGGATATGTAAGAACAGCAACAACTACTGACCAACGCTTCTTTACCTGACATCAGGATTTCTTCAttt                                                                                 <  1:26780/100‑1 (MQ=255)
                  ttGGGAAGCGTTGGGTAAGGTAGAATGCCCGGATATGTAAGAACAGCAACAACTACTGACCAACGCTTCTTTACCTGACATCAGGATTTCTTCATTTTCa                                                                              <  2:553414/100‑1 (MQ=255)
                  ttGGGAAGCGTTGGGTAAGGTAGAATGCCCGGATATGTAAGAACAGCAACAACTACTGACCAACGCTTCTTTACCTGACATCAGGATTTCTTCATTTTCa                                                                              <  1:870870/100‑1 (MQ=255)
                  ttGGGAAGCGTTGGGTAAGGTAGAATGCCCGGATATGTAAGAACAGCAACAACTACTGACCAACGCTTCTTTACCTGACATCAGGATTTCTTCATTTTCa                                                                              >  2:172855/1‑100 (MQ=255)
                   tGGGAAGCGTTGGGTAAGGTAGAATGCCCGGATATGTAAGAACAGCAACAACTACTGACCAACGCTTCTTTACCTGACATCAGGATTTCTTCATTTTCAc                                                                             >  1:243281/1‑100 (MQ=255)
                    gggAAGCGTTGGGTAAGGTAGAATGCCCGGATATGTAAGAACAGCAACAACTACTGACCAACGCTTCTTTACCTGACATCAGGATTTCTTCATTTTCAcc                                                                            >  1:1177903/1‑100 (MQ=255)
                     ggAAGCGTTGGGTAAGGTAGAATGCCCGGATATGTAAGAACAGCAACAACTACTGACCAACGCTTCTTTACCTGACATCAGGATTTCTTCATTTTCACCt                                                                           <  2:547948/100‑1 (MQ=255)
                      gAAGCGTTGGGTAAGGTAGAATGCCCGGATATGTAAGAACAGCAACAACTACTGACCAACGCTTCTTTACCTGACATCAGGATTTCTTCATTTTCACCtt                                                                          >  2:191892/1‑100 (MQ=255)
                       aaGCGTTGGGTAAGGTAGAATGCCCGGATATGTAAGAACAGCAACAACTACTGACCAACGCTTCTTTACCTGACATCAGGATTTCTTCATTTTCACCttt                                                                         <  2:407436/100‑1 (MQ=255)
                        aGCGTTGGGTAAGGTAGAATGCCCGGATATGTAAGAACAGCAACAACTACTGACCAACGCTTCTTTACCTGACATCAGGATTTCTTCATTTTCACCTTTg                                                                        <  2:681605/100‑1 (MQ=255)
                        aGCGTTGGGTAAGGTAGAATGCCCGGATATGTAAGAACAGCAACAACTACTGACCAACGCTTCTTTACCTGACATCAGGATTTCTTCATTTTCACCTTTg                                                                        <  1:837299/100‑1 (MQ=255)
                         gCGTTGGGTAAGGTAGAATGCCCGGATATGTAAGAACAGCAACAACTACTGACCAACGCTTCTTTACCTGACATCAGGATTTCTTCATTTTCACCTTTGc                                                                       <  2:135734/100‑1 (MQ=255)
                         gCGTTGGGTAAGGTAGAATGCCCGGATATGTAAGAACAGCAACAACTACTGACCAACGCTTCTTTACCTGACATCAGGATTTCTTCATTTTCACCTTTGc                                                                       <  2:466825/100‑1 (MQ=255)
                         gCGTTGGGTAAGGTAGAATGCCCGGATATGTAAGAACAGCAACAACTACTGACCAACGCTTCTTTACCTGACATCAGGATTTCTTCATTTTCACCTTTGc                                                                       <  1:424816/100‑1 (MQ=255)
                          cGTTGGGTAAGGTAGAATGCCCGGATATGTAAGAACAGCAACAACTACTGACCAACGCTTCTTTACCTGACATCAGGATTTCTTCATTTTCACCTTTGCa                                                                      >  2:531689/1‑100 (MQ=255)
                          cGTTGGGTAAGGTAGAATGCCCGGATATGTAAGAACAGCAACAACTACTGACCAACGCTTCTTTACCTGACATCAGGATTTCTTCATTTTCACCTTTGCa                                                                      >  2:1729387/1‑100 (MQ=255)
                          cGTTGGGTAAGGTAGAATGCCCGGATATGTAAGAACAGCAACAACTACTGACCAACGCTTCTTTACCTGACATCAGGATTTCTTCATTTTCACCTTTGCa                                                                      >  1:1050999/1‑100 (MQ=255)
                           gTTGGGTAAGGTAGAATGCCCGGATATGTAAGAACAGCAACAACTACTGACCAACGCTTCTTTACCTGATATCAGGATTTCTTCATTTTCACCTTTGCAg                                                                     <  2:1175058/100‑1 (MQ=255)
                           gTTGGGTAAGGTAGAATGCCCGGATATGTAAGAACAGCAACAACTACTGACCAACGCTTCTTTACCTGACATCAGGATTTCTTCATTTTCACCTTTGCAg                                                                     <  2:1842905/100‑1 (MQ=255)
                              gggTAAGGTAGAATGCCCGGATATGTAAGAACAGCAACAACTACTGACCAACGCTTCTTTACCTGACATCAGGATTTCTTCATTTTCACCTTTGCAGTaa                                                                  >  1:1502952/1‑100 (MQ=255)
                               ggTAAGGTAGAATGCCCGGATATGTAAGAACAGCAACAACTACTGACCAACGCTTCTTTACCTGATATCAGGATTTCTTCATTTTCACCTTTGCAGTAAc                                                                 <  2:507866/100‑1 (MQ=255)
                               ggTAAGGTAGAATGCCCGGATATGTAAGAACAGCAACAACTACTGACCAACGCTTCTTTACCTGACATCAGGATTTCTTCATTTTCACCTTTGCAGTAAc                                                                 <  2:1177905/100‑1 (MQ=255)
                               ggTAAGGTAGAATGCCCGGATATGTAAGAACAGCAACAACTACTGACCAACGCTTCTTTACCTGACATCAGGATTTCTTCATTTTCACCTTTGCAGTAAc                                                                 >  1:537789/1‑100 (MQ=255)
                               ggTAAGGTAGAATGCCCGGATATGTAAGAACAGCAACAACTACTGACCAACGCTTCTTTACCTGACATCAGGATTTCTTCATTTTCACCTTTGCAGTAAc                                                                 >  2:1501034/1‑100 (MQ=255)
                                gTAAGGTAGAATGCCCGGATATGTAAGAACAGCAACAACTACTGACCAACGCTTCTTTACCTGACATCAGGATTTCTTCATTTTCACCTTTGCAGTAAcc                                                                >  2:1699448/1‑100 (MQ=255)
                                  aaGGTAGAATGCCCGGATATGTAAGAACAGCAACAACTACTGACCAACGCTTCTTTACCTGACATCAGGATTTCTTCATTTTCACCTTTGCAGTAACcaa                                                              >  1:61641/1‑100 (MQ=255)
                                  aaGGTAGAATGCCCGGATATGTAAGAACAGCAACAACTACTGACCAACGCTTCTTTACCTGACATCAGGATTTCTTCATTTTCACCTTTGCAGTAACcaa                                                              >  1:1517161/1‑100 (MQ=255)
                                    ggTAGAATGCCCGGATATGTAAGAACAGCAACAACTACTGACCAACGCTTCTTTACCTGACATCAGGATTTCTTCATTTTCACCTTTGCAGTAACcaaca                                                            >  2:85289/1‑100 (MQ=255)
                                    ggTAGAATGCCCGGATATGTAAGAACAGCAACAACTACTGACCAACGCTTCTTTACCTGACATCAGGATTTCTTCATTTTCACCTTTGCAGTAACcaaca                                                            >  2:1783713/1‑100 (MQ=255)
                                    ggTAGAATGCCCGGATATGTAAGAACAGCAACAACTACTGACCAACGCTTCTTTACCTGACATCAGGATTTCTTCATTTTCACCTTTGCAGTAACcaaca                                                            >  2:85193/1‑100 (MQ=255)
                                    ggTAGAATGCCCGGATATGTAAGAACAGCAACAACTACTGACCAACGCTTCTTTACCTGACATCAGGATTTCTTCATTTTCACCTTTGCAGTAACcaaca                                                            <  1:133009/100‑1 (MQ=255)
                                     gTAGAATGCCCGGATATGTAAGAACAGCAACAACTACTGACCAACGCTTCTTTACCTGACATCAGGATTTCTTCATTTTCACCTTTGCAGTAACcaacaa                                                           <  1:741642/100‑1 (MQ=255)
                                      tAGAATGCCCGGATATGTAAGAACAGCAACAACTACTGACCAACGCTTCTTTACCTGATATCAGGATTTCTTCATTTTCACCTTTGCAGTAACCAACAAt                                                          >  2:904201/1‑100 (MQ=255)
                                      tAGAATGCCCGGATATGTAAGAACAGCAACAACTACTGACCAACGCTTCTTTACCTGATATCAGGATTTCTTCATTTTCACCTTTGCAGTAACCAACAAt                                                          >  1:599586/1‑100 (MQ=255)
                                        gAATGCCCGGATATGTAAGAACAGCAACAACTACTGACCAACGCTTCTTTACCTGACATCAGGATTTCTTCATTTTCACCTTTGCAGTAACCAACAATGc                                                        <  1:191891/100‑1 (MQ=255)
                                         aaTGCCCGGATATGTAAGAACAGCAACAACTACTGACCAACGCTTCTTTACCTGACATCAGGATTTCTTCATTTTCACCTTTGCAGTAACCAACAATGCa                                                       >  1:1163047/1‑100 (MQ=255)
                                          aTGCCCGGATATGTAAGAACAGCAACAACTACTGACCAACGCTTCTTTACCTGATATCAGGATTTCTTCATTTTCACCTTTGCAGTAACCAACAATGCAg                                                      <  1:875276/100‑1 (MQ=255)
                                          aTGCCCGGATATGTAAGAACAGCAACAACTACTGACCAACGCTTCTTTACCTGACATCAGGATTTCTTCATTTTCACCTTTGCAGTAACCAACAATGCAg                                                      <  2:1278854/100‑1 (MQ=255)
                                            gCCCGGATATGTAAGAACAGCAACAACTACTGACCAACGCTTCTTTACCTGACATCAGGATTTCTTCATTTTCACCTTTGCAGTAACCAACAATGCAGTc                                                    >  2:1007318/1‑100 (MQ=255)
                                            gCCCGGATATGTAAGAACAGCAACAACTACTGACCAACGCTTCTTTACCTGACATCAGGATTTCTTCATTTTCACCTTTGCAGTAACCAACAATGCAGTc                                                    <  2:1004148/100‑1 (MQ=255)
                                            gCCCGGATATGTAAGAACAGCAACAACTACTGACCAACGCTTCTTTACCTGACATCAGGATTTCTTCATTTTCACCTTTGCAGTAACCAACAATGCAGTc                                                    <  1:199322/100‑1 (MQ=255)
                                               cGGATATGTAAGAACAGCAACAACTACTGACCAACGCTTCTTTACCTGACATCAGGATTTCTTCATTTTCACCTTTGCAGTAACCAACAATGCAGTCAGc                                                 <  2:1239054/100‑1 (MQ=255)
                                               cGGATATGTAAGAACAGCAACAACTACTGACCAACGCTTCTTTACCTGACATCAGGATTTCTTCATTTTCACCTTTGCAGTAACCAACAATGCAGTCAGc                                                 <  2:338368/100‑1 (MQ=255)
                                                ggATATGTAAGAACAGCAACAACTACTGACCAACGCTTCTTTACCTGACATCAGGATTTCTTCATTTTCACCTTTGCAGTAACCAACAATGCAGTCAGCa                                                <  1:1286566/100‑1 (MQ=255)
                                                 gATATGTAAGAACAGCAACAACTACTGACCAACGCTTCTTTACCTGACATCAGGATTTCTTCATTTTCACCTTTGCAGTAACCAACAATGCAGTCAGCaa                                               <  2:230389/100‑1 (MQ=255)
                                                  atatGTAAGAACAGCAACAACTACTGACCAACGCTTCTTTACCTGACATCAGGATTTCTTCATTTTCACCTTTGCAGTAACCAACAATGCAGTCAGCAAc                                              <  2:130726/100‑1 (MQ=255)
                                                    atGTAAGAACAGCAACAACTACTGACCAACGCTTCTTTACCTGACATCAGGATTTCTTCATTTTCACCTTTGCAGTAACCAACAATGCAGTCAGCAACAt                                            >  2:890587/1‑100 (MQ=255)
                                                     tGTAAGAACAGCAACAACTACTGACCAACGCTTCTTTACCTGACATCAGGATTTCTTCATTTTCACCTTTGCAGTAACCAACAATGCAGTCATCAACATc                                           <  1:738295/100‑1 (MQ=255)
                                                         aGAACAGCAACAACTACTGACCAACGCTTCTTTACCTGACATCAGGATTTCTTCATTTTCACCTTTGCAGTAACCAACAATGCAGTCAGCAACATCAAtg                                       <  1:1058033/100‑1 (MQ=255)
                                                         aGAACAGCAACAACTACTGACCAACGCTTCTTTACCTGACATCAGGATTTCTTCATTTTCACCTTTGCAGTAACCAACAATGCAGTCAGCAACATCAAt                                        >  1:519373/1‑99 (MQ=255)
                                                             cAGCAACAACTACTGACCAACGCTTCTTTACCTGACATCAGGATTTCTTCATTTTCACCTTTGCAGTAACCAACAATGCAGTCAGCAACATCAATGTGcc                                   <  1:1673301/100‑1 (MQ=255)
                                                             cAGCAACAACTACTGACCAACGCTTCTTTACCTGACATCAGGATTTCTTCATTTTCACCTTTGCAGTAACCAACAATGCAGTCAGCAACATCAATGTGcc                                   <  1:948430/100‑1 (MQ=255)
                                                             cAGCAACAACTACTGACCAACGCTTCTTTACCTGACATCAGGATTTCTTCATTTTCACCTTTGCAGTAACCAACAATGCAGTCAGCAACATCAATGTGcc                                   <  2:341526/100‑1 (MQ=255)
                                                              aGCAACAACTACTGACCAACGCTTCTTTACCTGACATCAGGATTTCTTCATTTTCACCTTTGCAGTAACCAACAATGCAGTCAGCAACATCAATGTGCCg                                  <  1:1615148/100‑1 (MQ=255)
                                                               gCAACAACTACTGACCAACGCTTCTTTACCTGATATCAGGATTTCTTCATTTTCACCTTTGCAGTAACCAACAATGCAGTCAGCAACATCAATGTGCCgg                                 <  1:904200/100‑1 (MQ=255)
                                                               gCAACAACTACTGACCAACGCTTCTTTACCTGACATCAGGATTTCTTCATTTTCACCTTTGCAGTAACCAACAATGCAGTCAGCAACATCAATGTGCCgg                                 >  2:1098436/1‑100 (MQ=255)
                                                               gCAACAACTACTGACCAACGCTTCTTTACCTGACATCAGGATTTCTTCATTTTCACCTTTGCAGTAACCAACAATCCAGTCAGCAACATCAATGTGcc                                   >  2:1820565/1‑98 (MQ=255)
                                                                 aacaacTACTGACCAACGCTTCTTTACCTGACATCAGGATTTCTTCATTTTCACCTTTGCAGTAACCAACAATGCAGTCAGCAACATCAATGTGCCGGac                               <  1:1668835/100‑1 (MQ=255)
                                                                    agttactGACCATCGCTTCTTTTCCTGTCATCAGGTTTTCTTCATTTTCACCTTTGCAGTAACCAATATTGCAGTCAGCAACATCAATGTGCCg                                  <  1:1157846/91‑1 (MQ=255)
                                                                    aacTACTGACCAACGCTTCTTTACCTGACATCAGGATTTCTTCATTTTCACCTTTGCAGTAACCAACAATGCAGTCAGCAACATCAATGTGCCGGacaac                            <  1:85289/100‑1 (MQ=255)
                                                                    aacTACTGACCAACGCTTCTTTACCTGACATCAGGATTTCTTCATTTTCACCTTTGCAGTAACCAACAATGCAGTCAGCAACATCAATGTGCCGGacaac                            <  2:413822/100‑1 (MQ=255)
                                                                    aacTACTGACCAACGCTTCTTTACCTGACATCAGGATTTCTTCATTTTCACCTTTGCAGTAACCAACAATGCAGTCAGCAACATCAATGTGCCGGacaac                            <  1:85193/100‑1 (MQ=255)
                                                                     actactGACCAACGCTTCTTTACCTGACATCAGGATTTCTTCATTTTCACCTTTGCAGTAACCAACAATGCAGTCAGCAACATCAATGTGCCGGacaaca                           <  2:967420/100‑1 (MQ=255)
                                                                      ctactGACCAACGCTTCTTTACCTGACATCAGGATTTCTTCATTTTCACCTTTGCAGTAACCAACAATGCAGTCAGCAACATCAATGTGCCGGACAACAt                          <  2:1062402/100‑1 (MQ=255)
                                                                        actGACCAACGCTTCTTTACCTGACATCAGGATTTCTTCATTTTCACCTTTGCAGTAACCAACAATGCAGTCAGCAACATCAATGTGCCGGACAACATCa                        >  1:1548196/1‑100 (MQ=255)
                                                                         ctGACCAACGCTTCTTTACCTGACATCAGGATTTCTTCATTTTCACCTTTGCAGTAACCAACAATGCAGTCAGCAACATCAATGTGCCGGACAACATCaa                       >  1:275350/1‑100 (MQ=255)
                                                                           gACCAACGCTTCTTTACCTGATATCAGGATTTCTTCATTTTCACCTTTGCAGTAACCAACAATGCAGTCAGCAACATCAATGTGCCGGACAACATCAACg                     <  1:470320/100‑1 (MQ=255)
                                                                           gACCAACGCTTCTTTACCTGACATCAGGATTTCTTCATTTTCACCTTTGCAGTAACCAACAATGCAGTCAGCAACATCAATGTGCCGGACAACATCAACg                     >  2:866088/1‑100 (MQ=255)
                                                                           gACCAACGCTTCTTTACCTGACATCAGGATTTCTTCATTTTCACCTTTGCAGTAACCAACAATGCAGTCAGCAACATCAATGTGCCGGACAACATCAACg                     >  2:493693/1‑100 (MQ=255)
                                                                           gACCAACGCTTCTTTACCTGACATCAGGATTTCTTCATTTTCACCTTTGCAGTAACCAACAATGCAGTCAGCAACATCAATGTGCCGGACAACATCAACg                     >  1:1932952/1‑100 (MQ=255)
                                                                           gACCAACGCTTCTTTACCTGACATCAGGATTTCTTCATTTTCACCTTTGCAGTAACCAACAATGCAGTCAGCAACATCAATGTGCCGGACAACATCAACg                     >  2:1462623/1‑100 (MQ=255)
                                                                           gACCAACGCTTCTTTACCTGACATCAGGATTTCTTCATTTTCACCTTTGCAGTAACCAACAATGCAGTCAGCAACATCAATGTGCCGGACAACATCAACg                     <  2:113020/100‑1 (MQ=255)
                                                                           gACCAACGCTTCTTTACCTGACATCAGGATTTCTTCATTTTCACCTTTGCAGTAACCAACAATGCAGTCAGCAACATCAATGTGCCGGACAACATCAACg                     >  2:554961/1‑100 (MQ=255)
                                                                            aCCAACGCTTCTTTACCTGACATCAGGATTTCTTCATTTTCACCTTTGCAGTAACCAACAATGCAGTCAGCAACATCAATGTGCCGGACAACATCAACg                     <  1:1521546/99‑1 (MQ=255)
                                                                             ccAACGCTTCTTTACCTGACATCAGGATTTCTTCATTTTCACCTTTGCAGTAACCAACAATGCAGTCAGCAACATCAATGTGCCGGACAACATCAACGGc                   <  2:61641/100‑1 (MQ=255)
                                                                               aaCGCTTCTTTACCTGACATCAGGATTTCTTCATTTTCACCTTTGCAGTAACCAACAATGCAGTCAGCAACATCAATGTGCCGGACAACATCAACGGCGa                 >  2:345567/1‑100 (MQ=255)
                                                                                 cgatcttttACCTGACATCAGGATTTCTTCATTTTCACCTTTGCAGTAACCAACAATGCAGTCAGCAACATCAATGTGCCGGACAACATCAACGGCGAAg               <  2:1834777/94‑1 (MQ=255)
                                                                                 cGCTTCTTTACCTGACATCAGGATTTCTTCATTTTCACCTTTGCAGTAACCAACAATGCAGTCAGCAACATCAATGTGCCGGACAACATCAACGGCGAAg               <  1:1699437/100‑1 (MQ=255)
                                                                                      cttTACCTGACATCAGGATTTCTTCATTTTCACCTTTGCAGTAACCAACAATGCAGTCAGCAACATCAATGTGCCGGACAACATCAACGGCGAAGTCAGa          >  1:1585213/1‑100 (MQ=255)
                                                                                       tttACCTGACATCAGGATTTCTTCATTTTCACCTTTGCAGTAACCAACAATGCAGTCAGCAACATCAATGTGCCGGACAACATCAACGGCGAAgagatcg         >  1:1834765/1‑94 (MQ=255)
                                                                                       tttACCTGACATCAGGATTTCTTCATTTTCACCTTTGCAGTAACCAACAATGCAGTCAGCAACATCAATGTGCCGGACAACATCAACGGCGAAGTCAGAc         <  2:299679/100‑1 (MQ=255)
                                                                                           ccTGACATCAGGATTTCTTCATTTTCACCTTTGCAGTAACCAACAATGCAGTCAGCAACATCAATGTGCCGGACAACATCAACGGCGAAGTCAGACCAat     <  1:1918995/100‑1 (MQ=255)
                                                                                           ccTGACATCAGGATTTCTTCATTTTCACCTTTGCAGGAACCAACAATGCAGTCAGCAACATCAATGTGCCGGACAACATCAACGGCGAAGTCAGACCAat     >  2:1374426/1‑100 (MQ=255)
                                                                                              gATATCAGGATTTCTTCATTTTCACCTTTGCAGTAACCAACAATGCAGTCAGCAACATCAATGTGCCGGACAACATCAACGGCGAAGTCAGACCAAtatt  >  2:1925828/1‑100 (MQ=255)
                                                                                                |                                                                                                 
TAATAATCAGAGCAATTATTGGGAAGCGTTGGGTAAGGTAGAATGCCCGGATATGTAAGAACAGCAACAACTACTGACCAACGCTTCTTTACCTGATATCAGGATTTCTTCATTTTCACCTTTGCAGTAACCAACAATGCAGTCAGCAACATCAATGTGCCGGACAACATCAACGGCGAAGTCAGACCAATATT  >  NC_000913/3840451‑3840644

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