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.6% 938.6 / NA 242intergenic (+119/‑21)metC/yghBcystathionine beta‑lyase/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (106/135);  minor base A (1/0);  total (107/135)
Fisher's exact test for biased strand distribution p-value = 4.42e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 245 total aligned reads displayed.

CATTTTATTGTGAAAAGTCTTAAATTGTTGCGTCCGGGATCAAGGCGTCCCGGACGATTCAGGAGTACAATAGGCAGATAAAGGCTTAAACGCTGTTCCACAGGAAAGTCCATGGCTGTTATTCAAGATATCATCGCTGCGCTCTGGCAACACGACTTTGCCGCGCTGGCGGATCCTCAT  >  NC_000913/3151467‑3151646
                                                                                          |                                                                                         
 aTTTTATTGTGAAAAGTCTTAAATTGTTGCGTCCGGGATCAAGGCGTCCCGGACGATTCAGGAGTACAATAGGCAGATAAAGGCTTAAATGCt                                                                                        >  1:5789742/1‑93 (MQ=255)
  ttttATTGTGAAAAGTCTTAAATTGTTGCGTCCGGGATCAAGGCGTCCCGGACGATTCAGGAGTACAATAGGCAGATAAAGGCTTAAATGCTg                                                                                       <  2:5326557/93‑1 (MQ=255)
   tttATTGTGAAAAGTCTTAAATTGTTGCGTCCGGGATCAAGGCGTCCCGGACGATTCAGGAGTACAATAGGCAGATAAAGGCTTAAATGCTGt                                                                                      >  1:4068682/1‑93 (MQ=255)
    ttATTGTGAAAAGTCTTAAATTGTTGCGTCCGGGATCAAGGCGTCCCGGACGATTCAGGAGTACAATAGGCAGATAAAGGCTTAAATGCTGtt                                                                                     >  1:5118328/1‑93 (MQ=255)
     tATTGTGAAAAGTCTTAAATTGTTGCGTCCGGGATCAAGGCGTCCCGGACGATTCAGGAGTACAATAGGCAGATAAAGGCTTAAATGCTGTTc                                                                                    >  2:3712410/1‑93 (MQ=255)
      aTTGTGAAAAGTCTTAAATTGTTGCGTCCGGGATCAAGGCGTCCCGGACGATTCAGGAGTACAATAGGCAGATAAAGGCTTAAATGCTGTTcc                                                                                   <  1:1083002/93‑1 (MQ=255)
       ttGTGAAAAGTCTTAAATTGTTGCGTCCGGGATCAAGGCGTCCCGGACGATTCAGGAGTACAATAGGCAGATAAAGGCTTAAATGCTGTTCca                                                                                  >  2:1868524/1‑93 (MQ=255)
       ttGTGAAAAGTCTTAAATTGTTGCGTCCGGGATCAAGGCGTCCCGGACGATTCAGGAGTACAATAGGCAGATAAAGGCTTAAATGCTGTTCca                                                                                  >  1:262984/1‑93 (MQ=255)
         gtgAAAAGTCTTAAATTGTTGCGTCCGGGATCAAGGCGTCCCGGACGATTCAGGAGTACAATAGGCAGATAAAGGCTTAAATGCTGTTCcaca                                                                                >  1:3894223/1‑93 (MQ=255)
          tgAAAAGTCTTAAATTGTTGCGTCCGGGATCAAGGCGTCCCGGACGATTCAGGAGTACAATAGGCAGATAAAGGCTTAAATGCTGTTCCACAg                                                                               <  1:3508749/93‑1 (MQ=255)
           gAAAAGTCTTAAATTGTTGCGTCCGGGATCAAGGCGTCCCGGACGATTCAGGAGTACAATAGGCAGATAAAGGCTTAAATGCTGTTCCACAgg                                                                              >  2:5510438/1‑93 (MQ=255)
              aaGTCTTAAATTGTTGCGTCCGGGATCAAGGCGTCCCGGACGATTCAGGAGTACAATAGGCAGATAAAGGCTTAAATGCTGTTCCACAGGcaa                                                                           >  1:3050746/1‑93 (MQ=255)
              aaGTCTTAAATTGTTGCGTCCGGGATCAAGGCGTCCCGGACGATTCAGGAGTACAATAGGCAGATAAAGGCTTAAATGCTGTTCCACAGGaaa                                                                           <  1:3188791/93‑1 (MQ=255)
              aaGTCTTAAATTGTTGCGTCCGGGATCAAGGAGTCCCGGACGATTCAGGAGTACAATAGGCAGATAAAGGCTTAAATGCTGTTCCACAGGaaa                                                                           >  2:3887994/1‑93 (MQ=255)
               aGTCTTAAATTGTTGCGTCCGGGATCAAGGCGTCCCGGACGATTCAGGAGTACAATAGGCAGATAAAGGCTTAAATGCTGTTCCACAGGAAAg                                                                          >  2:2649038/1‑93 (MQ=255)
                   ttAAATTGTTGCGTCCGGGATCAAGGCGTCCCGGACGATTCAGGAGTACAATAGGCAGATAAAGGCTTAAATGCTGTTCCACAGGAAAGTCCa                                                                      >  1:4170258/1‑93 (MQ=255)
                   ttAAATTGTTGCGTCCGGGATCAAGGCGTCCCGGACGATTCAGGAGTACAATAGGCAGATAAAGGCTTAAATGCTGTTCCACAGGAAAGTCCa                                                                      >  1:3667971/1‑93 (MQ=255)
                    tAAATTGTTGCGTCCGGGATCAAGGCGTCCCGGACGATTCAGGAGTACAATAGGCAGATAAAGGCTTAAATGCTGTTCCACAGGAAAGTCCAt                                                                     <  1:4850891/93‑1 (MQ=255)
                    tAAATTGTTGCGTCCGGGATCAAGGCGTCCCGGACGATTCAGGAGTACAATAGGCAGATAAAGGCTTAAATGCTGTTCCACAGGAAAGTCCAt                                                                     <  1:1324156/93‑1 (MQ=255)
                    tAAATTGTTGCGTCCGGGATCAAGGCGTCCCGGACGATTCAGGAGTACAATAGGCAGATAAAGGCTTAAATGCTGTTCCACAGGAAAGTCCAt                                                                     >  1:1286105/1‑93 (MQ=255)
                     aaaTTGTTGCGTCCGGGATCAAGGCGTCCCGGACGATTCAGGAGTACAATAGGCAGATAAAGGCTTAAATGCTGTTCCACAGGAAAGTCCATg                                                                    <  2:6308495/93‑1 (MQ=255)
                       aTTGTTGCGTCCGGGATCAAGGCGTCCCGGACGATTCAGGAGTACAATAGGCAGATAAAGGCTTAAATGCTGTTCCACAGGAAAGTCCATGGc                                                                  >  2:3555053/1‑93 (MQ=255)
                       aTTGTTGCGTCCGGGATCAAGGCGTCCCGGACGATTCAGGAGTACAATAGGCAGATAAAGGCTTAAATGCTGTTCCACAGGAAAGTCCATGGc                                                                  >  2:4359971/1‑93 (MQ=255)
                         tgttgCGTCCGGGATCAAGGCGTCCCGGACGATTCAGGAGTACAATAGGCAGATAAAGGCTTAAATGCTGTTCCACAGGAAAGTCCATGGCTg                                                                >  1:2943710/1‑93 (MQ=255)
                          gttgCGTCCGGGATCAAGGCGTCCCGGACGATTCAGGAGTACAATAGGCAGATAAAGGCTTAAATGCTGTTCCACAGGAAAGTCCATGGCTGt                                                               <  2:1038770/93‑1 (MQ=255)
                           ttgCGTCCGGGATCAAGGCGTCCCGGACGATTCAGGAGTACAATAGGCAGATAAAGGCTTAAATGCTGTTCCACAGGAAAGTCCATGGCTGtt                                                              <  1:2071766/93‑1 (MQ=255)
                            tgCGTCCGGGATCAAGGCGTCCCGGACGATTCAGGCGTACAATAGGCAGATACAGGCTTCAATGCTGTTCCACAGGAAAGTCCATGGCTGTTa                                                             >  1:2371147/1‑93 (MQ=255)
                             gCGTCCGGGATCAAGGCGTCCCGGACGATTCAGGAGTACAATAGGCAGATAAAGGCTTAAATGCTGTTCCACAGGAAAGTCCATGGCTGTTAt                                                            <  1:4733789/93‑1 (MQ=255)
                                tCCGGGATCAAGGCGTCCCGGACGATTCAGGAGTACAATAGGCAGATAAAGGCTTAAATGCTGTTCCACAGGAAAGTCCATGGCTGTTATTCa                                                         >  1:6028281/1‑93 (MQ=255)
                                 ccGGGATCAAGGCGTCCCGGACGATTCAGGAGTACAATAGGCAGATAAAGGCTTAAATGCTGTTCCACAGGAAAGTCCATGGCTGTTATTCaa                                                        <  1:864523/93‑1 (MQ=255)
                                   gggATCAAGGCGTCCCGGACGATTCAGGAGTACAATAGGCAGATAAAGGCTTAAATGCTGTTCCACAGGAAAGTCCATGGCTGTTATTCAAGa                                                      >  1:2695126/1‑93 (MQ=255)
                                    ggATCAAGGCGTCCCGGACGATTCAGGAGTACAATAGGCAGATAAAGGCTTAAATGCTGTTCCACAGGAAAGTCCATGGCTGTTATTCAAGat                                                     >  1:1146849/1‑93 (MQ=255)
                                    ggATCAAGGCGTCCCGGACGATTCAGGAGTACAATAGGCAGATAAAGGCTTAAATGCTGTTCCACAGGAAAGTCCATGGCTGTTATTCAAGat                                                     >  2:4996353/1‑93 (MQ=255)
                                    ggATCAAGGCGTCCCGGACGATTCAGGAGTACAATAGGCAGATAAAGGCTTAAATGCTGTTCCACAGGAAAGTCCATGGCTGTTATTCAAGat                                                     >  1:348016/1‑93 (MQ=255)
                                     gATCAAGGCGTCCCGGACGATTCAGGAGTACAATAGGCAGATAAAGGCTTAAATGCTGTTCCACAGGAAAGTCCATGGCTGTTATTCAAGata                                                    <  1:5765112/93‑1 (MQ=255)
                                       tCAGGGCGGCCCGGACGAGTAAAGAAAAAAAGAGGCAAATAAAGGCTTAAATGCTGTTCCACAGGAAAGTCCATGGCTGTTATTCAAGATatc                                                  <  2:6423464/93‑1 (MQ=255)
                                       tCAAGGCGTCCCGGACGATTCAGGAGTACAATAGGCAGATAAAGGCTTAAATGCTGTTCCACAGGAAAGTCCATGGCTGTTATTCAAGATatc                                                  >  2:4221592/1‑93 (MQ=255)
                                        cAAGGCGTCCCGGACGATTCAGGAGTACAATAGGCAGATAAAGGCTTAAATGCTGTTCCACAGGAAAGTCCATGGCTGTTATTCAAGATatca                                                 <  2:1774864/93‑1 (MQ=255)
                                        cAAGGCGTCCCGGACGATTCAGGAGTACAATAGGCAGATAAAGGCTTAAATGCTGTTCCACAGGAAAGTCCATGGCTGTTATTCAAGATatca                                                 >  2:6098167/1‑93 (MQ=255)
                                          aGGCGTCCCGGACGATTCAGGAGTACAATAGGCAGATAAAGGCTTAAATGCTGTTCCACAGGAAAGTCCATGGCTGTTATTCAAGATatcatc                                               <  1:4937540/93‑1 (MQ=255)
                                          aGGCGTCCCGGACGATTCAGGAGTACAATAGGCAGATAAAGGCTTAAATGCTGTTCCACAGGAAAGTCCATGGCTGTTATTCAAAATatcatc                                               <  1:6310184/93‑1 (MQ=255)
                                           ggCGTCCCGGACGATTCAGGAGTACAATAGGCAGATAAAGGCTTAAATGCTGTTCCACAGGAAAGTCCATGGCTGTTATTCAAGATATCATCg                                              <  2:4609968/93‑1 (MQ=255)
                                            gCGTCACGGACGATTCAGGAGTACAATAGGCAGATAAAGGCTTAAATGCTGTTCCACAGGAAAGTCCATGGCTGTTATTCAAGATATCATCGc                                             <  1:5198889/93‑1 (MQ=255)
                                               tCCCGGACGATTCAGGAGTACAATAGGCAGATAAAGGCTTAAATGCTGTTCCACAGGAAAGTCCATGGCTGTTATTCAAGATATCATCGCTgc                                          >  1:540573/1‑93 (MQ=255)
                                               tCCCGGACGATTCAGGAGTACAATAGGCAGATAAAGGCTTAAATGCTGTTCCACAGGAAAGTCCATGGCTGTTATTCAAGATATCATCGCTgc                                          <  2:2035792/93‑1 (MQ=255)
                                                cccGGACGATTCAGGAGTACAATAGGCAGATAAAGGCTTAAATGCTGTTCCACAGGAAAGTCCATGGCTGTTATTCAAGATATCATCGCTgcg                                         <  1:5679520/93‑1 (MQ=255)
                                                 ccGGACGATTCAGGAGTACAATAGGCAGATAAAGGCTTAAATGCTGTTCCACAGGAAAGTCCATGGCTGTTATTCAAGATATCATCGCTgcgc                                        >  2:5658480/1‑93 (MQ=255)
                                                 ccGGACGATTCAGGAGTACAATAGGCAGATAAAGGCTTAAATGCTGTTCCACAGGAAAGTCCATGGCTGTTATTCAAGATATCATCGCTgcgc                                        <  1:103376/93‑1 (MQ=255)
                                                  cGGACGATTCAGGAGTACAATAGGCAGATAAAGGCTTAAATGCTGTTCCACAGGAAAGTCCATGGCTGTTATTCAAGATATCATCGCTGCGct                                       >  2:2119626/1‑93 (MQ=255)
                                                     gcGATTCAGGAGTACAATAGGCAGATAAAGGCTTAAATGCTGTTCCACAGGAAAGTCCATGGCTGTTATTCAAGATATCATCGCTGCGCTCTg                                    <  1:3935019/92‑1 (MQ=255)
                                                     aCGATTCAGGAGTACAATAGGCAGATAAAGGCTTAAATGCTGTTCCACAGGAAAGTCCATGGCTGTTATTCAAGATATCATCGCTGCGCTCTg                                    >  2:5078658/1‑93 (MQ=255)
                                                     aCGATTCAGGAGTACAATAGGCAGATAAAGGCTTAAATGCTGTTCCACAGGAAAGTCCATGGCTGTTATTCAAGATATCATCGCTGCGCTCTg                                    <  1:883311/93‑1 (MQ=255)
                                                     aCGATTCAGGAGTACAATAGGCAGATAAAGGCTTAAATGCTGTTCCACAGGAAAGTCCATGGCTGTTATTCAAGATATCATCGCTGCGCTCTg                                    >  1:6260012/1‑93 (MQ=255)
                                                       gATTCAGGAGTACAATAGGCAGATAAAGGCTTAAATGCTGTTCCACAGGAAAGTCCATGGCTGTTATTCAAGATATCATCGCTGCGCTCTGGc                                  <  2:2941949/93‑1 (MQ=255)
                                                       gATTCAGGAGTACAATAGGCAGATAAAGGCTTAAATGCTGTTCCACAGGAAAGTCCATGGCTGTTATTCAAGATATCATCGCTGCGCTCTGGc                                  >  2:6029722/1‑93 (MQ=255)
                                                       gATTCAGGAGTACAATAGGCAGATAAAGGCTTAAATGCTGTTCCACAGGAAAGTCCATGGCTGTTATTCAAGATATCATCGCTGCGCTCTGGc                                  >  1:4244328/1‑93 (MQ=255)
                                                          tCAGGAGTACAATAGGCAGATAAAGGCTTAAATGCTGTTCCACAGGAAAGTCCATGGCTGTTATTCAAGATATCATCGCTGCGCTCTGGCAac                               >  1:3299241/1‑93 (MQ=255)
                                                           cAGGAGTACAATAGGCAGATAAAGGCTTAAATGCTGTTCCACAGGAAAGTCCATGGCTGTTATTCAAGATATCATCGCTGCGCTCTGGCAaca                              <  2:6258947/93‑1 (MQ=255)
                                                            aGGAGTACAATAGGCAGATAAAGGCTTAAATGCTGTTCCACAGGAAAGTCCATGGCTGTTATTCAAGATATCATCGCTGCGCTCTGGCAacac                             <  1:6005099/93‑1 (MQ=255)
                                                             ggAGTACAATAGGCAGATAAAGGCTTAAATGCTGTTCCACAGGAAAGTCCATGGCTGTTATTCAAGATATCATCGCTGCGCTCTGGCAACACg                            >  1:1758474/1‑93 (MQ=255)
                                                                 tACAATAGGCAGATAAAGGCTTAAATGCTGTTCCACAGGAAAGTCCATGGCTGTTATTCAAGATATCATCGCTGCGCTCTGGCAACACGACtt                        <  2:5564681/93‑1 (MQ=255)
                                                                 tACAATAGGCAGATAAAGGCTTAAATGCTGTTCCACAGGAAAGTCCATGGCTGTTATTCAAGATATCATCGCTGCGCTCTGGCAACACGACtt                        <  2:2591709/93‑1 (MQ=255)
                                                                 tACAATAGGCAGATAAAGGCTTAAATGCTGTTCCACAGGAAAGTCCATGGCTGTTATTCAAGATATCATCGCTGCGCTCTGGCAACACGACtt                        >  1:4705024/1‑93 (MQ=255)
                                                                   cAATAGGCAGATAAAGGCTTAAATGCTGTTCCACAGGAAAGTCCATGGCTGTTATTCAAGATATCATCGCTGCGCTCTGGCAACACGACTTTg                      >  1:4619179/1‑93 (MQ=255)
                                                                   cAATAGGCAGATAAAGGCTTAAATGCTGTTCCACAGGAAAGTCCATGGCTGTTATTCAAGATATCATCGCTGCGCTCTGGCAACACGACTTTg                      <  2:1290427/93‑1 (MQ=255)
                                                                    aaTAGGCAGATAAAGGCTTAAATGCTGTTCCACAGGAAAGTCCATGGCTGTTATTCAAGATATCATCGCTGCGCTCTGGCAACACGACTTTGc                     <  2:5753775/93‑1 (MQ=255)
                                                                    aaTAGGCAGATAAAGGCTTAAATGCTGTTCCACAGGAAAGTCCATGGCTGTTATTCAAGATATCATCGCTGCGCTCTGGCAACACGACTTTGc                     >  2:4471152/1‑93 (MQ=255)
                                                                     aTAGGCAGATAAAGGCTTAAATGCTGTTCCACAGGAAAGTCCATGGCTGTTATTCAAGATATCATCGCTGCGCTCTGGCAACACGACTTTGcc                    <  2:1506119/93‑1 (MQ=255)
                                                                      tAGGCAGATAAAGGCTTAAATGCTGTTCCACAGGAAAGTCCATGGCTGTTATTCAAGATATCATCGCTGCGCTCTGGCAACACGACTTTGCcg                   <  2:3614850/93‑1 (MQ=255)
                                                                       aGGCAGATAAAGGCTTAAATGCTGTTCCACAGGAAAGTCCATGGCTGTTATTCAAGATATCATCGCTGCGCTCTGGCAACACGACTTTGCcgc                  <  2:886814/93‑1 (MQ=255)
                                                                       aGGCAGATAAAGGCTTAAATGCTGTTCCACAGGAAAGTCCATGGCTGTTATTCAAGATATCATCGCTGCGCTCTGGCAACACGACTTTGCcgc                  >  1:658997/1‑93 (MQ=255)
                                                                       aGGCAGATAAAGGCTTAAATGCTGTTCCACAGGAAAGTCCATGGCTGTTATTCAAGATATCATCGCTGCGCTCTGGCAACACGACTTTGCcgc                  <  2:3996307/93‑1 (MQ=255)
                                                                       aGGCAGATAAAGGCTTAAAAGCTGTTCCACAGGAAAGTCCATGGCTGTTATTCAAGATATCATCGCTGCGCTCTTGCAACACGACTTTGCcgc                  >  1:6320233/1‑93 (MQ=255)
                                                                        ggCAGATAAAGGCTTAAATGCTGTTCCACAGGAAAGTCCATGGCTGTTATTCAAGATATCATCGCTGCGCTCTGGCAACACGACTTTGCcgcg                 >  1:4342972/1‑93 (MQ=255)
                                                                        ggCAGATAAAGGCTTAAATGCTGTTCCACAGGAAAGTCCATGGCTGTTATTCAAGATATCATCGCTGCGCTCTGGCAACACGACTTTGCcgcg                 >  1:2268405/1‑93 (MQ=255)
                                                                        ggCAGATAAAGGCTTAAATGCTGTTCCACAGGAAAGTCCATGGCTGTTATTCAAGATATCATCGCTGCGCTCTGGCAACACGACTTTGCcgcg                 <  2:4458183/93‑1 (MQ=255)
                                                                         gCAGATAAAGGCTTAAATGCTGTTCCACAGGAAAGTCCATGGCTGTTATTCAAGATATCATCGCTGCGCTCTGGCAACACGACTTTGCCgcgc                >  2:2303938/1‑93 (MQ=255)
                                                                         gCAGATAAAGGCTTAAATGCTGTTCCACAGGAAAGTCCATGGCTGTTATTCAAGATATCATCGCTGCGCTCTGGCAACACGACTTTGCCgcgc                <  2:5410434/93‑1 (MQ=255)
                                                                         gCAGATAAAGGCTTAAATGCTGTTCCACAGGAAAGTCCATGGCTGTTATTCAAGATATCATCGCTGCGCTCTGGCAACACGACTTTGCCgcgc                <  2:1679753/93‑1 (MQ=255)
                                                                          cAGATAAAGGCTTAAATGCTGTTCCACAGGAAAGTCCATGGCTGTTATTCAAGATATCATCGCTGCGCTCTGGCAACACGACTTTGCCGCGCt               <  1:551698/93‑1 (MQ=255)
                                                                           aGATAAAGGCTTAAATGCTGTTCCACAGGAAAGTCCATGGCTGTTATTCAAGATATCATCGCTGCGCTCTGGCAACACGACTTTGCCGCGCTg              >  2:1021181/1‑93 (MQ=255)
                                                                            gATAAAGGCTTAAATGCTGTTCCACAGGAAAGTCCATGGCTGTTATTCAAGATATCATCGCTGCGCTCTGGCAACACGACTTTGCCGCGCTgg             >  1:1900410/1‑93 (MQ=255)
                                                                            gATAAAGGCTTAAATGCTGTTCCACAGGAAAGTCCATGGCTGTTATTCAAGATATCATCGCTGCGCTCTGGCAACACGACTTTGCCGCGCTgg             >  1:4506100/1‑93 (MQ=255)
                                                                            gATAAAGGCTTAAATGCTGTTCCACAGGAAAGTCCATGGCTGTTATTCAAGATATCATCGCTGCGCTCTGGCAACACGACTTTGCCGCGCTgg             <  2:4770731/93‑1 (MQ=255)
                                                                            gATAAAGGCTTAAATGCTGTTCCACAGGAAAGTCCATGGCTGTTATTCAAGATATCATCGCTGCGCTCTGGCAACACGACTTTGCCGCGCTgg             <  2:306061/93‑1 (MQ=255)
                                                                            gATAAAGGCTTAAATGCTGTTCCACAGGAAAGTCCATGGCTGTTATTCAAGATATCATCGCTGCGCTCTGGCAACACGACTTTGCCGCGCTgg             <  2:5291504/93‑1 (MQ=255)
                                                                             aTAAAGGCTTAAATGCTGTTCCACAGGAAAGTCCATGGCTGTTATTCAAGATATCATCGCTGCGCTCTGGCAACACGACTTTGCCGCGCTGGc            <  2:790605/93‑1 (MQ=255)
                                                                              tAAAGGCTTAAATGCTGTTCCACAGGAAAGTCCATGGCTGTTATTCAAGATATCATCGCTGCGCTCTGGCAACACGACTTTGCCGCGCTGGCg           >  1:5276728/1‑93 (MQ=255)
                                                                              tAAAGGCTTAAATGCTGTTCCACAGGAAAGTCCATGGCTGTTATTCAAGATATCATCGCTGCGCTCTGGCAACACGACTTTGCCGCGCTGGCg           <  2:591615/93‑1 (MQ=255)
                                                                               aaaGGCTTAAATGCTGTTCCACAGGAAAGTCCATGGCTGTTATTCAAGATATCATCGCTGCGCTCTGGCAACACGACTTTGCCGCGCTGGCgg          >  1:1630761/1‑93 (MQ=255)
                                                                               aaaGGCTTAAATGCTGTTCCACAGGAAAGTCCATGGCTGTTATTCAAGATATCATCGCTGCGCTCTGGCAACACGACTTTGCCGCGCTGGCgg          <  1:6222923/93‑1 (MQ=255)
                                                                               aaaGGCTTAAATGCTGTTCCACAGGAAAGTCCATGGCTGTTATTCAAGATATCATCGCTGCGCTCTGGCAACACGACTCTGCCGCGCTGGCgg          <  2:239992/93‑1 (MQ=255)
                                                                                aaGGCTTAAATGCTGTTCCACAGGAAAGTCCATGGCTGTTATTCAAGATATCATCGCTGCGCTCTGGCAACACGACTTTGCCGCGCTGGCGGa         >  2:3471415/1‑93 (MQ=255)
                                                                                 aGGCTTAAATGCTGTTCCACAGGAAAGTCCATGGCTGTTATTCAAGATATCATCGCTGCGCTCTGGCAACACGACTTTGCCGCGCTGGCGGAt        <  1:1452337/93‑1 (MQ=255)
                                                                                 aGGCTTAAATGCTGTTCCACAGGAAAGTCCATGGCTGTTATTCAAGATATCATCGCTGCGCTCTGGCAACACGACTTTGCCGCGCTGGCGGAt        <  2:1342589/93‑1 (MQ=255)
                                                                                 aGGCTTAAATGCTGTTCCACAGGAAAGTCCATGGCTGTTATTCAAGATATCATCGCTGCGCTCTGGCAACACGACTTTGCCGCGCTGGCGGAt        <  2:1994324/93‑1 (MQ=255)
                                                                                    cTTAAATGCTGTTCCACAGGAAAGTCCATGGCTGTTATTCAAGATATCATCGCTGCGCTCTGGCAACACGACTTTGCCGCGCTGGCGGATCCt     <  1:2500615/93‑1 (MQ=255)
                                                                                     ttAAATGCTGTTCCACAGGAAAGTCCATGGCTGTTATTCAAGATATCATCGCTGCGCTCTGGCAACACGACTTTGCCGCGCTGGCGGATCCTc    >  2:292004/1‑93 (MQ=255)
                                                                                       aaaTGCTGTTCCACAGGAAAGTCCATGGCTGTTATTCAAGATATCATCGCTGCGCTCTGGCAACACGACTTTGCCGCGCTGGCGGATCCTCat  <  2:464620/93‑1 (MQ=255)
                                                                                       aaaTGCTGTTCCACAGGAAAGTCCATGGCTGTTATTCAAGATATCATCGCTGCGCTCTGGCAACACGACTTTGCCGCGCTGGCGGATCCTCat  <  1:2557294/93‑1 (MQ=255)
                                                                                          |                                                                                         
CATTTTATTGTGAAAAGTCTTAAATTGTTGCGTCCGGGATCAAGGCGTCCCGGACGATTCAGGAGTACAATAGGCAGATAAAGGCTTAAACGCTGTTCCACAGGAAAGTCCATGGCTGTTATTCAAGATATCATCGCTGCGCTCTGGCAACACGACTTTGCCGCGCTGGCGGATCCTCAT  >  NC_000913/3151467‑3151646

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