Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,150,581 T→C Y111H (TAT→CAT)  metC → cystathionine beta‑lyase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,150,5810TC99.7% 1300.2 / ‑4.8 376Y111H (TAT→CAT) metCcystathionine beta‑lyase
Reads supporting (aligned to +/- strand):  ref base T (0/1);  new base C (188/187);  total (188/188)
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.43e-01

Only 100 of 383 total aligned reads displayed.

GTGCTATTTCCCTGCGGGGCGGCAGCGGTTGCTAATTCCATTCTTGCTTTTATCGAACAGGGCGATCATGTGTTGATGACCAACACCGCCTATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGCAAACTGGGCGTAACGACATCATGGTTTGATCCGCTGATTGGTGCCGATATC  >  NC_000913/3150491‑3150670
                                                                                          |                                                                                         
gTGCTATTTCCCTGCGGGGCGGCAGCGGTTGCTAATTCCATTCTTGCTTTTATCGAACAGGGCGATCATGTGTTGATGACCAACACCGCCCAt                                                                                         >  1:3762817/1‑93 (MQ=255)
 tGCTATTTCCCTGCGGGGCGGCAGCGGTTGCTAATTCCATTCTTGCTTTTATCGAACAGGGCGATCATGTGTTGATGACCAACACCGCCCATg                                                                                        <  1:4856317/93‑1 (MQ=255)
    tATTTCCCTGCGGGGCGGCAGCGGTTGCTAATTCCATTCTTGCTTTTATCGAACAGGGCGATCATGTGTTGATGACCAACACCGCCCATGAAc                                                                                     <  1:5963821/93‑1 (MQ=255)
    tATTTCCCTGCGGGGCGGCAGCGGTTGCTAATTCCATTCTTGCTTTTATCGAACAGGGCGATCATGTGTTGATGACCAACACCGCCCATGAAc                                                                                     >  1:717379/1‑93 (MQ=255)
    tATTTCCCTGCGGGGCGGCAGCGGTTGCTAATTCCATTCTTGCTTTTATCGAACAGGGCGATCATGTGTTGATGACCAACACCGCCCATGAAc                                                                                     <  2:1577604/93‑1 (MQ=255)
     aTTTCCCTGCGGGGCGGCAGCGGTTGCTAATTCCATTCTTGCTTTTATCGAACAGGGCGATCATGTGTTGATGACCAACACCGCCCATGAAcc                                                                                    >  1:8208736/1‑93 (MQ=255)
     aTTTCCCTGCGGGGCGGCAGCGGTTGCTAATTCCATTCTTGCTTTTATCGAACAGGGCGATCATGTGTTGATGACCAACACCGCCCATGAAcc                                                                                    <  2:3168826/93‑1 (MQ=255)
       ttCCCTGCGGGGCGGCAGCGGTTGCTAATTCCATTCTTGCTTTTATCGAACAGGGCGATCATGTGTTGATGACCAACACCGCCCATGAACCGa                                                                                  <  2:3691190/93‑1 (MQ=255)
        tCCCTGCGGGGCGGCAGCGGTTGCTAATTCCATTCTTGCTTTTATCGAACAGGGCGATCATGTGTTGATGACCAACACCGCCCATGAACCGAg                                                                                 >  2:26254/1‑93 (MQ=255)
         cccTGCGGGGCGGCAGCGGTTGCTAATTCCATTCTTGCTTTTATCTAACAGGGCGATCATGTGTTGATGACCAACACCGCCCATGAACCGAGt                                                                                >  1:7115243/1‑93 (MQ=255)
         cccTGCGGGGCGGCAGCGGTTGCTAATTCCATTCTTGCTTTTATCGAACAGGGCGATCATGTGTTGATGACCAACACCGCCCATGAACCGAGt                                                                                >  1:43600/1‑93 (MQ=255)
          ccTGCGGGGCGGCAGCGGTTGCTAATTCCATTCTTGCTTTTATCGAACAGGGCGATCATGTGTTGATGACCAACACCGCCCATGAACCGAGTc                                                                               >  1:8075354/1‑93 (MQ=255)
              cGGGGCGGCAGCGGTTGCTCATTCCATTCTTGCTTTTATCGAACAGGGCGATCATGTGTTGATGACCAACACCGCCCATGAACCGAGTCAGGa                                                                           <  1:6219550/93‑1 (MQ=255)
              cGGGGCGGCAGCGGTTGCTAATTCCATTCTTGCTTTTATCGAACAGGGCGATCATGTGTTGATGACCAACACCGCCCATGAACCGAGTCAGGa                                                                           <  2:6982381/93‑1 (MQ=255)
              cGGGGCGGCAGAGGTTGCTAATTCCATTCTTGCTTTTATCGAACAGGGCGATCATGTGTTGATGACCAACACCGCCCATGAACCGAGTCAGGa                                                                           <  1:4467929/93‑1 (MQ=255)
                gggCGGCAGCGGTTGCTAATTCCATTCTTGCTTTTATCGAACAGGGCGATCATGTGTTGATGACCAACACCGCCCATGAACCGAGTCAGGAtt                                                                         <  2:2105621/93‑1 (MQ=255)
                 ggcggcAGCGGTTGCTAATTCCATTCTTGCTTTTATCGAACAGGGCGATCATGTGTTGATGACCAACACCGCCCATGAACCGAGTCAGGAttt                                                                        <  2:5908942/93‑1 (MQ=255)
                  gcggcAGCGGTTGCTAATTCCATTCTTGCTTTTATCGAACAGGGCGATCATGTGTTGATGACCAACACCGCCCATGAACCGAGTCAGGATTTc                                                                       <  1:1157915/93‑1 (MQ=255)
                     gcAGCGGTTGCTAATTCCATTCTTGCTTTTATCGAACAGGGCGATCATGTGTTGATGACCAACACCGCCCATGAACCGAGTCAGGATTTCTGt                                                                    >  1:3661750/1‑93 (MQ=255)
                       cgCGGTTGCTAATTCCTATCTTGCTTTTATCGAACAGGGCGCTCATGAGTTGATGCCCAACACCGCCCATGAACCGAGTCAAGATTTCTGTAg                                                                  <  2:7787458/92‑1 (MQ=255)
                        gCGGTTGCTAATTCCATTCTTGCTTTTATCGAACAGGGCGATCATGTGTTGATGACCAACACCGCCCATGAACCGAGTCAGGATTTCTGTAGc                                                                 <  1:3131639/93‑1 (MQ=255)
                         cGGTTGCTAATTCCATTCTTGCTTTTATCGAACAGGGCGATCATGTGTTGATGACCAACACCGCCCATGAACCGAGTCAGGATTTCTGTAGCa                                                                >  2:1520649/1‑93 (MQ=255)
                         cGGTTGCTAATTCCATTCTTGCTTTTATCGAACAGGGCGATCATGTGTTGATGACCAACACCGCCCATGAACCGAGTCAGGATTTCTGTAGCa                                                                >  2:2181227/1‑93 (MQ=255)
                          ggTTGCTAATTCCATTCTTGCTTTTATCGAACAGGGCGATCATGTGTTGATGACCAACACCGCCCATGAACCGAGTCAGGATTTCTGTAGCaa                                                               >  2:6231374/1‑93 (MQ=255)
                            ttGCTAATTCCATTCTTGCTTTTATCGAACAGGGCGATCATGTGTTGATGACCAACACCGCCCATGAACCGAGTCAGGATTTCTGTAGCaaaa                                                             <  2:2881166/93‑1 (MQ=255)
                            gtGCTAATTCCATTCTTGCTTTTATCGAACAGGGCGATCATGTGTTGATGACCAACACCGCCCATGAACCGAGTCAGGATTTCTGTAGCaaaa                                                             <  2:3408502/92‑1 (MQ=255)
                             tGCTAATTCCATTCTTGCTTTTATCGAACAGGGCGATCATGTGTTGATGACCAACACCGCCCATGAACCGAGTCAGGATTTCTGTAGCAAAAt                                                            >  1:4972055/1‑93 (MQ=255)
                             tGCTAATTCCATTCTTGCTTTTATCGAACAGGGCGATCATGTGTTGATGACCAACACCGCCCATGAACCGAGTCAGGATTTCTGTAGCAAAAt                                                            >  1:7849161/1‑93 (MQ=255)
                             tGCTAATTCCATTCTTGCTTTTATCGAACAGGGCGATCATGTGTTGATGACCAACACCGCCCATGAACCGAGTCAGGATTTCTGTAGCAAAAt                                                            >  1:1486147/1‑93 (MQ=255)
                              gCTAATTCCATTCTTGCTTTTATCGAACAGGGCGATCATGTGTTGATGACCAACACCGCCCATGAACCGAGTCAGGATTTCTGTAGCAAAATc                                                           >  1:1000793/1‑93 (MQ=255)
                               cTAATTCCATTCTTGCTTTTATCGAACAGGGCGATCATGTGTTGATGACCAACACCGCCCATGAACCGAGTCAGGATTTCTGTAGCAAAATcc                                                          <  1:41888/93‑1 (MQ=255)
                               cTAATTCCATTCTTGCTTTTATCGAACAGGGCGATCATGTGTTGATGACCAACACCGCCCATGAACCGAGTCAGGATTTCTGTAGCAAAATcc                                                          <  2:4989359/93‑1 (MQ=255)
                                tAATTCCATTCTTGCTTTTATCGAACAGGGCGATCATGTGTTGATGACCAACACCGCCCATGCACCGAGTCAGGATTTCTGTAGCAAAATCCt                                                         >  2:540027/1‑93 (MQ=255)
                                tAATTCCATTCTTGCTTTTATCGAACAGGGCGATCATGTGTTGATGACCAACACCGCCCATGAACCGAGTCAGGATTTCTGTAGCAACATCCt                                                         >  2:1995962/1‑93 (MQ=255)
                                 aaTTCCATTCTTGCTTTTATCGAACAGGGCGATCATGTGTTGATGACCAACACCGCCCATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTc                                                        <  2:3535255/93‑1 (MQ=255)
                                  aTTCCATTCTTGCTTTTATCGAACAGGGCGATCATGTGTTGATGACCAACACCGCCCATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCa                                                       <  1:3004942/93‑1 (MQ=255)
                                   ttCCATTCTTGCTTTTATCGAACAGGGCGATCATGTGTTGATGACCAACACCGCCCATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAg                                                      <  1:3566868/93‑1 (MQ=255)
                                   ttCCATTCTTGCTTTTATCGAACAGGGCGATCATGTGTTGATGACCAACACCGCCCATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAg                                                      >  2:2691138/1‑93 (MQ=255)
                                    tCCATTCTTGCTTTTATCGAACAGGGCGATCATGTGTTGATGACCAACACCGCCCATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGc                                                     >  2:8021747/1‑93 (MQ=255)
                                     ccATTCTTGCTTTTATCGAACAGGGCGATCATGTGTTGATGACCAACACCGCCCATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGCa                                                    >  1:2753423/1‑93 (MQ=255)
                                      cATTCTTGCTTTTATCGAACAGGGCGATCATGTGTTGATGACCAACACCGCCCATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGCaa                                                   >  1:5288415/1‑93 (MQ=255)
                                        ttCTTGCTTTTATCGAACAGGGCGATCATGTGTTGATGACCAACACCGCCCATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGCAAAc                                                 >  1:7453709/1‑93 (MQ=255)
                                          cTTGCTTTTATCGAACAGGGCGATCATGTGTTGATGACCAACACCGCCCATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGCAAACTg                                               >  2:5199670/1‑93 (MQ=255)
                                          cTTGCTTTTATCGAACAGGGCGATCATGTGTTGATGACCAACACCGCCCATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGCAAACTg                                               <  1:2711978/93‑1 (MQ=255)
                                           ttGCTTTTATCGAACAGGGCGATCATGTGTTGATGACCAACACCGCCCATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGCAAACTgg                                              <  1:6496429/93‑1 (MQ=255)
                                           ttGCTTTTATCGAACAGGGCGATCATGTGTTGATGACCAACACCGCCCATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGCAAACTgg                                              >  2:7700001/1‑93 (MQ=255)
                                             gCTTTTATCGAACAGGGCGATCATGTGTTGATGACCAACACCGCCCATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGCAAACTGGGc                                            <  1:4995128/93‑1 (MQ=255)
                                             gCTTTTATCGAACAGGGCGATCATGTGTTGATGACCAACACCGCCCATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGCAAACTGGGc                                            >  1:6624206/1‑93 (MQ=255)
                                               ttttATCGAACAGGGCGATCATGTGTTGATGACCAACACCGCCCATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGCAAACTGGGCGt                                          >  2:8224350/1‑93 (MQ=255)
                                               ttttATCGAACAGGGCGATCATGTGTTGATGACCAACACCGCCCATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGCAAACTGGGCGt                                          <  1:1235675/93‑1 (MQ=255)
                                               ttttATCGAACAGGGCGATCATGTGTTGATGACCAACACCGCCCATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGCAAACTGGGCGt                                          >  1:2374204/1‑93 (MQ=255)
                                                 ttATCGAACAGGGCGATCATGTGTTGATGACCAACACCGCCCATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGCAAACTGGGCGTaa                                        <  1:5400182/93‑1 (MQ=255)
                                                 ttATCGAACAGGGCGATCATGTGTTGATGACCAACACCGCCCATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGCAAACTGGGCGTaa                                        >  2:6070031/1‑93 (MQ=255)
                                                  tATCGAACAGGGCGATCATGTGTTGATGACCAACACCGCCCATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGCAAACTGGGCGTAAc                                       <  2:8735664/93‑1 (MQ=255)
                                                  tATCGAACAGGGCGATCATGTGTTGATGACCAACACCGCCCATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGCAAACTGGGCGTAAc                                       >  1:934343/1‑93 (MQ=255)
                                                  tATCGAACAGGGCGATCATGTGTTGATGACCAACACCGCCCATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGCAAACTGGGCGTAAc                                       <  1:1670163/93‑1 (MQ=255)
                                                  tATCGAACAGGGCGATCATGTGTTGATGACCAACACCGCCCATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGCAAACTGGGCGTAAc                                       >  2:4909302/1‑93 (MQ=255)
                                                   aTCGAACAGGGCGATCATGTGTTGATGACCAACACCGCCCATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGCAAACTGGGCGTAACg                                      >  2:5882186/1‑93 (MQ=255)
                                                      gAACAGGGCGATCATGTGTTGATGACCAACACCGCCCATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGCAAACTGGGCGTAACGAca                                   >  2:1094315/1‑93 (MQ=255)
                                                      gAACAGGGCGATCATGTGTTGATGACCAACACCGCCCATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGCAAACTGGGCGTAACGAca                                   >  1:7923218/1‑93 (MQ=255)
                                                      gAACAGGGCGATCATGTGTTGATGACCAACACCGCCCATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGCAAACTGGGCGTAACGAca                                   >  1:8438325/1‑93 (MQ=255)
                                                       aaCAGGGCGATCATGTGTTGATGACCAACACCGCCCATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGCAAACTGGGCGTAACGAcat                                  <  2:4823057/93‑1 (MQ=255)
                                                        aCAGGGCGATCATGTGTTGATGACCAACACCGCCCATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGCAAACTGGGCGTAACGAcatc                                 <  2:2531052/93‑1 (MQ=255)
                                                          aGGGCGATCATGTGTTGATGACCAACACCGCCCATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGCAAACTGGGCGTAACGAcatcat                               <  1:3302433/93‑1 (MQ=255)
                                                          aGGGCGATCATGTGTTGATGACCAACACCGCCCATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGCAAACTGGGCGTAACGAcatcat                               <  1:1806165/93‑1 (MQ=255)
                                                           gggCGATCATGTGTTGATGACCAACACCGCCCATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGCAAACTGGGCGTAACGACATCATg                              >  2:1238622/1‑93 (MQ=255)
                                                            ggCGATCATGTGTTGATGACCAACACCGCCCATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGCAAACTGGGCGTAACGACATCATgg                             >  2:4133474/1‑93 (MQ=255)
                                                             gCGATCATGTGTTGATGACCAACACCGCCCATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGCAAACTGGGCGTAACGACATCATGGt                            <  2:4559646/93‑1 (MQ=255)
                                                              cGATCATGTGTTGATGACCAACACCGCCCATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGCAAACTGGGCGTAACGACATCATGGtt                           <  1:6856457/93‑1 (MQ=255)
                                                               gATCATGTGTTGATGACCAACACCGCCCATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGCAAACTGGGCGTAACGACATCATGGttt                          <  1:3986666/93‑1 (MQ=255)
                                                                aTCATGTGTTGATGACCAACACCGCCCATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGCAAACTGGGCGTAACGACATCATGGTTTg                         <  1:2006244/93‑1 (MQ=255)
                                                                 tCATGTGTTGATGACCAACACCGCCCATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGCAAACTGGGCGTAACGACATCATGGTTTGa                        <  2:7619806/93‑1 (MQ=255)
                                                                 tCATGTGTTGATGACCAACACCGCCCATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGCAAACTGGGCGTAACGACATCATGGTTTGa                        >  1:5512399/1‑93 (MQ=255)
                                                                 tCATGTGTTGATGACCAACACCGCCCATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGCAAACTGGGCGTAACGACATCATGGTTTGa                        <  2:6484606/93‑1 (MQ=255)
                                                                 tCATGTGTTGATGACCAACACCGCCCATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGCAAACTGGGCGTAACGACATCATGGTTTGa                        <  2:3755611/93‑1 (MQ=255)
                                                                  cATGTGTTGATGACCAACACCGCCCATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGCAAACTGGGCGTAACGACATCATGGTTTGAt                       >  2:3080580/1‑93 (MQ=255)
                                                                  cATGTGTTGATGACCAACACCGCCCATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGCAAACTGGGCGTAACGACATCATGGTTTGAt                       >  2:673840/1‑93 (MQ=255)
                                                                  cATGTGTTGATGACCAACACCGCCCATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGCAAACTGGGCGTAACGACATCATGGTTTGAt                       <  1:7562673/93‑1 (MQ=255)
                                                                  cATGTGTTAATGACCAACACCGCCCATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGCAAACTGGGCGTAACGACATCATGGTTTGAt                       <  2:430793/93‑1 (MQ=255)
                                                                   aTGTGTTGATGACCAACACCGCCCATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGCAAACTGGGCGTAACGACATCATGGTTTGATc                      >  2:5699695/1‑93 (MQ=255)
                                                                   aTGTGTTGATGACCAACACCGCCCATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGCAAACTGGGCGTAACGACATCATGGTTTGATc                      <  1:2267480/93‑1 (MQ=255)
                                                                    tgtgTTGATGACCAACACCGCCCATTAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGCAAACTGGGCGTAACGACATCATGGTTTGATcc                     <  2:1922318/93‑1 (MQ=255)
                                                                    tgtgTTGATGACCAACACCGCCCATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGCAAACTGGGCGTAACGACATCATGGTTTGATcc                     >  1:509237/1‑93 (MQ=255)
                                                                      tgtTGATGACCAACACCGCCCATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGCAAACTGGGCGTAACGACATCATGGTTTGATCCGc                   >  2:657718/1‑93 (MQ=255)
                                                                        ttGATGACCAACACCGCCCATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGCAAACTGGGCGTAACGACATCATGGTTTGATCCGCTg                 >  2:6876891/1‑93 (MQ=255)
                                                                             gaCCAACACCGCCCATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGCAAACTGGGCGTAACGACATCATGGTTTGATCCGCTGATTgg            >  2:8617935/1‑93 (MQ=255)
                                                                             gaCCAACACCGCCCATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGCAAACTGGGCGTAACGACATCATGGTTTGATCCGCTGATTgg            <  1:466005/93‑1 (MQ=255)
                                                                               ccAACACCGCCCATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGCAAACTGGGCGTAACGACATCATGGTTTGATCCGCTGATTGGTg          <  2:5463647/93‑1 (MQ=255)
                                                                               ccAACACCGCCCATGAACCGAGACAGGATTTCTGTAGCAAAATCCTCAGCAAACTGGGCGTAACGACATCATGGTTTGATCCGCTGATTGGTg          >  2:5528432/1‑93 (MQ=255)
                                                                                  acacCGCCCATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGCAAACTGGGCGTAACGACATCATGGTTTGATCCGCTGATTGGTGCCg       >  2:1413648/1‑93 (MQ=255)
                                                                                   caccacCCATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGCAAACTGGGCGTAACGACATCATGGTTTGATCCGCTGATTGGTGCCGa      <  2:1736947/93‑1 (MQ=255)
                                                                                   cacCGCCCATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGCAAACTGGGCGTAACGACATCATGGTTTGATCCGCTGATTGGTGCCGa      >  2:2971103/1‑93 (MQ=255)
                                                                                   cacCGCCCATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGCAAACTGGGCGTAACGACATCATGGTTTGATCCGCTGATTGGTGCCGa      >  2:222524/1‑93 (MQ=255)
                                                                                   cacCGCCCAGGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGCAAACTGGGCGTAACGACATCATGGTTTGATCCGCTGATTGGTGCCGa      >  1:8681921/1‑93 (MQ=255)
                                                                                    acCGCCCATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGCAAACTGGGCGTAACGACATCATGGTTTGATCCGCTGATTGGTGCCGat     >  2:6905930/1‑93 (MQ=255)
                                                                                     ccGCCCATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGCAAACTGGGCGTAACGACATCATGGTTTGATCCGCTGATTGGTGCCGata    <  2:7370828/93‑1 (MQ=255)
                                                                                     ccGCCCATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGCAAACTGGGCGTAACGACATCATGGTTTGATCCGCTGATTGGTGCCGata    <  2:713625/93‑1 (MQ=255)
                                                                                     ccGCCCATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGCAAACTGGGCGTAACGACATCATGGTTTGATCCGCTGATTGGTGCCGata    <  1:8537830/93‑1 (MQ=255)
                                                                                      cGCCCATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGCAAACTGGGCGTAACGACATCATGGTTTGATCCGCTGATTGGTGCCGatat   <  2:127588/93‑1 (MQ=255)
                                                                                       gCCCATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGCAAACTGGGCGTAACGACATCATGGTTTGATCCGCTGATTGGTGCCGATATc  <  2:2422190/93‑1 (MQ=255)
                                                                                          |                                                                                         
GTGCTATTTCCCTGCGGGGCGGCAGCGGTTGCTAATTCCATTCTTGCTTTTATCGAACAGGGCGATCATGTGTTGATGACCAACACCGCCTATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGCAAACTGGGCGTAACGACATCATGGTTTGATCCGCTGATTGGTGCCGATATC  >  NC_000913/3150491‑3150670

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