Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,304,813 (C)7→6 100% coding (1005/1245 nt) mtr ← tryptophan transporter of high affinity

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,304,8070C.96.0% 765.4 / 18.4 176coding (1011/1245 nt)mtrtryptophan transporter of high affinity
Reads supporting (aligned to +/- strand):  ref base C (3/2);  new base . (83/86);  total (86/90)
Fisher's exact test for biased strand distribution p-value = 6.80e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.61e-01

Only 100 of 191 total aligned reads displayed.

CGATGCACGGGCTAACAGCGCCGGAACAATTGCCGCCCAGATGGTAGCCGCTAAACCAGCATAACCAATGGCGTACAGGAATCCGTTCGGGAACAACAGCCCCCCCACAACTGGCGGGGCAAAGGTCAGCAATGCCGTTTTCAAGCGGCCCACAGCCGAGTCGTCGAAACCAAACAGATCTGCCAGATAGTCAAA  >  NC_000913/3304708‑3304902
                                                                                                   |                                                                                               
 gATGCACGGGCTAACAGCGCCGGAACAATTGCCGCCCAGATGGTAGCCGCTAAACCAGCATAACCAATGGCGTACAGGAATCCGTTCGGGAACAACAGcc                                                                                                <  2:290924/100‑1 (MQ=255)
   tGCACGGGCTAACAGCGCCGGAACAATTGCCGCCCAGATGGTAGCCGCTAAACCAGCATAACCAATGGCGTACAGGAATCCGTTCGGGAACAACAGcccc                                                                                              >  1:533525/1‑100 (MQ=255)
     cACGGGCTAACAGCGCCGGAACAATTGCTGCCTAGATTGTAGCCGCTAAACCAGCATAACCAATGGCGTACAGGAATCCGTTCGGGAACAACAGcccccc                                                                                            <  1:1019466/100‑1 (MQ=255)
     cACGGGCTAACAGCGCCGGAACAATTGCCGCCCAGATGGTAGCCGCTAAACCAGCATAACCAATGGCGTACAGGAATCCGTTCGGGAACAACAGcccccc                                                                                            <  2:2968087/100‑1 (MQ=255)
     cACGGGCTAACAGCGCCGGAACAATTGCCGCCCAGATGGTAGCCGCTAAACCAGCATAACCAATGGCGTACAGGAATCCGTTCGGGAACAACAGcccccc                                                                                            >  1:2644435/1‑100 (MQ=255)
     cACGGGCTAACAGCGCCGGAACAATTGCCGCCCAGATGGTAGCCGCTAAACAAACATAAACAATAACGTACAGGAATCCGTTCGGGAACAACAGcccccc                                                                                            >  1:31255/1‑100 (MQ=255)
      aCGGGCTAACAGCGCCGGAACAATTGCCGCCCAGATGGTAGCCGCTAAACCAGCATAACCAATGGCGTACAGGAATCCGTTCGGGAACAACAGcccccca                                                                                           >  2:3427552/1‑99 (MQ=255)
        gggCTAACAGCGCCGGAACAATTGCCGCCCAGATGGTAGCCGCTAAACCAGCATAACCAATGGCGTACAGGAATCCGTTCGGGAACAACAGccccccaca                                                                                         <  2:338510/100‑4 (MQ=255)
         ggCTAACAGCGCCGGAACAATTGCCGCCCAGATGGTAGCCGCTAAACCAGCATAACCAATGGCGTACAGGAATCCGTTCGGGAACAACAGccccccacaa                                                                                        >  2:3096408/1‑96 (MQ=255)
            tAACAGCGCCGGAACAATTGCCGCCCAGATGGTAGCCGCTAAACCAGCATAACCAAAGGCGTACAGGAATCCGTTCGGGAACAACAG‑CCCCCCACAACTg                                                                                    >  1:1490956/1‑100 (MQ=255)
             aaCAGCGCCGGAACAATTGCCGCCCAGATGGTAGCCGCTAAACCAGCATAACCAATGGCGTACAGGAATCCGTTCGGGAACAACAG‑CCCCCCACAACTgg                                                                                   >  1:1369096/1‑100 (MQ=255)
             aaCAGCGCCGGAACAATTGCCGCCAAAATGGTAGCCGCTAAACCAGCATAACCAATGGCGTACAGGCATCCGTCCGGGAACAACAG‑CCCCCCACCACTgg                                                                                   >  2:3072428/1‑100 (MQ=255)
               cAGCGCCGGAACAATTGCCGCCCAGATGGTAGCCGCTAAACCAGCATAACCAATGGCGTACAGGAATCCGTTCGGGAACAACAG‑CCCCCCACAACTGGCg                                                                                 >  2:2205148/1‑100 (MQ=255)
                aGCGCCGGAACAATTGCCGCCCAGATGGTAGCCGCTAAACCAGCATAACCAATGGCGTACAGGAATCCGTTCGGGAACAACAG‑CCCCCCACAACTGGCgg                                                                                >  1:2508351/1‑100 (MQ=255)
                 gcgcCGGAACAATTGCCGCCCAGATGGTAGCCGCTAAACCAGCATAACCAATGGCGTACAGGAATCCGTTCGGGAACAACAG‑CCCCCCACAACTGGc                                                                                  >  2:7922/1‑97 (MQ=255)
                 gcgcCGGAACAATTGCCGCCCAGATGGTAGCCGCTAAACCAGCATAACCAATGGCGTACAGGAATCCGTTCGGGAACAACAG‑CCCCCCACAACTGGCggg                                                                               <  1:1639002/100‑1 (MQ=255)
                 gcgcCGGAACAATTGCCGCCCAGATGGTAGCCGCTAAACCAGCATAACCAATGGCGTACAGGAATCCGTTCGGGAACAACAG‑CCCCCCACAACTGGCggg                                                                               <  2:11292/100‑1 (MQ=255)
                  cgcCGGAACAATTGCCGCCCAGATGGTAGCCGCTAAACCAGCATAACCAATGGCGTACAGGAATCCGTTCGGGAACAACAG‑CCCCCCACAACTGGCgggg                                                                              >  2:3490328/1‑100 (MQ=255)
                   gcCGGAACAATTGCCGCCCAGATGGTAGCCGCTAAACCAGCATAACCAATGGCGTACAGGAATCCGTTCGGGAACAACAG‑CCCCCCACAACTGGCGGGGc                                                                             >  2:2017578/1‑100 (MQ=255)
                   gcCGGAACAATTGCCGCCCAGATGGTAGCCGCTAAACCAGCATAACCAATGGCGTACAGGAATCCGTTCGGGAACAACAG‑CCCCCCACAACTGGCGGGGc                                                                             >  2:879400/1‑100 (MQ=255)
                   gcCGGAACAATTGCCGCCCAGATGGTAGCCGCTAAACCAGCATAACCAATGGCGTACAGGAATCCGTTCGGGAACAACAG‑CCCCCCACAACTGGCGGGGc                                                                             <  1:1758944/100‑1 (MQ=255)
                    ccGGAACAATTGCCGCCCAGATGGTAGCCGCTAAACCAGCATAACCAATGGCGTACAGGAATCCGTTCGGGAACAACAG‑CCCCCCACAACTGGCGGGGc                                                                             <  1:1404217/99‑1 (MQ=255)
                    ccGGAACAATTGCCGCCCAGATGGTAGCCGCTAAACCAGCATAACCAATGGCGTACAGGAATCCGTTCGGGAACAACAG‑CCCCCCACAACTGGCGGGGCa                                                                            <  2:1660474/100‑1 (MQ=255)
                      ggAACAATTGCCGCCCAGATGGTAGCCGCTAAACCAGCATAACCAATGGCGTACAGGAATCCGTTCGGGAACAACAG‑CCCCCCACAACTGGCGGGGCaaa                                                                          >  2:247262/1‑100 (MQ=255)
                       gAACAATTGCCGCCCAGATGGTAGCCGCTAAACCAGCATAACCAATGGCGTACAGGAATCCGTTCGGGAACAACAG‑CCCCCCACAACTGGCGGGGCAAAg                                                                         >  1:3584734/1‑100 (MQ=255)
                        aaCAATTGCCGCCCAGATGGTAGCCGCTAAACCAGCATAACCAATGGCGTACAGGAATCCGTTCGGGAACAACAG‑CCCCCCACAACTGGCGGGGCAAAgg                                                                        <  2:1074457/100‑1 (MQ=255)
                        aaCAATTGCCGCCCAGATGGTAGCCGCTAAACCAGCATAACCAATGGCGTACAGGAATCCGTTCGGGAACAACAG‑CCCCCCACAACTGGCGGGGCAAAgg                                                                        >  1:1207163/1‑100 (MQ=255)
                        aaCAATTGCCGCCCAGATGGTAGCCGCTAAACCAGCATAACCAATGGCGTACAGGAATCCGTTCGGGAACAACAG‑CCCCCCACAACTGGCGGGGCAAAgg                                                                        <  1:2033303/100‑1 (MQ=255)
                          cAATTGCCGCCCAGATGGTAGCCGCTAAACCAGCATAACCAATGGCGTACAGGAATCCGTTCGGGAACAACAG‑CCCCCCACAACTGGCGGGGCAAAGGTc                                                                      >  1:2407217/1‑100 (MQ=255)
                          cAATTGCCGCCCAGATGGTAGCCGCTAAACCAGCATAACCAATGGCGTACAGGAATCCGTTCGAGAACAACAG‑CCCCCCACAACTGGCGGGGCAAAGGTc                                                                      >  1:1074572/1‑100 (MQ=255)
                            aTTGCCGCCTAGATGGTAGCCGCTAAACCAGCATAACCAATGGCGTACAGGAATCCGTTCGGGAACAACAG‑CCCCCCACAACTGGCGGGGCAAAGGTCAg                                                                    <  1:3926506/100‑1 (MQ=255)
                            aTTGCCGCCCAGATGGTAGCCGCTAAACCAGCATAACCAATGGCGTACAGGAATCCGTTCGGGAACAACAG‑CCCCCCACAACTGGCGGGGCAAAGGTCAg                                                                    <  1:3812136/100‑1 (MQ=255)
                             ttGCCGCCCAGATGGTAGCCGCTAAACCAGCATAACCAATGGCGTACAGGAATCCGTTCGGGAACAACAG‑CCCCCCACAACTGGCGGGGCAAAGGTCAGc                                                                   <  1:3903220/100‑1 (MQ=255)
                              tGCCGCCCAGATGGTAGCCGCTAAACCAGCATAACCAATGGCGTACAGGAATCCGTTCGGGAACAACAG‑CCCCCCACAACTGGCGGGGCAAAGGTCAGCa                                                                  >  1:3489563/1‑100 (MQ=255)
                               gccgccCAGATGGTAGCCGCTAAACCAGCATAACCAATGGCGTACAGGAATCCGTTCGGGAACAACAG‑CCCCCCACAACTGGCGGGGCAAAGGTCAGCaa                                                                 >  1:108603/1‑100 (MQ=255)
                               gccgccCAGATGGTAGCCGCTAAACCAGCATAACCAATGGCGTACAGGAATCCGTTCGGGAACAACAG‑CCCCCCACAACTGGCGGGGCAAAGGTCAGCaa                                                                 >  2:1689960/1‑100 (MQ=255)
                               gccgccCAGATGGTAGCCGCTAAACCAGCATAACCAATGGCGTACAGGAATCCGTTCGGGAACAACAG‑CCCCCCACAACTGGCGGGGCAAAGGTCAGCaa                                                                 <  1:3646036/100‑1 (MQ=255)
                                ccgccCAGATGGTAGCCGCTAAACCAGCATAACCAATGGCGTACAGGAATCCGTTCGGGAACAACAG‑CCCCCCACAACTGGCGGGGCAAAGGTCAGCAAt                                                                >  2:3058462/1‑100 (MQ=255)
                                ccgccCAGATGGTAGCCGCTAAACCAGCATAACCAATGGCGTACAGGAATCCGTTCGGGAACAACAG‑CCCCCCACAACTGGCGGAGCAAAGGTCAGCAAt                                                                >  1:3780515/1‑100 (MQ=255)
                                  gccCAGATGGTAGCCGCTAAACCAGCATAACCAATGGCGTACAGGAATCCTTTCGGGAACAACAG‑CCCCCCACAACTGGCGGGGCAAAGGTCAGCAATGc                                                              <  1:1969779/100‑1 (MQ=255)
                                   cccAGATGGTAGCCGCTAAACCAGCATAACCAATGGCGTACAGGAATCCGTTCGGGAACAACAG‑CCCCCCACAACTAGCGGGGCAAAGGTCAGCAATGcc                                                             >  1:2126384/1‑100 (MQ=255)
                                    ccAGATGGTAGCCGCTAAACCAGCATAATCTATGGCGTACAGGTATCCGTTTGGGAACAACAG‑CCCCCCACAACTTGCGGTGCAAAGGTCAGCAATGCCg                                                            <  2:3320599/100‑1 (MQ=255)
                                     cAGATGGTAGCCGCTAAACCAGCATAACCAATGGCGTACAGGAATCCGTTCGGGAACAACAG‑CCCCCCACAACTGGCGGGGCAAAGGTCAGCAATGCCGt                                                           <  2:2631946/100‑1 (MQ=255)
                                        aTGGTAGCCGCTAAACCAGCATAACCAATGGCGTACAGGAATCCGTTCGGGAACAACAGCCCCCCCACAACTGGCGGGGCAAAGGTCAGCAATGCCGttt                                                         <  1:2955095/100‑1 (MQ=255)
                                        aTGGTAGCCGCTAAACCAGCATAACCAATGGCGTACAGGAATCCGTTCGGGAACAACAG‑CCCCCCACAACTGGCGGGGCAAAGGTCAGCAATGCCGtttt                                                        <  1:2369346/100‑1 (MQ=255)
                                           gTAGCCGCTAAACCAGCATAACCAATGGCGTACAGGAATCCGTTCGGGAACAACAG‑CCCCCCACAACTGGCGGGGCAAAGGTCAGCAATGCCGTTTTCaa                                                     <  2:488161/100‑1 (MQ=255)
                                            tAGCCGCTAAACCAGCATAACCAATGGCGTACAGGAATCCGTTCGGGAACAACAG‑CCCCCCACAACTGGCGGGGCAAAGGTCAGCAATGCCGTTTTCAAg                                                    <  2:2705456/100‑1 (MQ=255)
                                                cGCTAAACCAGCATAACCAATGGCGTACAGGAATCCGTTCGGGAACAACAG‑CCCCCCACAACTGGCGGGGCAAAGGTCAGCAATGCCGTTTTCAAGCGGc                                                >  1:848658/1‑100 (MQ=255)
                                                 gCTAAACCAGCATAACCAATGTCGTACAGTAATCCGTTCGGGAACAACAG‑CCCCCCACAACTGGCGGGGCAAAGGTCAGCAATGCCGTTTTCAAGCGGcc                                               <  1:684641/100‑1 (MQ=255)
                                                  cTAAACCAGCATAACCAATGGCGTACAGGAATCCGTTCGGGAACAACAG‑CCCCCCACAACTGGCGGGGCAAAGGTCAGCAATGCCGTTTTCAAGCGGccc                                              <  2:533518/100‑1 (MQ=255)
                                                   tAAACCAGCATAACCAATGGCGTACAGGAATCCGTTCGGGAACAACAG‑CCCCCCACAACTGGCGGGGCAAAGGTCAGCAATGCCGTTTTCAAGCGGCCca                                             <  1:1859920/100‑1 (MQ=255)
                                                   tAAACCAGCATAACCAATGGCGTACAGGAATCCGTTCGGGAACAACAG‑CCCCCCACAACTGGCGGGGCAAAGGTCAGCAATGCCGTTTTCAAGCGGCCca                                             >  2:3225227/1‑100 (MQ=255)
                                                   tAAACCAGCATAACCAATGGCGTACAGGAATCCGTTCGGGAACAACAG‑CCCCCCACAACTGGCGGGGCAAAGGTCAGCAATGCCGTTTTCAAGCGGCCca                                             >  2:1493414/1‑100 (MQ=255)
                                                   tAAACCAGCATAACCAATGGCGTACAGGAATCCGTTCCGGAACAACAG‑CCCCCCACAACTGGCGGGGCAAAGGTCAGCAATGCCGTTTTCAAGCGGCCca                                             >  2:2082550/1‑100 (MQ=255)
                                                      atcaccATAACCAATGCCGTAAATGCTTCCGTTCATGAACAACA‑TCCCCCCACACCGGGCAGCGCAAATCTCATCAATGCCGTTTTCAATCGGcc                                               <  2:3773420/90‑1 (MQ=255)
                                                      aCCAGCATAACCAATGGCGTACAGGAATCCGTTCGGGAACAACAG‑CCCCCCACAACTGGCGGGGCAAAGGTCAGCAATGCCGTTTTCAAGCGGCCCACAg                                          <  2:1920479/100‑1 (MQ=255)
                                                       ccAGCATAACCTATGGTGTATAGGAATCCGTTTGGGAACAACAG‑CCCCCCACAACTGGCGGGGCAAAGGTCAGCAATGCCGTTTTCAAGCGGCCCACAGc                                         <  1:1560104/100‑1 (MQ=255)
                                                       ccAGCATAACCAATGGCGTACAGGAATCCGTTCGGGAACAACAG‑CCCCCCACAACTGGCGGGGCAAAGGTCAGCAATGCCGTTTTCAAGCGGCCCACAGc                                         >  1:733262/1‑100 (MQ=255)
                                                         aGCATAACCAATGGCGTACAGGAATCCGTTCGCGAACAATAG‑CCCCCCACAACTGGCGGGGCAAAGGTCAGCAATGCCGTTTTCAAGCGGCCCACAGCCg                                       <  2:2753601/100‑1 (MQ=255)
                                                          gCATAACCAATGGCGTACAGGAATCCGTTCGGGAACAACAG‑CCCCCCACAACTGGCGGGGCAAAGGTCAACAATGCCGTTTTCAAGCGGCCCACAGCCGa                                      >  1:1599805/1‑100 (MQ=255)
                                                            aTAACCAATGGCGTACAGGAATCCGTTCGGGAACAACAGCCCCCCC‑CAACTGGCGGGGCAAAGGTCAGCAATGCCGTTTTCAAGCGGCCCACAGCCGAgt                                    <  2:451633/100‑1 (MQ=255)
                                                            aTAACCAATGGCGTACAGGAATCCGTTCGGGAACAACAG‑CCCCCCACAACTGGCGGGGCAAAGGTCAGCAATGCCGTTTTCAAGCGGCCCACAGCCGAgt                                    <  2:3465146/100‑1 (MQ=255)
                                                            aTAACCAATGGCGTACAGGAATCCGTTCGGGAACAACAG‑CCCCCCACAACTGGCGGGGCAAAGGTCAGCAATGCCGTCTTCAAGCGGCCCACAGCCg                                       >  2:3717387/1‑97 (MQ=255)
                                                             tAACCAATGGCGTACAGGAATCCGTTCGGGAACAACAG‑CCCCCCACAACTGGCGGGGCAAAGGTCAGCAATGCCGTTTTCAAGCGGCCCACAGCCGAgtc                                   >  1:3237818/1‑100 (MQ=255)
                                                             tAACCAATGGCGTACAGGAATCCGTTCGGGAACAACAG‑CCCCCCACAACTGGCGGAGCAAAGGTCAGCAAAACCGTTTTCAACCGACCCACAGCCGAgtc                                   >  1:911883/1‑100 (MQ=255)
                                                                ccAATGGCGTACAGGAATCCGTTCGGGAACAACAG‑CCCCCCACAACTGGCGGGGCAAAGGTCAGCAATGCCGTTTTCAAGCGGCCCACAGCCGAgtcgtc                                >  1:3838758/1‑100 (MQ=255)
                                                                   aTGGCGTACAGGAATCCGTGCGGGAACAACAG‑CCCCCCACAACTGGCGGGGCAAAGGTCAGCAATGCCGTTTTCAAGCGGCCCACAGCCGAGTCGTCGaa                             >  2:3727196/1‑100 (MQ=255)
                                                                    tGGCTTACAGGAATCCGTTCGGGAATAACAG‑CCCCCCACAACTGGCGGGGCAAAGGTCAGCAATGCCGTTTTCAAGCGGCCCACAGCCGAGtcgtcg                               <  1:3552155/97‑1 (MQ=255)
                                                                    tGGCGTACAGGAATCCGTTCGGGAACAACAG‑CCCCCCACAACTGGCGGGGCAAAGGTCAGCAATGCCGTTTTCAAGCGGCCCACAGCCGAGTCGTCGaaa                            <  1:3848538/100‑1 (MQ=255)
                                                                      gCGTACAGGAATCCGTGCGGGAACAACAG‑CCCCCCACAACTGGCGGGGCAAAGGTCAGCAATGCCGTTTTCAAGCGGCCCACAGCCGAGTCGTCGAAAc                           <  2:2999185/99‑1 (MQ=255)
                                                                       cGTACAGGAATCCGTTCGGGAACAACAG‑CCCCCCACAACTGGCGGGGCAAAGGTCAGCAATGCCGTTTTCAAGCGGCCCACAGCCGAGTCGTCGAAACCa                         >  2:1757911/1‑100 (MQ=255)
                                                                          aCAGGAATCCGTTCGGGAACAACAG‑CCCCCCACAACTGGCGGGGCAAAGGTCAGCAATGCCGTTTTCAAGCGGCCCACAGCCGAGTCGTCGAAACCAAAc                      <  1:247267/100‑1 (MQ=255)
                                                                          aCAGGAATCCGTTCGGGAACAACAG‑CCCCCCACAACTGGCGGGGCAAAGGTCAGCAAGGCCGTTTTCAAGCGGCCCACAGCCGAGTCGTCGAAACCAAAc                      <  1:143935/100‑1 (MQ=255)
                                                                           cAGGAATCCGTTCGGGAACAACAG‑CCCCCCACAACTGGCGGGGCAAAGGTCAGCAATGCCGTTTTCAAGCGGCCCACAGCCGAGTCGTCGAAACCaaaaa                     >  1:3272542/1‑98 (MQ=255)
                                                                             ggAATCCGTTCGGGAACAACAG‑CCCCCCACAACTGGCGGGGCAAAGGTCAGCAATGCCGTTTTCAAGCGGCCCACAGCCGAGTCGTCGAAACCAAACAg                    >  2:2493329/1‑99 (MQ=255)
                                                                              gAATCCGTTCGGGAACAACAG‑CCCCCCACAACTGGCGGGGCAAAGGTCAGCAATGCCGTTTTCAAGCGGCCCACAGCCGAGTCGTCGAAACCAAACAg                    <  2:284934/98‑1 (MQ=255)
                                                                               atctcGTTCGGGAACAACAG‑CCCCCCACAACTGGCGGGGCAAAGGTCAGCAATGCCGTTTTCAAGCGGCCCACAGCCGAGTCGTCGAAACCAAACAGATc                 <  1:794279/96‑1 (MQ=255)
                                                                                aTCCGTTCGGGAACAACAGCCCCCCCACAACTGGCGGGGCAAAGGTCAGCAATGCCGTTTTCAAGCGGCCCACAGCCGAGTCGTCGAAACCAAACAGATc                 >  2:2237932/1‑100 (MQ=255)
                                                                                aTCCGTTCGGGAACAACAG‑CCCCCCACAACTGGCGGGGCAAAGGTCAGCAATGCCGTTTTCAAGCGGCCCACAGCCGAGTCGTCGAAACCAAACAGATCt                <  1:2874423/100‑1 (MQ=255)
                                                                                aTCCGTTCGGGAACAACAG‑CCCCCCACAACTGGCGGGGCAAAGGTCAGCAATGCCGTTTTCAAGCGGCCCACAGCCGAGTCGTCGAAACCAAACAGATCt                >  2:3636279/1‑100 (MQ=255)
                                                                                  ccGTTCGGGAACAACAG‑CCCCCCACAACTGGCGGGGCAAAGGTCAGCAATGCCGTTTTCAAGCGGCCCACAGCCGAGTCGTCGAAACCAAACAGATCTGc              <  1:1114329/100‑1 (MQ=255)
                                                                                   cGTTCGGGAACAACAG‑CCCCCCACAACTGGCGGGGCAAAGGTCAGCAATGCCGTTTTCAAGCGGCCCACAGCCGAGTCGTCGAAACCAAACAGATCTGcc             >  2:232446/1‑100 (MQ=255)
                                                                                   cGTTCGGGAACAACAG‑CCCCCCACAACTGGCGGGGCAAAGGTCAGCAATGCCGTTTTCAAGCGGCCCACAGCCGAGTCGTCGAAACCAAACAGATCTGc              >  1:393378/1‑99 (MQ=255)
                                                                                   cGTTCGGGAACAACAG‑CCCCCCACAACTGGCGGGGCAAAGGTCAGCAATGCCGTTTTCAAGCGGCCCACAGCCGAGTCGTCGAAACCAAACACATCTg               >  2:1706466/1‑98 (MQ=255)
                                                                                    gTTCGGGAACAACAG‑CCCCCCACAACTGGCGGGGCAAAGGTCAGCAATGCCGTTTTCAAGCGGCCCACAGCCGAGTCGTCGAAACCAAACAGATCTGCCa            >  1:3883934/1‑100 (MQ=255)
                                                                                    gTTCGGGAACAACAG‑CCCCCCACAACTGGCGGGGCAAAGGTCAGCAATGCCGTTTTCAAGCGGCCCACAGCCGAGTCGTCGAAACCAAACAGATCTGCCa            >  1:999210/1‑100 (MQ=255)
                                                                                     ttCGGGAACAACAG‑CCCCCCACAACTGGCGGGGCAAAGGTCAGCAATGCCGTTTTCAAGCGGCCCACAGCCGAGTCGTCGAAACCAAACAGATCTGCCAg           >  2:1286695/1‑100 (MQ=255)
                                                                                       cGGGAACAACAG‑CCCCCCACAACTGGCGGGGCAAAGGTCAGCAATGCCGTTTTCAAGCGGCCCACAGCCGAGTCGTCGAAACCAAACAGATCTGCCa            >  2:102868/1‑97 (MQ=255)
                                                                                       cGGGAACAACAG‑CCCCCCACAACTGGCGGGGCAAAGGTCAGCAATGCCGTTTTCAAGCGGCCCACAGCCGAGTCGTCGAAACCAAACAGATCTGCCAGAt         <  2:3489524/100‑1 (MQ=255)
                                                                                       cGGGAACAACAG‑CCCCCCACAACTGGCGGGGCAAAGGTCAGCAATGCCGTTTTCAAGCGGCCCACAGCCGAGTCGTCGAAACCAAACAGATCTGCCAGAt         <  1:2266403/100‑1 (MQ=255)
                                                                                        gggAACAACAG‑CCCCCCACAACTGGCGGGGCAAAGGTCAGCAATGCCGTTTTCAAGCGGCCCACAGCCGAGTCGTCGAAACCAAACAGATCTGCCAGATa        >  1:1561394/1‑100 (MQ=255)
                                                                                        gggAACAACAG‑CCCCCCACAACTGGCGGGGCAAAGGTCAGCAATGCCGTTTTCAAGCGGCCCACAGCCGAGTCGTCGAAACCAAACAGATCTGCCAGATa        >  2:1567392/1‑100 (MQ=255)
                                                                                         ggAACAACAG‑CCCCCCACAACTGGCGGGGCAAAGGTCAGCAATGCCGTTTTCAAGCGGCCCACAGCCGAGTCGTCGAAACCAAACAGATCTGCCAGATAg       >  1:3060163/1‑100 (MQ=255)
                                                                                           aacaacAG‑CCCCCCACAACTGGCGTGGCAAAGGTCAGCAATGCCGTTTTCAAGCGGCCCACAGCCGAGTCGTCGAAACCAAACAGATCTGCCAGATAGTc     <  1:475974/100‑1 (MQ=255)
                                                                                           aacaacAG‑CCCCCCACAACTGGCGGGGCAAAGGTCAGCAATGCCGTTTTCAAGCGGCCCACAGCCGAGTCGTCGAAACCAAACAGATCTGCCAGATAGTc     <  1:2852204/100‑1 (MQ=255)
                                                                                           aacaacAG‑CCCCCCACAACTGGCGGGGAAAAGGTCAGCAATGCCGTTTTCAAGCGGCCCACAGCCGAGTCGTCGAAACCAAACAGATCTGCCAGATAGTc     <  2:2769456/100‑1 (MQ=255)
                                                                                             taacaG‑CCCCCCACAACTGGCGGGGCAAAGGTCAGCAATGCCGTTTTCAAGCGGCCCACAGCCGAGTCGTCGAAACCAAACAGATCTGCCAGATAGTcaa   >  2:568849/2‑100 (MQ=255)
                                                                                             caacaG‑CCCCCCACAACTGGCGGGGCAAAGGTCAGCAATGCCGTTTTCAAGCGGCCCACAGCCGAGTCGTCGAAACCAAACAGATCTGCCAGATAGTcaa   <  1:3498835/100‑1 (MQ=255)
                                                                                             caacaG‑CCCCCCACAACTGGCGGGGCAAAGGTCAGCAATGCCGTTTTCAAGCGGCCCACAGCCGAGTCGTCGAAAACAAACAGATCTGCCAGATAGTcaa   >  1:478061/1‑100 (MQ=255)
                                                                                             caacaG‑CCCCCCACAACTGGCGGGGAAAAGGTCAGCAATGCCGTTTTCAAGCGTCCCACACCCGAGTCGTCGAAACCAAACAGATCTGCCAGATAGTcaa   >  2:2966333/1‑100 (MQ=255)
                                                                                                   |                                                                                               
CGATGCACGGGCTAACAGCGCCGGAACAATTGCCGCCCAGATGGTAGCCGCTAAACCAGCATAACCAATGGCGTACAGGAATCCGTTCGGGAACAACAGCCCCCCCACAACTGGCGGGGCAAAGGTCAGCAATGCCGTTTTCAAGCGGCCCACAGCCGAGTCGTCGAAACCAAACAGATCTGCCAGATAGTCAAA  >  NC_000913/3304708‑3304902

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