Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,784,448 T→C intergenic (‑19/+81) ypjB ← / ← ypjC DUF5508 domain‑containing protein YpjB/DUF5507 domain‑containing protein YpjC

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,784,4480TC99.8% 1748.4 / ‑4.6 408intergenic (‑19/+81)ypjB/ypjCDUF5508 domain‑containing protein YpjB/DUF5507 domain‑containing protein YpjC
Reads supporting (aligned to +/- strand):  ref base T (0/1);  new base C (189/218);  total (189/219)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.07e-01

Only 100 of 416 total aligned reads displayed.

ATATCTATGGGGTGTTCCTTTGACTCTGACAAAGATGAAATCTCATCTATTTTGTTTTCATATGAATTACGTGATTCCATACAGACATTTGGCGGCGTTTCTAAAATAACACTACGCGTACTACTTGGCTTAAGTAAACCAACATGAAAATCACTTTTTCTTATATTATCGAAAAGGTTCTATTCATTTCTTTTAG  >  NC_000913/2784350‑2784545
                                                                                                  |                                                                                                 
 tatCTATGGGGTGTTCCTTTGACTCTGACAAAGATGAAATCTCATCTATTTTGTTTTCATATGAATTACGTGATTCCATACAGACATTTGGCGGCGTctct                                                                                                <  2:5161183/101‑1 (MQ=255)
 tatCTATGGGGTGTTCCTTTGACTCTGACAAAGATGAAATCTCATCTATTTTGTTTTCATATGAATTACGTGATTCCATACAGACATTTGGCGGCGTctct                                                                                                >  1:5161975/1‑101 (MQ=255)
 tatCTATGGGGTGTTCCTTTGACTCTGACAAAGATGAAATCTCATCTATTTTGTTTTCATATGAATTACGTGATTCCATACAGACATTTGGCGGCGTctct                                                                                                >  2:5161973/1‑101 (MQ=255)
   tCTATGGGGTGTTCCTTTGACTCTGACAAAGATGAAATCTCATCTATTTTGTTTTCATATGAATTACGTGATTCCATACAGACATTTGGCGGCGTCTCTaa                                                                                              <  1:5161331/101‑1 (MQ=255)
     tATGGGGTGTTCCTTTGACTCTGACAAAGATGAAATCTCATCTATTTTGTTTTCATATGAATTACGTGATTCCATACAGACATTTGGCGGCGTCTCTaaaa                                                                                            >  1:5161990/1‑101 (MQ=255)
     tATGGGGTGTTCCTTTGACTCTGACAAAGATGAAATCTCATCTATTTTGTTTTCATATGAATTACGTGATTCCATACAGACATTTGGCGGCGTCTCTaaaa                                                                                            >  1:5161983/1‑101 (MQ=255)
     tATGGGGTGTTCCTTTGACTCTGACAAAGATGAAATCTCATCTATTTTGTTTTCATATGAATTACGTGATTCCATACAGACATTTGGCGGCGTCTCTaaaa                                                                                            >  2:5161985/1‑101 (MQ=255)
      aTGGGGTGTTCCTTTGACTCTGACAAAGATGAAATCTCATCTATTTTGTTTTCATATGAATTACGTGATTCCATACAGACATTTGGCGGCGTCTCTAAAAt                                                                                           >  2:5161992/1‑101 (MQ=255)
       tGGGGTGTTCCTTTGACTCTGACAAAGATGAAATCTCATCTATTTTGTTTTCATATGAATTACGTGATTCCATACAGACATTTGGCGGCGTCTCTAAAATa                                                                                          <  2:5161534/101‑1 (MQ=255)
        ggggTGTTCCTTTGACTCTGACAAAGATGAAATCTCATCTATTTTGTTTTCATATGAATTACGTGATTCCATACAGACATTTGGCGGCGTCTCTAAAATaa                                                                                         >  2:5161999/1‑101 (MQ=255)
        ggggTGTTCCTTTGACTCTGACAAAGATGAAATCTCATCTATTTTGTTTTCATATGAATTACGTGATTCCATACAGACATTTGGCGGCGTCTCTAAAATaa                                                                                         >  1:5162001/1‑101 (MQ=255)
           gtgtTCCTTTGACTCTGACAAAGATGAAATCTCATCTATTTTGTTTTCATATGAATTACGTGATTCCATACAGACATTTGGCGGCGTCTCTAAAATAacac                                                                                      <  2:5161311/101‑1 (MQ=255)
            tgtTCCTTTGACTCTGACAAAGATGAAATCTCATCTATTTTGTTTTCATATGAATTACGTGATTCCATACAGACATTTGGCGGCGTCTCTAAAATAACACt                                                                                     >  1:5162005/1‑101 (MQ=255)
             gtTCCTTTGACTCTGACAAAGATGAAATCTCATCTATTTTGTTTTCATATGAATTACGTGATTCCATACAGACATTTGGCGGCGTCTCTAAAATAACACTa                                                                                    >  2:5162009/1‑101 (MQ=255)
              ttCCTTTGACTCTGACAAAGATGAAATCTCATCTATTTTGTTTTCATATGAATTACGTGATTCCATACAGACATTTGGCGGCGTCTCTAAAATAACACTAc                                                                                   <  1:5161491/101‑1 (MQ=255)
               tCCTTTGACTCTGACAAAGATGAAATCTCATCTATTTTGTTTTCATATGAATTACGTGATTCCATACAGACATTTGGCGGCGTCTCTAATATAACACTAcg                                                                                  >  2:5162014/1‑101 (MQ=255)
                ccTTTGACTCTGACAAAGATGAAATCTCATCTATTTTGTTTTCATATGAATTACGTGATTCCATACAGACATTTGGCGGCGTCTCTAAAATAACACTAcgc                                                                                 <  1:5161062/101‑1 (MQ=255)
                  tttGACTCTGACAAAGATGAAATCTCATCTATTTTGTTTTCATATGAATTACGTGATTCCATACAGACATTTGGCGGCGTCTCTAAAATAACACTACGCGt                                                                               <  2:5161204/101‑1 (MQ=255)
                   ttGACTCTGACAAAGATGAAATCTCATCTATTTTGTTTTCATATGAATTACGTGATTCCATACAGACATTTGGCGGCGTCTCTAAAATAACACTACGCGta                                                                              <  2:5161390/101‑1 (MQ=255)
                    tGACTCTGACAAAGATGAAATCTCATCTATTTTGTTTTCATATGAATTACGTGATTCCATACAGACATTTGGCGGCGTCTCTAAAATAACACTACGCGtac                                                                             >  1:5162017/1‑101 (MQ=255)
                    tGACTCTGACAAAGATGAAATCTCATCTATTTTGTTTTCATATGAATTACGTGATTCCATACAGACATTTGGCGGCGTCTCTAAAATAACACTACGCGtac                                                                             <  2:5161427/101‑1 (MQ=255)
                     gACTCTGACAAAGATGAAATCTCATCTATTTTGTTTTCATATGAATTACGTGATTCCATACAGACATTTGGCGGCGTCTCTAAAATAACACTACGCGtact                                                                            <  1:5161450/101‑1 (MQ=255)
                      aCTCTGACAAAGATGAAATCTCATCTATTTTGTTTTCATATGAATTACGTGATTCCATACAGACATTTGGCGGCGTCTCTAAAATAACACTACGCGtacta                                                                           >  2:5162022/1‑101 (MQ=255)
                        tctGACAAAGATGAAATCTCATCTATTTTGTTTTCATATGAATTACGTGATTCCATACAGACATTTGGCGGCGTCTCTAAAATAACACTACGCGTactact                                                                         >  2:5162027/1‑101 (MQ=255)
                        tctGACAAAGATGAAATCTCATCTATTTTGTTTTCATATGAATTACGTGATTCCATACAGACATTTGGCGGCGTCTCTAAAATAACACTACGCGTactact                                                                         >  1:5162029/1‑101 (MQ=255)
                         ctGACAAAGATGAAATCTCATCTATTTTGTTTTCATATGAATTACGTGATTCCATACAGACATTTGGCGGCGTCTCTAAAATAACACTACGCGTACTACtt                                                                        <  2:5161365/101‑1 (MQ=255)
                          tGACAAAGATGAAATCTCATCTATTTTGTTTTCATATGAATTACGTGATTCCATACAGACATTTGGCGGCGTCTCTAAAATAACACTACGCGTACTACTTg                                                                       >  1:5162034/1‑101 (MQ=255)
                             cAAAGATGAAATCTCATCTATTTTGTTTTCATATGAATTACGTGATTCCATACAGACATTTGGCGGCGTCTCTAAAATAACACTACGCGTACTACTTGGCt                                                                    >  2:5162040/1‑101 (MQ=255)
                             cAAAGATGAAATCTCATCTATTTTGTTTTCATATGAATTACGTGATTCCATACAGACATTTGGCGGCGTCTCTAAAATAACACTACGCGTACTACTTGGCt                                                                    <  2:5161664/101‑1 (MQ=255)
                              aaaGATGAAATCTCATCTATTTTGTTTTCATATGAATTACGTGATTCCATACAGACATTTGGCGGCGTTTCTAAAATAACACTACGCGTACTACTTGGCtt                                                                   <  1:5161172/101‑1 (MQ=255)
                              aaaGATGAAATCTCATCTATTTTGTTTTCATATGAATTACGTGATTCCATACAGACATTTGGCGGCGTCTCTAAAATAACACTACGCGTACTACTTGGCtt                                                                   <  1:5161352/101‑1 (MQ=255)
                              aaaGATGAAATCTCATCTATTTTGTTTTCATATGAATTACGTGATTCCATACAGACATTTGGCGGCGTCTCTAAAATAACACTACGCGTACTACTTGGCtt                                                                   >  1:5162043/1‑101 (MQ=255)
                               aaGATGAAATCTCATCTATTTTGTTTTCATATGAATTACGTGATTCCATACAGACATTTGGCGGCGTCTCTAAAATAACACTACGCGTACTACTTGGCTTa                                                                  <  1:5161308/101‑1 (MQ=255)
                               aaGATGAAATCTCATCTATTTTGTTTTCATATGAATTACGTGATTCCATACAGACATTTGGCGGCGTCTCTAAAATAACACTACGCGTACTACTTGGCTTa                                                                  >  2:5162048/1‑101 (MQ=255)
                                aGATGAAATCTCATCTATTTTGTTTTCATATGAATTACGTGATTCCATACAGACATTTGGCGGCGTCTCTAAAATAACACTACGCGTACTACTTGGCTTaa                                                                 >  1:5162049/1‑101 (MQ=255)
                                    gAAATCTCATCTATTTTGTTTTCATATGAATTACGTGATTCCATACAGACATTTGGCGGCGTCTCTAAAATAACACTACGCGTACTACTTGGCTTAAGTaa                                                             >  1:5162055/1‑101 (MQ=255)
                                    gAAATCTCATCTATTTTGTTTTCATATGAATTACGTGATTCCATACAGACATTTGGCGGCGTCTCTAAAATAACACTACGCGTACTACTTGGCTTAAGTaa                                                             <  1:5160976/101‑1 (MQ=255)
                                     aaaTCTCATCTATTTTGTTTTCATATGAATTACGTGATTCCATACAGACATTTGGCGGCGTCTCTAAAATAACACTACGCGTACTACTTGGCTTAAGTaaa                                                            >  2:5162057/1‑101 (MQ=255)
                                      aaTCTCATCTATTTTGTTTTCATATGAATTACGTGATTCCATACAGACATTTGGCGGCGTCTCTAAAATAACACTACGCGTACTACTTGGCTTAAGTAAAc                                                           <  2:5161342/101‑1 (MQ=255)
                                      aaTCTCATCTATTTTGTTTTCATATGAATTACGTGATTCCATACAGACATTTGGCGGCGTCTCTAAAATAACACTACGCGTACTACTTGGCTTAAGTAAAc                                                           <  2:5161478/101‑1 (MQ=255)
                                       aTCTCATCTATTTTGTTTTCATATGAATTACGTGATTCCATACAGACATTTGGCGGCGTCTCTAAAATAACACTACGCGTACTACTTGGCTTAAGTAAAcc                                                          >  2:5162061/1‑101 (MQ=255)
                                        tctcATCTATTTTGTTTTCATATGAATTACGTGATTCCATACAGACATTTGGCGGCGTCTCTAAAATAACACTACGCGTACTACTTGGCTTAAGTAAACCa                                                         <  2:5161286/101‑1 (MQ=255)
                                          tcATCTATTTTGTTTTCATATGAATTACGTGATTCCATACAGACATTTGGCGGCGTCTCTAAAATAACACTACGCGTACTACTTGGCTTAAGTAAACCAAc                                                       <  2:5161254/101‑1 (MQ=255)
                                           cATCTATTTTGTTTTCATATGAATTACGTGATTCCATACAGACATTTGGCGGCGTCTCTAAAATAACACTACGCGTACTACTTGGCTTAAGTAAACCAACa                                                      >  1:5162069/1‑101 (MQ=255)
                                           cATCTATTTTGTTTTCATATGAATTACGTGATTCCATACAGACATTTGGCGGCGTCTCTAAAATAACACTACGCGTACTACTTGGCTTAAGTAAACCAACa                                                      <  2:5161383/101‑1 (MQ=255)
                                           cATCTATTTTGTTTTCATATGAATTACGTGATTCCATACAGACATTTGGCGGCGTCTCTAAAATAACACTACGCGTACTACTTGGCTTAAGTAAACCAACa                                                      <  2:5161350/101‑1 (MQ=255)
                                             tCTATTTTGTTTTCATATGAATTACGTGATTCCATACAGACATTTGGCGGCGTCTCTAAAATAACACTACGCGTACTACTTGGCTTAAGTAAACCAACATg                                                    >  2:5162073/1‑101 (MQ=255)
                                              cTATTTTGTTTTCATATGAATTACGTGATTCCATACAGACATTTGGCGGCGTCTCTAAAATAACACTACGCGTACTACTTGGCTTAAGTAAACCAACATGa                                                   <  2:5161457/101‑1 (MQ=255)
                                              cTATTTTGTTTTCATATGAATTACGTGATTCCATACAGACATTTGGCGGCGTCTCTAAAATAACACTACGCGTACTACTTGGCTTAAGTAAACCAACATGa                                                   >  1:5162076/1‑101 (MQ=255)
                                               tATTTTGTTTTCATATGAATTACGTGATTCCATACAGACATTTGGCGGCGTCTCTAAAATAACACTACGCGTACTACTTGGCTTAAGTAAACCAACATGaa                                                  <  1:5161284/101‑1 (MQ=255)
                                               tATTTTGTTTTCATATGAATTACGTGATTCCATACAGACATTTGGCGGCGTCTCTAAAATAACACTACGCGTACTACTTGGCTTAAGTAAACCAACATGaa                                                  <  2:5161732/101‑1 (MQ=255)
                                                 ttttGTTTTCATATGAATTACGTGATTCCATACAGACATTTGGCGGCGTCTCTAAAATAACACTACGCGTACTACTTGGCTTAAGTAAACCAACATGaaaa                                                >  2:5162083/1‑101 (MQ=255)
                                                  tttGTTTTCATATGAATTACGTGATTCCATACAGACATTTGGCGGCGTCTCTAAAATAACACTACGCGTACTACTTGGCTTAAGTAAACCAACATGAAAAt                                               >  1:5162085/1‑101 (MQ=255)
                                                   ttGTTTTCATATGAATTACGTGATTCCATACAGACATTTGGCGGCGTCTCTAAAATAACACTACGCGTACTACTTGGCTTAAGTAAACCAACATGAAAATc                                              <  1:5161423/101‑1 (MQ=255)
                                                   ttGTTTTCATATGAATTACGTGATTCCATACAGACATTTGGCGGCGTCTCTAAAATAACACTACGCGTACTACTTGGCTTAAGTAAACCAACATGAAAATc                                              <  1:5161380/101‑1 (MQ=255)
                                                   ttGTTTTCATATGAATTACGTGATTCCATACAGACATTTGGCGGCGTCTCTAAAATAACACTACGCGTACTACTTGGCTTAAGTAAACCAACATGAAAATc                                              <  1:5161675/101‑1 (MQ=255)
                                                    tGTTTTCATATGAATTACGTGATTCCATACAGACATTTGGCGGCGTCTCTAAAATAACACTACGCGTACTACTTGGCTTAAGTAAACCAACATGAAAATCa                                             >  2:5162089/1‑101 (MQ=255)
                                                     gTTTTCATATGAATTACGTGATTCCATACAGACATTTGGCGGCGTCTCTAAAATAACACTACGCGTACTACTTGGCTTAAGTAAACCAACATGAAAATCAc                                            <  2:5161468/101‑1 (MQ=255)
                                                       tttCATATGAATTACGTGATTCCATACAGACATTTGGCGGCGTCTCTAAAATAACACTACGCGTACTACTTGGCTTAAGTAAACCAACATGAAAATCACtt                                          >  2:5162095/1‑101 (MQ=255)
                                                       tttCATATGAATTACGTGATTCCATACAGACATTTGGCGGCGTCTCTAAAATAACACTACGCGTACTACTTGGCTTAAGTAAACCAACATGAAAATCACtt                                          <  1:5161404/101‑1 (MQ=255)
                                                       tttCATATGAATTACGTGATTCCATACAGACATTTGGCGGCGTCTCTAAAATAACACTACGCGTACTACTTGGCTTAAGTAAACCAACATGAAAATCACtt                                          <  2:5160909/101‑1 (MQ=255)
                                                        ttCATATGAATTACGTGATTCCATACAGACATTTGGCGGCGTCTCTAAAATAACACTACGCGTACTACTTGGCTTAAGTAAACCAACATGAAAATCACttt                                         <  2:5161357/101‑1 (MQ=255)
                                                         tCATATGAATTACGTGATTCCATACAGACATTTGGCGGCGTCTCTAAAATAACACTACGCGTACTACTTGGCTTAAGTAAACCAACATGAAAATCACtttt                                        >  1:5162103/1‑101 (MQ=255)
                                                         tCATATGAATTACGTGATTCCATACAGACATTTGGCGGCGTCTCTAAAATAACACTACGCGTACTACTTGGCTTAAGTAAACCAACATGAAAATCACtttt                                        >  2:5162099/1‑101 (MQ=255)
                                                            tatGAATTACGTGATTCCATACAGACATTTGGCGGCGTCTCTAAAATAACACTACGCGTACTACTTGGCTTAAGTAAACCAACATGAAAATCACTTTTTCt                                     <  1:5161515/101‑1 (MQ=255)
                                                               gAATTACGTGATTCCATACAGACATTTGGCGGCGTCTCTAAAATAACACTACGCGTACTACTTGGCTTAAGTAAACCAACATGAAAATCACTTTTTCTtat                                  >  2:5162109/1‑101 (MQ=255)
                                                                  ttACGTGATTCCATACAGACATTTGGCGGCGTCTCTAAAATAACACTACGCGTACTACTTGGCTTAAGTAAACCAACATGAAAATCACTTTTTCTTAtatt                               >  2:5162115/1‑101 (MQ=255)
                                                                  ttACGTGATTCCATACAGACATTTGGCGGCGTCTCTAAAATAACACTACGCGTACTACTTGGCTTAAGTAAACCAACATGAAAATCACTTTTTCTTAtatt                               >  1:5162114/1‑101 (MQ=255)
                                                                    aCGTGATTCCATACAGACATTTGGCGGCGTCTCTAAAATAACACTACGCGTACTACTTGGCTTAAGTAAACCAACATGAAAATCACTTTTTCTTAtattat                             >  2:5162121/1‑101 (MQ=255)
                                                                     cGTGATTCCATACAGACATTTGGCGGCGTCTCTAAAATAACACTACGCGTACTACTTGGCTTAAGTAAACCAACATGAAAATCACTTTTTCTTATATTATc                            <  1:5161540/101‑1 (MQ=255)
                                                                     cGTGATTCCATACAGACATTTGGCGGCGTCTCTAAAATAACACTACGCGTACTACTTGGCTTAAGTAAACCAACATGAAAATCACTTTTTCTTATATTATc                            <  1:5161564/101‑1 (MQ=255)
                                                                      gTGATTCCATACAGACATTTGGCGGCGTCTCTAAAATAACACTACGCGTACTACTTGGCTTAAGTAAACCAACATGAAAATCACTTTTTCTTATATTATCg                           >  2:5162126/1‑101 (MQ=255)
                                                                       tGATTCCATACAGACATTTGGCGGCGTCTCTAAAATAACACTACGCGTACTACTTGGCTTAAGTAAACCAACATGAAAATCACTTTTTCTTATATTATCGa                          >  1:5162128/1‑101 (MQ=255)
                                                                         aTTCCATACAGACATTTGGCGGCGTCTCTAAAATAACACTACGCGTACTACTTGGCTTAAGTAAACCAACATGAAAATCACTTTTTCTTATATTATCGaaa                        <  1:5161253/101‑1 (MQ=255)
                                                                         aTTCCATACAGACATTTGGCGGCGTCTCTAAAATAACACTACGCGTACTACTTGGCTTAAGTAAACCAACATGAAAATCACTTTTTCTTATATTATCGaaa                        >  2:5162133/1‑101 (MQ=255)
                                                                         aTTCCATACAGACATTTGGCGGCGTCTCTAAAATAACACTACGCGTACTACTTGGCTTAAGTAAACCAACATGAAAATCACTTTTTCTTATATTATCGaaa                        <  2:5161398/101‑1 (MQ=255)
                                                                         aTTCCATACAGACATTTGGCGGCGTCTCTAAAATAACACTACGCGTACTACTTGGCTTAAGTAAACCAACATGAAAATCACTTTTTCTTATATTATCGaaa                        <  1:5161742/101‑1 (MQ=255)
                                                                          ttCCATACAGACATTTGGCGGCGTCTCTAAAATAACACTACGCGTACTACTTGGCTTAAGTAAACCAACATGAAAATCACTTTTTCTTATATTATCGaaaa                       <  1:5161119/101‑1 (MQ=255)
                                                                           tCCATACAGACATTTGGCGGCGTCTCTAAAATAACACTACGCGTACTACTTGGCTTAAGTAAACCAACATGAAAATCACTTTTTCTTATATTATCGAAAAg                      <  1:5161934/101‑1 (MQ=255)
                                                                             cATACAGACATTTGGCGGCGTCTCTAAAATAACACTACGCGTACTACTTGGCTTAAGTAAACCAACATGAAAATCACTTTTTCTTATATTATCGAAAAGGt                    <  2:5162067/101‑1 (MQ=255)
                                                                             cATACAGACATTTGGCGGCGTCTCTAAAATAACACTACGCGTACTACTTGGCTTAAGTAAACCAACATGAAAATCACTTTTTCTTATATTATCGAAAAGGt                    >  2:5162138/1‑101 (MQ=255)
                                                                              aTACAGACATTTGGCGGCGTCTCTAAAATAACACTACGCGTACTACTTGGCTTAAGTAAACCAACATGAAAATCACTTTTTCTTATATTATCGAAAAGGtt                   <  2:5161117/101‑1 (MQ=255)
                                                                               tACAGACATTTGGCGGCGTCTCTAAAATAACACTACGCGTACTACTTGGCTTAAGTAAACCAACATGAAAATCACTTTTTCTTATATTATCGAAAAGGTTc                  >  1:5162140/1‑101 (MQ=255)
                                                                               tACAGACATTTGGCGGCGTCTCTAAAATAACACTACGCGTACTACTTGGCTTAAGTAAACCAACATGAAAATCACTTTTTCTTATATTATCGAAAAGGTTc                  <  2:5161511/101‑1 (MQ=255)
                                                                                aCAGACATTTGGCGGCGTCTCTAAAATAACACTACGCGTACTACTTGGCTTAAGTAAACCAACATGAAAATCACTTTTTCTTATATTATCGAAAAGGTTCt                 <  2:5161579/101‑1 (MQ=255)
                                                                                 cAGACATTTGGCGGCGTCTCTAAAATAACACTACGCGTACTACTTGGCTTAAGTAAACCAACATGAAAATCACTTTTTCTTATATTATCGAAAAGGTTCTa                >  1:5162146/1‑101 (MQ=255)
                                                                                 cAGACATTTGGCGGCGTCTCTAAAATAACACTACGCGTACTACTTGGCTTAAGTAAACCAACATGAAAATCACTTTTTCTTATATTATCGAAAAGGTTCTa                <  2:5161093/101‑1 (MQ=255)
                                                                                  aGACATTTGGCGGCGTCTCTAAAATAACACTACGCGTACTACTTGGCTTAAGTAAACCAACATGAAAATCACTTTTTCTTATATTATCGAAAAGGTTCTAt               <  2:5161599/101‑1 (MQ=255)
                                                                                   gACATTTGGCGGCGTCTCTAAAATAACACTACGCGTACTACTTGGCTTAAGTAAACCAACATGAAAATCACTTTTTCTTATATTATCGAAAAGGTTCTAtt              <  1:5161602/101‑1 (MQ=255)
                                                                                   gACATTTGGCGGCGTCTCTAAAATAACACTACGCGTACTACTTGGCTTAAGTAAACCAACATGAAAATCACTTTTTCTTATATTATCGAAAAGGTTCTAtt              <  2:5161818/101‑1 (MQ=255)
                                                                                    aCATTTGGCGGCGTCTCTAAAATAACACTACGCGTACTACTTGGCTTAAGTAAACCAACATGAAAATCACTTTTTCTTATATTATCGAAAAGGTTCTATTc             <  1:5161527/101‑1 (MQ=255)
                                                                                    aCATTTGGCGGCGTCTCTAAAATAACACTACGCGTACTACTTGGCTTAAGTAAACCAACATGAAAATCACTTTTTCTTATATTATCGAAAAGGTTCTATTc             >  2:5162151/1‑101 (MQ=255)
                                                                                      aTTTGGCGGCGTCTCTAAAATAACACTACGCGTACTACTTGGCTTAAGTAAACCAACATGAAAATCACTTTTTCTTATATTATCGAAAAGGTTCTATTCAt           >  1:5162155/1‑101 (MQ=255)
                                                                                        ttGGCGGCGTCTCTAAAATAACACTACGCGTACTACTTGGCTTAAGTAAACCAACATGAAAATCACTTTTTCTTATATTATCGAAAAGGTTCTATTCAttt         >  2:5162158/1‑101 (MQ=255)
                                                                                        ttGGCGGCGTCTCTAAAATAACACTACGCGTACTACTTGGCTTAAGTAAACCAACATGAAAATCACTTTTTCTTATATTATCGAAAAGGTTCTATTCAttt         <  1:5161684/101‑1 (MQ=255)
                                                                                        ttGGCGGCGTCTCTAAAATAACACTACGCGTACTACTTGGCTTAAGTAAACCAACATGAAAATCACTTTTTCTTATATTATCGAAAAGGTTCTATTCAttt         <  2:5161443/101‑1 (MQ=255)
                                                                                         tGGCGGCGTCTCTAAAATAACACTACGCGTACTACTTGGCTTAAGTAAACCAACATGAAAATCACTTTTTCTTATATTATCGAAAAGGTTCTATTCATTTc        <  1:5161643/101‑1 (MQ=255)
                                                                                           gcggcgTCTCTAAAATAACACTACGCGTACTACTTGGCTTAAGTAAACCAACATGAAAATCACTTTTTCTTATATTATCGAAAAGGTTCTATTCATTTCtt      <  2:5161322/101‑1 (MQ=255)
                                                                                            cggcgTCTCTAAAATAACACTACGCGTACTACTTGGCTTAAGTAAACCAACATGAAAATCACTTTTTCTTATATTATCGAAAAGGTTCTATTCATTTCttt     <  2:5161630/101‑1 (MQ=255)
                                                                                             ggcgTCTCTAAAATAACACTACGCGTACTACTTGGCTTAAGTAAACCAACATGAAAATCACTTTTTCTTATATTATCGAAAAGGTTCTATTCATTTCtttt    >  1:5162164/1‑101 (MQ=255)
                                                                                                  |                                                                                                 
ATATCTATGGGGTGTTCCTTTGACTCTGACAAAGATGAAATCTCATCTATTTTGTTTTCATATGAATTACGTGATTCCATACAGACATTTGGCGGCGTTTCTAAAATAACACTACGCGTACTACTTGGCTTAAGTAAACCAACATGAAAATCACTTTTTCTTATATTATCGAAAAGGTTCTATTCATTTCTTTTAG  >  NC_000913/2784350‑2784545

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