Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,151,557 C→T intergenic (+119/‑21) metC → / → yghB cystathionine beta‑lyase/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,151,5570CT99.2% 1333.1 / ‑2.6 378intergenic (+119/‑21)metC/yghBcystathionine beta‑lyase/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (2/1);  new base T (190/185);  total (192/186)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.17e-01

Only 100 of 386 total aligned reads displayed.

CATTTTATTGTGAAAAGTCTTAAATTGTTGCGTCCGGGATCAAGGCGTCCCGGACGATTCAGGAGTACAATAGGCAGATAAAGGCTTAAACGCTGTTCCACAGGAAAGTCCATGGCTGTTATTCAAGATATCATCGCTGCGCTCTGGCAACACGACTTTGCCGCGCTGGCGGATCCTCAT  >  NC_000913/3151467‑3151646
                                                                                          |                                                                                         
cATTTTATTGTGAAAAGTCTTAAATTGTTGCGTCCGGGATCAAGGCGTCCCGGACGATTCAGGAGTACAATAGGCAGATAAAGGCTTAAATGc                                                                                         >  1:3626217/1‑93 (MQ=255)
cATTTTATTGTGAAAAGTCTTAAATTGTTGCGTCCGGGATCAAGGCGTCCCGGACGATTCAGGAGTACAATAGGCAGATAAAGGCTTAAATGc                                                                                         >  1:5080892/1‑93 (MQ=255)
cATTTTATTGTGAAAAGTCTTAAATTGTTGCGTCCGGGATCAAGGCGTCCCGGACGATTCAGGAGTACAATAGGCAGATAAAGGCTTAAATGc                                                                                         >  1:1690504/1‑93 (MQ=255)
  ttttATTGTGAAAAGTCTTAAATTGTTGCGTCCGGGATCAAGGCGTCCCGGACGATTCAGGAGTACAATAGGCAGATAAAGGCTTAAATGCTg                                                                                       >  1:4149688/1‑93 (MQ=255)
  ttttATTGTGAAAAGTCTTAAATTGTTGCGTCCGGGATCAAGGCGTCCCGGACGATTCAGGAGTACAATAGGCAGATAAAGGCTTAAATGCTg                                                                                       >  2:4512089/1‑93 (MQ=255)
   tttATTGTGAAAAGTCTTAAATTGTTGCGTCCGGGGTCAAGGCGTCCCGGACGATTCAGGAGTACAATAGGCAGATAAAGGCTTAAATGCTGt                                                                                      >  2:681797/1‑93 (MQ=255)
   tttATTGTGAAAAGTCTTAAATTGTTGCGTCCGGGATCAAGGCGTCCCGGACGATTCAGGAGTACAATAGGCAGATAAAGGCTTAAATGCTGt                                                                                      <  1:3373577/93‑1 (MQ=255)
    ttATTGTGAAAAGTCTTAAATTGTTGCGTCCGGGATCAAGGCGTCCCGGACGATTCAGGAGTACAATAGGCAGATAAAGGCTTAAATGCTGtt                                                                                     <  2:2989809/93‑1 (MQ=255)
     tATTGTGAAAAGTCTTAAATTGTTGCGTCCGGGATCAAGGCGTCCCGGACGATTCAGGAGTACAATAGGCAGATAAAGGCTTAAATGCTGTTc                                                                                    <  2:899755/93‑1 (MQ=255)
       ttGTGAAAAGTCTTAAATTGTTGCGTCCGGGATCAAGGCGTCCCGGACGATTCAGGAGTACAATAGGCAGATAAAGGCTTAAATGCTGTTCca                                                                                  <  2:154765/93‑1 (MQ=255)
        tgtgAAAAGTCTTAAATTGTTGCGTCCGGGATCAAGGCGTCCCGGACGATTCAGGAGTACAATAGGCAGATAAAGGCTTAAATGCTGTTCcac                                                                                 <  2:2535321/93‑1 (MQ=255)
         gtgAAAAGTCTTAAATTGTTGCGTCCGGGATCAAGGCGTCCCGGACGATTCAGGAGTACAATAGGCAGATAAAGGCTTAAATGCTGTTCcaca                                                                                >  2:7105154/1‑93 (MQ=255)
         gtgAAAAGTCTTAAATTGTTGCGTCCGGGATCAAGGCGTCCCGGACGATTCAGGAGTACAATAGGCAGATAAAGGCTTAAATGCTGTTCcaca                                                                                >  1:1431433/1‑93 (MQ=255)
          tgAAAAGTCTTAAATTGTTGCGTCCGGGATCAAGGCGTCCCGGACGATTCAGGAGTACAATAGGCAGATAAAGGCTTAAATGCTGTTCCACAg                                                                               >  2:7362781/1‑93 (MQ=255)
          tgAAAAGTCTTAAATTGTTGCGTCCGGGATCAAGGCGTCCCGGACGATTCAGGAGTACAATAGGCAGATAAAGGCTTAAATGCTGTTCCACAg                                                                               >  1:1246411/1‑93 (MQ=255)
           gAAAAGTCTTAAATTGTTGCGTCCGGGATCAAGGCGTCCCGGACGATTCAGGAGTACAATAGGCAGATAAAGGCTTAAATGCTGTTCCACAgg                                                                              >  2:5033384/1‑93 (MQ=255)
           gAAAAGTCTTAAATTGTTGCGTCCGGGATCAAGGCGTCCCGGACGATTCAGGAGTACAATAGGCAGATAAAGGCTTAAATGCTGTTCCACAgg                                                                              >  2:339085/1‑93 (MQ=255)
           gAAAAGTCTTAAATTGTTGCGTCCGGGATCAAGGCGTCCCGGACGATTCAGGAGTACAATAGGCAGATAAAGGCTTAAATGCTGTTCCACAgg                                                                              >  1:3751435/1‑93 (MQ=255)
            aaaaGTCTTAAATTGTTGCGTCCGGGATCAAGGCGTCCCGGACGATTCAGGAGTACAATAGGCAGATAAAGGCTTAAATGCTGTTCCACAGGa                                                                             >  1:7287298/1‑93 (MQ=255)
             aaaGTCTTAAATTGTTGCGTCCGGGATCAAGGCGTCCCGGACGATTCAGGAGTACAATAGGCAGATAAAGGCTTAAATGCTGTTCCACAGGaa                                                                            >  1:8646021/1‑93 (MQ=255)
             aaaGTCTTAAATTGTTGCGTCCGGGATCAAGGCGTCCCGGACGATTCAGGAGTACAATAGGCAGATAAAGGCTTAAATGCTGTTCCACAGGaa                                                                            <  1:6959899/93‑1 (MQ=255)
             aaaGTCTTAAATTGTTGCGTCCGGGATCAAGGCGTCCCGGACGATTCAGGAGTACAATAGGCAGATAAAGGCTTAAATGCTGTTCCACAGGaa                                                                            >  1:4164549/1‑93 (MQ=255)
              aaGTCTTAAATTGTTGCGTCCGGGATCAAGGCGTCCCGGACGATTCAGGAGTACAATAGGCAGATAAAGGCTTAAATGCTGTTCCACAGGaaa                                                                           >  1:1125448/1‑93 (MQ=255)
               aGTCTTAAATTGTTGCGTCCGGGATCAAGGCGTCCCGGACGATTCAGGAGTACAATAGGCAGATAAAGGCTTAAATGCTGTTCCACAGGAAAg                                                                          <  2:3538498/93‑1 (MQ=255)
                gTCTTAAATTGTTGCGTCCGGGGTCAAGGCGTCCCCGACGGTTTAGGAGTACAATAGGCAGATAAAGGCTTAAATGCTGGTACACAGGAAAGt                                                                         >  2:4709268/1‑93 (MQ=255)
                gTCTTAAATTGTTGCGTCCGGGATCAAGGCGTCCCGGACGATTCAGGAGTACAATAGGCAGATAAAGGCTTAAATGCTGTTCCACAGGAAAGt                                                                         >  2:7475832/1‑93 (MQ=255)
                 tCTTAAATTGTTGCGTCCGGGATCAAGGCGTCCCGGACGATTCAGGAGTACAATAGGCAGATAAAGGCTTAAATGCTGTTCCACAGGAAAGTc                                                                        >  1:8415576/1‑93 (MQ=255)
                 tCTTAAATTGTTGCGTACGGGATCAAGGCGTCCCGGACGATTCAGGAGTACAATAGGCAGATAAAGGCTTAAATGCTGTTCCACAGGAAAGTc                                                                        <  1:7760712/93‑1 (MQ=255)
                   ttAAATTGTTGCGTCCGGGATCAAGGCGTCCCGGACGATTCAGGAGTACAATAGGCAGATAAAGGCTTAAATGCTGTTCCACAGGAAAGTCCa                                                                      >  2:5162285/1‑93 (MQ=255)
                   ttAAATTGTTGCGTCCGGGATCAAGGCGTCCCGGACGATTCAGGAGTACAATAGGCAGATAAAGGCTTAAATGCTGTTCCACAGGAAAGTCCa                                                                      <  1:2344463/93‑1 (MQ=255)
                    tAAATTGTTGCGTCCGGGATCAAGGCGTCCCGGACGATTCAGGAGTACAATAGGCAGATAAAGGCTTAAATGCTGTTCCACAGGAAAGTCCAt                                                                     <  1:8092673/93‑1 (MQ=255)
                     aaaTTGTTGCGTCCGGGATCAAGGCGTCCCGGACGATTCAGGAGTACAATAGGCAGATAAAGGCTTAAATGCTGTTCCACAGGAAAGTCCATg                                                                    <  2:6562528/93‑1 (MQ=255)
                     aaaTTGTTGCGTCCGGGATCAAGGCGTCCCGGACGATTCAGGAGTACAATAGGCAGATAAAGGCTTAAATGCTGTTCCACAGGAAAGTCCATg                                                                    <  2:2717554/93‑1 (MQ=255)
                     aaaTTGTTGCGTCCGGGATCAAGGCGTCCCGGACGATTCAGGAGTACAATAGGCAGATAAAGGCTTAAATGCTGTTCCACAGGAAAGTCCATg                                                                    >  1:6654366/1‑93 (MQ=255)
                     aaaTTGTTGCGTCCGGGATCAAGGCGTCCCGGACGATTCAGGAGTACAATAGGCAGATAAAGGCTTAAATGCTGTTCCACAGGAAAGTCCATg                                                                    <  1:4230024/93‑1 (MQ=255)
                       aTTGTTGCGTCCGGGATCAAGGCGTCCCGGACGATTCAGGAGTACAATAGGCAGATAAAGGCTTAAATGCTGTTCCACAGGAAAGTCCATGGc                                                                  >  2:4096316/1‑93 (MQ=255)
                        ttgttgCGGTCGGGATCAAGCCGTCCCGGACGATTCAGGAGTACAATAGTCAGATAAAGGCTTAAATTCTGTTCCACAGGAAAGTCCATGGCt                                                                 <  2:7498285/93‑1 (MQ=255)
                         tgttgCGTCCGGGATCAAGGCGTCCCGGACGATTCAGGAGTACAATAGGCAGATAAAGGCTTAAATGCTGTTCCACAGGAAAGTCCATGGCTg                                                                <  1:1034559/93‑1 (MQ=255)
                         tgttgCGTCCGGGATCAAGGCGTCCCGGACGATTCAGGAGTACAATAGGCAGATAAAGGCTTAAATGCTGTTCCACAGGAAAGTCCATGGCTg                                                                <  2:3146306/93‑1 (MQ=255)
                         tgttgCGTCCGGGATCAAGGCGTCCCGGACGATTCAGGAGTACAATAGGCAGATAAAGGCTTAAATGCTGTTCCACAGGAAAGTCCATGGCTg                                                                <  2:2850013/93‑1 (MQ=255)
                         tgttgCGTCCGGGATCAAGGCGTCCCGGACGATTCAGGAGTACAATAGGCAGATAAAGGCTTAAATGCTGTTCCACAGGAAAGTCCATGGCTg                                                                <  1:8331337/93‑1 (MQ=255)
                          gttgCGTCCGGGATCAAGGCGTCCCGGACGATTCAGGAGTACAATAGGCAGATAAAGGCTTAAATGCTGTTCCACAGGAAAGTCCATGGCTGt                                                               <  1:2672989/93‑1 (MQ=255)
                              cGTCCGGGATCAAGGCGTCCCGGACGATTCAGGAGTACAATAGGCAGATAAAGGCTTAAATGCTGTTCCACAGGAAAGTCCATGGCTGTTAtt                                                           >  1:4344536/1‑93 (MQ=255)
                                tCCGGGATCAAGGCGTCCCGGACGATTCAGGAGTACAATAGGCAGATAAAGGCTTAAATGCTGTTCCACAGGAAAGTCCATGGCTGTTATTCa                                                         >  2:5623558/1‑93 (MQ=255)
                                 ccGGGATCAAGGCGTCCCGGACGATTCAGGAGTACAATAGGCAGATAAAGGCTTAAATGCTGTTCCACAGGAAAGTCCATGGCTGTTATTCaa                                                        <  1:5730592/93‑1 (MQ=255)
                                 ccGGGATCAAGGCGTCCCGGACGATTCAGGAGTACAATAGGCAGATAAAGGCTTAAATGCTGTTCCACAGGAAAGTCCATGGCTGTTATTCaa                                                        >  2:3946999/1‑93 (MQ=255)
                                 ccGGGATCAAGGCGTCCCGGACGATTCAGGAGTACAATAGGCAGATAAAGGCTTAAATGCTGTTCCACAGGAAAGTCCATGGCTGCTATTCaa                                                        <  2:5323223/93‑1 (MQ=255)
                                   gggATCAAGGCGTCCCGGACGATTCAGGAGTACAATAGGCAGATAAAGGCTTAAATGCTGTTCCACAGGAAAGTCCATGGCTGTTATTCAAGa                                                      <  2:8621943/93‑1 (MQ=255)
                                   gggATCAAGGCGTCCCGGACGATTCAGGAGTACAATAGGCAGATAAAGGCTTAAATGCTGTTCCACAGGAAAGTCCATGGCTGTTATTCAAGa                                                      <  2:6237369/93‑1 (MQ=255)
                                   gggATCAAGGCGTCCCGGACGATTCAGGAGTACAATAGGCAGATAAAGGCTTAAATGCTGTTCCACAGGAAAGTCCATGGCTGTTATTCAAGa                                                      >  1:555277/1‑93 (MQ=255)
                                    ggATCAAGGCGTCCCGGACGATTCAGGAGTACAATAGGCAGATAAAGGCTTAAATGCTGTTCCACAGGAAAGTCCATGGCTGTTATTCAAGat                                                     >  1:2074481/1‑93 (MQ=255)
                                     gATCAAGGCGTCCCGGACGATTCAGGAGTACAATAGGCAGATAAAGGCTTAAATGCTGTTCCACAGGAAAGTCCATGGCTGTTATTCAAGata                                                    <  2:2253370/93‑1 (MQ=255)
                                     gATCAAGGCGTCCCGGACGATTCAGGAGTACAATAGGCAGATAAAGGCTTAAATGCTGTTCCACAGGAAAGTCCATGGCTGTTATTCAAGata                                                    >  2:1310546/1‑93 (MQ=255)
                                      aTCAAGGCGTCCCGGACGATTCAGGAGTACAATAGGCAGATAAAGGCTTAAATGCTGTTCCACAGGAAAGTCCATGGCTGTTATTCAAGatat                                                   >  2:8066780/1‑93 (MQ=255)
                                      aTCAAGGCGTCCCGGACGATTCAGGAGTACAATAGGCAGATAAAGGCTTAAATGCTGTTCCACAGGAAAGTCCATGGCTGTTATTCAAGatat                                                   <  1:600415/93‑1 (MQ=255)
                                       tCAAGGCGTCCCGGACGATTCAGGAGTACAATAGGCAGATAAAGGCTTAAATGCTGTTCCACAGGAAAGTCCATGGCTGTTATTCAAGATatc                                                  >  2:4871170/1‑93 (MQ=255)
                                       tCAAGGCGTCCCGGACGATTCAGGAGTACAATAGGCAGATAAAGGCTTAAATGCTGTTCCACAGGAAAGTCCATGGCTGTTATTCAAGATatc                                                  <  2:8499435/93‑1 (MQ=255)
                                        cAAGGCGTCCCGGACGATTCAGGAGTACAATAGGCAGATAAAGGCTTAAATGCTGTTCCACAGGAAAGTCCATGGCTGTTATTCAAGATatca                                                 <  1:8283958/93‑1 (MQ=255)
                                         aaGGCGTCCCGGACGATTCAGGAGTACAATAGGCAGATAAAGGCTTAAATGCTGTTCCACAGGAAAGTCCATGGCTGTTATTCAAGATatcat                                                >  1:3540085/1‑93 (MQ=255)
                                         aaGGCGTCCCGGACGATTCAGGAGTACAATAGGCAGATAAAGGCTTAAATGCTGTTCCACAGGAAAGTCCATGGCTGTTATTCAAGATatcat                                                <  1:4519712/93‑1 (MQ=255)
                                          aGGCGTCCCGGACGATTCAGGAGTACAATAGGCAGATAAAGGCTTAAATGCTGTTCCACAGGAAAGTCCATGGCTGTTATTCAAGATatcatc                                               >  2:2304305/1‑93 (MQ=255)
                                          aGGCGTCCCGGACGATTCAGGAGTACAATAGGCAGATAAAGGCTTAAATGCTGTTCCACAGGAAAGTCCATGGCTGTTATTCAAGATatcatc                                               >  1:7861494/1‑93 (MQ=255)
                                            gCGTCCCGGACGATTCAGGAGTACAATAGGCAGATAAAGGCTTAAATGCTGTTCCACAGGAAAGTCCATGGCTGTTATTCAAGATATCATCGc                                             <  1:2532181/93‑1 (MQ=255)
                                               tCCCGTACGATTCAGGAGTACAATAGGCAGATAAAGGCTTAAATGCTGTTCCACAGGAAAGTCCATGGCTGTTATTCAAGATATCATCGCTgc                                          <  2:700128/93‑1 (MQ=255)
                                               tCCCGGACGCTTCAGGAGTAAAATAGGCAGATAAAGGCTTAAATGCTGTTCCACAGGACAGTACATGGCTGTTATCAAAGATATAATGGTTgc                                          >  1:3072300/1‑93 (MQ=255)
                                               tCCCGGACGATTCAGGAGTACAATAGGCAGATAAAGGCTTAAATGCTGTTCCACAGGAAAGTCCATGGCTGTTATTCAAGATATCATCGCTgc                                          >  1:4710000/1‑93 (MQ=255)
                                               tCCCGGACGATTCAGGAGTACAATAGGCAGATAAAGGCTTAAATGCTGTTCCACAGGAAAGTCCATGGCTGTTATTCAAGATATCATCGCTgc                                          >  1:3228891/1‑93 (MQ=255)
                                               tCCCGGACGATTCAGGAGGACAATAGGCAGATAAAGGCTTAAATGCTGTTCCACAGGAAAGTCCATGGCTGTTATTCAAGATATAATCGCTgc                                          >  2:8413451/1‑93 (MQ=255)
                                                 ccGGACGATTCAGGAGTACAATAGGCAGATAAAGGCTTAAATGCTGTTCCACAGGAAAGTCCATGGCTGTTATTCAAGATATCATCGCTgcgc                                        <  1:5957635/93‑1 (MQ=255)
                                                 ccGGACGATTCAGGAGTACAATAGGCAGATAAAGGCTTAAATGCTGTTCCACAGGAAAGTCCATGGCTGTTATTCAAGATATCATCGCTgcgc                                        >  1:2414396/1‑93 (MQ=255)
                                                   ggACGATTCAGGAGTACAATAGGCAGATAAAGGCTTAAATGCTGTTCCACAGGAAAGTCCATGGCTGTTATTCAAGATATCATCGCTGCGctc                                      <  1:2944556/93‑1 (MQ=255)
                                                      cGATTCAGGAGTACAATAGGCAGATAAAGGCTTAAATGCTGTTCCACAGGAAAGTCCATGGCTGTTATTCAAGATATCATCGCTGCGCTCTgg                                   >  2:7267305/1‑93 (MQ=255)
                                                       gATTCAGGAGTACAATAGGCAGATAAAGGCTTAAATGCTGTTCCACAGGAAAGTCCATGGCTGTTATTCAAGATATCATCGCTGCGCTCTGGc                                  >  1:1517991/1‑93 (MQ=255)
                                                         ttCAGGAGTACAATAGGCAGATAAAGGCTTAAATGCTGTTCCACAGGAAAGTCCATGGCTGTTATTCAAGATATCATCGCTGCGCTCTGGCac                                >  1:7454790/1‑92 (MQ=255)
                                                           cAGGAGTACAATAGGCAGATAAAGGCTTAAATGCTGTTCCACAGGAAAGTCCATGGCTGTTATTCAAGATATCATCGCTGCGCTCTGGCAaca                              <  1:6343733/93‑1 (MQ=255)
                                                            aGGAGTACAATAGGCAGATAAAGGCTTAAATGCTGTTCCACAGGAAAGTCCATGGCTGTTATTCAAGATATCATCGCTGCGCTCTGGCAacac                             <  1:8538113/93‑1 (MQ=255)
                                                             ggAGTACAATAGGCAGATAAAGGCTTAAATGCTGTTCCACAGGAAAGTCCATGGCTGTTATTCAAGATATCATCGCTGCGCTCTGGCAACACg                            >  2:5363143/1‑93 (MQ=255)
                                                              gAGTACTATAGGCAGCTAAGGGCTTAAATGCTGTTCCACAGGAAAGTACATGGCTGTTATTCAAGATATCATCGCTGCGCTCTGGCAACACGa                           <  2:7829770/93‑1 (MQ=255)
                                                               aGTACAATAGGCAGATAAAGGCTTAAATGCTGTTCCACATGAAAGTCCATGGCTGTTATTCAAGATATCATCGCTGCGCTCTGGCAACACGAc                          >  1:4934529/1‑93 (MQ=255)
                                                               aGTACAATAGGCAGATAAAGGCTTAAATGCTGTTCCACAGGAAAGTCCATGGCTGTTATTCAAGATATCATCGCTGCGCTCTGGCAACACGAc                          <  2:1090123/93‑1 (MQ=255)
                                                                   cAATAGGCAGATAAAGGCTTAAATGCTGTTCCACAGGAAAGTCCATGGCTGTTATTCAAGATATCATCGCTGCGCTCTGGCAACACGACTTTg                      <  2:574319/93‑1 (MQ=255)
                                                                     aTAGGCAGATAAAGGCTTAAATGCTGTTCCACAGGAAAGTCCATGGCTGTTATTCAAGATATCATCGCTGCGCTCTGGCAACACGACTTTGcc                    >  1:1913448/1‑93 (MQ=255)
                                                                       aGGCAGATAAAGGCTTAAATGCTGTTCCACAGGAAAGTCCATGGCTGTTATTCAAGATATCATCGCTGCGCTCTGGCAACACGACTTTGCcgc                  <  1:6744931/93‑1 (MQ=255)
                                                                        ggCAGATAAAGGCTTAAATGCTGTTCCACAGGAAAGTCCATGGCTGTTATTCAAGATATCATCGCTGCGCTCTGGCAACACGACTTTGCcgcg                 >  1:5500712/1‑93 (MQ=255)
                                                                         gCAGATAAAGGCTTAAATGCTGTTCCACAGGAAAGTCCATGGCTGTTATTCAAGATATCATCGCTGCGCTCTGGCAACACGACTTTGCCgcgc                >  1:4978040/1‑93 (MQ=255)
                                                                          cAGATAAAGGCTTAAATGCTGTTCCACAGGAAAGTCCATGGCTGTTATTCAAGATATCATCGCTGCGCTCTGGCAACACGACTTTGCCGCGCt               <  2:2872548/93‑1 (MQ=255)
                                                                          cAGATAAAGGCTTAAATGCTGTTCCACAGGAAAGTCCATGGCTGTTATTCAAGATATCATCGCTGCGCTCTGGCAACACGACTTTGCCGCGCt               <  1:7082870/93‑1 (MQ=255)
                                                                             aTAAAGGCTTAAATGCTGTTCCACAGGAAAGTCCATGGCTGTTATTCAAGATATCATCGCTGCGCTCTGGCAACACGACTTTGCCGCGCTGGc            >  2:3683920/1‑93 (MQ=255)
                                                                             aTAAAGGCTTAAATGCTGTTCCACAGGAAAGTCCATGGCTGTTATTCAAGATATCATCGCTGCGCTCTGGCAACACGACTTTGCCGCGCTGGc            >  1:5827032/1‑93 (MQ=255)
                                                                              tAAAGGCTTAAATGCTGTTCCACAGGAAAGTCCATGGCTGTTATTCAAGATATCATCGCTGCGCTCTGGCAACACGACTTTGCCGCGCTGGCg           >  2:6134475/1‑93 (MQ=255)
                                                                                 aGGCTTAAATGCTGTTCCACAGGAAAGTCCATGGCTGTTATTCAAGATATCATCGCTGCGCTCTGGCAACACGACTTTGCCGCGCTGGCGGAt        >  2:1686334/1‑93 (MQ=255)
                                                                                  ggCTTAAATGCTGTTCCACAGGAAAGTCCATGGCTGTTATTCAAGATATCATCGCTGCGCTCTGGCAACACGACTTTGCCGCGCTGGCGGATc       <  1:5355650/93‑1 (MQ=255)
                                                                                  ggCTTAAATGCTGTTCCACAGGAAAGTCCATGGCTGTTATTCAAGATATCATCGCTGCGCTCTGGCAACACGACTTTGCCGCGCTGGCGGATc       >  2:7607978/1‑93 (MQ=255)
                                                                                   gCTTAAATGCTGTTCCACAGGAAAGTCCATGGCTGTTATTCAAGATATCATCGCTGCGCTCTGGCAACACGACTTTGCCGCGCTGGCGGATcc      >  1:1819607/1‑93 (MQ=255)
                                                                                   gCTTAAATGCTGTTCCACAGGAAAGTCCATGGCTGTTATTCAAGATATCATCGCTGCGCTCTGGCAACACGACTTTGCCGCGCTGGCGGATcc      <  1:7185792/93‑1 (MQ=255)
                                                                                     ttAAATGCTGTTCCACAGGAAAGTCCATGGCTGTTATTCAAGATATCATCGCTGCGCTCTGGCAACACGACTTTGCCGCGCTGGCGGATCCTc    >  1:8812048/1‑93 (MQ=255)
                                                                                     ttAAATGCTGTTCCACAGGAAAGTCCATGGCTGTTATTCAAGATATCATCGCTGCGCTCTGGCAACACGACTTTGCCGCGCTGGCGGATCCTc    >  2:1848621/1‑93 (MQ=255)
                                                                                      tAAATGCTGTTCCACAGGAAAGTCCATGGCTGTTATTCAAGATATCATCGCTGCGCTCTGGCAACACGACTTTGCCGCGCTGGCGGATCCTCa   <  1:2799010/93‑1 (MQ=255)
                                                                                      tAAATGCTGTTCCACAGGAAAGTCCATGGCTGTTATTCAAGATATCATCGCTGCGCTCTGGCAACACGACTTTGCCGCGCTGGCGGATCCTCa   >  1:7598843/1‑93 (MQ=255)
                                                                                      tAAATGCTGTTCCACAGGAAAGTCCATGGCTGTTATTCAAGATATCATCGCTGCGCTCTGGCAACACGACTTTCCCGCGCTGGCGGATCCTCa   <  2:4272441/93‑1 (MQ=255)
                                                                                          |                                                                                         
CATTTTATTGTGAAAAGTCTTAAATTGTTGCGTCCGGGATCAAGGCGTCCCGGACGATTCAGGAGTACAATAGGCAGATAAAGGCTTAAACGCTGTTCCACAGGAAAGTCCATGGCTGTTATTCAAGATATCATCGCTGCGCTCTGGCAACACGACTTTGCCGCGCTGGCGGATCCTCAT  >  NC_000913/3151467‑3151646

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