Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,577,071 A→G *114Q (TAA→CAA)  yjiW ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,577,0710AG99.6% 866.3 / NA 240*114Q (TAA→CAA) yjiWhypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base G (116/123);  minor base T (1/0);  total (117/123)
Fisher's exact test for biased strand distribution p-value = 4.88e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 257 total aligned reads displayed.

AAATCAGACTGGTAAATATTTACTGAGGGGGAAAGTTTCTATTGAGTCAGTGGAAGGCTCCCGGTGGTTAACCGGGAGTAAACGCTGTTACGCGACTTTCTGTTTACCGGCAATCACTCCAATAAACGCCTGCACCTGCTTTTGTTTACGCGCCGACAGTTTGCACACCTGGCGTAGCG  >  NC_000913/4576982‑4577160
                                                                                         |                                                                                         
aaaTCAGACTGGTAAATATTTACTGAGGGGGAAAGTTTCTATTGAGTCAGTGGAAGGCTCCCGGTGGTTAACCGGGAGTAAACGCTGTTgcgc                                                                                        >  2:6574755/1‑93 (MQ=255)
aaaTCAGACTGGTAAATATTTACTGAGGGGGAAAGTTTCTATTGAGTCAGTGGAAGGCTCCCGGTGGTTAACCGGGAGTAAACGCTGTTgcgc                                                                                        >  1:3903867/1‑93 (MQ=255)
  aTCAGACTGGTAAATATTTACTGAGGGGGAAAGTTTCTATTGAGTCAGTGGAAGGCTCCCGGTGGTTAACCGGGAGTAAACGCTGTTGCGCGa                                                                                      >  2:1728732/1‑93 (MQ=255)
      gACTGGTAAATATTTACTGAGGGGGAAAGTTTCTATTGAGTCAGTGGAAGGGTCCCGGTGGTTAACCGGGAGGAAACGCTGTTGCGCGACttt                                                                                  >  1:1727563/1‑93 (MQ=255)
      gACTGGTAAATATTTACTGAGGGGGAAAGTTTCTATTGAGTCAGTGGAAGGCTCCCGGTGGTTAACCGGGAGTAAACGCTGTTGCGCGACttt                                                                                  >  1:5014810/1‑93 (MQ=255)
      gACTGGTAAATATTTACTGAGGGGGAAAGTTTCTATTGAGTCAGTGGAAGGCTCCCGGTGGTTAAACCAGGGCAAACGCTGTTGCGCGAGttt                                                                                  >  2:3846266/1‑93 (MQ=255)
       aCTGGTAAATATTTACTGAGGGGGAAAGTTTCTATTGAGTCAGTGGAAGGCTCCCGGTGGTTAACCGGGAGTAAACGCTGTTGCGCGACTTTc                                                                                 <  1:3380204/93‑1 (MQ=255)
       aCTGGTAAATATTTACTGAGGGGGAAAGTTTCTATTGAGTCAGTGGAAGGCTCCCGGTGGTTAACCGGGAGTAAACGCTGTTGCGCGACTTTc                                                                                 <  2:6472050/93‑1 (MQ=255)
       aCTGGTAAATATTTACTGAGGGGGAAAGTTTCTATTGAGTCAGTGGAAGGCTCCCGGTGGTTAACCGGGAGTAAACGCTGTTGCGCGACTTTc                                                                                 >  2:3150128/1‑93 (MQ=255)
       aCTGGTAAATATTTACTGAGGGGGAAAGTTTCTATTGAGTCAGTGGAAGGCTCCCGGTGGTTAACCGGGAGTAAACGCTGTTGCGCGACTTTc                                                                                 <  1:5386348/93‑1 (MQ=255)
        cTGGTAAATATTTACTGAGGGGGAAAGTTTCTATTGAGTCAGTGGAAGGCTCCCGGTGGTTAACCGGGAGTAAACGCTGTTGCGCGACTTTCt                                                                                <  2:3605350/93‑1 (MQ=255)
        cTGGTAAATATTTACTGAGGGGGAAAGTTTCTATTGAGTCAGTGGAAGGCTCCCGGTGGTTAACCGGGAGTAAACGCTGTTGCGCGACTTTCt                                                                                <  2:6182160/93‑1 (MQ=255)
         tGGTAAATATTTACTGAGGGGGAAAGTTTCTATTGAGTCAGTGGAAGGCTCCCGGTGGTTAACCGGGAGTAAACGCTGTTGCGCGACTTTCTg                                                                               >  1:1665692/1‑93 (MQ=255)
         tGGTAAATATTTACTGAGGGGGAAAGTTTCTATTGAGTCAGTGGAAGGCTCCCGGTGGTTAACCGGGAGTAAACGCTGTTGCGCGACTTTCTg                                                                               <  2:201169/93‑1 (MQ=255)
           gTAAATATTTACTGAGGGGGAAAGTTTCTATTGAGTCAGTGGAAGGCTGCCGGTGGGTAACCGGGGGGAAACGCCGCTGCGCGACTTTTTgct                                                                             >  1:42411/1‑91 (MQ=255)
           gTAAATATTTACTGAGGGGGAAAGTTTCTATTGAGTCAGTGGAAGGCTCCCGGTGGTTAACCGGGAGTAAACGCTGTTGCGCGACTTTCTGtt                                                                             <  1:6075975/93‑1 (MQ=255)
           gTAAATATTTACTGAGGGGGAAAGGTTCTATTGAGTCAGTGGAAGGCTCCCGGTGGTTAACCGGGAGTAAACGCTGTTGCGCGACTTTCTGtt                                                                             <  1:6523462/93‑1 (MQ=255)
            tAAATATTTACTGAGGGGGAAAGTTTCTATTGAGTCAGTGGAAGGCTCCCGGTGGTTAACCGGGAGTAAACGCTGTTGCGCGACTTTCTGttt                                                                            <  1:123836/93‑1 (MQ=255)
            tAAATATTTACTGAGGGGGAAAGTTTCTATTGAGTCAGTGGAAGGCTCCCGGTGGTTAACCGGGAGTAAACGCTGTTGCGCGACTTTCTGttt                                                                            <  2:4072567/93‑1 (MQ=255)
              aaTATTTACTGAGGGGGAAAGTTTCTATTGAGTCAGTGGAAGGCTCCCGGTGGTTAACCGGGAGTAAACGCTGTTGCGCGACTTTCTGTTTAc                                                                          <  2:2707768/93‑1 (MQ=255)
              aaTATTTACTGAGGGGGAAAGTTTCTATTGAGTCAGTGGAAGGCTCCCGGTGGTTAACCGGGAGTAAACGCTGTTGCGCGACTTTCTGTTTAc                                                                          <  2:3962573/93‑1 (MQ=255)
              aaTATTTACTGAGGGGGAAAGTTTCTATTGAGTCAGTGGAAGGCTCCCGGTGGTTAACCGGGAGTAAACGCTGTTGCGCGACTTTCTGTTTAc                                                                          >  1:5616550/1‑93 (MQ=255)
               atatTTACTGAGGGGGAAAGTTTCTATTGAGTCAGTGGAAGGCTCCCGGTGGTTAACCGGGAGTAAACGCTGTTGCGCGACTTTCTGTTTAcc                                                                         <  1:4254176/93‑1 (MQ=255)
                tatTTACTGAGGGGGAAAGTTTCTATTGAGTCAGTGGAAGGCTCCCGGTGGTTAACCGGGAGTAAACGCTGTTGCGCGACTTTCTGTTTACCg                                                                        >  1:4594237/1‑93 (MQ=255)
                  tttACTGAGGGGGAAAGTTTCTATTGAGTCAGTGGAAGGCTCCCGGTGGTTAACCGGGAGTAAACGCTGTTGCGCGACTTTCTGTTTACCGGc                                                                      <  2:604716/93‑1 (MQ=255)
                  tttACTGAGGGGGAAAGTTTCTATTGAGTCAGTGGAAGGCTCCCGGTGGTTAACCGGGAGTAAACGCTGTTGCGCGACTTTCTGTTTACCGGc                                                                      <  1:5284474/93‑1 (MQ=255)
                     aCTGAGGGGGAAAGTTTCTATTGAGTCAGTGGAAGGCTCCCGGTGGTTAACCGGGAGTAAACGCTGTTGCGCGACTTTTTGTTTACCGGCAAt                                                                   >  1:75910/1‑93 (MQ=255)
                     aCTGAGGGGGAAAGTTTCTATTGAGTCAGTGGAAGGCTCCCGGTGGTTAAACGGGAGTAAACGCTGTTGCGCGACTTTCTGTTTACCGGCAAt                                                                   >  2:2996662/1‑93 (MQ=255)
                          gggggAAAGTTTCTATTGAGTCAGTGGAAGGCTCCCGGTGGTTAACCGGGAGTAAACGCTGTTGCGCGACTTTCTGTTTACCGGCAATCACTc                                                              <  2:1181437/93‑1 (MQ=255)
                          gggggAAAGTTTCTATTGAGTCAGTGGAAGGCTCCCGGTGGTTAACCGGGAGTAAACGCTGTTGCGCGACTTTCTGTTTACCGGCAATCACTc                                                              >  1:791738/1‑93 (MQ=255)
                          gggggAAAGTTTCTATTGAGTCAGTGGAAGGCTCCCGGTGGTTAACCGGGAGTAAACGCTGTTGCGCGACTTTCTGTTTACCGGCAATCACTc                                                              >  2:5356518/1‑93 (MQ=255)
                           ggggAAAGTTTCTATTGAGTCAGTGGAAGGCTCCCGGTGGTTAACCGGGAGTAAACGCTGTTGCGCGACTTTCTGTTTACCGGCAATCACTcc                                                             <  1:6329834/93‑1 (MQ=255)
                            gggAAAGTTTCTATTGAGTCAGTGGAAGGCTCCCGGTGGTTAACCGGGAGTAAACGCTGTTGCGCGACTTTCTGTTTACCGGCAATCACTCCa                                                            >  1:348624/1‑93 (MQ=255)
                            gggAAAGTTTCTATTGAGTCAGTGGAAGGCTCCCGGTGGTTAACCGGGAGTAAACGCTGTTGCGCGACTTTCTGTTTACCGGCAATCACTCCa                                                            >  2:5800666/1‑93 (MQ=255)
                             ggAAAGTTTCTATTGAGTCAGTGGAAGGCTCCCGGTGGTTAACCGGGAGTAAACGCTGTTGCGCGACTTTCTGTTTACCGGCAATCACTCCaa                                                           >  2:2751110/1‑93 (MQ=255)
                             ggAAAGTTTCTATTGAGTCAGTGGAAGGCTCCCGGTGGTTAACCGGGAGTAAACGCTGTTGCGCGACTTTCTGTTTACCGGCAATCACTCCaa                                                           >  2:2209765/1‑93 (MQ=255)
                              gAAAGTTTCTATTGAGTCAGTGGAAGGCTCCCGGTGGTTAACCGGGAGTAAACGCTGTTGCGCGACTTTCTGTTTACCGGCAATCACTCCAAt                                                          <  1:5135691/93‑1 (MQ=255)
                               aaaGTTTCTATTGAGTCAGTGGAAGGCTCCCGGTGGTTAACCGGGAGTAAACGCTGTTGCGCGACTTTCTGTTTACCGGCAATCACTCCAATa                                                         >  2:4886687/1‑93 (MQ=255)
                                aaGTTTCTATTGAGTCAGTGGAAGGCTCCCGGTGGTTAACCGGGAGTAAACGCTGTTGCGCGACTTTCTGTTTACCGGCAATCACTCCAATaa                                                        <  1:3767418/93‑1 (MQ=255)
                                 aGTTTCTATTGAGTCAGTGGAAGGCTCCCGGTGGTTAACCGGGAGTAAACGCTGTTGCGCGACTTTCTGTTTACCGGCAATCACTCCAATaaa                                                       <  2:322443/93‑1 (MQ=255)
                                 aGTTTCTATTGAGTCAGTGGAAGGCTCCCGGTGGTTAACCGGGAGTAAACGCTGTTGCGCGACTTTCTGTTTACCGGCAATCACTCCAATaaa                                                       <  2:3357685/93‑1 (MQ=255)
                                 aGTTTCTATTGAGTCAGTGGAAGGCTCCCGGTGGTTAACCGGGAGTAAACGCTGTTGCGCGACTTTCTGTTTACCGGCAATCACTCCAATaaa                                                       >  2:3477045/1‑93 (MQ=255)
                                  gTTTCTATTGAGTCAGTGGAAGGCTCCCGGTGGTTAACCGGGAGTAAACGCTGTTGCGCGACTTTCTGTTTACCGGCAATCACTCCAATAAAc                                                      <  1:6116731/93‑1 (MQ=255)
                                  gTTTCTATTGAGTCAGTGGAAGGCTCCCGGTGGTTAACCGGGAGTAAACGCTGTTGCGCGACTTTCTGTTTACCGGCAATCACTCCAATAAAc                                                      <  2:1526063/93‑1 (MQ=255)
                                   tttCTATTGAGTCAGTGGAAGGCTCCCGGTGGTTAACCGGGAGTAAACGCTGTTGCGCGACTTTCTGTTTACCGGCAATCACTCCAATAAACg                                                     <  2:5106470/93‑1 (MQ=255)
                                    ttCTATTGAGTCAGTGGAAGGCTCCCGGTGGTTAACCGGGAGTAAACGCTGTTGCGCGACTTTCTGTTTACCGGCAATCACTCCAATAAACGc                                                    <  1:2847154/93‑1 (MQ=255)
                                    ttCTATTGAGTCAGTGGAAGGCTCCCGGTGGTTAACCGGGAGTAAACGCTGTTGCGCGACTTTCTGTTTACCGGCAATCACTCCAATAAACGc                                                    <  2:6355306/93‑1 (MQ=255)
                                     tCTATTGAGTCAGTGGAAGGCTCCCGGTGGTTAACCGGGAGTAAACGCTGTTGCGCGACTTTCTGTTTACCGGCAATCACTCCAATAAACGcc                                                   <  1:3984292/93‑1 (MQ=255)
                                      cTATTGAGGCAGTGGAAGGCTCCCGGTGGTTAACCGGGAGTAAACGCTGTTGCGCGACTTTCTGTTTACCGGCAATCACTCCAATAAACGCCt                                                  >  1:6467732/1‑93 (MQ=255)
                                       tATTGAGTCAGTGGAAGGCTCCCGGTGGTTAAACGGGAGTAAACGCTGTTTCGCGACTTTCTGTTTACCGGCAATCACTCCAATAAACGCCTg                                                 >  1:5756997/1‑93 (MQ=255)
                                        aTTGAGTCAGTGGAAGGCTCCCGGTGGTTAACCGGGAGGAAACGCTGTTGCGCGACTTTCTGTTTACCGGCAATCACTCCAATAAACGCCTGc                                                >  1:5860676/1‑93 (MQ=255)
                                           gAGTCAGTGGAAGGCTCCCGGTGGTTAACCGGGAGTAAACGCTGTTGCGCGACTTTCTGTTTACCGGCAATCACTCCAATAAACGCCTGCAcc                                             <  1:2431592/93‑1 (MQ=255)
                                            aGTCAGTGGAAGGCTCCCGGTGGTTAACCGGGAGTAAACGCTGTTGCGCGACTTTCTGTTTACCGGCAATCACTCCAATAAACGCCTGCACCt                                            >  2:5652291/1‑93 (MQ=255)
                                            aGTCAGTGGAAGGCTCCCGGTGGTTAACCGGGAGTAAACGCTGTTGCGCGACTTTCTGTTTACCGGCAATCACTCCAATAAACGCCTGCACCt                                            <  1:1457622/93‑1 (MQ=255)
                                             gTCAGTGGAAGGCTCCCGGTGGTTAACCGGGAGTAAACGCTGTTGCGCGACTTTCTGTTTACCGGCAATCACTCCAATAAACGCCTGCACCTg                                           >  1:2240458/1‑93 (MQ=255)
                                             gTCAGTGGAAGGCTCCCGGTGGTTAACCGGGAGTAAACGCTGTTGCGCGACTTTCTGTTTACCGGCAATCACTCCAATAAACGCCTGCACCTg                                           >  1:2177127/1‑93 (MQ=255)
                                             gTCAGTGGAAGGCTCCCGGTGGTTAACCGGGAGTAAACGCTGTTGCGCGACTTTCTGTTTACCGGCAATCACTCCAATAAACGCCTGCACCTg                                           >  1:1801009/1‑93 (MQ=255)
                                              tCAGTGGAAGGCTCCCGGTGGTTAACCGGGAGTAAACGCTGTTGCGCGACTTTCTGTTTACCGGCAATCACTCCAATAAACGCCTGCACCTGc                                          >  1:3637182/1‑93 (MQ=255)
                                               cAGTGGAAGGCTCCCGGTGGTTAACCGGGAGTAAACGCTGTTGCGCGACTTTCTGTTTACCGGCAATCACTCCAATAAACGCCTGCACCTGCt                                         <  1:4095960/93‑1 (MQ=255)
                                                 gTGGAAGGCTCCCGGTGGTTAACCGGGAGTAAACGCTGTTGCGCGACTTTCTGTTTACCGGCAATCACTCCAATAAACGCCTGCACCTGCttt                                       >  1:4836779/1‑93 (MQ=255)
                                                  tGGAAGGCTCCCGGTGGTTAACCGGGAGTAAACGCTGTTGCGCGACTTTCTGTTTACCGGCAATCACTCCAATAAACGCCTGCACCTGCtttt                                      >  2:5975743/1‑93 (MQ=255)
                                                  tGGAAGGCTCCCGGTGGTTAACCGGGAGTAAACGCTGTTGCGCGACTTTCTGTTTACCGGCAATCACTCCAATAAACGCCTGCACCTGCtttt                                      <  2:5168559/93‑1 (MQ=255)
                                                   ggAAGGCTCCCGGTGGTTAACCGGGAGTAAACGCTGTTGCGCGACTTTCTGTTTACCGGCAATCACTCCAATAAACGCCTGCACCTGCTTTTg                                     <  1:6264081/93‑1 (MQ=255)
                                                    gAAGGCTCCCGGTGGTTAACCGGGAGTAAACGCTGTTGCGCGACTTTCTGTTTACCGGCAATCACTCCAATAAACGCCTGCACCTGCTTTTGt                                    >  2:5503508/1‑93 (MQ=255)
                                                    gAAGGCTCCCGGTGGTTAACCGGGAGTAAACGCTGTTGCGCGACTTTCTGTTTACCGGCAATCACTCCAATAAACGCCTGCACCTGCTTTTGt                                    >  1:5951730/1‑93 (MQ=255)
                                                      aGGCTCCCGGTGGTTAACCGGGAGTAAACGCTGTTGCGCGACTTTCTGTTTACCGGCAATCACTCCAATAAACGCCTGCACCTGCTTTTGttt                                  >  2:267074/1‑93 (MQ=255)
                                                      aGGCTCCCGGTGGTTAACCGGGAGTAAACGCTGTTGCGCGACTTTCTGTTTACCGGCAATCACTCCAATAAACGCCTGCACCTGCTTTTGttt                                  <  1:2574277/93‑1 (MQ=255)
                                                       ggCTCCCGGTGGTTAACCGGGAGTAAACGCTGTTGCGCGACTTTCTGTTTACCGGCAATCACTCCAATAAACGCCTGCACCTGCTTTTGTTTa                                 <  2:5979/93‑1 (MQ=255)
                                                       ggCTCCCGGTGGTTAACCGGGAGTAAACGCTGTTGCGCGACTTTCTGTTTACCGGCAATCACTCCAATAAACGCCTGCACCTGCTTTTGTTTa                                 >  2:1497938/1‑93 (MQ=255)
                                                        gCTCACGGTGGTTAACCGGGAGTAAACGCTGTTGCGCGACTTTCTGTTTACCGGCAATCACTCCAATAAACGCCTGCACCTGCTTTTGTTTAc                                <  2:457270/93‑1 (MQ=255)
                                                         cTCCCGGTGGTTAACCGGGAGTAAACGCTGTTGCGCGACTTTCTGTTTACCGGCAATCACTCCAATAAACGCCTGCACCTGCTTTTGTTTAcg                               <  1:4186442/93‑1 (MQ=255)
                                                          tCCCGGTGGTTAACCGGGAGTAAACGCTGTTGCGCGACTTTCTGTTTACCGGCAATCACTCCAATAAACGCCTGCACCTGCTTTTGTTTAcgc                              >  2:6439041/1‑93 (MQ=255)
                                                          tCCCGGTGGTTAACCGGGAGTAAACGCTGTTGCGCGACTTTCTGTTTACCGGCAATCACTCCAATAAACGCCTGCACCTGCTTTTGTTTAcgc                              >  2:6261035/1‑93 (MQ=255)
                                                           cccGGTGGTTAACCGGGAGTAAACGCTGTTGCGCGACTTTCTGTTTACCGGCAATCACTCCAATAAACGCCTGCACCTGCTTTTGTTTAcgcg                             >  1:2043829/1‑93 (MQ=255)
                                                            ccGGTGGTTAACCGGGAGTAAACGCTGTTGCGCGACTTTCTGTTTACCGGCAATCACTCCAATAAACGCCTGCACCTGCTTTTGTTTACgcgc                            >  1:1976522/1‑93 (MQ=255)
                                                             cGGTGGTTAACCGGGAGTAAACGCTGTTGCGCGACTTTCTGTTTACCGGCAATCACTCCAATAAACGCCTGCACCTGCTTTTGTTTACGCGcc                           >  1:6576861/1‑93 (MQ=255)
                                                             cGGTGGTTAACCGGGAGTAAACGCTGTTGCGCGACTTTCTGTTTACCGGCAATCACTCCAACAAACGCCTGCACCTGCTTTTGTTTACGCGcc                           <  2:3275410/93‑1 (MQ=255)
                                                              ggtggtTAACCGGGAGTAAACGCTGTTGCGCGACTTTCTGTTTACCGGCAATCACTCCAATAAACGCCTGCACCTGCTTTTGTTTACGCGCCg                          >  2:5873828/1‑93 (MQ=255)
                                                              ggtggtTAACCGGGAGTAAACGCTGTTGCGCGACTTTCTGTTTACCGGCAATCACTCCAATAAACGCCTGCACCTGCTTTTGTTTACGCGCCg                          >  2:3787101/1‑93 (MQ=255)
                                                              ggtggtTAACCGGGAGTAAACGCTGTTGCGCGACTTTCTGTTTACCGGCAATCACTCCAATAAACGCCTGCACCTGCTTTTGTTTACGCGCCg                          >  2:4964873/1‑93 (MQ=255)
                                                              ggtggtTAACCGGGAGTAAACGCTGTTGCGCGACTTTCTGTTTACCGGCAATCACTCCAATAAACGCCTGCACCTGCTTTTGTTTACGCGCCg                          <  2:5309169/93‑1 (MQ=255)
                                                               gtggtTAACCGGGAGTAAACGCTGTTGCGCGACTTTCTGTTTACCGGCAATCACTCCAATAAACGCCTGCACCTGCTTTTGTTTACGCGCCGa                         <  2:2818169/93‑1 (MQ=255)
                                                                 ggtTAACCGGGAGTAAACGCTGTTGCGCGACTTTCTGTTTACCTGCAATCACTCCAATAAACGCCTGCACCTGCTTTTGTTTACGCGCCGACa                       >  1:3338401/1‑93 (MQ=255)
                                                                  gtTAACCGGGAGTAAACGCTGTTGCGCGACTTTCTGTTTACCGGCAATCACTCCAATAAACGCCTGCACCTGCTTTTGTTTACGCGCCGACAg                      >  1:3128963/1‑93 (MQ=255)
                                                                  gtTAACCGGGAGTAAAAGCTGTAGCGCGACTTTCTGTTTACCGGCAATCACTCCACTAACCGCCTGCACCTACTTTTGTTTACGCGCCGACAg                      <  1:6640822/93‑1 (MQ=255)
                                                                   ttAACCGGGAGTAAACGCTGTTGCGCGACTTTCTGTTTACCGGCAATCACTCCAATAAACGCCTGCACCTGCTTTTGTTTACGCGCCGACAGt                     <  2:828930/93‑1 (MQ=255)
                                                                      aCCGGGAGTAAACGCTGTTGCGCGACTTTCTGTTTACCGGCAATCACTCCAATAAACGCCTGCACCTGCTTTTGTTTACGCGCCGACAGTTTg                  >  1:6216617/1‑93 (MQ=255)
                                                                       ccGGGAGTAAACGCTGTTGCGCGACTTTCTGTTTACCGGCAATCACTCCAATAAACGCCTGCACCTGCTTTTGTTTACGCGCCGACAGTTTGc                 >  2:4456476/1‑93 (MQ=255)
                                                                         gggAGTAAACGCTGTTGCGCGACTTTCTGTTTACCGGCAATCACTCCAATAAACGCCTGCACCTGCTTTTGTTTACGCGCCGACAGTTTGcac               >  2:2351044/1‑93 (MQ=255)
                                                                         gggAGTAAACGCTGTTGCGCGACTTTCTGTTTACCGGCAATCACTCCAATAAACGCCTGCACCTGCTTTTGTTTACGCGCCGACAGTTTGcac               >  1:1268017/1‑93 (MQ=255)
                                                                            aGTAAACGCTGTTGCGCGACTTTCTGTTTACCGGCAATCACTCCAATAAACGCCTGCACCTGCTTTTGTTTACGCGCCGACAGCTTGCACAcc            <  2:1677936/93‑1 (MQ=255)
                                                                             gTAAACGCTGTTGCGCGACTTTCTGTTTACCGGCAATCACTCCAATAAACGCCTGCACCTGCTTTTGTTTACGCGCCGACAGTTTGCACACCt           <  1:4745143/93‑1 (MQ=255)
                                                                              tAAACGCTGTTGCGCGACTTTCTGTTTACCGGCAATCACTCCAATAAACGCCTGCACCTGCTTTTGTTTACGCGCCGACAGTTTGCACACCTg          >  2:2518753/1‑93 (MQ=255)
                                                                                aaCGCTGTTGCGCGACTTTCTGTTTACCGGCAATCACTCCAATAAACGCCTGCACCTGCTTTTGTTTACGCGCCGACAGTTTGCACACCTGGc        >  2:4206339/1‑93 (MQ=255)
                                                                                 aCGCTGTTGCGCGACTTTCTGTTTACCGGCAATCACTCCAATAAACGCCTGCACCTGCTTTTGTTTACGCGCCGACAGTTTGCACACCTGGCg       <  1:1441718/93‑1 (MQ=255)
                                                                                   gCTGTTGCGCGACTTTCTGTTTACCGGCAATCACTCCAATAAACGCCTGCACCTGCTTTTGTTTACGCGCCGACAGTTTGCACACCTGGCGTa     >  1:5708542/1‑93 (MQ=255)
                                                                                    cTGTTGCGCGACTTTCTGTTTACCGGCAATCACTCCAATAAACGCCTGCACCTGCTTTTGTTTACGCGCCGACAGTTTGCACACCTGGCGTAg    <  2:5550379/93‑1 (MQ=255)
                                                                                    cTGTTGCGCGACTTTCTGTTTACCGGCAATCACTCCAATAAACGCCTGCACCTGCTTTTGTTTACGCGCCGACAGTTTGCACACCTGGCGTAg    >  1:2456899/1‑93 (MQ=255)
                                                                                    cTGTTGCGCGACTTTCTGTTTACCGGCAATCACTACAATAAACGCCTGCACCTGCTTTTGTTTACGCGCCGACAGTTTGCACACCTGGCGTAg    >  2:2628364/1‑93 (MQ=255)
                                                                                     tGTTGCGCGACTTTCTGTTTACCGGCAATCACTCCAATAAACGCCTGCACCTGCTTTTGTTTACGCGCCGACAGTTTGCACACCTGGCGTAGc   <  2:4657507/93‑1 (MQ=255)
                                                                                         |                                                                                         
AAATCAGACTGGTAAATATTTACTGAGGGGGAAAGTTTCTATTGAGTCAGTGGAAGGCTCCCGGTGGTTAACCGGGAGTAAACGCTGTTACGCGACTTTCTGTTTACCGGCAATCACTCCAATAAACGCCTGCACCTGCTTTTGTTTACGCGCCGACAGTTTGCACACCTGGCGTAGCG  >  NC_000913/4576982‑4577160

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