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.6% 904.4 / NA 278S550G (AGT→GGT) yicFhypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (0/0);  major base C (147/130);  minor base A (1/0);  total (148/130)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 285 total aligned reads displayed.

TACGCCTTAAGCTGGGATTACCGGTCTTTTACTACAGCCACGAAGGACACTAAGGTTCAAAACCTGTGATCTGCTGGGCAGCCAGCCAACTGCCCAGCTTCTTGATTTGCGCATTTTCCTTCCATTCAATAACCTGTCTGGCGCGTCCCGATCCAGTCCCCGGCAGCTGCTGCCAGAACT  >  NC_000913/3817066‑3817245
                                                                                          |                                                                                         
tACGCTTTAAGCTGAGATTACCTGTCTTTTCCTACAGCTACGAAGGACACTAAGGTTCAAAACCTGTGATCTGCTGGGCAGCCAGCCAACCGc                                                                                         <  1:3134770/93‑1 (MQ=255)
    ccTTAAGCTGGGATTACCGGTCTTTTACTACAGCCACGAAGGACACTAAGGTTCAAAACCTGTGATCTGCTGGGCAGCCAGCCAACCGCCCAg                                                                                     >  1:3094299/1‑93 (MQ=255)
      ttAAGCTGGGATTACCGGTCTTTTACTACAGCCACGAAGGACACTAAGGTTCAAAACCTGTGATCTGCTGGGCAGCCAGCCAACCGCCCAGct                                                                                   >  1:6574026/1‑93 (MQ=255)
      ttAAGCTGGGATTACCGGTCTTTTACTACAGCCACGAAGGACACTAAGGTTCAAAACCTGTGATCTGCTGGGCAGCCAGCCAACCGCCCAGct                                                                                   <  1:2507603/93‑1 (MQ=255)
      ttAAGCTGGGATTACCGGTCTTTTACTACAGCCACGAAGGACACTAAGGTTCAAAACCTGAGATCTGCTGGGCAGCCAGCCAACCGCCCAGct                                                                                   >  1:5300751/1‑93 (MQ=255)
       taggCTGCGGTTACCGGTCTTTTACTACAGCGACGAAGGACACTAAGGTTCAAAACCTGTGATCTGCTGGGCAGCCAGCCAACCGCCCAGctt                                                                                  <  1:942005/90‑1 (MQ=255)
        aaGCTGGGATTACCGGTCTTTTACTACAGCCACGAAGGACACTAAGGTTCAAAACCTGTGATCTGCTGGGCAGCCAGCCAACCGCCCAGcttc                                                                                 <  1:5956546/93‑1 (MQ=255)
         aGCTGGGATTACCGGTCTTTTACTACAGCCACGAAGGACACTAAGGTTCAAAACCTGTGATCTGCTGGGCAGCCAGCCAACCGCCCAGcttct                                                                                >  2:637641/1‑93 (MQ=255)
         aGCTGGGATTACCGGTCTTTTACTACAGCCACGAAGGACACTAAGGTTCAAAACCTGTGATCTGCTGGGCAGCCAGCCAACCGCCCAGcttct                                                                                >  2:3460799/1‑93 (MQ=255)
         aGCTGGGATTACCGGTCTTTTACTACAGCCACGAAGGACACTAAGGTTCAAAACCTGTGATCTGCTGGGCAGCCAGCCAACCGCCCAGcttct                                                                                >  1:3984879/1‑93 (MQ=255)
         aGCTGGGATTACCGGTCTTTTACTACAGCCACGAAGGACACTAAGGTTCAAAACCTGTGATCTGCTGGGCAGCCAGCCAACCGCCCAGcttct                                                                                <  1:1563862/93‑1 (MQ=255)
           cTGGGATTACCGGTCTTTTACTACAGCCACGAAGGACACTAAGGTTCAAAACCTGTGATCTGCTGGGCAGCCAGCCAACCGCCCAGCTTCTTg                                                                              >  1:376799/1‑93 (MQ=255)
            tGGGATTACCGGTCTTTTACTACAGCCACGAAGGACACTAAGGTTCAAAACCTGTGATCTGCTGGGCAGCCAGCCAACCGCCCAGCTTCTTGa                                                                             <  2:1377153/93‑1 (MQ=255)
             gggATTACCGGTCTTTTACTACAGCCACGAAGGACACTAAGGTTCAAAACCTGTGATCTGCTGGGCAGCCAGCCAACCGCCCAGCTTCTTGAt                                                                            <  1:3020445/93‑1 (MQ=255)
             gggATTACCGGTCTTTTACTACAGCCACGAAGGACACTAAGGTTCAAAACCTGTGATCTGCTGGGCAGCCAGCCAACCGCCCAGCTTCTTGAt                                                                            >  1:5753310/1‑93 (MQ=255)
              ggATTACCGGTCTTTTACTACAGCCACGAAGGACACTAAGGTTCAAAACCTGTGATCTGCTGGGCAGCCAGCCAACCGCCCAGCTTCTTGAtt                                                                           <  2:3377087/93‑1 (MQ=255)
              ggATTACCGGTCTTTTACTACAGCCACGAAGGACACTAAGGTTCAAAACCTGTGATCTGCTGGGCAGCCAGCCAACCGCCCAGCTTCTTGAtt                                                                           <  2:5537181/93‑1 (MQ=255)
              ggATTACCGGTCTTTTACTACAGCCACGAAGGACACTAAGGTTCAAAACCTGTGATCTGCTGGGCAGCCAGCCAACCGCCCAGCTTCTTGAtt                                                                           <  2:475592/93‑1 (MQ=255)
              ggATTACCGGTCTTTTACTACAGCCACGAAGGACACTAAGGTTCAAAACCTGTGATCTGCTGGGCAGCCAGCCAACCGCCCAGCTTCTTGAtt                                                                           <  2:5191430/93‑1 (MQ=255)
              ggATTACCGGTCTTTTACTACAGCCACGAAGGACACTAAGGTTCAAAACCTGTGATCTGCTGGGCAGCCAGCCAACCGCCCAGCTTCTTGAtt                                                                           >  2:4044541/1‑93 (MQ=255)
               gATTACCGGTCTTTTACTACAGCCACGAAGGACACTAAGGTTCAAAACCTGTGATCTGCTGGGCAGCCAGCCAACCGCCCAGCTTCTTGAttt                                                                          >  2:4945886/1‑93 (MQ=255)
               gATTACCGGTCTTTTACTACAGCCACGAAGGACACTAAGGTTCAAAACCTGTGATCTGCTGGGCAGCCAGCCAACCGCCCAGCTTCTTGAttt                                                                          <  1:6364723/93‑1 (MQ=255)
               gATTACCGGTCTTTTACTACAGCCACGAAGGACACTAAGGTTCAAAACCTGTGATCTGCTGGGCAGCCAGCCAACCGCCCAGCTTCTTGAttt                                                                          <  2:4365306/93‑1 (MQ=255)
                 ttACCGGTCTTTTACTACAGCCACGAAGGACACTAAGGTTCAAAACCTGTGATCTGCTGGGCAGCCAGCCAACCGCCCAGCTTCTTGATTTgc                                                                        >  2:63517/1‑93 (MQ=255)
                 ttACCGGTCTTTTACTACAGCCACGAAGGACACTAAGGTTCAAAACCTGTGATCTGCTGGGCAGCCAGCCAACCGCCCAGCTTCTTGATTTgc                                                                        <  1:6482208/93‑1 (MQ=255)
                 ttACCGGTCTTTTACTACAGCCACGAAGGACACTAAGGTTCAAAACCTGTGATCTGCTGGGCAGCCAGCCAACCGCCCAGCTTCTTGATTTgc                                                                        >  2:6216055/1‑93 (MQ=255)
                  tACCGGTCTTTTACTACAGCCACGAAGGACACTAAGGTTCAAAACCTGTGATCTGCTGGGCAGCCAGCCAACCGCCCAGCTTCTTGATTTgcg                                                                       >  2:6622040/1‑93 (MQ=255)
                  tACCGGTCTTTTACTACAGCCACGAAGGACACTAAGGTTCAAAACCTGTGATCTGCTGGGCAGCCAGCCAACCGCCCAGCTTCTTGATTTgcg                                                                       >  2:4206803/1‑93 (MQ=255)
                   aCCGGTCTTTTACTACAGCCACGAAGGACACTAAGGTTCAAAACCTGTGATCTGCTGGGCAGCCAGCCAACCGCCCAGCTTCTTGATTTgcgc                                                                      <  2:1136237/93‑1 (MQ=255)
                   aCCGGTCTTTTACTACAGCCACGAAGGACACTAAGGTTCAAAACCTGTGATCTGCTGGGCAGCCAGCCAACCGCCCAGCTTCTTGATTTgcgc                                                                      >  2:232723/1‑93 (MQ=255)
                   aCCGGTCTTTTACTACAGCCACGAAGGACACTAAGGTTCAAAACCTGTGATCTGCTGGGCAGCCAGCCAAACGCCCAGCTTCTTGATTTgcgc                                                                      >  2:5702332/1‑93 (MQ=255)
                    ccGGTCTTTTACTACAGCCACGAAGGACACTAAGGTTCAAAACCTGTGATCTGCTGGGCAGCCAGCCAACCGCCCAGCTTCTTGATTTGCGCa                                                                     <  1:6179744/93‑1 (MQ=255)
                     cGGTCTTTTACTACAGCCACGAAGGACACTAAGGTTCAAAACCTGTGATCTGCTGGGCAGCCAGCCAACCGCCCAGCTTCTTGATTTGCGCAt                                                                    >  2:2142961/1‑93 (MQ=255)
                       gTCTTTTACTACAGCCACGAAGGACACTAAGGTTCAAAACCTGTGATCTGCTGGGCAGCCAGCCAACCGCCCAGCTTCTTGATTTGCGCAttt                                                                  >  2:5019838/1‑93 (MQ=255)
                       gTCTTTTACTACAGCCACGAAGGACACTAAGGTTCAAAACCTGTGATCTGCTGGGCAGCCAGCCAACCGCCCAGCTTCTTGATTTGCGCAttt                                                                  <  1:2295578/93‑1 (MQ=255)
                        tCTTTTACTACAGCCACGAAGGACACTAAGGTTCAAAACCTGTGATCTGCTGGGCAGCCAGCCAACCGCCCAGCTTCTTGATTTGCGCAtttt                                                                 >  2:1713376/1‑93 (MQ=255)
                          ttttAGTACAGGCACGAAAGACACTAAGGTTCAAAACCTGTGATCTGCTGGGCAGCCAGCCAACCGCCCAGCTTCTTGCTTTGAGCCTTttct                                                               >  1:3425429/1‑92 (MQ=255)
                           tttACTACAGCCACGAAGGACACTAAGGTTCAAAACCTGTGATCTGCTGGGCAGCCAGCCAACCGCCCAGCTTCTTGATTTGCGCATTttcct                                                              >  1:4248945/1‑93 (MQ=255)
                             tactacAGCCACGAAGGACACTAAGGTTCAAAACCTGTGATCTGCTGGGCAGCCAGCCAACCGCCCAGCTTCTTGATTTGCGCATTttccttc                                                            <  1:6388738/93‑1 (MQ=255)
                               ctacAGCCACGAAGGACACTAAGGTTCAAAACCTGTGATCTGCTGGGCAGCCAGCCAACCGCCCAGCTTCTTGATTTGCGCGTTTTCCTTCCa                                                          >  1:1900969/1‑93 (MQ=255)
                               ctacAGCCACGAAGGACACTAAGGTTCAAAACCTGTGATCTGCTGGGCAGCCAGCCAACCGCCCAGCTTCTTGATTTGCGCATTTTCCTTCCa                                                          <  2:3726941/93‑1 (MQ=255)
                               ctacAGCCACGAAGGACACTAAGGTTCAAAACCTGTGATCTGCTGGGCAGCCAGCCAACCGCCCAGCTTCTTGATTTGCGCATTTTCCTTCCa                                                          <  1:1593694/93‑1 (MQ=255)
                                tacAGCCACGAAGGACACTAAGGTTCAAAACCTGTGATCTGCTGGGCAGCCAGCCAACCGCCCAGCTTCTTGATTTGCGCATTTTCCTTCCAt                                                         >  2:2368981/1‑93 (MQ=255)
                                  cAGCCACGAAGGACACTAAGGTTCAAAACCTGTGATCTGCTGGGCAGCCAGCCAACCGCCCAGCTTCTTGATTTGCGCATTTTCCTTCCATTc                                                       >  2:1618683/1‑93 (MQ=255)
                                  cAGCCACGAAGGACACTAAGGTTCAAAACCTGTGATCTGCTGGGCAGCCAGCCAACCCCCCAGCTTCTTGATTTGCGCATTTTACCTCCATTc                                                       >  2:5926669/1‑93 (MQ=255)
                                   aGCCACGAAGGACACTAAGGTTCAAAACCTGTGATCTGCTGGGCAGCCAGCCAACCGCCCAGCTTCTTGATTTGCGCATTTTCCTTCCATTCa                                                      >  2:4279647/1‑93 (MQ=255)
                                   aGCCACGAAGGACACTAAGGTTCAAAACCTGTGATCTGCTGGGCAGCCAGCCAACCGCCCAGCTTCTTGATTTGCGCATTTTCCTTCCATTCa                                                      >  1:2167535/1‑93 (MQ=255)
                                    gCCACGAAGGACACTAAGGTTCAAAACCTGTGATCTGCTGGGCAGCCAGCCAACCGCCCAGCTTCTTGATTTGCGCATTTTCCTTCCATTCaa                                                     >  1:4268312/1‑93 (MQ=255)
                                     ccACGAAGGACACTAAGGTTCAAAACCTGTGATCTGCTGGGCAGCCAGCCAACCGCCCAGCTTCTTGATTTGCGCATTTTCCTTCCATTCAAt                                                    >  2:1234097/1‑93 (MQ=255)
                                      cACGAAGGACACTAAGGTTCAAAACCTGTGATCTGCTGGGCAGCCAGCCAACCGCCCAGCTTCTTGATTTGCGCATTTTCCTTCCATTCAATa                                                   >  2:2990165/1‑93 (MQ=255)
                                      cACGAAGGACACTAAGGTTCAAAACCTGTGATCTGCTGGGCAGCCAGCCAACCGCCCAGCTTCTTGATTTGCGCATTTTCCTTCCATTCAATa                                                   <  1:5220114/93‑1 (MQ=255)
                                      cACGAAGGACACTAAGGTTCAAAACCTGTGATCTGCTGGGCAGCCAGCCAACCGCCCAGCTTCTTGATTTGCGCATTTTCCTTCCATTCAATa                                                   <  1:791898/93‑1 (MQ=255)
                                       aCGAAGGACACTAAGGTTCAAAACCTGTGATCTGCTGGGCAGCCAGCCAACCGCCCAGCTTCTTGATTTGCGCATTTTCCTTCCATTCAATaa                                                  >  1:3263857/1‑93 (MQ=255)
                                         gAAGGACACTAAGGTTCAAAACCTGTGATCTGCTGGGCAGCCAGCCAACCGCCCAGCTTCTTGATTTGCGCATTTTCCTTCCATTCAATAAcc                                                >  2:293847/1‑93 (MQ=255)
                                          aaGGACACTAAGGTTCAAAACCTGTGATCTGCTGGGCAGCCAGCCAACCGCCCAGCTTCTTGATTTGCGCATTTTCCTTCCATTCAATAACCt                                               >  2:3577089/1‑93 (MQ=255)
                                           aGGACACTAAGGTTCAAAACCTGTGATCTGCTGGGCAGCCAGCCAACCGCCCAGCTTCTTGATTTGCGCATTTTCCTTCCATTCAATAACCTg                                              >  1:4104613/1‑93 (MQ=255)
                                           aGGACACTAAGGTTCAAAACCTGTGATCTGCTGGGCAGCCAGCCAACCGCCCAGCTTCTTGATTTGCGCATTTTCCTTCCATTCAATAACCTg                                              >  2:1913702/1‑93 (MQ=255)
                                           aGGACACTAAGGTTCAAAACCTGTGATCTGCTGGGCAGCCAGCCAACCGCCCAGCTTCTTGATTTGCGCATTTTCCTTCCATTCAATAACCTg                                              >  2:5073406/1‑93 (MQ=255)
                                           aGGACACTAAGGTTCAAAACCTGTGATCTGCTGGGCAGCCAGCCAACCGCCCAGCTTCTTGATTTGCGCATTTTCCTTCCATTCAATAACCTg                                              >  2:5227892/1‑93 (MQ=255)
                                           aGGACACTAAGGTTCAAAACCTGTGATCTGCTGGGCAGCCAGCCAACCGCCCAGCTTCTTGATTTGCGCATTTTCCTTCCATTCAATAACCTg                                              <  2:4622951/93‑1 (MQ=255)
                                             gACACTAAGGTTCAAAACCTGTGATCTGCTGGGCAGCCAGCCAACCGCCCAGCTTCTTGATTTGCGCATTTTCCTTCCATTCAATAACCTGTc                                            >  1:5584778/1‑93 (MQ=255)
                                             gACACTAAGGTTCAAAACCTGTGATCTGCTGGGCAGCCAGCCAACCGCCCAGCTTCTTGATTTGCGCATTTTCCTTCCATTCAATAACCTGTc                                            <  1:866952/93‑1 (MQ=255)
                                              acacTAAGGTTCAAAACCTGTGATCTGCTGGGCAGCCAGCCAACCGCCCAGCTTCTTGATTTGCGCATTTTCCTTCCATTCAATAACCTGTCt                                           >  2:1995460/1‑93 (MQ=255)
                                              acacTAAGGTTCAAAACCTGTGATCTGCTGGGCAGCCAGCCAACCGCCCAGCTTCTTGATTTGCGCATTTTCCTTCCATTCAATAACCTGTCt                                           <  1:4650651/93‑1 (MQ=255)
                                                acTAAGGTTCAAAACCTGTGATCTGCTGGGCAGCCAGCCAACCGCCCAGCTTCTTGATTTGCGCATTTTCCTTCCATTCAATAACCTGTCTgg                                         <  1:1077661/93‑1 (MQ=255)
                                                acTAAGGTTCAAAACCTGTGATCTGCTGGGCAGCCAGCCAACCGCCCAGCTTCTTGATTTGCGCATTTTCCTTCCATTCAATAACCTGTCTgg                                         >  1:3489568/1‑93 (MQ=255)
                                                acTAAGGTTCAAAACCTGTGATCTGCTGGGCAGCCAGCCAACCGCCCAGCTTCTTGATTTGCGCATTTTCCTTCCATTCAATAAACTGTCTgg                                         >  2:4636601/1‑93 (MQ=255)
                                                  tAAGGTTCAAAACCTGTGATCTGCTGGGCAGCCAGCCAACCGCCCAGCTTCTTGATTTGCGCATTTTCCTTCCATTCAATAACCTGTCTGgcg                                       <  1:1493439/93‑1 (MQ=255)
                                                  tAAGGTTCAAAACCTGTGATCTGCTGGGCAGCCAGCCAACCGCCCAGCTTCTTGATTTGCGCATTTTCCTTCCATTCAATAACCTGTCTGgcg                                       >  1:132370/1‑93 (MQ=255)
                                                    aGGTTCAAAACCTGTGATCTGCTGGGCAGCCAGCCAACCGCCCAGCTTCTTGATTTGCGCATTTTCCTTCCATTCAATAACCTGTCTGGcgcg                                     >  1:1586056/1‑93 (MQ=255)
                                                    aGGTTCAAAACCTGTGATCTGCTGGGCAGCCAGCCAACCGCCCAGCTTCTTGATTTGCGCATTTTCCTTCCATTCAATAACCTGTCTGGcgcg                                     <  1:3813950/93‑1 (MQ=255)
                                                    aGGTTCAAAACCTGTGATCTGCTGGGCAGCCAGCCAACCGCCCAGCTTCTTGATTTGCGCATTTTCCTTCCATTCAATAACCTGTCTGGcgcg                                     >  1:5711105/1‑93 (MQ=255)
                                                      gTTCAAAACCTGTGATCTGCTGGGCAGCCAGCCAACCGCCCAGCTTCTTGATTTGCGCATTTTCCTTCCATTCAATAACCTGTCTGGCGCGTc                                   <  2:5630191/93‑1 (MQ=255)
                                                       ttCAAAACCTGTGATCGGCTGGGCAGCCAGCCAACCGCCCAGCTTCTTGATTTGCGCATTTTCCTTCCATTCAATAACCTGTCTGGCGCGTcc                                  <  1:1241527/93‑1 (MQ=255)
                                                           aaaCCTGTGATCTGCTGGGCAGCCAGCCAACCGCCCAGCTTCTTGATTTGCGCATTTTCCTTCCATTCAATAACCTGTCTGGCGCGTCCCGAt                              >  1:274818/1‑93 (MQ=255)
                                                            aaCCTGTGATCTGCTGGGCAGCCAGCCAACCGCCCAGCTTCTTGATTTGCGCATTTTCCTTCCATTCAATAACCTGTCTGGCGCGTCCCGATc                             <  1:4878967/93‑1 (MQ=255)
                                                            aaCCTGTGATCTGCTGGGCAGCCAGCCAACCGCCCAGCTTCTTGATTTGCGCATTTTCCTTCCATTCAATAACCTGTCTGGCGCGTCCCGATc                             >  2:6115374/1‑93 (MQ=255)
                                                              ccTGTGATCTGCTGGGCAGCCAGCCAACCGCCCAGCTTCTTGATTTGCGCATTTTCCTTCCATTCAATAACCTGTCTGGCGCGTCCCGATCCa                           <  2:3027975/93‑1 (MQ=255)
                                                               cTGTGATCTGCTGGGCAGCCAGCCAACCGCCCAGCTTCTTGATTTGCGCATTTTCCTTCCATTCAATAACCTGTCTGGCGCGTCCCGATCCAg                          >  2:1557428/1‑93 (MQ=255)
                                                               cTGTGATCTGCTGGGCAGCCAGCCAACCGCCCAGCTTCTTGATTTGCGCATTTTCCTTCCATTCAATAACCTGTCTGGCGCGTCCCGATCCAg                          >  1:550443/1‑93 (MQ=255)
                                                               cTGTGATCTGCTGGGCAGCCAGCCAACCGCCCAGCTTCTTGATTTGCGCATTTTCCTTCCATTCAATAACCTGTCTGGCGCGTCCCGATCCAg                          >  2:1790007/1‑93 (MQ=255)
                                                               cTGTGATCTGCTGGGCAGCCAGCCAACCGCCCAGCTTCTTGATTTGCGCATTTTCCTTCCATTCAATAACCTGTCTGGCGCGTCCCGATCCAg                          >  2:4699260/1‑93 (MQ=255)
                                                                tgtgATCTGCTGGGCAGCCAGCCAACCGCCCAGCTTCTTGATTTGCGCATTTTCCTTCCATTCAATAACCTGTCTGGCGCGTCCCGATCCAGt                         <  2:4493071/93‑1 (MQ=255)
                                                                 gtgATCTGCTGGGCAGCCAGCCAACCGCCCAGCTTCTTGATTTGCGCATTTTCCTTCCATTCAATAACCTGTCTGGCGCGTCCCGATGCAGTc                        >  2:2719427/1‑93 (MQ=255)
                                                                  tgATCTGCTGGGCAGCCAGCCAACCGCCCAGCTTCTTGATTTGCGCATTTTCCTTCCATTCAATAACCTGTCTGGCGCGTCCCGATCCAGTcc                       >  2:1083082/1‑93 (MQ=255)
                                                                   gATCTGCTGGGCAGCCAGCCAACCGCCCAGCTTCTTGATTTGCGCATTTTCCTTCCATTCAATAACCTGTCTGGCGCGTCCCGATCCAGTccc                      <  1:2604078/93‑1 (MQ=255)
                                                                     tCTGCTGGGCAGCCAGCCAACCGCCCAGCTTCTTGATTTGCGCATTTTCCTTCCATTCAATAACCTGTCTGGCGCGTCCCGATCCAGTCCCCg                    >  2:286762/1‑93 (MQ=255)
                                                                       tgctgGGCAGCCAGCCAACCGCCCAGCTTCTTGATTTGCGCATTTTCCTTCCATTCAATAACCTGTCTGGCGCGTCCCGATCCAGTCCCCGGc                  >  1:2401673/1‑93 (MQ=255)
                                                                       tgctgGGCAGCCAGCCAACCGCCCAGCTTCTTGATTTGCGCATTTTCCTTCCATTCAATAACCTGTCTGGCGCGTCCCGATCCAGTCCCCGGc                  >  2:3905479/1‑93 (MQ=255)
                                                                        gctgGGCAGCCAGCCAACCGCCCAGCTTCTTGATTTGCGCATTTTCCTTCCATTCAATAACCTGTCTGGCGCGTCCCGATCCAGTCCCCGGCa                 <  1:2825454/93‑1 (MQ=255)
                                                                          tgGGCAGCCAGCCACCCGCCCAGCTTCTTGATTTGCGCATTTTCCTTCCATTCAATAACCTGTCTGGCGCGTCCCGATCCAGTCCCCGGCAgc               <  1:6008848/93‑1 (MQ=255)
                                                                          tgGGCAGCCAGCCAACCGCCCAGCTTCTTGATTTGCGCATTTTCCTTCCATTCAATAACCTGTCTGGCGCGTCCCGATCCAGTCCCCGGCAgc               <  2:4834363/93‑1 (MQ=255)
                                                                          tgGGCAGCCAGCCAACCGCCCAGCTTCTTGATTTGCGCATTTTCCTTCCATTCAATAACCTGTCTGGCGCGTCCCGATCCAGTCCCCGGCAgc               >  1:6706765/1‑93 (MQ=255)
                                                                          tgGGCAGCCAGCCAACCGCCCAGCTTCTTGATTTGCGCATTTTCCTTCCATTCAATAACCTGTCTGGCGCGTCCCGATCCAGTCCCCGGCAgc               <  2:6662767/93‑1 (MQ=255)
                                                                              cagccagcCAACCGCCCAGCTTCTTGATTTGCGCATTTTCCTTCCATTCAATAACCTGTCTGGCGCGTCCCGATCCAGTCCCCGGCAgctgct           <  2:5432933/93‑1 (MQ=255)
                                                                              cagccagcCAACCGCCCAGCTTCTTGATTTGCGCATTTTCCTTCCATTCAATAACCTGTCTGGCGCGTACCGATCCAGTCCCCGGCAgctgct           >  1:4362151/1‑93 (MQ=255)
                                                                                gccagccaACCGCCCAGCTTCTTGATTTGCGCATTTTCCTTCCATTCAATAACCTGTCTGGCGCGTCCCGATCCAGTCCCCGGCAGCtgctgc         >  1:1151941/1‑93 (MQ=255)
                                                                                   agccaACCGCCCAGCTTCTTGATTTGCGCATTTTCCTTCCATTCAATAACCTGTCTGGCGCGTCCCGATCCAGTCCCCGTCCGCTGGTGCCAg      >  2:2416517/1‑93 (MQ=255)
                                                                                       aaCCGCCCAGCTTCTTGATTTGCGCATTTTCCTTCCATTCAATAACCTGTCTGGCGCGTCCCGATCCAGTCCCCGGCAGCTGCTGCCAGAACt  <  2:5649187/93‑1 (MQ=255)
                                                                                       aaCCGCCCAGCTTCTTGATTTGCGCATTTTCCTTCCATTCAATAACCTGTCTGGCGCGTCCCGATCCAGTCCCCGGCAGCTGCTGCCAGAACt  <  2:804695/93‑1 (MQ=255)
                                                                                          |                                                                                         
TACGCCTTAAGCTGGGATTACCGGTCTTTTACTACAGCCACGAAGGACACTAAGGTTCAAAACCTGTGATCTGCTGGGCAGCCAGCCAACTGCCCAGCTTCTTGATTTGCGCATTTTCCTTCCATTCAATAACCTGTCTGGCGCGTCCCGATCCAGTCCCCGGCAGCTGCTGCCAGAACT  >  NC_000913/3817066‑3817245

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