Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,817,156 T→C S550G (AGT→GGT)  yicF ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,817,1560TC99.8% 1516.7 / NA 466S550G (AGT→GGT) yicFhypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (1/0);  new base C (216/249);  total (217/249)
Fisher's exact test for biased strand distribution p-value = 4.66e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 479 total aligned reads displayed.

TACGCCTTAAGCTGGGATTACCGGTCTTTTACTACAGCCACGAAGGACACTAAGGTTCAAAACCTGTGATCTGCTGGGCAGCCAGCCAACTGCCCAGCTTCTTGATTTGCGCATTTTCCTTCCATTCAATAACCTGTCTGGCGCGTCCCGATCCAGTCCCCGGCAGCTGCTGCCAGAACT  >  NC_000913/3817066‑3817245
                                                                                          |                                                                                         
 aCGCCTTAAGCTGGGATTACCGGTCTTTTACTACAGCCACGAAGGACACTAAGGTTCAAAACCTGTGATCTGCTGGGCAGCCAGCCAACCGcc                                                                                        <  2:2707592/93‑1 (MQ=255)
  cGCCTTAAGCTGGGATTACCGGTCTTTTACTACAGCCACGAAGGACACTAAGGTTCAAAACCTGTGATCTGCTGGGCAGCCAGCCAACCGccc                                                                                       <  1:4349466/93‑1 (MQ=255)
      ttAAGCTGGGATTACCGGTCTTTTACTACAGCCACGAAGGACACTAAGGTTCAAAACCTGTGATCTGCTGGGCAGCCAGCCAACCGCCCAGct                                                                                   <  1:1906536/93‑1 (MQ=255)
        aaGCTGGGATTACCGGTCTTTTACTACAGCCACGAAGGACACTAAGGTTCAAAACCTGTGATCTGCTGGGCAGCCAGCCAACCGCCCAGcttc                                                                                 <  1:45355/93‑1 (MQ=255)
         aGCTGGGATTACCGGTCTTTTACTACAGCCACGAAGGACACTAAGGTTCAAAACCTGTGATCTGCTGGGCAGCCAGCCAACCGCCCAGcttct                                                                                <  2:5787024/93‑1 (MQ=255)
          gCTGGGATTACCGGTCTTTTACTACAGCCACGAAGGACACTAAGGTTCAAAACCTGTGATCTGCTGGGCAGCCAGCCAACCGCCCAGcttctt                                                                               >  2:8712163/1‑93 (MQ=255)
          gCTGGGATTACCGGTCTTTTACTACAGCCACGAAGGACACTAAGGTTCAAAACCTGTGATCTGCTGGGCAGCCAGCCAACCGCCCAGcttctt                                                                               <  2:5530294/93‑1 (MQ=255)
            tGGGATTACCGGTCTTTTACTACAGCCACGAAGGACACTAAGGTTCAAAACCTGTGATCTGCTGGGCAGCCAGCCAACCGCCCAGCTTCTTGa                                                                             >  2:1784805/1‑93 (MQ=255)
              ggATTACCGGTCTTTTACTACAGCCACGTAGGACACTAAGGTTCAAAACCTGTGATCTGCTGGGCAGCCAGCCAACCGCCCAGCTTCTTGAtt                                                                           <  2:7661689/93‑1 (MQ=255)
              ggATTACCGGTCTTTTACTACAGCCACGAAGGACACTAAGGTTCAAAACCTGTGATCTGCTGGGCAGCCAGCCAACCGCCCAGCTTCTTGAtt                                                                           <  1:3144351/93‑1 (MQ=255)
              ggATTACCGGTCTTTTACTACAGCCACGAAGGACACTAAGGTTCAAAACCTGTGATCTGCTGGGCAGCCAGCCAACCGCCCAGCTTCTTGAtt                                                                           >  2:9379037/1‑93 (MQ=255)
               gATTACCGGTCTTTTACTACAGCCACGAAGGACACTAAGGTTCAAAACCTGTGATCTGCTGGGCAGCCAGCCAACCGCCCAGCTTCTTGAttt                                                                          <  1:1401387/93‑1 (MQ=255)
               gATTACCGGTCTTTTACTACAGCCACGAAGGACACTAAGGTTCAAAACCTGTGATCTGCTGGGCAGCCAGCCAACCGCCCAGCTTCTTGAttt                                                                          >  2:1083084/1‑93 (MQ=255)
               gATTACCGGTCTTTTACTACAGCCACGAAGGACACTAAGGTTCAAAACCTGTGATCTGCTGGGCAGCCAGCCAACCGCCCAGCTTCTTGAttt                                                                          >  1:3829475/1‑93 (MQ=255)
               gATTACCGGTCTTTTACTACAGCCACGAAGGACACTAAGGTTCAAAACCTGTGATCTGCTGGGCAGCCAGCCAACCGCCCAGCTTCTTGAttt                                                                          <  1:6270240/93‑1 (MQ=255)
                aTTACCGGTCTTTTACTACAGCCACGAAGGACACTAAGGTTCAAAACCTGTGATCTGCTGGGCAGCCAGCCAACCGCCCAGCTTCTTGATTTg                                                                         >  2:8934554/1‑93 (MQ=255)
                 ttACCGGTCTTTTACTACAGCCACGAAGGACACTAAGGTTCAAAACCTGTGATCTGATGGGCAGCCAGCCAACCGCCCAGCTTCTTGATTTgc                                                                        <  2:2534900/93‑1 (MQ=255)
                  tACCGGTCTTTTACTACAGCCACGAAGGACACTAAGGTTCAAAACCTGTGATCTGCTGGGCAGCCAGCCAACCGCCCAGCTCCTTGATTTgcg                                                                       >  1:8508313/1‑93 (MQ=255)
                   aCCGGTCTTTTACTACAGCCACGAAGGACACTAAGGTTCAAAACCTGTGATCTGCTGGGCAGCCAGCCAACCGCCCAGCTTCTTGATTTgcgc                                                                      <  2:9849637/93‑1 (MQ=255)
                    ccGGTCTTTTACTACAGCCACGAAGGACACTAAGGTTCAAAACCTGTGATCTGCTGGGCAGCCAGCCAACCGCCCAGCTTCTTGATTTGCGCa                                                                     >  1:9471587/1‑93 (MQ=255)
                    ccGGTCTTTTACTACAGCCACGAAGGACACTAAGGTTCAAAACCTGTGATCTGCTGGGCAGCCAGCCAACCGCCCAGCTTCTTGATTTGCGCa                                                                     <  1:3439013/93‑1 (MQ=255)
                    ccGGTCTTTTACTACAGCCACGAAGGACACTAAGGTTCAAAACCTGTGATCTGCTGGGCAGCCAGCCAACCGCCCAGCTTCTTGATTTGCGCa                                                                     <  2:8262110/93‑1 (MQ=255)
                    ccGGTCTTTTACTACAGCCACGAAGGACACTAAGGTTCAAAACCTGTGATCTGCTGGGCAGCCAGCCAACCGCCCAGCTTCTTGATTTCCGCa                                                                     >  1:7837688/1‑93 (MQ=255)
                      ggTCTTTTACTACAGCCACGAAGGACACTAAGGTTCAAAACCTGTGATCTGCTGGGCAGCCAGCCAACCGCCCAGCTTCTTGATTTGCGCAtt                                                                   <  1:6692923/93‑1 (MQ=255)
                      ggTCTTTTACTACAGCCACGAAGGACACTAAGGTTCAAAACCTGTGATCTGCTGGGCAGCCAGCCAACCGCCCAGCTTCTTGATTTGCGCAtt                                                                   <  1:5768857/93‑1 (MQ=255)
                        tCTTTTACTACAGCCACGAAGGACACTAAGGTTCAAAACCTGTGATCTGCTGGGCAGCCAGCCAACCGCCCAGCTTCTTGATTTGCGCCtttt                                                                 >  2:7110665/1‑93 (MQ=255)
                        tCTTTTACTACAGCCACGAAGGACACTAAGGTTCAAAACCTGTGATCTGCTGGGCAGCCAGCCAACCGCCCAGCTTCTTGATTTGCGCAtttt                                                                 <  2:386025/93‑1 (MQ=255)
                          ttttACTACAGCCACGAAGGACACTAAGGTTCAAAACCTGTGATCTGCTGGGCAGCCAGCCAACCGCCCAGCTTCTTGATTTGCGCATTttcc                                                               >  1:9250486/1‑93 (MQ=255)
                          ttttACTACAGCCACGAAGGACACTAAGGTTCAAAACCTGTGATCTGCTGGGCAGCCAGCCAACCGCCCAGCTTCTTGATTTGCGCATTttcc                                                               <  1:3775203/93‑1 (MQ=255)
                           tttACTACAGCCACGAAGGACACTAAGGTTCAAAACCTGTGATCTGCTGGGCAGCCAGCCAACCGCCCAGCTTCTTGATTTGCGCATTttcct                                                              >  2:5682864/1‑93 (MQ=255)
                             tactacAGCCACGAAGGACACTAAGGTTCAAAACCTGTGATCTGCTGGGCAGCCAGCCAACCGCCCAGCTTCTTGATTTGCGCATTttccttc                                                            <  1:9713196/93‑1 (MQ=255)
                              actacAGCCACGAAGGACACTAAGGTTCAAAACCTGTGATCTGCTGGGCAGCCAGCCAACCGCCCAGCTTCTTGATTTGCGCATTttccttcc                                                           <  2:1312660/93‑1 (MQ=255)
                              actacAGCCACGAAGGACACTAAGGTTCAAAACCTGTGATCTGCTGAGCAGCCAGCCAACCGCCCAGCTTCTTGATTTGCGCATTttccttcc                                                           >  2:290283/1‑93 (MQ=255)
                               ctacAGCCACGAAGGACACTAAGGTTCAAAACCTGTGATCTGCTGGGCAGCCAGCCAACCGCCCAGCTTCTTGATTTGCGCATTTTCCTTCCa                                                          <  1:3623855/93‑1 (MQ=255)
                               ctacAGCCACGAAGGACACTAAGGTTCAAAACCTGTGATCTGCTGGGCAGCCAGCCAACCGCCCAGCTTCTTGATTTGCGCATTTTCCTTCCa                                                          <  2:5197450/93‑1 (MQ=255)
                               ctacAGCCACGAAGGACACTAAGGTTCAAAACCTGTGATCTGCTGGGCAGCCAGCCAACCGCCCAGCTTCTTGATTTGCGCATTTTCCTTCCa                                                          >  2:2133964/1‑93 (MQ=255)
                                 acAGCCACGAAGGACACTAAGGTTCAAAACCTGTGATCTGCTGGGCAGCCAGCCAACCGCCCAGCTTCTTGATTTGCGCATTTTCCTTCCAtt                                                        >  2:9113717/1‑93 (MQ=255)
                                   aGCCACGAAGGACACTAAGGTTCAAAACCTGTGATCTGCTGGGCAGCCAGCCAACCGCCCAGCTTCTTGATTTGCGCATTTTCCTTCCATTCa                                                      <  1:6975779/93‑1 (MQ=255)
                                   aGCCACGAAGGACACTAAGGTTCAAAACCTGTGATCTGCTGGGCAGCCAGCCAACCGCCCAGCTTCTTGATTTGCGCATTTTCCTTCCATTCa                                                      <  1:5645473/93‑1 (MQ=255)
                                    gCCACGAAGGACACTAAGGTTCAAAACCTGTGATCTGCTGGGCAGCCAGCCAACCGCCCAGCTTCTTGATTTGCGCATTTTCCTTCCATTCaa                                                     <  2:3983734/93‑1 (MQ=255)
                                    gCCACGAAGGACACTAAGGTTCAAAACCTGTGATCTGCTGGGCAGCCAGCCAACCGCCCAGCTTCTTGATTTGCGCATTTTCCTTCCATTCaa                                                     >  1:2511416/1‑93 (MQ=255)
                                     ccACGAAGGACACTAAGGTTCAAAACCTGTGATCTGCTGGGCAGCCAGCCAACCGCCCAGCTTCTTGATTTGCGCATTTTCCTTCCATTCAAt                                                    <  2:1514290/93‑1 (MQ=255)
                                     ccACGAAGGACACTAAGGTTCAAAACCTGTGATCTGCTGGGCAGCCAGCCAACCGCCCAGCTTCTTGATTTGCGCATTTTCCTTCCATTCAAt                                                    <  2:4169008/93‑1 (MQ=255)
                                       aCGAAGGACACTAAGGTTCAAAACCTGTGATCTGCTGGGCAGCCAGCCAACCGCCCAGCTTCTTGATTTGCGCATTTTCCTTCCATTCAATaa                                                  >  1:7219228/1‑93 (MQ=255)
                                       aCGAAGGACACTAAGGTTCAAAACCTGTGATCTGCTGGGCAGCCAGCCAACCGCCCAGCTTCTTGATTTGCGCATTTTCCTTCCATTCAATaa                                                  <  1:8140937/93‑1 (MQ=255)
                                       aCGAAGGACACTAAGGTTCAAAACCTGTGATCTGCTGGGCAGCCAGCCAACCGCCCAGCTTCTTGATTTGCGCATTTTCCTTCCATTCAATaa                                                  >  1:2019688/1‑93 (MQ=255)
                                        cGAAGGACACTAAGGTTCAAAACCTGTGATCTGCTGGGCAGCCAGCCAACCGCCCAGCTTCTTGATTTGCGCATTTTCCTTCCATTCAATAAc                                                 <  2:4625678/93‑1 (MQ=255)
                                        cGAAGGACACTAAGGTTCAAAACCTGTGATCTGCTGGGCAGCCAGCCAACCGCCCAGCTTCTTGATTTGCGCATTTTCCTTCCATTCAATAAc                                                 >  2:7794359/1‑93 (MQ=255)
                                         gAAGGACACTAAGGTTCAAAACCTGTGATCTGCTGGGCAGCCAGCCAACCGCCCAGCTTCTTGATTTGCGCATTTTCCTTCCATTCAATAAcc                                                >  2:808521/1‑93 (MQ=255)
                                          aaGGACACTAAGGTTCAAAACCTGTGATCTGCTGGGCAGCCAGCCAACCGCCCAGCTTCTTGATTTGCGCATTTTCCTTCCATTCAATAATCt                                               >  1:8547353/1‑93 (MQ=255)
                                          aaGGACACTAAGGTTCAAAACCTGTGATCTGCTGGGCAGCCAGCCAACCGCCCAGCTTCTTGATTTGCGCATTTTCCTTCCATTCAATAACCt                                               >  1:9090405/1‑93 (MQ=255)
                                          aaGGACACTAAGGTTCAAAACCTGTGATCTGCTGGGCAGCCAGCCAACCGCCCAGCTTCTTGATTTGCGCATTTTCCTTCCATTCAATAACCt                                               <  1:1835756/93‑1 (MQ=255)
                                           aGGACACTAAGGTTCAAAACCTGTGATCTGCTGGGCAGCCAGCCAACCGCCCAGCTTCTTGATTTGCGCATTTTCCTTCCATTCAATAACCTg                                              <  2:1737429/93‑1 (MQ=255)
                                           aGGAAACTAAGGTTCAGAACCTGTGATCTGCTGCGGAGCCAGCCAACCGCTCAGCTTCTTGATTTTGGCATTTTCCTTCCATTCAATaaaccg                                              >  2:2028506/1‑89 (MQ=255)
                                            ggACACTAAGGTTCAAAACCTGTGATCTGCTGGGCAGCCAGCCAACCGCCCAGCTTCTTGATTTGCGCATTTTCCTTCCATTCAATAACCTGt                                             >  2:6813372/1‑93 (MQ=255)
                                            ggACACTAAGGTTCAAAACCTGTGATCTGCTGGGCAGCCAGCCAACCGCCCAGCTTCTTGATTTGCGCATTTTCCTTCCATTCAATAACCTGt                                             <  1:485074/93‑1 (MQ=255)
                                            ggACACTAAGGTTCAAAACCTGTGATCTGCTGGGCAGCCAGCCAACCGCCCAGCTTCTTGATTTGCGCATTTTCCTTCCATTCAATAACCTGt                                             >  2:3207289/1‑93 (MQ=255)
                                            ggACACTAAGGTTCAAAACCTGTGATCTGCTGGGCAGCCAGCCAACCGCCCAGCTTCTTGATTTGCGCATTTTCCTTCCATTCAATAACCTGt                                             <  1:8737863/93‑1 (MQ=255)
                                             gACACTAAGGTTCAAAACCTGTGATCTGCTGGGCAGCCAGCCAACCGCCCAGCTTCTTGATTTGCGCATTTTCCTTCCATTCAATAACCTGTc                                            <  2:8765734/93‑1 (MQ=255)
                                             gACACTAAGGTTCAAAACCTGTGATCTGCTGGGCAGCCAGCCAACCGCCCAGCTTCTTGATTTGCGCATTTTCCTTCCATTCAATAACCTGTc                                            >  1:755084/1‑93 (MQ=255)
                                              acacTAAGGTTCAAAACCTGTGATCTGCTGGGCAGCCAGCCAACCGCCCAGCTTCTTGATTTGCGCATTTTCCTTCCATTCAATAACCTGTCt                                           <  2:9642044/93‑1 (MQ=255)
                                                  tAAGGTTCAAAACCTGTGATCTGCTGGGCAGCCAGCCAACCGCCCAGCTTCTTGATTTGCGCATTTTCCTTCCATTCAATAACCTGTCTGgcg                                       >  1:8018894/1‑93 (MQ=255)
                                                      gTTCAAAACCTGTGATCTGCTGGGCAGCCAGCCAACCGCCCAGCTTCTTGATTTGCGCATTTTCCTTCCATTCAATAACCTGTCTGGCGCGTc                                   >  2:7582413/1‑93 (MQ=255)
                                                      gTTCAAAACCTGTGATCTGCTGGGCAGCCAGCCAACCGCCCAGCTTCTTGATTTGCGCATTTTCCTTCCATTCAATAACCTGTCTGGCGCGTc                                   <  1:6400746/93‑1 (MQ=255)
                                                      gTTCAAAACCTGTGATCTGCTGGGCAGCCAGCCAACCGCCCAGCTTCTTGATTTGCGCATTTTCCTTCCATTCAATAACCTGTCTGGCGCGTc                                   >  1:6101957/1‑93 (MQ=255)
                                                        tCAAAACCTGTGATCTGCTGGGCAGCCAGCCAACCGCCCAGCTTCTTGATTTGCGCATTTTCCTTCCATTCAATAACCTGTCTCGCGCCTccc                                 >  2:4556635/1‑93 (MQ=255)
                                                          aaaaCCTGTGATCTGCTGGGCAGCCAGCCAACCGCCCAGCTTCTTGATTTGCGCATTTTCCTTCCATTCAATAACCTGTCTGGCGCGTCCCGa                               >  2:2346429/1‑93 (MQ=255)
                                                            aaCCTGTGATCTGCTGGGCAGCCAGCCAACCGCCCAGCTTCTTGATTTGCGCATTTTCCTTCCATTCAATAACCTGTCTGGCGCGTCCCGATc                             >  2:1573077/1‑93 (MQ=255)
                                                            aaCCTGTGATCTGCTGGGCAGCCAGCCAACCGCCCAGCTTCTTGATTTGCGCATTTTCCTTCCATTCAATAACCTGTCTGGCGCGTCCCGATc                             <  2:850136/93‑1 (MQ=255)
                                                            aaCCTGTGATCTGCTGGGCAGCCAGCCAACCGCCCAGCTTCTTGATTTGCGCATTTTCCTTCCATTCAATAACCTGTCTGGCGCGTCCCGATc                             <  2:1215135/93‑1 (MQ=255)
                                                             aCCTGTGATCTGCTGGGCAGCCAGCCAACCGCCCAGCTTCTTGATTTGCGCATTTTCCTTCCATTCAATAACCTGTCTGGCGCGTCCCGATcc                            <  1:6133073/93‑1 (MQ=255)
                                                             aCCTGTGATCTGCTGGGCAGCCAGCCAACCGCCCAGCTTCTTGATTTGCGCATTTTCCTTCCATTCAATAACCTGTCTGGCGCGTCCCGATcc                            >  1:5417144/1‑93 (MQ=255)
                                                             aCCTGTGATCTGCTGGGCAGCCAGCCAACCGCCCAGCTTCTTGATTTGCGCATTTTCCTTCCATTCAATAACCTGTCTGGCGCGTCCCGATcc                            <  1:459563/93‑1 (MQ=255)
                                                              ccTGTGATCTGCTGGGCAGCCACCCAACCGCCCAGCTTCTTGATTTGCGCATTTTCCTTCCATTCAATAACCTGTCTGGCGCGTCCCGATCCa                           >  2:7390230/1‑93 (MQ=255)
                                                               cTGTGATCTGCTGGGCAGCCAGCCAACCGCCCAGCTTCTTGATTTGCGCATTTTCCTTCCATTCAATAACCTGTCTGGCGCGTCCCGATCCAg                          <  2:5978674/93‑1 (MQ=255)
                                                               cTGTGATCTGCTGGGCAGCCAGCCAACCGCCCAGCTTCTTGATTTGCGCATTTTCCTTCCATTCAATAACCTGTCTGGCGCGTCCCGATCCAg                          >  1:2340051/1‑93 (MQ=255)
                                                               cTGTGATCTGCTGGGCAGCCAGCCAACCGCCCAGCTTCTTGATTTGCGCATTTTCCTTCCATTCAATAACCTGTCTGGCGCGTCCCGATCCAg                          <  1:1533478/93‑1 (MQ=255)
                                                                tgtgATCTGCTGGGCAGCCAGCCAACCGCCCAGCTTCTTGATTTGCGCATTTTCCTTCCATTCAATAACCTGTCTGGCGCGTCCCGATCCAgc                         >  1:2195516/1‑92 (MQ=255)
                                                                tgtgATCTGCTGGGCAGCCAGCCAACCGCCCAGCTTCTTGATTTGCGCATTTTCCTTCCATTCAATAACCTGTCTGGCGCGTCCCGATCCAGt                         >  2:3623617/1‑93 (MQ=255)
                                                                 gtgATCTGCTGGGCAGCCAGCCAACCGCCCAGCTTCTTGATTTGCGCATTTTCCTTCCATTCAATAACCTGTCTGGCGCGTCCCGATCCAGTc                        >  2:8412630/1‑93 (MQ=255)
                                                                 gtgATCTGCTGGGCAGCCAGCCAACCGCCCAGCTTCTTGATTTGCGCATTTTCCTTCCATTCAATAACCTGTCTGGCGCGTCCCGATCCAGTc                        >  2:6300071/1‑93 (MQ=255)
                                                                  tgATCTGCTGGGCAGCCAGCCAACCGCCCAGCTTCTTGATTTGCGCATTTTCCTTCCATTCAATAACCTGTCTGGCGCGTCCCGATCCAGTcc                       <  1:8382122/93‑1 (MQ=255)
                                                                    aTCTGCTGGGCAGCCAGCCAACCGCCCAGCTTCTTGATTTGCGCATTTTCCTTCCATTCAATAACCTGTCTGGCGCGTCCCGATCCAGTcccc                     <  1:7725750/93‑1 (MQ=255)
                                                                     tCTGCTGGGCAGCCAGCCAACCGCCCAGCTTCTTGATTTGCGCATTTTCCTTCCATTCAATAACCTCTCTGGCGCGTCACGACCCAGTCCCCg                    >  2:6148780/1‑93 (MQ=255)
                                                                      ctgctgGGCAGCCAGCCAACCGCCCAGCTTCTTGATTTGCGCATTTTCCTTCCATTCAATAACCTGTCTGGCGCGTCCCGATCCAGTCCCCgg                   <  1:5065869/93‑1 (MQ=255)
                                                                       tgctgGGCAGCCAGCCAACCGCCCAGCTTCTTGATTTGCGCATTTTCCTTCCATTCAATAACCTGTCTGGCGCGTCCCGATCCAGTCCCCGGc                  >  1:2796555/1‑93 (MQ=255)
                                                                        gctgGGCAGCCAGCCAACCGCCCAGCTTCTTGATTTGCGCATTTTCCTTCCATTCAATAACCTGTCTGGCGCGTCCCGATCCAGTCCCCCGCa                 >  1:8900644/1‑93 (MQ=255)
                                                                           gggCAGCCAGCCAACCGCCCAGCTTCTTGATTTGCGCATTTTCCTTCCATTCAATAACCTGTCTGGCGCGTCCCGATCCAGTCCCCGGCAgct              >  2:3547228/1‑93 (MQ=255)
                                                                                gccagccaACCGCCCAGCTTCTTGATTTGCGCATTTTCCTTCCATTCAATAACCTGTCTGGCGCGTCCCGATCCAGTCCCCGGCAGCtgctgc         >  1:5289190/1‑93 (MQ=255)
                                                                                gccagccaACCGCCCAGCTTCTTGATTTGCGCATTTTCCTTCCATTCAATAACCTGTCTGGCGCGTCCCGATCCAGTCCCCGGCAGCtgctgc         >  2:8877306/1‑93 (MQ=255)
                                                                                 ccagccaACCGCCCAGCTTCTTGATTTGCGCATTTTCCTTCCATTCAATAACCTGTCTGGCGCGTCCCGATCCAGTCCCCGGCAGCTGCTGcc        <  1:7916967/93‑1 (MQ=255)
                                                                                 ccagccaACCGCCCAGCTTCTTGATTTGCGCATTTTCCTTCCATTCAATAACCTGTCTGGCGCGTCCCGATCCAGTCCCCGGCAGCTGCTGcc        <  1:10011415/93‑1 (MQ=255)
                                                                                  cagccaACCGCCCAGCTTCTTGATTTGCGCATTTTCCTTCCATTCAATAACCTGTCTGGCGCGTCCCGATCCAGTCCCCGGCAGCTGCTGCCa       <  1:4126920/93‑1 (MQ=255)
                                                                                    gccaACCGCCCAGCTTCTTGATTTGCGCATTTTCCTTCCATTCAATAACCTGTCTGGCGCGTCCCGATCCAGTCCCCGGCAGCTGCTGCCAGa     >  2:4933387/1‑93 (MQ=255)
                                                                                     ccaACCGCCCAGCTTCTTGATTTGCGCATTTTCCTTCCATTCAATAACCTGTCTGGCGCGTCCCGATCCAGTCCCCGGCAGCTGCTGCCAGaa    <  1:2889430/93‑1 (MQ=255)
                                                                                      caACCGCCCAGCTTCTTGATTTGCGCATTTTCCTTCCATTCAATAACCTGTCTGGCGCGTCCCGATCCAGTCCCCGGCAGCTGCTGCCAGAAc   <  1:9582991/93‑1 (MQ=255)
                                                                                      caACCGCCCAGCTTCTTGATTTGCGCATTTTCCTTCCATTCAATAACCTGTCTGGCGCGTCCCGATCCAGTCCCCGGCAGCTGCTGCCAGAAc   <  2:7215073/93‑1 (MQ=255)
                                                                                      caACCGCCCAGCTTCTTGATTTGCGCATTTTCCTTCCATTCAATAACCTGTCTGGCGCGTCCCGATCCAGTCCCCGGCAGCTGCTGCCAGAAc   <  1:6864874/93‑1 (MQ=255)
                                                                                      caACCGCCCAGCTTCTTGATTTGCGCATTTTCCTTCCATTCAATAACCTGTCTGGCGCGTCCCGATCCAGTCCCCGGCAGCTGCTGCCAGAAc   <  1:4473259/93‑1 (MQ=255)
                                                                                       aaCCGCCCAGCTTCTTGATTTGCGCATTTTCCTTCCATTCAATAACCTGTCTGGCGCGTCCCGATCCAGTCCCCGGCAGCTGCTGCCAGAACt  >  1:1140667/1‑93 (MQ=255)
                                                                                          |                                                                                         
TACGCCTTAAGCTGGGATTACCGGTCTTTTACTACAGCCACGAAGGACACTAAGGTTCAAAACCTGTGATCTGCTGGGCAGCCAGCCAACTGCCCAGCTTCTTGATTTGCGCATTTTCCTTCCATTCAATAACCTGTCTGGCGCGTCCCGATCCAGTCCCCGGCAGCTGCTGCCAGAACT  >  NC_000913/3817066‑3817245

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