Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 985,236 A→G Y324H (TAC→CAC)  ompF ← outer membrane protein 1a

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913985,2360AG82.4% 475.1 / 101.9 205Y324H (TAC→CAC) ompFouter membrane protein 1a
Reads supporting (aligned to +/- strand):  ref base A (19/17);  new base G (79/90);  total (98/107)
Fisher's exact test for biased strand distribution p-value = 5.83e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.74e-02

Only 100 of 215 total aligned reads displayed.

CGGTGTCGTCTGAACCTACGCCCAGTTTGTTGTCAGAATCGATCTGGTTGATGATGTAGTCAACATAGGTGGACATGTTTTTGTTGAAGTAGTAGGTTGCGCCCACTTCAAAGTAGTTCACCAGATCAACATCACCGATACCTTCTACGTCTTTCGCTTTAGATTTGGTGTAAGCGATG  >  NC_000913/985146‑985324
                                                                                          |                                                                                        
 ggTGTCGTCTGAACCTACGCCCAGTTTGTTGTCAGAATCGATCTGGTTGATGATGTAGTCAACATAGGTGGACATGTTTTTGTTGAAgtggta                                                                                       <  1:5431153/93‑1 (MQ=255)
  gtgtCGTCTGAACCTACGCCCAGTTTGTTGTCAGAATCGATCTGGTTGATGATGTAGTCAACATAGGTGGACATGTTTTTGTTGAAGtagtag                                                                                      >  1:3994966/1‑93 (MQ=255)
    gtcgtcTGAACCTACGCCCAGTTTGTTGTCAGAATCGATCTGGTTGATGATGTAGTCAACATAGGTGGACATGTTTTTGTTGAAGTGGTAGGt                                                                                    >  2:1018456/1‑93 (MQ=255)
    gtcgtcTGAACCTACGCCCAGTTTGTTGTCAGAATCGATCTGGTTGATGATGTAGTCAACATAGGTGGACATGTTTTTGTTGAAGTGGTAGGt                                                                                    <  1:2340349/93‑1 (MQ=255)
     tcgtcTGAACCTACGCCCAGTTTGTTGTCAGAATCGATCTGGTTGATGATGTAGTCAACATAGGTGGACATGTTTTTGTTGAAGTGGTAGGtt                                                                                   >  2:1524948/1‑93 (MQ=255)
     tcgtcTGAACCTACGCCCAGTTTGTTGTCAGAATCGATCTGGTTGATGATGTAGTCAACATAGGTGGACATGTTTTTGTTGAAGTGGTAGGtt                                                                                   >  2:1987172/1‑93 (MQ=255)
     tcgtcTGAACCTACGCCCAGTTTGTTGTCAGAATCGATCTGGTTGATGATGTAGTCAACATAGGTGGACATGTTTTTGTTGAAGTGGTAGGtg                                                                                   >  2:3981268/1‑92 (MQ=255)
      cgtcTGAACCTACGCCCAGTTTGTTGTCAGAATCGATCTGGTTGATGATGTAGTCAACATAGGTGGACATGTTTTTGTTGAAGTGGTAGGTTg                                                                                  >  2:5136189/1‑93 (MQ=255)
        tcTGAACCTACGCCCAGTTTGTTGTCAGAATCGATCTGGTTGATGATGTAGTCAACATAGGTGGACATGTTTTTGTTGAAGTGGTAGGTTgcg                                                                                >  1:960417/1‑93 (MQ=255)
        tcTGAACCTACGCCCAGTTTGTTGTCAGAATCGATCTGGTTGATGATGTAGTCAACATAGGTGGACATGTTTTTGTTGAAGTGGTAGGTTgcg                                                                                >  1:2720169/1‑93 (MQ=255)
        tcTGAACCTACGCCCAGTTTGTTGTCAGAATCGATCTGGTTGATGATGTAGTCAACATAGGTGGACATGTTTTTGTTGAAGTGGTAGGTTgcg                                                                                >  1:2403919/1‑93 (MQ=255)
          tGAACCTACGCCCAGTTTGTTGTCAGAATCGATCTGGTTGATGATGTAGTCAACATAGGTGGACATGTTTTTGTTGAAGTAGTAGGTTGCGcc                                                                              >  1:6040962/1‑93 (MQ=255)
            aaCCTACGCCCAGTTTGTTGTCAGAATCGATCTGGTTGATGATGTAGTCAACATAGGTGGACATGTTTTTGTTGAAGTGGTAGGTTGCGCCCa                                                                            >  2:2305294/1‑93 (MQ=255)
             aCCTACGCCCAGTTTGTTGTCAGAATCGATCTGGTTGATGATGTAGTCAACATAGGTGGACATGTTTTTGTTGAAGTGGTAGGTTGCGCCCAc                                                                           >  1:1972207/1‑93 (MQ=255)
              ccTACGCCCAGTTTGTTGTCAGAATCGATCTGGTTGATGATGTAGTCAACATAGGTGGACATGTTTTTGTTGAAGTGGTAGGTTGCGCCCACt                                                                          <  1:6329016/93‑1 (MQ=255)
              ccTACGCCCAGTTTGTTGTCAGAATCGATCTGGTTGATGATGTAGTCAACATAGGTGGACATGTTTTTGTTGAAGTAGTAGGTTGCGCCCACt                                                                          <  1:2267042/93‑1 (MQ=255)
               cTACGCCCAGTTTGTTGTCAGAATCGATCTGGTTGATGATGTAGTCAACATAGGTGGACATGTTTTTGTTGAAGTGGTAGGTTGCGCCCACtt                                                                         <  1:3910260/93‑1 (MQ=255)
                tACGCCCAGTTTGTTGTCAGAATCGATCTGGTTGATGATGTAGTCAACATAGGTGGACATGTTTTTGTTGAAGTGGTAGGTTGCGCCCACTTc                                                                        <  1:573366/93‑1 (MQ=255)
                tACGCCCAGTTTGTTGTCAGAATCGATCTGGTTGATGATGTAGTCAACATAGGTGGACATGTTTTTGTTGAAGTGGTAGGTTGCGCCCACTTc                                                                        >  2:4447086/1‑93 (MQ=255)
                tACGCCCAGTTTGTTGTCAGAATCGATCTGGTTGATGATGTAGTCAACATAGGTGGACATGTTTTTGTTGAAGTGGTAGGTTGCGCCCACTTc                                                                        <  1:690932/93‑1 (MQ=255)
                tACGCCCAGTTTGTTGTCAGAATCGATCTGGTTGATGATGTAGTCAACATAGGTGGACATGTTTTTGTTGAAGTGGTAGGTTGCGCCCACTTc                                                                        <  1:6340668/93‑1 (MQ=255)
                 aCGCCCAGTTTGTTGTCAGAATCGATCTGGTTGATGATGTAGTCAACATAGGTGGACATGTTTTTGTTGAAGTGGTAGGTTGCGCCCACTTCa                                                                       <  1:2828694/93‑1 (MQ=255)
                 aCGCCCAGTTTGTTGTCAGAATCGATCTGGTTGATGATGTAGTCAACATAGGTGGACATGTTTTTGTTGAAGTGGTAGGTTGCGCCCACTTCa                                                                       <  1:5211648/93‑1 (MQ=255)
                  cGCCCAGTTTGTTGTCAGAATCGATCTGGTTGATGATGTAGTCAACATAGGTGGACATGTTTTTGTTGAAGTGGTAGGTTGCGCCCACTTCaa                                                                      <  1:1869703/93‑1 (MQ=255)
                   gCCCAGTTTGTTGTCAGAATCGATCTGGTTGATGATGTAGTCAACATAGGTGGACATGTTTTTGTTGAAGTGGTAGGTTGCGCCCACTTCaaa                                                                     <  1:4597278/93‑1 (MQ=255)
                     ccAGTTTGTTGTCAGAATCGATCTGGTTGATGATGTAGTCAACATAGGTGGACATGTTTTTGTTGAAGTAGTAGGTTGCGCCCACTTCAAagt                                                                   >  1:2206651/1‑93 (MQ=255)
                      cAGTTTGTTGTCAGAATCGATTTGGTTGATGATGTAGTCAACATAGGTGGACATGTTTTTGTTGAAGTAGTAGGTTGCGCCCACTTCAAagta                                                                  >  2:1080699/1‑93 (MQ=255)
                      cAGTTTGTTGTCAGAATCGATCTGGTTGATGATGTAGTCAACATAGGTGGACATGTTTTTGTTGAAGTGGTAGGTTGCGCCCACTTCAAagta                                                                  <  1:4805774/93‑1 (MQ=255)
                      cAGTTTGTTGTCAGAATCGATCTGGTTGATGATGTAGTCAACATAGGTGGACATGTTTTTGTTGAAGTGGTAGGTTGCGCCCACTTCAAagta                                                                  >  1:4702640/1‑93 (MQ=255)
                        gTTTGTTGTCAGAATCGATCTGGTTGATGATGTAGTCAACATAGGTGGACATGTTTTTGTTGAAGTGGTAGGTTGCGCCCACTTCAAagtagt                                                                <  1:3436868/93‑1 (MQ=255)
                        gTTTGTTGTCAGAATCGATCTGGTTGATGATGTAGTCAACATAGGTGGACATGTTTTTGTTGAAGTGGTAGGTTGCGCCCACTTCAAagtagt                                                                >  1:4427348/1‑93 (MQ=255)
                        gTTTGTTGTCAGAATCGATCTGGTTGATGATGTAGTCAACATAGGTGGACATGTTTTTGTTGAAGTGGTAGGTTGCGCCCACTTCAAagtagt                                                                >  2:2650052/1‑93 (MQ=255)
                        gTTTGTTGTCAGAATCGATCTGGTTGATGATGTAGTCAACATAGGTGGACATGTTTTTGTTGAAGTAGTAGGTTGCGCCCACTTCAAagtagt                                                                >  2:2187992/1‑93 (MQ=255)
                           tgttgtCAGAATCGATCTGGTTGATGATGTAGTCAACATAGGTGGACATGTTTTTGTTGAAGTGGTAGGTTGCGCCCACTTCAAAGTAGTTCa                                                             <  2:6416356/93‑1 (MQ=255)
                           tgttgtCAGAATCGATCTGGTTGATGATGTAGTCAACATAGGTGGACATGTTTTTGTTGAAGTGGTAGGTTGCGCCCACTTCAAAGTAGTTCa                                                             >  1:4349127/1‑93 (MQ=255)
                            gttgtCAGAATCGATCTGGTTGATGATGTAGTCAACATAGGTGGACATGTTTTTGTTGAAGTGGTAGGTTGCGCCCACTTCAAAGTAGTTCAc                                                            <  2:3325671/93‑1 (MQ=255)
                            gttgtCAGAATCGATCTGGTTGATGATGTAGTCAACATAGGTGGACATGTTTTTGTTGAAGTGGTAGGTTGCGCCCACTTCAAAGTAGTTCAc                                                            >  1:2679097/1‑93 (MQ=255)
                             ttgtCAGAATCGATCTGGTTGATGATGTAGTCAACATAGGTGGACATGTTTTTGTTGAAGTGGTAGGTTGCGCCCACTTCAAAGTAGTTCAcc                                                           <  1:316208/93‑1 (MQ=255)
                               gtCAGAATCGATCTGGTTGATGATGTAGTCAACATAGGTGGACATGTTTTTGTTGAAGTGGTAGGTTGCGCCCACTTCAAAGTAGTTCACCAg                                                         <  2:1640773/93‑1 (MQ=255)
                               gtCAGAATCGATCTGGTTGATGATGTAGTCAACATAGGTGGACATGTTTTTGTTGAAGTGGTAGGTTGCGCCCACTTCAAAGTAGTTCACCAg                                                         >  2:4231795/1‑93 (MQ=255)
                                  aGAATCGATCTGGTTGATGATGTAGTCAACATAGGTGGACATGTTTTTGTTGAAGTGGTAGGTTGCGCCCACTTCAAAGTAGTTCACCAGATc                                                      <  2:5439274/93‑1 (MQ=255)
                                  aGAATCGATCTGGTTGATGATGTAGTCAACATAGGTGGACATGTTTTTGTTGAAGTGGTAGGTTGCGCCCACTTCAAAGTAGTTCACCAGATc                                                      >  1:6174895/1‑93 (MQ=255)
                                   gAATCGATCTGGTTGATGATGTAGTCAACATAGGTGGACATGTTTTTGTTGAAGTGGTAGGTTGCGCCCACTTCAAAGTAGTTCACCAGATCa                                                     <  2:6633766/93‑1 (MQ=255)
                                    aaTCGATCTGGTTGATGATGTAGTCAACATAGGTGGACATGTTTTTGTTGAAGTGGTAGGTTGCGCCCACTTCAAAGTAGTTCACCAGATCaa                                                    >  1:3510927/1‑93 (MQ=255)
                                    aaTCGATCTGGTTGATGATGTAGTCAACATAGGTGGACATGTTTTTGTTGAAGTGGTAGGTTGCGCCCACTTCAAAGTAGTTCACCAGATCaa                                                    <  1:1688317/93‑1 (MQ=255)
                                     aTCGATCTGGTTGATGATGTAGTCAACATAGGTGGACATGTTTTTGTTGAAGTGGTAGGTTGCGCCCACTTCAAAGTAGTTCACCAGATCAAc                                                   <  1:1923848/93‑1 (MQ=255)
                                     aTCGATCTGGTTGATGATGTAGTCAACATAGGTGGACATGTTTTTGTTGAAGTAGTAGGTTGCGCCCACTTCAAAGTAGTTCACCAGATCAAc                                                   <  1:3091649/93‑1 (MQ=255)
                                     aTCGATCTGGTTGATGATGTAGTCAACATAGGTGGACATGTTTTTGTTGAAGTAGTAGGTTGCGCCCACTTCAAAGTAGTTCACCAGATCAAc                                                   >  1:6229125/1‑93 (MQ=255)
                                      tCGATCTGGTTGATGATGTAGTCAACATAGGTGGACATGTTTTTGTTGAAGTGGTAGGTTGCGCCCACTTCAAAGTAGTTCACCAGATCAACa                                                  >  2:38240/1‑93 (MQ=255)
                                       cGATCTGGTTGATGATGTAGTCAACATAGGTGGACATGTTTTTGTTGAAGTAGTAGGTTGCGCCCACTTCAAAGTAGTTCACCAGATCAACAt                                                 <  1:1746080/93‑1 (MQ=255)
                                        gATCTGGTTGATGATGTAGTCAACATAGGTGGACATGTTTTTGTTGAAGTGGTAGGTTGCGCCCACTTCAAAGTAGTTCACCAGATCAACATc                                                >  2:5080644/1‑93 (MQ=255)
                                        gATCTGGTTGATGATGTAGTCAACATAGGTGGACATGTTTTTGTTGAAGTGGTAGGTTGCGCCCACTTCAAAGTAGTTCACCAGATCAACATc                                                <  1:6543651/93‑1 (MQ=255)
                                         aTCTGGTTGATGATGTAGTCAACATAGGTGGACATGTTTTTGTTGAAGTGGTAGGTTGCGCCCACTTCAAAGTAGTTCACCAGATCAACATCa                                               <  1:5971467/93‑1 (MQ=255)
                                         aTCTGGTTGATGATGTAGTCAACATAGGTGGACATGTTTTTGTTGAAGTAGTAGGTTGCGCCCACTTCAAAGTAGTTCACCAGATCAACATCa                                               >  1:1619524/1‑93 (MQ=255)
                                           cTGGTTGATGATGTAGTCAACATAGGTGGACATGTTTTTGTTGAAGTGGTAGGTTGCGCCCACTTCAAAGTAGTTCACCAGATCAACATCAcc                                             <  1:4272076/93‑1 (MQ=255)
                                            tGGTTGATGATGTAGTCAACATAGGTGGACATGTTTTTGTTGAAGTAGTAGGTTGCGCCCACTTCAAAGTAGTTCACCAGATCAACATCACCg                                            >  2:2885075/1‑93 (MQ=255)
                                             ggTTGATGATGTAGTCAACATAGGTGGACATGTTTTTGTTGAAGTGGTAGGTTGCGCCCACTTCAAAGTAGTTCACCAGATCAACATCACCGa                                           <  2:3050928/93‑1 (MQ=255)
                                              gTTGATGATGTAGTCAACATAGGTGGACATGTTTTTGTTGAAGTGGTAGGTTGCGCCCACTTCACAGTAGTTCACCTGATCAACATCACCGAt                                          >  1:1236538/1‑93 (MQ=255)
                                                 gatgatGTAGTCAACATAGGTGGACATGTTTTTGTTGAAGTGGTAGGTTGCGCCCACTTCAAAGTAGTTCACCAGATCAACATCACCGATAcc                                       >  2:585483/1‑93 (MQ=255)
                                                 gatgatGTAGTCAACATAGGTGGACATGTTTTTGTTGAAGTGGTAGGTTGCGCCCACTTCAAAGTAGTTCACCAGATCAACATCACCGATAcc                                       >  2:5227712/1‑93 (MQ=255)
                                                 gatgatGTAGTCAACATAGGTGGACATGTTTTTGTTGAAGTGGTAGGTTGCGCCCACTTCAAAGTAGTTCACCAGATCAACATCACCGATAcc                                       >  2:3013011/1‑93 (MQ=255)
                                                   tgatgTAGTCAACATAGGTGGACATGTTTTTGTTGAAGTGGTAGGTTGCGCCCACTTCAAAGTAGTTCACCAGATCAACATCACCGATACCtt                                     >  2:4046339/1‑93 (MQ=255)
                                                   tgatgTAGTCAACATAGGTGGACATGTTTTTGTTGAAGTGGTAGGTTGCGCCCACTTCAAAGTAGGTCACCAGATCAACATCACCGATACCtt                                     >  1:5707796/1‑93 (MQ=255)
                                                     atgTAGTCAACATAGGTGGACATGTTTTTGTTGAAGTGGTAGGTTGCGCCCACTTCAAAGTAGTTCACCAGATCAACATCACCGATACCTTCt                                   <  1:739500/93‑1 (MQ=255)
                                                         aGTCAACATAGGTGGACATGTTTTTGTTGAAGTGGTAGGTTGCGCCCACTTCAAAGTAGTTCACCAGATCAACATCACCGATACCTTCTACGt                               <  2:2701388/93‑1 (MQ=255)
                                                          gTCAACATAGGTGGACATGTTTTTGTTGAAGTGGTAGGTTGCGCCCACTTCAAAGTAGTTCACCAGATCAACATCACCGATACCTTCTACGTc                              >  1:702353/1‑93 (MQ=255)
                                                          gTCAACATAGGTGGACATGTTTTTGTTGAAGTGGTAGGTTGCGCCCACTTCAAAGTAGTTCACCAGATCAACATCACCGATACCTTCTACGTc                              >  2:4692546/1‑93 (MQ=255)
                                                           tCAACATAGGTGGACATGTTTTTGTTGAAGTGGTAGGTTGCGCCCACTTCAAAGTAGTTCACCAGATCAACATCACCGATACCTTCTACGTCt                             <  1:4237005/93‑1 (MQ=255)
                                                             aaCATAGGTGGACATGTTTTTGTTGAAGTGGTAGGTTGCGCCCACTTCAAAGTAGTTCACCAGATCAACATCACCGATACCTTCTACGTCttt                           <  1:1756630/93‑1 (MQ=255)
                                                               cATAGGTGGACATGTTTTTGTTGAAGTGGTAGGTTGCGCCCACTTCAAAGTAGTTCACCAGATCAACATCACCGATACCTTCTACGTCTTTCg                         <  1:2462159/93‑1 (MQ=255)
                                                                aTAGGTGGACATGTTTTTGTTGAAGTGGTAGGTTGCGCCCACTTCAAAGTAGTTCACCAGATCAACATCACCGATACCTTCTACGTCTTTCGc                        <  2:864239/93‑1 (MQ=255)
                                                                aTAGGTGGACATGTTTTTGTTGAAGTGGTAGGTTGCGCCCACTTCAAAGTAGTTCACCAGATCAACATCACCGATACCTTCTACGTCTTTCGc                        <  1:2028391/93‑1 (MQ=255)
                                                                 tAGGTGGACATGTTTTTGTTGAAGTGGTAGGTTGCGCCCACTTCAAAGTAGTTCACCAGATCAACATCACCGATACCTTCTACGTCTTTCGCt                       <  1:330830/93‑1 (MQ=255)
                                                                 tAGGTGGACATGTTTTTGTTGAAGTGGTAGGTTGCGCCCACTTCAAAGTAGTTCACCAGATCAACATCACCGATACCTTCTACGTCTTTCGCt                       >  2:4096924/1‑93 (MQ=255)
                                                                   ggTGGACATGTTTTTGTTGAAGTAGTAGGTTGCGCCCACTTCAAAGTAGTTCACCAGATCAACATCACCGATACCTTCTACGTCTTTCGCttt                     <  2:1276621/93‑1 (MQ=255)
                                                                   ggTGGACATGTTTTTGTTGAAGTAGTAGGTTGCGCCCACTTCAAAGTAGTTCACCAGATCAACATCACCGATACCTTCTACGTCTTTCGCttt                     <  2:4602311/93‑1 (MQ=255)
                                                                    gTGGACATGTTTTTGTTGAAGTGGTAGGTTGCGCCCACTTCAAAGTAGTTCACCAGATCAACATCACCGATACCTTCTACGTCTTTCGCTTTa                    <  2:697796/93‑1 (MQ=255)
                                                                    gTGGACATGTTTTTGTTGAAGTGGTAGGTTGCGCCCACTTCAAAGTAGTTCACCAGATCAACATCACCGATACCTTCTACGTCTTTCGCTTTa                    <  1:4212194/93‑1 (MQ=255)
                                                                    gTGGACATGTTTTTGTTGAAGTAGTAGGTTGCGCCCACTTCAAAGTAGTTCACCAGATCAACATCACCGATACCTTCTACGTCTTTCGCTTTa                    >  1:4154966/1‑93 (MQ=255)
                                                                     tGGACATGTTTTTGTTGAAGTGGTAGGTTGCGCCCACTTCAAAGTAGTTCACCAGATCAACATCACCGATACCTTCTACGTCTTTCGCTTTAg                   <  1:2113418/93‑1 (MQ=255)
                                                                     tGGACATGTTTTTGTTGAAGTGGTAGGGTGCGCCCACTTCAAAGTAATTCACCAGATCAACATCACCGATACCTTCTACGTCTTTCGCTTTAg                   >  2:1900162/1‑93 (MQ=255)
                                                                      ggACATGTTTTTGATGAAGTGGTAGGTTGCGCCCACTTCAAAGTAGTTCACCAGATCAACATCACCGATACCTTCTACGTCTTTCGCTTTAGa                  >  1:6533193/1‑93 (MQ=255)
                                                                       gACATGTTTTTGTTGAAGTGGTAGGTTGCGCCCACTTCAAAGTAGTTCACCAGATCAACATCACCGATACCTTCTACGTCTTTCGCTTTAGAt                 <  2:5369547/93‑1 (MQ=255)
                                                                       gACATGTTTTTGTTGAAGTAGTAGGTTGCGCCCACTTCAAAGTAGTTCACCAGATCAACATCACCGATACCTTCTACGTCTTTCGCTTTAGAt                 >  2:1893685/1‑93 (MQ=255)
                                                                        aCATGTTTTTGTTGAAGTGGTAGGTTGCGCCCACTTCAAAGTAGTTCACCAGATCAACATCACCGATACCTTCTACGTCTTTCGCTTTAGAtt                <  1:5107970/93‑1 (MQ=255)
                                                                         cATGTTTTTGTTGAAGTGGTAGGTTGCGCCCACTTCAAAGTAGTTCACCAGATCAACATCACCGATACCTTCTACGTCTTTCGCTTTAGAttt               >  2:3171059/1‑93 (MQ=255)
                                                                         cATGTTTTTGTTGAAGTGGTAGGTTGCGCCCACTTCAAAGTAGTTCACCAGATCAACATCACCGATACCTTCTACGTCTTTCGCTTTAGAttt               >  1:2057573/1‑93 (MQ=255)
                                                                          aTGTTTTTGTTGAAGTGGTAGGTTGCGCCCACTTCAAAGTAGTTCACCAGATCAACATCACCGATACCTTCTACGTCTTTCGCTTTAGATTTg              <  2:3411822/93‑1 (MQ=255)
                                                                           tGTTTTTGTTGAAGTGGTAGGTTGCGCCCACTTCAAAGTAGTTCACCAGATCAACATCACCGATACCTTCTACGTCTTTCGCTTTAGATTTgg             >  1:1042820/1‑93 (MQ=255)
                                                                           tGTTTTTGTTGAAGTGGTAGGTTGCGCCCACTTCAAAGTAGTTCACCAGATCAACATCACCGATACCTTCTACGTCTTTCGCTTTAGATTTgg             >  1:6411506/1‑93 (MQ=255)
                                                                            gTTTTTGTTGAAGTGGTAGGTTGCGCCCACTTCAAAGTAGTTCACCAGATCAACATCACCGATACCTTCTACGTCTTTCGCTTTAGATTTGgt            <  1:4886607/93‑1 (MQ=255)
                                                                            gTTTTTGTTGAAGTGGTAGGTTGCGCCCACTTCAAAGTAGTTCACCAGATCAACATCACCGATACCTTCTACGTCTTTCGCTTTAGATTTGgt            <  2:5631021/93‑1 (MQ=255)
                                                                               tttGTTGAAGTGGTAGGTTGCGCCCACTTCAAAGTAGTTCACCAGATCAACATCACCGATACCTTCTACGTCTTTCGCTTTAGATTTGGTGTa         <  2:5573995/93‑1 (MQ=255)
                                                                               tttGTTGAAGTGGTAGGTTGCGCCCACTTCAAAGTAGTTCACCAGATCAACATCACCGATACCTTCTACGTCTTTCGCTTTAGATTTGGTGTa         <  2:5861171/93‑1 (MQ=255)
                                                                                   ttgAAGTAGTAGGTTGCGCCCACTTCAAAGTAGTTCACCAGATCAACATCACCGATACCTTCTACGTCTTTCGCTTTAGATTTGGTGTAAGCg     <  1:43344/93‑1 (MQ=255)
                                                                                     gAAGTGGTAGGTTGCGCCCACTTCAAAGTAGTTCACCAGATCAACATCACCGATACCTTCTACGTCTTTCGCTTTAGATTTGGTGTAAGCGAt   <  1:4518794/93‑1 (MQ=255)
                                                                                     gAAGTGGTAGGTTGCGCCCACTTCAAAGTAGTTCACCAGATCAACATCACCGATACCTTCTACGTCTTTCGCTTTAGATTTGGTGTAAGCGAt   <  2:4199827/93‑1 (MQ=255)
                                                                                      aaGTGGTAGGTTGCGCCCACTTCAAAGTAGTTCACCAGATCAACATCACCGATACCTTCTACGTCTTTCGCTTTAGATTTGGTGTAAGCGATg  >  2:1448551/1‑93 (MQ=255)
                                                                                      aaGTGGTAGGTTGCGCCCACTTCAAAGTAGTTCACCAGATCAACATCACCGATACCTTCTACGTCTTTCGCTTTAGATTTGGTGTAAGCGATg  <  2:3710100/93‑1 (MQ=255)
                                                                                      aaGTAGTAGGTTGCGCCCACTTCAAAGTAGTTCACCAGATCAACATCACCGATACCTTCTACGTCTTTCGCTTTAGATTTGGTGTAAGCGATg  <  1:2443882/93‑1 (MQ=255)
                                                                                          |                                                                                        
CGGTGTCGTCTGAACCTACGCCCAGTTTGTTGTCAGAATCGATCTGGTTGATGATGTAGTCAACATAGGTGGACATGTTTTTGTTGAAGTAGTAGGTTGCGCCCACTTCAAAGTAGTTCACCAGATCAACATCACCGATACCTTCTACGTCTTTCGCTTTAGATTTGGTGTAAGCGATG  >  NC_000913/985146‑985324

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