Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 195,014 A→G I38V (ATT→GTT)  yaeS → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913195,0140AG98.5% 603.1 / ‑1.0 193I38V (ATT→GTT) yaeShypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (1/2);  new base G (97/93);  total (98/95)
Fisher's exact test for biased strand distribution p-value = 6.17e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.38e-01

Only 100 of 199 total aligned reads displayed.

CCACTTAGCGAAAAATTGCCAGCGCATGGCTGCCGTCATGTTGCGATCATTATGGACGGCAATGGCCGCTGGGCAAAAAAGCAAGGGAAGATTCGTGCCTTTGGGCATAAAGCCGGGGCAAAATCCGTCCGCCGGGCTGTCTCTTTTGCGGCCAACAACGGTATTGAGGCGTTAACGCTG  >  NC_000913/194924‑195103
                                                                                          |                                                                                         
ccACTTAGCGAAAAATTGCCAGCGCATGGCTGCCGTCATGTTGCGATCATTATGGACGGCAATGGCCGCTGGGCAAAAAAGCAAGGGAAGGtt                                                                                         <  2:2179254/93‑1 (MQ=255)
ccACTTAGCGAAAAATTGCCAGCGCATGGCTGCCGTCATGTTGCGATCATTATGGACGGCAATGGCCGCTGGGCAAAAAAGCAAGGGAAGGtt                                                                                         <  2:1006315/93‑1 (MQ=255)
 cACTTAGCGAAAAATTGCCAGCGCATGGCTGCCGTCATGTTGCGATCATTATGGACGGCAATGGCCGCTGGGCAAAAAAGCAAGGGAAGGTTc                                                                                        >  2:4570934/1‑93 (MQ=255)
 cACTTAGCGAAAAATTGCCAGCGCATGGCTGCCGTCATGTTGCGATCATTATGGACGGCAATGGCCGCTGGGCAAAAAAGAAAGGGAAGGTTc                                                                                        >  2:1547947/1‑93 (MQ=255)
  aCTTAGCGAAAAATTGCCAGCGCATGGCTGCCGTCATGTTGCGATCATTATGGACGGCAATGGCCGCTGGGCAAAAAAGCAAGGGAAGGTTCg                                                                                       >  1:3065420/1‑93 (MQ=255)
  aCTTAGCGAAAAATTGCCAGCGCATGGCTGCCGTCATGTTGCGATCATTATGGACGGCAATGGCCGCTGGGCAAAAAAGCAAGGGAAGGTTCg                                                                                       <  1:2724075/93‑1 (MQ=255)
  aCTTAGCGAAAAATTGCCAGCGCATGGCTGCCGTCATGTTGCGATCATTATGGACGGCAATGGCCGCTGGGCAAAAAAGCAAGGGAAGGTTCg                                                                                       <  1:3857522/93‑1 (MQ=255)
   cTTAGCGAAAAATTGCCAGCGCATGGCTGCCGTCATGTTGCGATCATTATGGACGGCAATGGCCGCTGGGCAAAAAAGCAAGGGAAGGTTCGt                                                                                      <  2:3858390/93‑1 (MQ=255)
   cTTAGCGAAAAATTGCCAGCGCATGGCTGCCGTCATGTTGCGATCATTATGGACGGCAATGGCCGCTGGGCAAAAAAGCAAGGGAAGGTTCGt                                                                                      >  2:3072462/1‑93 (MQ=255)
    ttAGCGAAAAATTGCCAGCGCATGGCTGCCGTCATGTTGCGATCATTATGGACGGCAATGGCCGCTGGGCAAATAAGTAAGGGAAGGCTCGTg                                                                                     >  1:4213584/1‑93 (MQ=255)
      aGCGAAAAATTGCCAGCGCATGGCTGCCGTCATGTTGCGATCATTATGGACGGCAATGGCCGCTGGGCAAAAAAGCAAGGGAAGGTTCGTGcc                                                                                   <  2:834372/93‑1 (MQ=255)
       gCGAAAAATTGCCAGCGCATGGCTGCCGTCATGTTGCGATCATTATGGACGGCAATGGCCGCTGGGCAAAAAAGCAAGGGAAGGTTCGTGCCt                                                                                  <  1:4180810/93‑1 (MQ=255)
        cGAAAAATTGCCAGCGCATGGCTGCCGTCATGTTGCGATCATTATGGACGGCAATGGCCGCTGGGCAAAAAAGCAAGGGAAGGTTCGTGCCtt                                                                                 <  1:4699430/93‑1 (MQ=255)
         gAAAAATTGCCAGCGCATGGCTGCCGTCATGTTGCGATCATTATGGACGGCAATGGCCGCTGGGCAAAAAAGCAAGGGAAGGTTCGTGCCttt                                                                                <  2:2888192/93‑1 (MQ=255)
         gAAAAATTGCCAGCGCATGGCTGCCGTCATGTTGCGATCATTATGGACGGCAATGGCCGCTGGGCAAAAAAGCAAGGGAAGGTTCGTGCCttt                                                                                >  1:4012073/1‑93 (MQ=255)
          aaaaaTTGCCAGCGCATGGCTGCCGTCATGTTGCGATCATTATGGACGGCAATGGCCGCTGGGCAAAAAAGCAAGGGAAGGTTCGTGCCTTTg                                                                               <  2:1725223/93‑1 (MQ=255)
          aaaaaTTGCCAGCGCATGGCTGCCGTCATGTTGCGATCATTATGGACGGCAATGGCCGCTGGGCAAAAAAGCAAGGGAAGGTTAGTGCCTTTg                                                                               >  1:3085054/1‑93 (MQ=255)
           aaaaTTGCCAGCGCATGGCTGCCGTCATGTTGCGATCATTATGGACGGCAATGGCCGCTGGGCAAAAAAGCAAGGGAAGGTTCGTGCCTTTgg                                                                              >  1:2365382/1‑93 (MQ=255)
           aaaaTTGCCAGCGCATGGCTGCCGTCATGTTGCGATCATTATGGACGGCAATGGCCGCTGGGCAAAAAAGCAAGGGAAGGTTCGTGCCTTTgg                                                                              <  1:3606643/93‑1 (MQ=255)
           aaaaTTGCCAGCGCATGGCTGCCGTCATGTTGCGATCATTATGGACGGCAATGGCCGCTGGGCAAAAAAGCAAGGGAAGATTCGTGCCTTTgg                                                                              <  1:3922386/93‑1 (MQ=255)
           aaaaTTGCAAGCGCAGGGCTGCCGTCATGTTGCGATCATGATGGACGGCAATGGCCGCTGGGCAAAAAAGCAAGGGAAGGTTCGTGCCTTTgg                                                                              <  1:2596063/93‑1 (MQ=255)
                  ccAGCGCATGGCTGCCGTCATGTTGCGATCATTATGGACGGCAATGGCCGCTGGGCAAAAAAGCAAGGGAAGGTTCGTGCCTTTGGGCATaaa                                                                       <  1:863774/93‑1 (MQ=255)
                     gcgcATGGCTGCCGTCATGTTGCGATCATTATGGACGGCAATGGCCGCTGGGCAAAAAAGCAAGGGAAGGTTCGTGCCTTTGGGCATAAAGcc                                                                    <  1:3552309/93‑1 (MQ=255)
                     gcgcATGGCTGCCGTCAGGTTGCGATCATTATGGACGGCAATGGCCGCTGTGCAAAAAAGCAGGGGCAGGTTCGTGCCTTTGGGCATAAAGcc                                                                    >  1:1803787/1‑93 (MQ=255)
                        cATGGCTGCCGTCATGTTGCGATCATTATGGACGGCAATGGCCGCTGGGCAAAAAAGCAAGGGAAGGTTCGTGCCTTTGGGCATAAAGCCggg                                                                 >  2:2446673/1‑93 (MQ=255)
                        cATGGCTGCCGTCATGTTGCGATCATTATGGACGGCAATGGCCGCTGGGCAAAAAAGCAAGGGAAGGTTCGTGCCTTTGGGCATAAAGCCggg                                                                 >  1:2925146/1‑93 (MQ=255)
                        cATGGCTGCCGTCATGTTGCGATCATTATGGACGGCAATGGCCGCTGGGCAAAAAAGCAAGGGAAGGTTCGTGCCTTTGGGCATAAAGCCggg                                                                 >  1:3466296/1‑93 (MQ=255)
                        cATGGCTGCCGTCATGTTGCGATCATTAAGGACGGCAATGGCCGCTGGGCAAAAAAGCAAGGGAAGGTTCGTGCCTTTGGGCATAAAGCCtgg                                                                 >  2:162111/1‑93 (MQ=255)
                        cATGCCTGCCGTCATGTTGCGATCATTATGGACGGCAATGGCCGCTGGGCAAAAAAGCAAGGGAAGGTTCGTGCCTTTGGGCATAAAGCCggg                                                                 <  1:3635336/93‑1 (MQ=255)
                         aTGGCTGCGGTCATGTTGCGATCATTATGGACGGCAATGGCCGCTGGGCAAAAAAGCAAGGGAAGGTTCGTGCCTTTGGGCATAAAGCCgggg                                                                <  2:4532879/93‑1 (MQ=255)
                            gCTGCCGTCATGTTGCGATCATTATGGACGGCAATGGCCGCTGGGCAAAAAAGCAAGGGAAGGTTCGTGCCTTTGGGCATAAAGCCGGGGCaa                                                             <  2:3793971/93‑1 (MQ=255)
                             cTGCCGTCATGTTGCGATCATTATGGACGGCAATGGCCGCTGGGCAAAAAAGCAAGGGAAGGTTCGTGCCTTTGGGCATAAAGCCGGGGCaaa                                                            >  2:2283314/1‑93 (MQ=255)
                              tGCCGTCATGTTGCGATCATTATGGACGGCAATGGCCGCTGGGCAAAAAAGCAAGGGAAGGTTCGTGCCTTTGGGCATAAAGCCGGGGCaaaa                                                           >  1:1454737/1‑93 (MQ=255)
                              tGCCGTCATGTTGCGATCATTATGGACGGCAATGGCCGCTGGGCAAAAAAGCAAGGGAAGGTTCGTGCCTTTGGGCATAAAGCCGGGGCaaaa                                                           >  1:1663193/1‑93 (MQ=255)
                                   tCATGTTGCGATCATTATGGACGGCAATGGCCGCTGGGCAAAAAAGCAAGGGAAGGTTCGTGCCTTTGGGCATAAAGCCGGGGCAAAAtccgt                                                      >  1:4507495/1‑93 (MQ=255)
                                    cATGTTGCGATCATTATGGACGGCAATGGCCGCTGGGCAAAAAAGCAAGGGAAGGTTCGTGCCTTTGGGCATAAAGCCGGGGCAAAAtccgtc                                                     <  2:1590923/93‑1 (MQ=255)
                                    cATGTTGCGATCATTATGGACGGCAATGGCCGCTGGGCAAAAAAGCAAGGGAAGGTTCGTGCCTTTGGGCATAAAGCCGGGGCAAAAtccgtc                                                     <  2:3316221/93‑1 (MQ=255)
                                    cATGTTGCGATCATTATGGACGGCAATGGCCGCTGGGCAAAAAAGCAAGGGAAGGTTCGTGCCTTTGGGCATAAAGCCGGGGCAAAAtccgtc                                                     >  1:4302461/1‑93 (MQ=255)
                                     aTGTTGCGATCATTATGGACGGCAATGGCCGCTGGGCAAAAAAGCAAGGGAAGGTTCGTGCCTTTGGGCATAAAGCCGGGGCAAAAtccgtcc                                                    <  2:4030961/93‑1 (MQ=255)
                                     aTGTTGCGATCATTATGGACGGCAATGGCCGCTGGGCAAAAAAGCAAGGGAAGGGTCGTGCCTTTGGGCATAAAGCCGGGGCAAAAtccgtcc                                                    >  2:2648377/1‑93 (MQ=255)
                                      tGTTGCGATCATTATGGACGGCAATGGCCGCTGGGCAAAAAAGCAAGGGAAGGTTCGTGCCTTTGGGCATAAAGCCGGGGCAAAAtccgtccg                                                   >  2:1200481/1‑93 (MQ=255)
                                       gTTGCGATCATTATGGACGGCAATGGCCGCTGGGCAAAAAAGCAAGGGAAGGTTCGTGCCTTTGGGCATAAAGCCGGGGCAAAATCCGTccgc                                                  >  1:4547301/1‑93 (MQ=255)
                                          gCGATCATTATGGACGGCAATGGCCGCTGGGCAAAAAAGCAAGGGAAGGTTCGTGCCTTTGGGCATAAAGCCGGGGCAACATCCGTCCGCCgg                                               >  2:520782/1‑93 (MQ=255)
                                          gCGATCATTATGGACGGCAATGGCCGCTGGGCAAAAAAGCAAGGGAAGGTTCGTGCCTTTGGGCATAAAGCCGGGGCAAAATCCGTCCGCCgg                                               >  2:2482699/1‑93 (MQ=255)
                                            gATCATTATGGACGGCAATGGCCGCTGGGCAAAAAAGCAAGGGAAGGTTCGTGCCTTTGGGCATAAAGCCGGGGCAAAATCCGTCCGCCGGGc                                             <  1:3771018/93‑1 (MQ=255)
                                             aTCATTATGGACGGCAATGGCCGCTGGGCAAAAAAGCAAGGGAAGGTTCGTGCCTTTGGGCATAAAGCCGGGGCAAAATCCGTCCGCCGGGCt                                            <  1:700068/93‑1 (MQ=255)
                                             aTCATTATGGACGGCAATGGCCGCTGGGCAAAAAAGCAAGGGAAGGTTCGTGCCTTTGGGCATAAAGCCGGGGCAAAATCCGTCCGCCGGGCt                                            >  1:668833/1‑93 (MQ=255)
                                              tCATTATGGACGGCAATGGCCGCTGGGCAAAAAAGCAAGGGAAGGTTCGTGCCTTTGGGCATAAAGCCGGGGCAAAATCCGTCCGCCGGGCTg                                           <  1:2448010/93‑1 (MQ=255)
                                              tCATTATGGACGGCAATGGCCGCTGGGCAAAAAAGCAAGGGAAGGTTCGTGCCTTTGGGCATAAAGCCGGGGCAAAATCCGTCCGCCGGGCTg                                           >  2:1766855/1‑93 (MQ=255)
                                              tCATTATGGACGGCAATGGCCGCTGGGCAAAAAAGCAAGGGAAGGTTCGTGCCTTTGGGCATAAAGCCGGGGCAAAATCCGTCCGCCGGGCTg                                           >  1:3259686/1‑93 (MQ=255)
                                                aTTATGGACGGCAATGGCCGCTGGGCAAAAAAGCAAGGGAAGGTTCGTGCCTTTGGGCATAAAGCCGGGGCAAAATCCGTCCGCCGGGCTGtc                                         >  1:3965403/1‑93 (MQ=255)
                                                aTTATGGACGGCAATGGCCGCTGGGCAAAAAAGCAAGGGAAGGTTCGTGCCTTTGGGCATAAAGCCGGGGCAAAATCCGTCCGCCGGGCTGtc                                         <  1:1852718/93‑1 (MQ=255)
                                                     ggACGGCAATGGCCGCTGGGCAAAAAAGCAAGGGAAGGTTCGTGCCTTTGGGCATAAAGCCGGGGCAAAATCCGTCCGCCGGGCTGTCTCttt                                    <  2:740214/93‑1 (MQ=255)
                                                     ggACGGCAATGGCCGCTGGGCAAAAAAGCAAGGGAAGGTTCGTGCCTTTGGGCATAAAGCCGGGGCAAAATCCGTCCGCCGGGCTGTCTCttt                                    <  2:4145353/93‑1 (MQ=255)
                                                      gACGGCAATGGCCGCTGGGCAAAAAAGCAAGGGAAGGTTCGTGCCTTTGGGCATAAAGCCGGGGCAAAATCCGTCCGCCGGGCTGTCTCtttt                                   >  1:1443133/1‑93 (MQ=255)
                                                         ggCAATGGCCGCTGGGCAAAAAAGCAAGGGAAGGTTCGTGCCTTTGGGCATAAAGCCGGGGCAAAATCCGTCCGCCGGGCTGTCTCTTTTGcc                                >  1:2829629/1‑92 (MQ=255)
                                                         ggCAATGGCCGCTGGGCAAAAAAGCAAGGGAAGGTTCGTGCCTTTGGGCATAAAGCCGGGGCAAAATCCGTCCGCCGGGCTGTCTCTTTTGCg                                >  2:2763199/1‑93 (MQ=255)
                                                          gCAATGGCCGCTGGGCAAAAAAGCAAGGGAAGGTTCGTGCCTTTGGGCATAAAGCCGGGGCAAAATCCGTCCGCCGGGCTGTCTCTTTTGCgg                               >  2:4679871/1‑93 (MQ=255)
                                                           cAATGGCCGCTGGGCAAAAAAGCAAGGGAAGGTTCGTGCCTTTGGGCATAAAGCCGGGGCAAAATCCGTCCGCCGGGCTGTCTCTTTTGCGGc                              <  2:1106741/93‑1 (MQ=255)
                                                            aaTGGCCGCTGGGCAAAAAAGCAAGGGAAGGTTCGTGCCTTTGGGCATAAAGCCGGGGCAAAATCCGTCCGCCGGGCTGTCTCTTTTGCGGcc                             >  1:2005475/1‑93 (MQ=255)
                                                            aaTGGCCGCGGGGCAAAAAAGCAAGGGAAGGTTCGTGCCTTTGGGCATAAAGCCGGGGCAAAATCCGTCCGCCGGGCTGTCTCTTTTGCGGcc                             <  1:2058321/93‑1 (MQ=255)
                                                             aTGGCCGCTGGGCAAAAAAGCAAGGGAAGGTTCGTGCCTTTGGGCATAAAGCCGGGGCAAAATCCGTCCGCCGGGCTGTCTCTTTTGCGGCca                            >  2:3954367/1‑93 (MQ=255)
                                                              tGGCCGCTGGGCAAAAAAGCAAGGGAAGGTTCGTGCCTTTGGGCATAAAGCCGGGGCAAAATCCGTCCGCCGGGCTGTCTCTTTTGCGGCcaa                           >  2:1466623/1‑93 (MQ=255)
                                                              tGGCCGCTGGGCAAAAAAGCAAGGGAAGGTTCGTGCCCTTGGGCATAAAGCCGGGGCAAAATCCGTCCGCCGGGCTGTCTCTTTTGCGGCcaa                           >  2:3166230/1‑93 (MQ=255)
                                                               ggCCGCTGGGCAAAAAAGCAAGGGAAGGTTCGTGCCTTTGGGCATAAAGCCGGGGCAAAATCCGTCCGCCGGGCTGTCTCTTTTGCGGCcaac                          <  2:3557461/93‑1 (MQ=255)
                                                                gCCGCTGGGCAAAAAAGCAAGGGAAGGTTCGTGCCTTTGGGCATAAAGCCGGGGCAAAATCCGTCCGCCGGGCTGTCTCTTTTGCGGCcaaca                         >  1:655653/1‑93 (MQ=255)
                                                                  cGCTGGGCAAAAAAGCAAGGGAAGGTTCGTGCCTTTGGGCATAAAGCCGGGGCAAAATCCGTCCGCCGGGCTGTCTCTTTTGCGGCCaacaac                       >  1:155354/1‑93 (MQ=255)
                                                                  cGCTGGGCAAAAAAGCAAGGGAAGGTTCGTGCCTTTGGGCATAAAGCCGGGGCAAAATCCGTCCGCCGGGCTGTCTCTTTTGCGGCCaacaac                       >  1:155048/1‑93 (MQ=255)
                                                                  cGCTGGGCAAAAAAGCAAGGGAAGGTTCGTGCCTTTGGGCATAAAGCCGGGGCAAAATCCGTCCGCCGGGCTGTCTCTTTTGCGGCCaacaac                       <  1:4321367/93‑1 (MQ=255)
                                                                   gCTGGGCAAAAAAGCAAGGGAAGGTTCGTGCCTTTGGGCATAAAGCCGGGGCAAAATCCGTCCGCCGGGCTGTCTCTTTTGCGGCCAACAACg                      <  2:3676401/93‑1 (MQ=255)
                                                                   gCTGGGCAAAAAAGCAAGGGAAGGTTCGTGCCTTTGGGCATAAAGCCGGGGCAAAATCCGTCCGCCGGGCTGTCTCTTTTGCGGCCAACAACg                      <  1:3290576/93‑1 (MQ=255)
                                                                    cTGGGCAAAAAAGCAAGGGAAGGTTCGTGCCTTTGGGCATAAAGCCGGGGCAAAATCCGTCCGCCGGGCTGTCTCTTTTGCGGCCAACAACgg                     <  2:473976/93‑1 (MQ=255)
                                                                    cTGGGCAAAAAAGCAAGGGAAGGTTCGTGCCTTTGGGCATAAAGCCGGGGCAAAATCCGTCCGCCGGGCTGTCTCTTTTGCGGCCAACAACgg                     <  2:4516181/93‑1 (MQ=255)
                                                                     tGGGCAAAAAAGCAAGGGAAGGTTCGTGCCTTTGGGCATAAAGCCGGGGCAAAATCCGTCCGCCGGGCTGTCTCTTTTGCGGCCAACAACGGt                    >  1:2646616/1‑93 (MQ=255)
                                                                     tGGGCAAAAAAGCAAGGGAAGGTTCGTGCCTTTGGGCATAAAGCCGGGGCAAAATCCGTCCGCCGGGCTGTCTCTTTTGCGGCCAACAACGGt                    <  2:4341621/93‑1 (MQ=255)
                                                                      gggCAAAAAAGCAAGGGAAGGTTCGTGCCTTTGGGCATAAAGCCGGGGCAAAATCCGTCCGCCGGGCTGTCTCTTTTGCGGCCAACAACGGTa                   <  2:3399380/93‑1 (MQ=255)
                                                                      gggCAAAAAAGCAAGGGAAGGTTCGTGCCTTTGGGCATAAAGCCGGGGCAAAATCCGTCCGCCGGGCTGTCTCTTTTGCGGCCAACAACGGTa                   >  1:1039086/1‑93 (MQ=255)
                                                                      gggCAAAAAAGCAAGGGAAGGTTCGTGCCTTTGGGCATAAAGCCGGGGCAAAATCCGTCCGCCGGGCTGTCTCTTTTGCGGCCAACAACGGTa                   >  1:3340090/1‑93 (MQ=255)
                                                                      gggCAAAAAAGCAAGGGAAGGTTCGTGCCTTTGGGCATAAAGCCGGGGCAAAATCCGTCCGCCGGGCTGTCTCTTTTGCGGCCAACAACGGTa                   >  1:4048448/1‑93 (MQ=255)
                                                                       ggCAAAAAAGCAAGGGAAGGTTCGTGCCTTTGGGCATAAAGCCGGGGCAAAATCCGTCCGCCGGGCTGTCTCTTTTGCGGCCAACAACGGTAt                  <  1:2246940/93‑1 (MQ=255)
                                                                          aaaaaaGCAAGGGAAGGTTCGTGCCTTTGGGCATAAAGCCGGGGCAAAATCCGTCCGCCGGGCTGTCTCTTTTGCGGCCAACAACGGTATTGa               <  2:941423/93‑1 (MQ=255)
                                                                          aaaaaaGCAAGGGAAGGTTCGTGCCTTTGGGCATAAAGCCGGGGCAAAATCCGTCCGCCGGGCTGTCTCTTTTGCGGCCAACAACGGTATTGa               <  1:1248045/93‑1 (MQ=255)
                                                                          aaaaaaGCAAGGGAAGGTTCGTGCCTTTGGGCATAAAGCCGGGGCAAAATCCGTCCGCCGGGCTGTCTCTTTTGCGGCCAACAACGGTATTGa               <  1:3403624/93‑1 (MQ=255)
                                                                            aaaaGCAAGGGAAGGTTCGTGCCTTTGGGCATAAAGCCGGGGCAAAATCCGTCCGCCGGGCTGTCTCTTTTGCGGCCAACAACGGTATTGAgg             >  2:2707133/1‑93 (MQ=255)
                                                                            aaaaGCAAGGGAAGGTTCGTGCCTTTGGGCATAAAGCCGGGGCAAAATCCGTCCGCCGGGCTGTCTCTTTTGCGGCCAACAACGGTATTGAgg             <  2:1334474/93‑1 (MQ=255)
                                                                            aaaaGCAAGGGAAGGTTCGTGCCTTTGGGCATAAAGCCGGGGCAAAATCCGTCCGCCGGGCTGTCTCTTTTGCGGCCAACAACGGTATTGAgg             <  2:4450507/93‑1 (MQ=255)
                                                                             aaaGCAAGGGAAGGTTCGTGCCTTTGGGCATAAAGCCGGGGCAAAATCCGTCCGCCGGGCTGTCTCTTTTGCGGCCAACAACGGTATTGAGGc            >  2:683333/1‑93 (MQ=255)
                                                                             aaaGCAAGGGAAGGTTCGTGCCTTTGGGCATAAAGCCGGGGCAAAATCCGTCCGCCGGGCTGTCTCTTTTGCGGCCAACAACGGTATTGAGGc            >  2:4183576/1‑93 (MQ=255)
                                                                             aaaGCAAGGGAAGGTTCGTGCCTTTGGGCATAAAGCCGGGGCAAAATCCGTCCGCCGGGCTGTCTCTTTTGCGGCCAACAACGGTATTGAGGc            >  1:775866/1‑93 (MQ=255)
                                                                                gCAAGGGAAGGTTCGTGCCTTTGGGCATAAAGCCGGGGCAAAATCCGTCCGCCGGGCTGTCTCTTTTGCGGCCAACAACGGTATTGAGGCGtt         >  2:4075941/1‑93 (MQ=255)
                                                                                  aaGGGAAGGTTCGTGCCTTTGGGCATAAAGCCGGGGCAAAATCCGTCCGCCGGGCTGTCTCTTTTGCGGCCAACAACGGTATTGAGGCGTTaa       >  1:1027249/1‑93 (MQ=255)
                                                                                  aaGGGAAGGTTCGTGCCTTTGGGCATAAAGCCGGGGCAAAATCCGTCCGCCGGGCTGTCTCTTTTGCGGCCAACAACGGTATTGAGGCGTTaa       >  1:1154820/1‑93 (MQ=255)
                                                                                  aaGGGAAGGTTCGTGCCTTTGGGCATAAAGCCGGGGCAAAATCCGTCCGCCGGGCTGTCTCTTTTGCGGCCAACAACGGTATTGAGGCGTTaa       >  2:1820409/1‑93 (MQ=255)
                                                                                  aaGGGAAGGTTCGTGCCTTTGGGCATAAAGCCGGGGCAAAATCCGTCCGCCGGGCTGTCTCTTTTGCGGCCAACAACGGTATTGAGGCGTTaa       >  1:1709265/1‑93 (MQ=255)
                                                                                   aGGGAAGGTTCGTGCCTTTGGGCATAAAGCCGGGGCAAAATCCGTCCGCCGGGCTGTCTCTTTTGCGGCCAACAACGGTATTGAGGCGTTAAc      <  1:2160058/93‑1 (MQ=255)
                                                                                   aGGGAAGGTTCGTACATTTGGAAATAAGGCCGGGGAAAAGTCCGTCCGCCGGGCTGTCTCTTTTGCGGCCAACAACGGTATTGAGGCGTTAAc      <  2:203076/93‑1 (MQ=255)
                                                                                    gggAAGGTTCGTGCCTTTGGGCATAAAGCCGGGGCAAAATCCGTCCGCCGGGCTGTCTCTTTTGCGGCTAACAACGGTATTGAGGCGTTAACg     <  2:2705253/93‑1 (MQ=255)
                                                                                    gggAAGATTCGTGCCTTTGGGCATAAAGCCGGGGCAAAATCCGTCCGCCGGGCTGTCTCTTTTGCGGCCAACAACGGTATTGAGGCGTTAACg     >  2:2346665/1‑93 (MQ=255)
                                                                                      gAAGGTTCGTGCCTTTGGGCATAAAGCCGGGGCAAAATCCGTCCGCCGGGCTGTCTCTTTTGCGGCCAACAACGGTATTGAGGCGTTAACGCt   <  1:3109823/93‑1 (MQ=255)
                                                                                       aaGGTTCGTGCCTTGGGGCATAAAGCCGGGGCAAAATCCGTCCGCCGGGCTGTCTCTTTTGCGGCCAACAACGGTATTGAGGCGTTAACGCTg  <  1:3390161/93‑1 (MQ=255)
                                                                                          |                                                                                         
CCACTTAGCGAAAAATTGCCAGCGCATGGCTGCCGTCATGTTGCGATCATTATGGACGGCAATGGCCGCTGGGCAAAAAAGCAAGGGAAGATTCGTGCCTTTGGGCATAAAGCCGGGGCAAAATCCGTCCGCCGGGCTGTCTCTTTTGCGGCCAACAACGGTATTGAGGCGTTAACGCTG  >  NC_000913/194924‑195103

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