Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 675,063 A→G A236A (GCT→GCC ybeQ ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913675,0630AG99.6% 712.6 / NA 221A236A (GCT→GCCybeQhypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (0/1);  new base G (112/108);  total (112/109)
Fisher's exact test for biased strand distribution p-value = 4.93e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 232 total aligned reads displayed.

CCATTCAAAGGCAGCCTGATAATCTTTCTCAACGCCTTGTCCGTGTTTATACATATAACCGATGTTGACGTAAGCGCCGACGCTCCCCTGAGCGGCAGCCTTCAGATACCAGTACATCGCTTGCTTATAGTCCTGGTCCACCCCTTCGCCAGCGTTATAATCCCATGCAATTTGAAACTG  >  NC_000913/674973‑675152
                                                                                          |                                                                                         
  aTTCAAAGGCAGCCTGATAATCTTTCTCAACGCCTTGTCCGTGTTTATACATATAACCGATGTTGACGTAAGCGCCGACGCTCCCCTGGGCgg                                                                                       >  2:503694/1‑93 (MQ=255)
  aTTCAAAGGCAGCCTGATAATCTTTCTCAACGCCTTGTCCGTGTTTATACATATAACCGATGTTGACGTAAGCGCCGACGCTCCCCTGGGCgg                                                                                       <  1:2186564/93‑1 (MQ=255)
   ttCAAAGGCAGCCTGATAATCTTTCTCAACGCCTTGTCCGTGTTTATACATATAACCGATGTTGACGTAAGCGCCGACGCTCCCCTGggcggc                                                                                      >  2:6358360/1‑93 (MQ=255)
   ttCAAAGGCAGCCTGATAATCTTTCTCAACGCCTTGTCCGTGTTTATACATATAACCGATGTTGACGTAAGCGCCGACGCTCCCCTGggcggc                                                                                      >  2:5893345/1‑93 (MQ=255)
      aaaGGCAGCCTGATAATCTTTCTCAACGCCTTGTCCGTGTTTATACATATAACCGATGTTGACGTAAACGCCGACGCTCCCCTGGGCCGCAGc                                                                                   >  1:820678/1‑93 (MQ=255)
       aaGGCAGCCTGATAATCTTTCTCAACGCCTTGTCCGTGTTTATACATATAACCGATGTTGACGTAAGCGCCGACGCTCCCCTGGGCGGCAGcc                                                                                  >  2:2202388/1‑93 (MQ=255)
       aaGGCAGCCTGATAATCTTTCTCAACGCCTTGTCCGTGTTTATACATATAACCGATGTTGACGTAAGCGCCGACGCTCCCCTGGGCGGCAGcc                                                                                  >  1:2449108/1‑93 (MQ=255)
        aGGCAGCCTGATAATCTTTCTCAACGCCTTGTCCGTGTTTATACATATAACCGATGTTGACGTAAGCGCCGACGCTCCCCTGGGCGGCAGCCt                                                                                 >  1:33676/1‑93 (MQ=255)
         ggCAGCCTGATAATCTTTCTCAACGCCTTGTCCGTGTTTATACATATAACCGATGTTGACGTAAGCGCCGACGCTCCCCTGGGCGGCAGCCtt                                                                                >  2:835107/1‑93 (MQ=255)
           cAGCCTGATAATCTTTCTCAACGCCTTGTCCGTGTTTATACATATAACCGATGTTGACGTAAGCGCCGACGCTCCCCTGGGCGGCAGCCTTCa                                                                              <  1:1259222/93‑1 (MQ=255)
           cAGCCTGATAATCTTTCTCAACGCCTTGTCCGTGTTTATACATATAACCGATGTTGACGTAAGCGCCGACGCTCCCCTGGGCGGCAGCCTTCa                                                                              <  1:1343851/93‑1 (MQ=255)
           cAGCCTGATAATCTTTCTCAACGCCTTGTCCGTGTTTATACATATAACCGATGTTGACGTAAGCGCCGACGCTCCCCTGGGCGGCAGCCTTCa                                                                              <  2:2538579/93‑1 (MQ=255)
            aGCCTGATAATCTTTCTCAACGCCTTGTCCGTGTTTATACATATAACCGATGTTGACGTAAGCGCCGACGCTCCCCTGGGCGGCAGCCTTCAg                                                                             <  2:3829589/93‑1 (MQ=255)
            aGCCTGATAATCTTTCTCAACGCCTTGTCCGTGTTTATACATATAACCGATGTTGACGTAAGCGCCGACGCTCCCCTGGGCGGCAGCCTTCAg                                                                             <  1:6725689/93‑1 (MQ=255)
            aGCCTGATAATCTTTCTCAACGCCTTGTCCGTGTTTATACATATAACCGATGTTGACGTAAGCGCCGACGCTCCCCTGGGCGGCAGCCTTCAg                                                                             <  1:4060409/93‑1 (MQ=255)
             gCCTGATAATCTTTCTCAACGCCTTGTCCGTGTTTATACATATAACCGATGTTGACGTAAGCGCCGACGCTCCCCTGGGCGGCAGCCTTCAGa                                                                            <  2:5302989/93‑1 (MQ=255)
               cTGATAATCTTTCTCAACGCCTTGTCCGTGTTTATACATATAACCGATGTTGACGTAAGCGCCGACGCTCCCCTGGGCGGCAGCCTTCAGATa                                                                          >  1:314125/1‑93 (MQ=255)
               cTGATAATCTTTCTCAACGCCCTGTCCGTGTTTATACATATAACCGATGTTGACGTAAGCGCCGACGCTCCCCTGGGCGGCAGCCTTCAGATa                                                                          >  2:3243062/1‑93 (MQ=255)
                tGATAATCTTTCTCAACGCCTTGTCCGTGTTTATACATATAACCGATGTTGACGTAAGCGCCGACGCTCCCCTGGGCGGCAGCCTTCAGATAc                                                                         >  2:2325181/1‑93 (MQ=255)
                 gATAATCTTTCTCAACGCCTTGTCCGTGTTTATACATATAACCGATGTTGACGTAAGCGCCGACGCTCCCCTGGGCGGCAGCCTTCAGATAcc                                                                        >  1:3836951/1‑93 (MQ=255)
                 gATAATCTTTCTCAACGCCTTGTCCGTGTTTATACATATAACCGATGTTGACGTAAGCGCCGACGCTCCCCTGGGCGGCAGCCTTCAGATAcc                                                                        >  2:843181/1‑93 (MQ=255)
                 gATAATCTTTCTCAACGCCTTGTCCGTGTTTATACATATAACCGATGTTGACGTAAGCGCCGACGCTCCCCTGGGCGGCAGCCTTCAGATAcc                                                                        >  2:1095105/1‑93 (MQ=255)
                  aTAATCTTTCTCAACGCCTTGTCCGTGTTTATACATATAACCGATGTTGACGTAAGCGCCGACGCTCCCCTGGGCGGCAGCCTTCAGATACCa                                                                       <  1:3423352/93‑1 (MQ=255)
                  aTAATCTTTCTCAACGCCTTGTCCGTGTTTATACATATAACCGATGTTGACGTAAGCGCCGACGCTCCCCTGGGCGGCAGCCTTCAGATACCa                                                                       >  2:1374622/1‑93 (MQ=255)
                   tAATCTTTCTCAACGCCTTGTCCGTGTTTATACATATAACCGATGTTGACGTAAGCGCCGACGCTCCCCTGGGCGGCAGCCTTCAGATACCAg                                                                      <  2:39866/93‑1 (MQ=255)
                      tCTTTCTCAACGCCTTGTCCGTGTTTATACATATAACCGATGTTGACGTAAGCGCCGACGCTCCCCTGGGCGGCAGCCTTCAGATACCAGTAc                                                                   <  1:5190596/93‑1 (MQ=255)
                         ttCTCAACGCCTTGTCCGTGTTTATACATATAACCGATGTTGACGTAAGCGCCGACGCTCCCCTGGGCGGCAGCCTTCAGATACCAGTACATc                                                                <  2:2029226/93‑1 (MQ=255)
                         ttCTCAACGCCTTGTCCGTGTTTATACATATAACCGATGTTGACGTAAGCGCCGACGCTCCCCCGGGCGGCAGCCTTCAGATACCAGTACATc                                                                >  1:5690281/1‑93 (MQ=255)
                            tcAACGCCTTGTCCGTGTTTATACATATAACCGATGTTGACGTAAGCGCCGACGCTCCCCTGGGCGGCAGCCTTCAGATACCAGTACATCgct                                                             >  1:4172130/1‑93 (MQ=255)
                             cAACGCCTTGTCCGTGTTTATACATATAACCGATGTTGACGTAAGCGCCGACGCTCCCCTGGGCGGCAGCCTTCAGATACCAGTACATCgctt                                                            <  1:4583778/93‑1 (MQ=255)
                             cAACGCCTTGTCCGTGTTTATACATATAACCGATGTTGACGTAAGCGCCGACGCTCCCCTGGGCGGCAGCCTTCAGATACCAGTACATCgctt                                                            >  2:479802/1‑93 (MQ=255)
                             cAACGCCTTGTCCGTGTTTATACATATAACCGATGTTGACGTAAGCGCCGACGCTCCCCTGGGCGGCAGCCTTCAGATACCAGTACATCgctt                                                            >  2:6693234/1‑93 (MQ=255)
                               aCGCCTTGTCCGTGTTTATACATATAACCGATGTTGACGTAAGCGCCGACGCTCCCCTGGGCGGCAGCCTTCAGATACCAGTACATCgcttgc                                                          >  1:5331639/1‑93 (MQ=255)
                                 gCCTTGTCCGTGTTTATACATATAACCGATGTTGACGTAAGCGCCGACGCTCCCCTGGGCGGCAGCCTTCAGATACCAGTACATCgcttgctt                                                        <  2:241459/93‑1 (MQ=255)
                                  ccTTGTCCGTGTTTATACATATAACCGATGTTGACGTAAGCGCCGACGCTCCCCTGGGCGGCAGCCTTCAGATACCAGTACATCGCTTGCTta                                                       >  1:4821509/1‑93 (MQ=255)
                                  ccTTGTCCGTGTTTATACATATAACCGATGTTGACGTAAGCGCCGACGCTCCCCTGGGCGGCAGCCTTCAGATACCAGTACATCGCTTGCTta                                                       <  1:169745/93‑1 (MQ=255)
                                  ccTTGTCCGTGTTTATACATATAACCGATGTTGACGTAAGCGCCGACGCTCCCCTGGGCGGCAGCCTTCAGATACCAGTACATCGCTTGCTta                                                       <  1:2484075/93‑1 (MQ=255)
                                   cTTGTCCGTGTTTATACATATAACCGATGTTGACGTAAGCGCCGACGCTCCCCTGGGCGGCAGCCTTCAGATACCAGTACATCGCTTGCTtat                                                      <  2:4315036/93‑1 (MQ=255)
                                   cTTGTCCGTGTTTATACATATAACCGATGTTGACGTAAGCGCCGACGCTCCCCTGGGCGGCAGCCTTCAGATACCAGTACATCGCTTGCTtat                                                      <  1:3298097/93‑1 (MQ=255)
                                    ttGTCCGTGTTTATACATATAACCGATGTTGACGTAAGCGCCGACGCTCCCCTGGGCGGCAGCCTTCAGATACCAGTACATCGCTTGCTtata                                                     >  1:1970524/1‑93 (MQ=255)
                                    ttGTCCGTGTTTATACATATAACCGATGTTGACGTAAGCGCCGACGCTCCCCTGGGCGGCAGCCTTCAGATACCAGTACATCGCTTGCTtata                                                     <  2:3259286/93‑1 (MQ=255)
                                    ttGTCCGTGTTTATACATATAACCGATGTTGACGTAAGCGCCGACGCTCCCCTGGGCGGCAGCCTTCAGATACCAGTACATCGCTTGCTtata                                                     <  2:4097672/93‑1 (MQ=255)
                                    ttGTCCGTGTTTATACATATAACCGATGTTGACGTAAGCGCCGACGCTCCCCTGGGCGGCAGCCTTCAGATACCAGTACATCGCTTGCTtata                                                     <  2:5507305/93‑1 (MQ=255)
                                         cGTGTTTATACATATAACCGATGTTGACGTAAGCGCCGACGCTCCCCTGGGCGGCAGCCTTCAGATACCAGTACATCGCTTGCTTATAGTCCt                                                <  1:4223716/93‑1 (MQ=255)
                                          gtgtTTATACATATAACCGATGTTGACGTAAGCGCCGACGCTCCCCTGGGCGGCAGCCTTCAGATACCAGTACATCGCTTGCTTATAGTCCTg                                               <  1:470331/93‑1 (MQ=255)
                                           tgtTTATACATATAACCGATGTTGACGTAAGCGCCGACGCTCCCCTGGGCGGCAGCCTTCAGATACCAGTACATCGCTTGCTTATAGTCCTgg                                              >  1:2832159/1‑93 (MQ=255)
                                           tgtTTATACATATAACCGATGTTGACGTAAGCGCCGACGCTCCCCTGGGCGGCAGCCTTCAGATACCAGTACATCGCTTGCTTATAGTCCTgg                                              >  1:2723543/1‑93 (MQ=255)
                                              ttATACATATAACCGATGTTGACGTAAGCGCCGACGCTCCCCTGGGCGGCAGCCTTCAGATACCAGTACATCGCTTGCTTATAGTCCTGGTcc                                           >  1:3739944/1‑93 (MQ=255)
                                               tataCATATAACCGATGTTGACGTAAGCGCCGACGCTCCCCTGGGCGGCAGCCTTCAGATACCAGTACATCGCTTGCTTATAGTCCTGGTCCa                                          >  1:2998785/1‑93 (MQ=255)
                                                ataCATATAACCGATGTTGACGTAAGCGCCGACGCTCCCCTGGGCGGCAGCCTTCAGATACCAGTACATCGCTTGCTTATAGTCCTGGTCCAc                                         >  1:2693180/1‑93 (MQ=255)
                                                ataCATATAACCGATGTTGACGTAAGCGCCGACGCTCCCCTGGGCGGCAGCCTTCAGATACCAGTACATCGCTTGCTTATAGTCCTGGTCCAc                                         >  2:1515378/1‑93 (MQ=255)
                                                 taCATATAACCGATGTTGACGTAAGCGCCGACGCCCCCCTGGGCGGCAGCCTTCAGATACCAGTACATCGCTTGCTTATAGTCCTGGTCCAcc                                        >  2:5487125/1‑93 (MQ=255)
                                                   cATATAACCGATGTTGACGTAAGCGCCGACGCTCCCCTGGGCGGCAGCCTTCAGATACCAGTACATCGCTTGCTTATAGTCCTGGTCCAcccc                                      <  2:4838072/93‑1 (MQ=255)
                                                    atatAACCGATGTTGACGTAAGCGCCGACGCTCCCCTGGGCGGCAGCCTTCAGATACCAGTACATCGCTTGCTTATAGTCCTGGTCCACCCCt                                     >  1:1097836/1‑93 (MQ=255)
                                                     tataACCGATGTTGACGTAAGCGCCGACGCTCCCCTGGGCGGCAGCCTTCAGATACCAGTACATCGCTTGCTTATAGTCCTGGTCCACCCCtt                                    <  1:3208952/93‑1 (MQ=255)
                                                     tataACCGATGTTGACGTAAGCGCCGACGCTCCCCTGGGCGGCAGCCTTCAGATACCAGTACATCGCTTGCTTATAGTCCTGGTCCACCCCtt                                    <  2:1803572/93‑1 (MQ=255)
                                                      ataACCGATGTTGACGTAAGCGCCGACGCTCCCCTGGGCGGCAGCCTTCAGATACCAGTACATCGCTTGCTTATAGTCCTGGTCCACCCCTTc                                   <  2:1768280/93‑1 (MQ=255)
                                                       taACCGATGTTGACGTAAGCGCCGACGCTCCCCTGGGCGGCAGCCTTCAGATACCAGTACATCGCTTGCTTATAGTCCTGGTCCACCCCTTCg                                  >  1:6229466/1‑93 (MQ=255)
                                                       taACCGATGTTGACGTAAGCGCCGACGCTCCCCTGGGCGGCAGCCTTCAGATACCAGTACATCGCTTGCTTATAGTCCTGGTCCACCCCTTCg                                  >  1:3499588/1‑93 (MQ=255)
                                                        aaCCGATGTTGACGTAAGCGCCGACGCTCCCCTGGGCGGCAGCCTTCAGATACCAGTACATCGCTTGCTTATAGTCCTGGTCCACCCCTTCGc                                 >  2:5864938/1‑93 (MQ=255)
                                                        aaCCGATGTTGACGTAAGCGCCGACGCTCCCCTGGGCGGCAGCCTTCAGATACCAGTACATCGCTTGCTTATAGTCCTGGTCCACCCCTTCGc                                 >  2:2675729/1‑93 (MQ=255)
                                                         aCCGATGTTGACGTAAGCGCCGACGCTCCCCTGGGCGGCAGCCTTCAGATACCAGTACATCGCTTGCTTATAGTCCTGGTCCACCCCTTCGcc                                >  2:2760182/1‑93 (MQ=255)
                                                         aCCGATGTTGACGTAAGCGCCGACGCTCCCCTGGGCGGCAGCCTTCAGATACCAGTACATCGCTTGCTTATAGTCCTGGTCCACCCCTTCGcc                                >  1:1739841/1‑93 (MQ=255)
                                                           cGATGTTGACGTAAGCGCCGACGCTCCCCTGGGCGGCAGCCTTCAGATACCAGTACATCGCTTGCTTATAGTCCTGGTCCACCCCTTCGCCAg                              >  1:3001060/1‑93 (MQ=255)
                                                           cGATGTTGACGTAAGCGCCGACGCTCCCCTGGGCGGCAGCCTTCAGATACCAGTACATCGCTTGCTTATAGTCCTGGTCCACCCCTTCGCCAg                              >  1:3678173/1‑93 (MQ=255)
                                                            gATGTTGACGTAAGCGCCGACGCTCCCCTGGGCGGCAGCCTTCAGATACCAGTACATCGCTTGCTTATAGTCCTGGTCCACCCCTTCGCCAGc                             >  1:1001468/1‑93 (MQ=255)
                                                            gATGTTGACGTAAGCGCCGACGCTCCCCTGGGCGGCAGCCTTCAGATACCAGTACATCGCTTGCTTATAGTCCTGGTCCACCCCTTCGCCAGc                             >  1:5248511/1‑93 (MQ=255)
                                                            gATGTTGACGTAAGCGCCGACGCTCCCCTGGGCGGCAGCCTTCAGATAACAGTACATCGCTTGCTTATAGTCCTGGTCCACCCCTTCGCCAGc                             >  1:5029114/1‑93 (MQ=255)
                                                              tGTTGACGTAAGCGCCGACGCTCCCCTGGGCGGCAGCCTTCAGATACCAGTACATCGCTTGCTTATAGTCCTGGTCCACCCCTTCGCCAGCGt                           >  1:3680865/1‑93 (MQ=255)
                                                              tGTTGACGTAAGCGCCGACGCTCCCCTGGGCGGCAGCCTTCAGATACCAGTACATCGCTTGCTTATAGTCCTGGTCCACCCCTTCGCCAGCGt                           <  1:389989/93‑1 (MQ=255)
                                                                  gACGTAAGCGCCGACGCTCCCCTGGGCGGCAGCCTTCAGATACCAGTACATCGCTTGCTTATAGTCCTGGTCCACCCCTTCGCCAGCGTtata                       <  2:3607095/93‑1 (MQ=255)
                                                                  gACGTAAGCGCCGACGCTCCCCTGGGCGGCAGCCTTCAGATACCAGTACATCGCTTGCTTATAGTCCTGGTCCACCCCTTCGCCAGCGTtata                       <  1:2546233/93‑1 (MQ=255)
                                                                  gACGTAAGCGCCGACGCTCCCCTGGGCGGCAGCCTTCAGATACCAGTACATCGCTTGCTTATAGTCCTGGTCCACCCCTTCGCCAGCGTtata                       >  2:510038/1‑93 (MQ=255)
                                                                  gACGTAAGCGCCGACGCTCCCCTGGGCGGCAGCCTTCAGATACCAGTACATCGCTTGCTTATAGTCCTGGTCCACCCCTTCGCCAGCGTtata                       <  2:6290603/93‑1 (MQ=255)
                                                                   aCGTAAGCGCCGACGCTCCCCTGGGCGGCAGCCTTCAGATACCAGTACATCGCTTGCTTATAGTCCTGGTCCACCCCTTCGCCAGCGTTATaa                      >  2:4328243/1‑93 (MQ=255)
                                                                    cGTAAGCGCCGACGCTCCCCTGGGCGGCAGCCTTCAGATACCAGTACATCGCTTGCTTATAGTCCTGGTCCACCCCTTCGCCAGCGTTATAAt                     >  2:5172603/1‑93 (MQ=255)
                                                                    agTAAGCGCCGACGCTCCCCTGGGCGGCAGCCTTCAGATACCAGTACATCGCTTGCTTATAGTCCTGGTCCACCCCTTCGCCAGCGTTATAAt                     <  2:3942778/92‑1 (MQ=255)
                                                                     gTAAGCGCCGACGCTCCCCTGGGCGGCAGCCTTCAGATACCAGTACATCGCTTGCTTATAGTCCTGGTCCACCCCTTCGCCAGCGTTATAATc                    <  1:2109837/93‑1 (MQ=255)
                                                                      tAAGCGCCGACGCTCCCCTGGGCGGCAGCCTTCAGATACCAGTACATCGCTTGCTTATAGTCCTGGTCCACCCCTTCGCCAGCGTTATAATcc                   >  1:1037846/1‑93 (MQ=255)
                                                                      tAAGCGCCGACGCTCCCCTGGGCGGCAGCCTTCAGATACCAGTACATCGCTTGCTTATAGTCCTGGTCCACCCCTTCGCCAGCGTTATAATcc                   <  1:1217125/93‑1 (MQ=255)
                                                                      tAAGCGCCGACGCTCCCCTGGGCGGCAGCCTTCAGATACCAGTACATCGCTTGCTTATAGTCCTGGTCCACCCCTTCGCCAGCGTTATAATcc                   >  1:1564467/1‑93 (MQ=255)
                                                                       aaGCGCCGACGCTCCCCTGGGCGGCAGCCTTCAGATACCAGTACATCGCTTGCTTATAGTCCTGGTCCACCCCTTCGCCAGCGTTATAATccc                  <  2:1453132/93‑1 (MQ=255)
                                                                       aaGCGCCGACGCTCCCCTGGGCGGCAGCCTTCAGATACCAGTACATCGCTTGCTTATAGTCCTGGTCCACCCCTTCGCCAGCGTTATAATccc                  >  1:4258924/1‑93 (MQ=255)
                                                                        aGCGCCGACGCTCCCCTGGGCGGCAGCCTTCAGATACCAGTACATCGCTTGCTTATAGTCCTGGTCCACCCCTTCGCCAGCGTTATAATCCCa                 <  1:3265056/93‑1 (MQ=255)
                                                                        aGCGCCGACGCTCCCCTGGGCGGCAGCCTTCAGATACCAGTACATCGCTTGCTTATAGTCCTGGTCCACCCCTTCGCCAGCGTTATAATCCCa                 >  2:2115148/1‑93 (MQ=255)
                                                                         gcgcCGACGCTCCCCTGGGCGGCAGCCTTCAGATACCAGTACATCGCTTGCTTATAGTCCTGGTCCACCCCTTCGCCAGCGTTATAATCCCAt                <  2:3367680/93‑1 (MQ=255)
                                                                         gcgcCGACGCTCCCCTGGGCGGCAGCCTTCAGATACCAGTACATCGCTTGCTTATAGTCCTGGTCCACCCCTTCGCCAGCGTTATAATCCCAt                <  1:2238328/93‑1 (MQ=255)
                                                                          cgcCGACGCTCCCCTGGGCGGCAGCCTTCAGATACCAGTACATCGCTTGCTTATAGTCCTGGTCCACCCCTTCGCCAGCGTTATAATCCCATg               <  2:3322533/93‑1 (MQ=255)
                                                                          cgcCGACGCTCCCCTGGGCGGCAGCCTTCAGATACCAGTACATCGCTTGCTTATAGTCCTGGTCCACCCCTTCGCCAGCGTTATAATCCCATg               <  1:5945306/93‑1 (MQ=255)
                                                                           gcCGACGCTCCCCTGGGCGGCAGCCTTCAGATACCAGTACATCGCTTGCTTATAGTCCTGGTCCACCCCTTCGCCAGCGTTATAATCCCATGc              >  1:4931877/1‑93 (MQ=255)
                                                                              gACGCTCCCCTGGGCGGCAGCCTTCAGATACCAGTACATCGCTTGCTTATAGTCCTGGTCCACCCCTTCGCCAGCGTTATAATCCCATGCAAt           <  1:6654243/93‑1 (MQ=255)
                                                                               aCGCTCCCCTGGGCGGCAGCCTTCAGATACCAGTACATCGCTTGCTTATAGTCCTGGTCCACCCCTTCGCCAGCGTTATAATCCCATGCAAtt          >  2:3445215/1‑93 (MQ=255)
                                                                                cGCTCCCCTGGGCGGCAGCCTTCAGATACCAGTACATCGCTTGCTTATAGTCCTGGTCCACCCCTTCGCCAGCGTTATAATCCCATGCAAttt         >  2:6229503/1‑93 (MQ=255)
                                                                                 gCTCCCCTGGGCGGCAGCCTTCAGATACCAGTACATCGCTTGCTTATAGTCCTGGTCCACCCCTTCGCCAGCGTTATAATCCCATGCAATTTg        <  2:5329778/93‑1 (MQ=255)
                                                                                 gCTCCCCTGGGCGGCAGCCTTCAGATACCAGTACATCGCTTGCTTATAGTCCTGGTCCACCCCTTCGCCAGCGTTATAATCCCATGCAATTTg        <  2:2995560/93‑1 (MQ=255)
                                                                                    ccccTGGGCGGCAGCCTTCAGATACCAGTACATCGCTTGCTTATAGTCCTGGTCCACCCCTTCGCCAGCGTTATAATCCCATGCAATTTGaaa     >  2:4297360/1‑93 (MQ=255)
                                                                                     cccTGGGCGGCAGCCTTCAGATACCAGTACATCGCTTGCTTATAGTCCTGGTCCACCCCTTCGCCAGCGTTATAATCCCATGCAATTTGAAAc    <  2:4392189/93‑1 (MQ=255)
                                                                                      ccTGGGCGGCAGCCTTCAGATACCAGTACATCGCTTGCTTATAGTCCTGGTCCACCCCTTCGCCAGCGTTATAATCCCATGCAATTTGAAACt   <  1:5889534/93‑1 (MQ=255)
                                                                                       cTGGGCGGCAGCCTTCAGATACCAGTACATCGCTTGCTTATAGTCCTGGTCCACCCCTTCGCCAGCGTTATAATCCCATGCAATTTGAAACTg  >  1:3513100/1‑93 (MQ=255)
                                                                                       cTGGGCGGCAGCCTTCAGATACCAGTACATCGCTTGCTTATAGTCCTGGTCCACCCCTTCGCCAGCGTTATAATCCCATGCAATTTGAAACTg  <  1:5208238/93‑1 (MQ=255)
                                                                                          |                                                                                         
CCATTCAAAGGCAGCCTGATAATCTTTCTCAACGCCTTGTCCGTGTTTATACATATAACCGATGTTGACGTAAGCGCCGACGCTCCCCTGAGCGGCAGCCTTCAGATACCAGTACATCGCTTGCTTATAGTCCTGGTCCACCCCTTCGCCAGCGTTATAATCCCATGCAATTTGAAACTG  >  NC_000913/674973‑675152

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: