Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,447,463 T→C N4S (AAC→AGC)  aroC ← chorismate synthase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,447,4630TC100.0% 1481.2 / NA 355N4S (AAC→AGC) aroCchorismate synthase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (125/230);  total (125/230)

Only 100 of 361 total aligned reads displayed.

TCAGCGGAATGCCTGGCGGAACACCATCGACGATGCAGCCGAGCGCCAGCCCGTGCGATTCGCCGAAGGTGGTTACGCGAAAGAGTTGTCCAATTGTGTTTCCAGCCATCACGGCTCCGTTATTGTTGTGTTTGCGTGTTTACTTAATCTTTATAAATCGCGAAATGTTCTCGTGCGGCAATAAGCTGCTCTTTGG  >  NC_000913/2447365‑2447560
                                                                                                  |                                                                                                 
 cAGCGGAATGCCTGGCGGAACACCATCGACGATGCAGCCGAGCGCCAGCCCGTGCGATTCGCCGAAGGTGGTTACGCGAAAGAGTTGTCCAATTGTGCTTc                                                                                                <  1:3984615/101‑1 (MQ=255)
  aGCGGAATGCCTGGCGGAACACCATCGACGATGCAGCCGAGCGCCAGCCCGTGCGATTCGCCGAAGGTGGTTACGCGAAAGAGTTGTCCAATTGTGCTTcc                                                                                               >  1:3985213/1‑101 (MQ=255)
   gCGGAATGCCTGGCGGAACACCATCGACGATGCAGCCGAGCGCCAGCCCGTGCGATTCGCCGAAGGTGGTTACGCGAAAGAGTTGTCCAATTGTGCTTCCa                                                                                              <  2:3984681/101‑1 (MQ=255)
    cGGAATGCCTGGCGGAACACCATCGACGATGCAGCCGAGCGCCAGCCCGTGCGATTCGCCGAAGGTGGTTACGCGAAAGAGTTGTCCAATTGTGCTTcc                                                                                               >  2:3985216/1‑99 (MQ=255)
    cGGAATGCCTGGCGGAACACCATCGACGATGCAGCCGAGCGCCAGCCCGTGCGATTCGCCGAAGGTGGTTACGCGAAAGAGTTGTCCAATTGTGCTTCCAg                                                                                             <  1:3984691/101‑1 (MQ=255)
         tGCCTGGCGGAACACCATCGACGATGCAGCCGAGCGCCAGCCCGTGCGATTCGCCGAAGGTGGTTACGCGAAAGAGTTGTCCAATTGTGCTTCCAGCCATc                                                                                        >  1:3985221/1‑101 (MQ=255)
           ccTGGCGGAACACCATCGACGATGCAGCCGAGCGCCAGCCCGTGCGATTCGCCGAAGGTGGTTACGCGAAAGAGTTGTCCAATTGTGCTTCCAGcc                                                                                           <  1:3984728/96‑1 (MQ=255)
           ccTGGCGGAACACCATCGACGATGCAGCCGAGCGCCAGCCCGTGCGATTCGCCGAAGGTGGTTACGCGAAAGAGTTGTCCAATTGTGCTTCCAGCCATCAc                                                                                      <  2:3984956/101‑1 (MQ=255)
           ccTGGCGGAACACCATCGACGATGCAGCCGAGCGCCAGCCCGTGCGATTCGCCGAAGGTGGTTACGCGAAAGAGTTGTCCAATTGTGCTTCCAGCCATCAc                                                                                      >  2:3985223/1‑101 (MQ=255)
            cTGGCGGAACACCATCGACGATGCAGCCGAGCGCCAGCCCGTGCGATTCGCCGAAGGTGGTTACGCGAAAGAGTTGTCCAATTGTGCTTCCAGCCATCACg                                                                                     <  2:3985187/101‑1 (MQ=255)
              ggCGGAACACCATCGACGATGCAGCCGAGCGCCAGCCCGTGCGATTCGCCGAAGGTGGTTACGCGAAAGAGTTGTCCAATTGTGCTTCCAGCCATCACGGc                                                                                   <  2:3984653/101‑1 (MQ=255)
                 ggAACACCATCGACGATGCAGCCGAGCGCCAGCCCGTGCGATTCGCCGAAGGTGGTTACGCGAAAGAGTTGTCCAATTGTGCTTCCAGCCATCACGGCTcc                                                                                >  2:3985228/1‑101 (MQ=255)
                    acacCATCGACGATGCAGCCGAGCGCCAGCCCGTGCGATTCGCCGAAGGTGGTTACGCGAAAGAGTTGTCCAATTGTGCTTCCAGCCATCACGGCTCCGtt                                                                             >  1:3985231/1‑101 (MQ=255)
                    acacCATCGACGATGCAGCCGAGCGCCAGCCCGTGCGATTCGCCGAAGGTGGTTACGCGAAAGAGTTGTCCAATTGTGCTTCCAGCCATCACGGCTCCGtt                                                                             <  1:3984766/101‑1 (MQ=255)
                      acCATCGACGATGCAGCCGAGCGCCAGCCCGTGCGATTCGCCGAAGGTGGTTACGCGAAAGAGTTGTCCAATTGTGCTTCCAGCCATCACGGCTCCGTTAt                                                                           >  2:3985237/1‑101 (MQ=255)
                       ccATCGACGATGCAGCCGAGCGCCAGCCCGTGCGATTCGCCGAAGGTGGTTACGCGAAAGAGTTGTCCAATTGTGCTTCCAGCCATCACGGCTCCGTTAtt                                                                          <  1:3984449/101‑1 (MQ=255)
                         aTCGACGATGCAGCCGAGCGCCAGCCCGTGCGATTCGCCGAAGGTGGTTACGCGAAAGAGTTGTCCAATTGTGCTTCCAGCCATCAc                                                                                      <  2:3984464/87‑1 (MQ=255)
                          tCGACGATGCAGCCGAGCGCCAGCCCGTGCGATTCGCCGAAGGTGGTTACGCGAAAGAGTTGTCCAATTGTGCTTCCAGCCATCACGGCTCCGTTAttgtt                                                                       >  1:3985241/1‑101 (MQ=255)
                          tCGACGATGCAGCCGAGCGCCAGCCCGTGCGATTCGCCGAAGGTGGTTACGCGAAAGAGTTGTCCAATTGTGCTTCCAGCCATCACGGCTCCGTTAttgtt                                                                       <  2:3984588/101‑1 (MQ=255)
                          tCGACGATGCAGCCGAGCGCCAGCCCGTGCGATTCGCCGAAGGTGGTTACGCGAAAGAGTTGTCCAATTGTGCTTCCAGCCATCACGGCTCCGTTAttgtt                                                                       >  2:3985242/1‑101 (MQ=255)
                           cgacgaTGCAGCCGAGCGCCAGCCCGTGCGATTCGCCGAAGGTGGTTACGCGAAAGAGTTGTCCAATTGTGCTTCCAGCCATCACGGCTCCGTTAttgttg                                                                      <  1:3985095/101‑1 (MQ=255)
                           cgacgaTGCAGCCGAGCGCCAGCCCGTGCGATTCGCCGAAGGTGGTTACGCGAAAGAGTTGTCCAATTGTGCTTCCAGCCATCACGGCTCCGTTAttgttg                                                                      <  1:3984669/101‑1 (MQ=255)
                              cgaTGCAGCCGAGCGCCAGCCCGTGCGATTCGCCGAAGGTGGTTACGCGAAAGAGTTGTCCAATTGTGCTTCCAGCCATCACGGCTCCGTTATTGTTgtgt                                                                   >  2:3985246/1‑101 (MQ=255)
                              cgaTGCAGCCGAGCGCCAGCCCGTGCGATTCGCCGAAGGTGGTTACGCGAAAGAGTTGTCCAATTGTGCTTCCAGCCATCACGGCTCCGTTATTGTTgtgt                                                                   <  1:3985169/101‑1 (MQ=255)
                               gaTGCAGCCGAGCGCCAGCCCGTGCGATTCGCCGAAGGTGGTTACGCGAAAGAGTTGTCCAATTGTGCTTCCAGCCATCACGGCTCCGTTAttgtt                                                                       <  1:3984914/96‑1 (MQ=255)
                                 tGCAGCCGAGCGCCAGCCCGTGCGATTCGCCGAAGGTGGTTACGCGAAAGAGTTGTCCAATTGTGCTTCCAGCCATCACGGCTCCGTTATTGTTGTGTTTg                                                                <  1:3984589/101‑1 (MQ=255)
                                    aGCCGAGCGCCAGCCCGTGCGATTCGCCGAAGGTTGTTACGCGAAAGAGTTGTCCAATTGTGCTTCCAGCCATCACGGCTCCGTTATTGTTGTGttt                                                                 <  1:3984649/97‑1 (MQ=255)
                                    aGCCGAGCGCCAGCCCGTGCGATTCGCCGAAGGTGGTTACGCGAAAGAGTTGTCCAATTGTGCTTCCAGCCATCACGGCTCCGTTATTGTTGTGTTTGCgt                                                             <  1:3984639/101‑1 (MQ=255)
                                     gCCGAGCGCCAGCCCGTGCGATTCGCCGAAGGTGGTTACGCGAAAGAGTTGTCCAATTGTGCTTCCAGCCATCACGGCTCCGTTATTGTTGTGTTTGCgtg                                                            <  2:3984752/101‑1 (MQ=255)
                                     gCCGAGCGCCAGCCCGTGCGATTCGCCGAAGGTGGTTACGCGAAAGAGTTGTCCAATTGTGCTTCCAGCCATCACGGCTCCGTTATTGTTGTGTTTGCgtg                                                            <  2:3984829/101‑1 (MQ=255)
                                      ccGAGCGCCAGCCCGTGCGATTCGCCGAAGGTGGTTACGCGAAAGAGTTGTCCAATTGTGCTTCCAGCCATCACGGCTCCGTTATTGTTGTGTTTGCgtgt                                                           <  1:3984680/101‑1 (MQ=255)
                                       cGAGCGCCAGCCCGTGCGATTCGCCGAAGGTGGTTACGCGAAAGAGTTGTCCAATTGTGCTTCCAGCCATCACGGCTCCGTTATTGTTGTGTTTGCgtg                                                            <  1:3984890/99‑1 (MQ=255)
                                       cGAGCGCCAGCCCGTGCGATTCGCCGAAGGTGGTTACGCGAAAGAGTTGTCCAATTGTGCTTCCAGCCATCACGGCTCCGTTATTGTTGTGTTTGCGTGtt                                                          <  2:3984600/101‑1 (MQ=255)
                                        gAGCGCCAGCCCGTGCGATTCGCCGAAGGTGGTTACGCGAAAGAGTTGTCCAATTGTGCTTCCAGCCATCACGGCTCCGTTATTGTTGTGTTTGCGTGttt                                                         >  1:3985253/1‑101 (MQ=255)
                                         aGCGCCAGCCCGTGCGATTCGCCGAAGGTGGTTACGCGAAAGAGTTGTCCAATTGTGCTTCCAGCCATCACGGCTCCGTTATTGTTGTGTTTGCGTGTTTa                                                        <  2:3984749/101‑1 (MQ=255)
                                            gcCAGCCCGTGCGATTCGCCGAAGGTGGTTACGCGAAAGAGTTGTCCAATTGTGCTTCCAGCCATCACGGCTCCGTTATTGTTGTGTTTGCGTGTTTACtt                                                     >  2:3985255/1‑101 (MQ=255)
                                             ccAGCCCGTGCGATTCGCCGAAGGTGGTTACGCGAAAGAGTTGTCCAATTGTGCTTCCAGCCATCACGGCTCCGTTATTGTTGTGTTTGCGTGTTTACTTa                                                    <  2:3984904/101‑1 (MQ=255)
                                             ccAGCCCGTGCGATTCGCCGAAGGTGGTTACGCGAAAGAGTTGTCCAATTGTGCTTCCAGCCATCACGGCTCCGTTATTGTTGTGTTTGCGTGTTTACTTa                                                    >  2:3985259/1‑101 (MQ=255)
                                              cAGCCCGTGCGATTCGCCGAAGGTGGTTACGCGAAAGAGTTGTCCAATTGTGCTTCCAGCCATCACGGCTCCGTTATTGTTGTGTTTGCGTGTTTACTTaa                                                   >  1:3985261/1‑101 (MQ=255)
                                              cAGCCCGTGCGATTCGCCGAAGGTGGTTACGCGAAAGAGTTGTCCAATTGTGCTTCCAGCCATCACGGCTCCGTTATTGTTGTGTTTGCGTGTTTACTTaa                                                   <  1:3984724/101‑1 (MQ=255)
                                                gCCCGTGCGATTCGCCGAAGGTGGTTACGCGAAAGAGTTGTCCAATTGTGCTTCCAGCCATCACGGCTCCGTTATTGTTGTGTTTGCGTGTTTACTTAATc                                                 <  2:3984786/101‑1 (MQ=255)
                                                 cccGTGCGATTCGCCGAAGGTGGTTACGCGAAAGAGTTGTCCAATTGTGCTTCCAGCCATCACGGCTCCGTTATTGTTGTGTTTGCGTGTTTACTTAATCt                                                <  1:3984794/101‑1 (MQ=255)
                                                 cccGTGCGATTCGCCGAAGGTGGTTACGCGAAAGAGTTGTCCAATTGTGCTTCCAGCCATCACGGCTCCGTTATTGTTGTGTTTGCGTGTTTACTTAATCt                                                <  1:3984981/101‑1 (MQ=255)
                                                 cccGTGCGATTCGCCGAAGGTGGTTACGCGAAAGAGTTGTCCAATTGTGCTTCCAGCCATCACGGCTCCGTTATTGTTGTGTTTGCGTGTTTACTTAATCt                                                >  1:3985265/1‑101 (MQ=255)
                                                  ccGTGCGATTCGCCGAAGGTGGTTACGCGAAAGAGTTGTCCAATTGTGCTTCCAGCCATCACGGCTCCGTTATTGTTGTGTTTGCGTGTTTACTTAATCtt                                               <  1:3984758/101‑1 (MQ=255)
                                                  ccGTGCGATTCGCCGAAGGTGGTTACGCGAAAGAGTTGTCCAATTGTGCTTCCAGCCATCACGGCTCCGTTATTGTTGTGTTTGCGTGTTTACTTAATCtt                                               <  2:3984819/101‑1 (MQ=255)
                                                    gTGCGATTCGCCGAAGGTGGTTACGCGAAAGAGTTGTCCAATTGTGCTTCCAGCCATCACGGCTCCGTTATTGTTGTGTTTGCGTGTTTACTTAATCTTta                                             <  2:3984851/101‑1 (MQ=255)
                                                     tGCGATTCGCCGAAGGTGGTTACGCGAAAGAGTTGTCCAATTGTGCTTCCAGCCATCACGGCTCCGTTATTGTTGTGTTTGCGTGTTTGCTTAATCTTtat                                            <  2:3984884/101‑1 (MQ=255)
                                                     tGCGATTCGCCGAAGGTGGTTACGCGAAAGAGTTGTCCAATTGTGCTTCCAGCCATCACGGCTCCGTTATTGTTGTGTTTGCGTGTTTACTTAATCTTtat                                            <  1:3984710/101‑1 (MQ=255)
                                                     tGCGATTCGCCGAAGGTGGTTACGCGAAAGAGTTGTCCAATTGTGCTTCCAGCCATCACGGCTCCGTTATTGTTGTGTTTGCGTGTTTACTTAATCTTtat                                            <  1:3984866/101‑1 (MQ=255)
                                                     tGCGATTCGCCGAAGGTGGTTACGCGAAAGAGTTGTCCAATTGTGCTTCCAGCCATCACGGCTCCGTTATTGTTGTGTTTGCGTGTTTACTTAATCTTtat                                            <  1:3984798/101‑1 (MQ=255)
                                                       cGATTCGCCGAAGGTGGTTACGCGAAAGAGTTGTCCAATTGTGCTTCCAGCCATCACGGCTCCGTTATTGTTGTGTTTGCGTGTTTACTTAATCTTTATaa                                          <  1:3984847/101‑1 (MQ=255)
                                                       cGATTCGCCGAAGGTGGTTACGCGAAAGAGTTGTCCAATTGTGCTTCCAGCCATCACGGCTCCGTTATTGTTGTGTTTGCGTGTTTACTTAATCTTTATaa                                          <  1:3985014/101‑1 (MQ=255)
                                                        gATTCGCCGAAGGTGGTTACGCGAAAGAGTTGTCCAATTGTGCTTCCAGCCATCACGGCTCCGTTATTGTTGTGTTTGCGTGTTTACTTAATCTTTATaaa                                         <  2:3984704/101‑1 (MQ=255)
                                                        gATTCGCCGAAGGTGGTTACGCGAAAGAGTTGTCCAATTGTGCTTCCAGCCATCACGGCTCCGTTATTGTTGTGTTTGCGTGTTTACTTAATCTTTATaaa                                         <  2:3984950/101‑1 (MQ=255)
                                                         aTTCGCCGAAGGTGGTTACGCGAAAGAGTTGTCCAATTGTGCTTCCAGCCATCACGGCTCCGTTATTGTTGTGTTTGCGTGTTTACTTAATCTTTATAAAt                                        <  2:3984988/101‑1 (MQ=255)
                                                         aTTCGCCGAAGGTGGTTACGCGAAAGAGTTGTCCAATTGTGCTTCCAGCCATCACGGCTCCGTTATTGTTGTGTTTGCGTGTTTACTTAATCTTTATAAAt                                        >  1:3985271/1‑101 (MQ=255)
                                                          ttCGCCGAAGGTGGTTACGCGAAAGAGTTGTCCAATTGTGCTTCCAGCCATCACGGCTCCGTTATTGTTGTGTTTGCGTGTTTACTTAATCTTTATAAATc                                       <  2:3984808/101‑1 (MQ=255)
                                                          ttCGCCGAAGGTGGTTACGCGAAAGAGTTGTCCAATTGTGCTTCCAGCCATCACGGCTCCGTTATTGTTGTGTTTGCGTGTTTACTTAATCTTTATAAATc                                       >  2:3985274/1‑101 (MQ=255)
                                                          ttCGCCGAAGGTGGTTACGCGAAAGAGTTGTCCAATTGTGCTTCCAGCCATCACGGCTCCGTTATTGTTGTGTTTGCGTGTTTACTTAATCTTTATAAATc                                       <  1:3984550/101‑1 (MQ=255)
                                                          ttCGCCGAAGGTGGTTACGCGAAAGAGTTGTCCAATTGTGCTTCCAGCCATCACGGCTCCGTTATTGTTGTGTTTGCGTGTTTACTTAATCTTTATAAATc                                       <  1:3984881/101‑1 (MQ=255)
                                                           tCGCCGAAGGTGGTTACGCGAAAGAGTTGTCCAATTGTGCTTCCAGCCATCACGGCTCCGTTATTGTTGTGTTTGCGTGTTTACTTAATCTTTATAAATcg                                      <  2:3984895/101‑1 (MQ=255)
                                                            cGCCGAAGGTGGTTACGCGAAAGAGTTGTCCAATTGTGCTTCCAGCCATCACGGCTCCGTTATTGTTGTGTTTGCGTGTTTACTTAATCTTTATAAATcgc                                     >  1:3985275/1‑101 (MQ=255)
                                                               cGAAGGTGGTTACGCGAAAGAGTTGTCCAATTGTGCTTCCAGCCATCACGGCTCCGTTATTGTTGTGTTTGCGTGTTTACTTAATCTTTATAAATCGCGaa                                  >  2:3985279/1‑101 (MQ=255)
                                                                gAAGGTGGTTACGCGAAAGAGTTGTCCAATTGTGCTTCCAGCCATCACGGCTCCGTTATTGTTGTGTTTGCGTGTTTACTTAATCTTTATAAATCGCGaaa                                 <  1:3984838/101‑1 (MQ=255)
                                                                 aaGGTGGTTACGCGAAAGAGTTGTCCAATTGTGCTTCCAGCCATCACGGCTCCGTTATTGTTGTGTTTGCGTGTTTACTTAATCTTTATAAATCGCGAAAt                                <  1:3984827/101‑1 (MQ=255)
                                                                 aaGGTGGTTACGCGAAAGAGTTGTCCAATTGTGCTTCCAGCCATCACGGCTCCGTTATTGTTGTGTTTGCGTGTTTACTTAATCTTTATAAATCGCGAAAt                                >  1:3985285/1‑101 (MQ=255)
                                                                  aGGTGGTTACGCGAAAGAGTTGTCCAATTGTGCTTCCAGCCATCACGGCTCCGTTAttgtt                                                                       >  2:3985286/1‑61 (MQ=255)
                                                                    gtggtTACGCGAAAGAGTTGTCCAATTGTGCTTCCAGCCATCACGGCTCCGTTATTGTTGTGTTTGCGTGTTTACTTAATCTTTATAAATCGCGAAATGtt                             <  1:3984811/101‑1 (MQ=255)
                                                                    gtggtTACGCGAAAGAGTTGTCCAATTGTGCTTCCAGCCATCACGGCTCCGTTATTGTTGTGTTTGCGTGTTTACTTAATCTTTATAAATCGCGAAATGtt                             <  1:3984859/101‑1 (MQ=255)
                                                                    gtggtTACGCGAAAGAGTTGTCCAATTGTGCTTCCAGCCATCACGGCTCCGTTATTGTTGTGTTTGCGTGTTTACTTAATCTTTATAAATCGCGAAATGtt                             <  2:3984982/101‑1 (MQ=255)
                                                                    gtggtTACGCGAAAGAGTTGTCCAATTGTGCTTCCAGCCATCACGGCTCCGTTATTGTTGTGTTTGCGTGTTTACTTAATCTTTATAAATCGCGAAATGtt                             <  2:3984921/101‑1 (MQ=255)
                                                                    gtggtTACGCGAAAGAGTTGTCCAATTGTGCTTCCAGCCATCACGGCTCCGTTATTGTTGTGTTTGCGTGTTTACTTAATCTTTATAAATCGCGAAATGtt                             <  1:3984737/101‑1 (MQ=255)
                                                                      ggtTACGCGAAAGAGTTGTCCAATTGTGCTTCCAGCCATCACGGCTCCGTTATTGTTGTGTTTGCGTGTTTACTTAATCTTTATAAATCGCGAAATGTtct                           <  2:3984714/101‑1 (MQ=255)
                                                                       gtTACGCGAAAGAGTTGTCCAATTGTGCTTCCAGCCATCACGGCTCCGTTATTGTTGTGTTTGCGTGTTTACTTAATCTTTATAAATCGCGAAATGTtatc                          <  1:3984901/101‑1 (MQ=255)
                                                                         tACGCGAAAGAGTTGTCCAATTGTGCTTCCAGCCATCACGGCTCCGTTATTGTTGTGTTTGCGTGTTTACTTAATCTTTATAAATCGCGAAATGTTCTCGt                        <  1:3984779/101‑1 (MQ=255)
                                                                         tACGCGAAAGAGTTGTCCAATTGTGCTTCCAGCCATCACGGCTCCGTTATTGTTGTGTTTGCGTGTTTACTTAATCTTTATAAATCGCGAAATGTTCTCGt                        <  1:3984938/101‑1 (MQ=255)
                                                                         tACGCGAAAGAGTTGTCCAATTGTGCTTCCAGCCATCACGGCTCCGTTATTGTTGTGTTTGCGTGTTTACTTAATCTTTATAAATCGCGAAATGTTCTCGt                        <  2:3985079/101‑1 (MQ=255)
                                                                         tACGCGAAAGAGTTGCCCAATTGTGCTTCCAGCCATCACGGCTCCGTTATTGTTGTGTTTGCGTGTTTACTTAATCTTTATAAATCGCGAAATGTtct                           <  2:3985056/98‑1 (MQ=255)
                                                                          aCGCGAAAGAGTTGTCCAATTGTGCTTCCAGCCATCACGGCTCCGTTATTGTTGTGTTTGCGTGTTTACTTAATCTTTATAAATCGCGAAATGTTCTCGTg                       >  1:3985293/1‑101 (MQ=255)
                                                                           cgcgAAAGAGTTGTCCAATTGTGCTTCCAGCCATCACGGCTCCGTTATTGTTGTGTTTGCGTGTTTACTTAATCTTTATAAATCGCGAAATGTTCTCGTGc                      >  2:3985294/1‑101 (MQ=255)
                                                                             cgAAAGAGTTGTCCAATTGTGCTTCCAGCCATCACGGCTCCGTTATTGTTGTGTTTGCGTGTTTACTTAATCTTTATAAATCGCGAAATGTTCTCGTGCgg                    >  1:3985299/1‑101 (MQ=255)
                                                                             cgAAAGAGTTGTCCAATTGTGCTTCCAGCCATCACGGCTCCGTTATTGTTGTGTTTGCGTGTTTACTTAATCTTTATAAATCGCGAAATGTTCTCGTGCgg                    >  2:3985300/1‑101 (MQ=255)
                                                                               aaaGAGTTGTCCAATTGTGCTTCCAGCCATCACGGCTCCGTTATTGTTGTGTTTGCGTGTTTACTTAATCTTTATAAATCGCGAAATGTTCTCGTGCGGc                   <  2:3984769/100‑1 (MQ=255)
                                                                               aaaGAGTTGTCCAATTGTGCTTCCAGCCATCACGGCTCCGTTATTGTTGTGTTTGCGTGTTTACTTAATCTTTATAAATCGCGAAATGTTCTCGTGCGGCa                  >  1:3985305/1‑101 (MQ=255)
                                                                               aaaGAGTTGTCCAATTGTGCTTCCAGCCATCACGGCTCCGTTATTGTTGTGTTTGCGTGTTTACTTAATCTTTATAAATCGCGAAATGTTCTCGTGCGGCa                  <  1:3984772/101‑1 (MQ=255)
                                                                                 agagTTGTCCAATTGTGCTTCCAGCCATCACGGCTCCGTTATTGTTGTGTTTGCGTGTTTACTTAATCTTTATAAATCGCGAAATGTTCTCGTGCGGCAAt                <  1:3984973/101‑1 (MQ=255)
                                                                                 agagTTGTCCAATTGTGCTTCCAGCCATCACGGCTCCGTTATTGTTGTGTTTGCGTGTTTACTTAATCTTTATAAATCGCGAAATGTTCTCGTGCGGCAAt                >  2:3985309/1‑101 (MQ=255)
                                                                                  gagTTGTCCAATTGTGCTTCCAGCCATCACGGCTCCGTTATTGTTGTGTTTGCGTGTTTACTTAATCTTTATAAATCGCGAAATGTTCTCGTGCGGCAATa               >  1:3985311/1‑101 (MQ=255)
                                                                                       gTCCAATTGTGCTTCCAGCCATCACGGCTCCGTTATTGTTGTGTTTGCGTGTTTACTTAATCTTTATAAATCGCGAAATGTTCTCGTGCGGCAATAAgctg          >  2:3985318/1‑101 (MQ=255)
                                                                                       gTCCAATTGTGCTTCCAGCCATCACGGCTCCGTTATTGTTGTGTTTGCGTGTTTACTTAATCTTTATAAATCGCGAAATGTTCTCGTGCGGCAATAAgctg          >  1:3985319/1‑101 (MQ=255)
                                                                                       gTCCAATTGTGCTTCCAGCCATCACGGCTCCGTTATTGTTGTGTTTGCGTGTTTACTTAATCTTTATAAATCGCGAAATGTTCTCGTGCGGCAATAAgctg          <  1:3985061/101‑1 (MQ=255)
                                                                                         ccAATTGTGCTTCCAGCCATCACGGCTCCGTTATTGTTGTGTTTGCgtgt                                                           >  2:3985326/1‑50 (MQ=255)
                                                                                         ccAATTGTGCTTCCAGCCATCACGGCTCCGTTATTGTTGTGTTTGCGTGTTTACTTAATCTTTATAAATCGCGAAATGTTCTCGTGCGGCAATAAgc            >  1:3985325/1‑97 (MQ=255)
                                                                                           aaTTGTGCTTCCAGCCATCACGGCTCCGTTATTGTTGTGTTTGCGTGTTTACTTAATCTTTATAAATCGCGAAATGTTCTCGTGCGGCAATAAGCTGctct      <  2:3985011/101‑1 (MQ=255)
                                                                                            aTTGTGCTTCCAGCCATCACGGCTCCGTTATTGTTGTGTTTGCGTGTTTACTTAATCTTTATAAATCGCGAAATGTTCTCGTGCGGCAATAAGCTGCTCtt     >  1:3985332/1‑101 (MQ=255)
                                                                                            aTTGTGCTTCCAGCCATCACGGCTCCGTTATTGTTGTGTTTGCGTGTTTACTTAATCTTTATAAATCGCGAAATGTTCTCGTGCGGCAATAAGCTGCTCtt     >  2:3985333/1‑101 (MQ=255)
                                                                                             ttGTGCTTCCAGCCATCACGGCTCCGTTATTGTTGTGTTTGCGTGTTTACTTAATCTTTATAAATCGCGAAATGTTCTCGTGCGGCAATAAGCTGCTCttt    <  2:3984931/101‑1 (MQ=255)
                                                                                              tgtgCTTCCAGCCATCACGGCTCCGTTATTGTTGTGTTTGCGTGTTTACTTAATCTTTATAAATCGCGAAATGTTCTCGTGCGGCAATAAGCTGCTCTTTg   <  2:3985000/101‑1 (MQ=255)
                                                                                              tgtgCTTCCAGCCATCACGGCTCCGTTATTGTTGTGTTTGCGTGTTTACTTAATCTTTATAAATCGCGAAATGTTCTCGTGCGGCAATAAGCTGCTCTTTg   >  1:3985338/1‑101 (MQ=255)
                                                                                                  |                                                                                                 
TCAGCGGAATGCCTGGCGGAACACCATCGACGATGCAGCCGAGCGCCAGCCCGTGCGATTCGCCGAAGGTGGTTACGCGAAAGAGTTGTCCAATTGTGTTTCCAGCCATCACGGCTCCGTTATTGTTGTGTTTGCGTGTTTACTTAATCTTTATAAATCGCGAAATGTTCTCGTGCGGCAATAAGCTGCTCTTTGG  >  NC_000913/2447365‑2447560

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