Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 370,440 A→G N314D (AAC→GAC)  mhpB → 2,3‑dihydroxyphenylpropionate 1,2‑dioxygenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913370,4400AG98.9% 875.0 / ‑3.6 273N314D (AAC→GAC) mhpB2,3‑dihydroxyphenylpropionate 1,2‑dioxygenase
Reads supporting (aligned to +/- strand):  ref base A (2/1);  new base G (108/162);  total (110/163)
Fisher's exact test for biased strand distribution p-value = 5.67e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 294 total aligned reads displayed.

GCGTTTGGCAACTGGCGTAGCGAAGGGCGTTATTACCGCCCAATCCCGGAGTGGATTGCCGGATTTGGCTCGTTAAGCGCCAGAACAGAGAACTGAATATGCAGGAGAAGATGATGAGTTATCAGCCACAAACCGAAGCCGCCACCAGCCGTTTTCTGAATGTAGAAGAAGCGGGTAAAA  >  NC_000913/370350‑370529
                                                                                          |                                                                                         
 cGTTTGGCAACTGGCGTAGCGAAGGGCGTTATTACCGCCCAATCCCGGAGTGGATTGCCGGATTTGGCTCGTTAAGCGCCAGAACAGAGGACt                                                                                        <  2:7602305/93‑1 (MQ=255)
 cGTTTGGCAACTGGCGTAGCGAAGGGCGTTATTACCGCCCAATCCCGGAGTGGATTGCCGGATTTGGCTCGTTAAGCGCCAGAACAGAGGACt                                                                                        <  1:1747412/93‑1 (MQ=255)
   tttGGCAACTGGCGTAGCGAAGGGCGTTATTACCGCCCAATCCCGGAGTGGATTGCCGGATTTGGCTCGTTAAGCGCCAGAACAGAGGACTGa                                                                                      >  2:4346294/1‑93 (MQ=255)
   tttGGCAACTGGCGTAGCGAAGGGCGTTATTACCGCCCAATCCCGGAGTGGATTGCCGGATTTGGCTCGTTAAGCGCCAGAACAGAGGACTGa                                                                                      >  2:6019432/1‑93 (MQ=255)
    ttGGCAACTGGCGTAGCGAAGGGCGTTATTACCGCCCAATCCCGGAGTGGATTGCCGGATTTGGCTCGTTAAGCGCCAGAACAGAGGACTGaa                                                                                     <  1:5227757/93‑1 (MQ=255)
    ttGGCAACTGGCGTAGCGAAGGGCGTTATTACCGCCCAATCCCGGAGTGGATTGCCGGATTTGGCTCGTTAAGCGCCAGAACAGAGGACTGaa                                                                                     <  1:7025411/93‑1 (MQ=255)
    ttGGCAACTGGCGTAGCGAAGGGCGTTATTACCGCCCAATCCCGGAGTGGATTGCCGGATTTGGCTCGTTAAGCGCCAGAACAGAGGACTGaa                                                                                     >  1:4175115/1‑93 (MQ=255)
    ttGGCAACTGGCGTAGCGAAGGGCGTTATTACCGCCCAATCCCGGAGTGGATTGCCGGATTTGGCTCCTTAAGCGCCATGACAGAGGACTGaa                                                                                     >  2:6830256/1‑93 (MQ=255)
      ggCAACTGGCGTAGCGAAGGGCGTTATTACCGCCCAATCCCGGAGTGGATTGCCGGATTTGGCTCGTTAAGCGCCAGAACAGAGGACTGAata                                                                                   <  1:775482/93‑1 (MQ=255)
      ggCAACTGGCGTAGCGAAGGGCGTTATTACCGCCCAATCCCGGAGTGGATTGCCGGATTTGGCTCGTTAAGCGCCAGAACAGAGGACTGAata                                                                                   >  1:7326869/1‑93 (MQ=255)
      ggCAACTGGCGTAGCGAAGGGCGTTATTACCGCCCAATCCCGGAGTGGATTGCCGGATTTGGCTCGTTAAGCGCCAGAACAGAGGACTGAata                                                                                   <  1:6428645/93‑1 (MQ=255)
      ggCAACTGGCGTAGCGAAGGGCGTTATTACCGCCCAATCCCGGAGTGGATTGCCGGATTTGGCTCGTTAAGCGCCAGAACAGAGGACTGAata                                                                                   <  1:6131922/93‑1 (MQ=255)
       gCAACTGGCGTAGCGAAGGGCGTTATTACCGCCCAATCCCGGAGTGGATTGCCGGATTTGGCTCGTTAAGCGCCAGAACAGAGGACTGAatat                                                                                  >  2:9385017/1‑93 (MQ=255)
         aaCTGGCGTAGCGAAGGGCGTTATTACCGCCCAATCCCGGAGTGGATTGCCGGATTTGGCTCGTTAAGCGCCAGAACAGAGGACTGAATATGc                                                                                >  1:9062386/1‑93 (MQ=255)
           cTGGCGTAGCGAAGGGCGTTATTACCGCCCAATCCCGGAGTGGATTGCCGGATTTGGCTCGTTAAGCGCCAGAACAGAGGACTGAATATGCAg                                                                              <  1:566931/93‑1 (MQ=255)
           cTGGCGTAGCGAAGGGCGTTATTACCGCCCAATCCCGGAGTGGATTGCCGGATTTGGCTCGTTAAGCGCCAGAACAGAGGACTGAATATGCAg                                                                              <  1:8505590/93‑1 (MQ=255)
           cTGGCGTAGCGAAGGGCGTTAATACCGCCCAATCCCGGAGTGGATTGCCGGATTTGGCTCGTTAAGCGCCAGAACAGAGGACTGAATATGCAg                                                                              >  2:9888640/1‑93 (MQ=255)
            tGGCGTAGCGAAGGGCGTTATTACCGCCCAATCCCGGAGTGGATTGCCGGATTTGGCTCGTTAAGCGCCAGAACAGAGGACTGAATATGCAgg                                                                             >  2:3551425/1‑93 (MQ=255)
            tGGCGTAGCGAAGGGCGTTATTACCGCCCAATCCCGGAGTGGATTGCCGGATTTGGCTCGTTAAGCGCCAGAACAGAGGACTGAATATGCAgg                                                                             <  1:515929/93‑1 (MQ=255)
             ggCGTAGCGAAGGGCGTTATTACCGCCCAATCCCGGAGTGGATTGCCGGATTTGGCTCGTTAAGCGCCAGAACAGAGGACTGAATATGCAGga                                                                            <  2:673303/93‑1 (MQ=255)
              gCGTAGCGAAGGGCGTTATTACCGCCCAATCCCGGAGTGGATTGCCGGATTTGGCTCGTTAAGCGCCAGAACAGAGGACTGAATATGCAGgag                                                                           >  1:9830995/1‑93 (MQ=255)
                 tAGCGAAGGGCGTTATTACCGCCCAATCCCGGAGTGGATTGCCGGATTTGGCTCGTTAAGCGCCAGAACAGAGGACTGAATATGCAGGagaag                                                                        <  1:5899857/93‑1 (MQ=255)
                  aGCGAAGGGCGTTATTACCGCCCAATCCCGGAGTGGATTGCCGGATTTGGCTCGTTAAGCGCCAGAACAGAGGACTGAATATGCAGGagaaga                                                                       <  2:5983093/93‑1 (MQ=255)
                   gCGAAGGGCGTTATTACCGCCCAATCCCGGAGTGGATTGCCGGATTTGGCTCGTTAAGCGCCAGAACAGAGGACTGAATATGCAGGAGAAgat                                                                      <  2:1101744/93‑1 (MQ=255)
                   gCGAAGGGCGTTATTACCGCCCAATCCCGGAGTGGATTGCCGGATTTGGCTCGTTAAGCGCCAGAACAGAGGACTGAATATGCAGGAGAAgat                                                                      <  2:9641788/93‑1 (MQ=255)
                     gAAGGGCGTTATTACCGCCCAATCCCGGAGTGGATTGCCGGATTTGGCTCGTTAAGCGCCAGAACAGAGGACTGAATATGCAGGAGAAgatga                                                                    <  2:6524947/93‑1 (MQ=255)
                     gAAGGGCGTTATTACCGCCCAATCCCGGAGTGGATTGCCGGATTTGGCTCGTTAAGCGCCAGAACAGAGGACTGAATATGCAGGAGAAgatga                                                                    <  1:4510896/93‑1 (MQ=255)
                     gAAGGGCGTTATTACCGCCCAATCCCGGAGTGGATTGCCGGATTTGGCTCGTTAAGCGCCAGAACAGAGGACTGAATATGCAGGAGAAgatga                                                                    <  2:747515/93‑1 (MQ=255)
                     gAAGGGCGTTATTACCGCCCAATCCCGGAGTGGATTGCCGGATTTGGCTCGTTAAGCGCCAGAACAGAGGACTGAATATGCAGGAGAAgatga                                                                    <  1:3873268/93‑1 (MQ=255)
                      aaGGGCGTTATTACCGCCCAATCCCGGAGTGGATTGCCGGATTTGGCTCGTTAAGCGCCAGAACAGAGGACTGAATATGCAGGAGAAgatgat                                                                   <  2:3658276/93‑1 (MQ=255)
                      aaGGGCGTTATTACCGCCCAATCCCGGAGTGGATTGCCGGATTTGGCTCGTTAAGCGCCAGAACAGAGGACTGAATATGCAGGAGAAgatgat                                                                   <  2:1664220/93‑1 (MQ=255)
                      aaGGGCGTTATTACCGCCCAATCCCGGAGTGGATTGCCGGATTTGGCTCGTTAAGCGCCAGAACAGAGGACTGAATATGCAGGAGAAgatgat                                                                   >  1:3667771/1‑93 (MQ=255)
                       aGGGCGTTATTACCGCCCAATCCCGGAGTGGATTGCCGGATTTGGCTCGTTAAGCGCCAGAACAGAGGACTGAATATGCAGGAGAAGatgatg                                                                  <  2:1647734/93‑1 (MQ=255)
                          gCGTTATTACCGCCCAATCCCGGAGTGGATTGCCGGATTTGGCTCGTTAAGCGCCAGAACAGAGGACTGAATATGCAGGAGAAGATGATGAGt                                                               >  1:9596484/1‑93 (MQ=255)
                            gTTATTACCGCCCAATCCCGGAGTGGATTGCCGGATTTGGCTCGTTAAGCGCCAGAACAGAGGACTGAATATGCAGGAGAAGATGATGAGTTa                                                             <  2:9509032/93‑1 (MQ=255)
                            gTTATTACCGCCCAATCCCGGAGTGGATTGCCGGATTTGGCTCGTTAAGCGCCAGAACAGAGGACTGAATATGCAGGAGAAGATGATGAGTTa                                                             >  1:8022852/1‑93 (MQ=255)
                             ttattaCCGCCCAATCCCGGAGTGGATTGCCGGATTTGGCTCGTTAAGCGCCAGAACAGAGGGCTGAATATGCAGGAGAAGATGATGAGTTAt                                                            >  2:8626560/1‑93 (MQ=255)
                              tattaCCGCCCAATCCCGGAGTGGATTGCCGGATTTGGCTCGTTAAGGGCCAGAACAGAGGACTGAATATGCAGGAGAAGGTGATGAGTTATc                                                           >  2:758436/1‑93 (MQ=255)
                               attaCCGCCCAATCCCGGAGTGGATTGCCGGATTTGGCTCGTTAAGCGCCAGAACAGAGGACTGAATATGCAGGAGAAGATGATGAGTTATCa                                                          >  1:543344/1‑93 (MQ=255)
                                ttaCCGCCCAATCCCGGAGTGGATTGCCGGATTTGGCTCGTTAAGCGCCAGAACAGAGGACTGAATATGCAGGAGAAGATGATGAGTTATCAg                                                         >  1:3808139/1‑93 (MQ=255)
                                ttaCCGCCCAATCCCGGAGTGGATTGCCGGATTTGGCTCGTTAAGCGCCAGAACAGAGGACTGAATATGCAGGAGAAGATGATGAGTTATCAg                                                         >  1:6537734/1‑93 (MQ=255)
                                    cGCCCAATCCCGGAGTGGATTGCCGGATTTGGCTCGTTAAGCGCCAGAACAGAGGACTGAATATGCAGGAGAAGATGATGAGTTATCAGCcac                                                     <  2:4014239/93‑1 (MQ=255)
                                    cGCCCAATCCCGGAGTGGATTGCCGGATTTGGCTCGTTAAGCGCCAGAACAGAGGACTGAATATGCAGGAGAAGATGATGAGTTATCAGCcac                                                     >  1:5043411/1‑93 (MQ=255)
                                    cGCCCAATCCCGGAGTGGATTGCCGGATGTGGCTCGTTAAGCGCCAGAACAGAGAACTGAATATGCAGGAGAAGATGATGAGTTATCAGCcac                                                     <  2:5087133/93‑1 (MQ=255)
                                     gCCCAATCCCGGAGTGGATTGCCGGATTTGGCTCGTTAAGCGCCAGAACAGAGGACTGAATATGCAGGAGAAGATGATGAGTTATCAGCcaca                                                    >  1:3049355/1‑93 (MQ=255)
                                        cAATCCCGGAGTGGATTGCCGGATTTGGCTCGTTAAGCGCCAGAACAGAGGACTGAATATGCAGGAGAAGATGATGAGTTATCAGCCACAAAc                                                 >  2:512729/1‑93 (MQ=255)
                                        cAATCCCGGAGTGGATTGCCGGATTTGGCTCGTTAAGCGCCAGAACAGAGGACTGAATATGCAGGAGAAGATGATGAGTTATCAGCCACAAAc                                                 <  1:5538839/93‑1 (MQ=255)
                                        cAATCCCGGAGTGGATTGCCGGATTTGGCTCGTTAAGCGCCAGAACAGAGGACTGAATATGCAGGAGAAGATGATGAGTTATCAGCCACAAAc                                                 >  1:1639709/1‑93 (MQ=255)
                                         caTCCCGGAGGTGATTTCCCGATTTGGCTCGTTAAGCGCCAGAACAGAGGACTGAATAGGCAGGAGAAGATGATGATTTATCAGCCACAAAcc                                                <  1:7639274/92‑1 (MQ=255)
                                           tCCCGGAGTGGATTGCCGGATTTGGCTCGTTAAGCGCCAGAACAGAGGACTGAATATGCAGGAGAAGATGATGAGTTATCAGCCACAAACCGa                                              <  1:1363415/93‑1 (MQ=255)
                                             ccGGAGTGGATTGCCGGATTTGGCTCGTTAAGCGCCAGAACAGAGGACTGAATATGCAGGAGAAGATGATGAGTTATCAGCCACAAACCGAAg                                            >  2:8131385/1‑93 (MQ=255)
                                              cGGAGTGGATTGCCGGATTTGGCTCGTTAGGCGCCAGAACAGAGGACTGAATATGCAGGAGAAGATGATGAGTTATCAGCCACAAACCGAAgc                                           >  2:6220149/1‑93 (MQ=255)
                                              cGGAGTGGATTGCCGGATTTGGCTCGTTAAGCGCCAGAACAGAGGACTGAATATGCAGGAGAAGATGATGAGTTATCAGCCACAAACCGAAgc                                           <  2:3489539/93‑1 (MQ=255)
                                              cGGAGTGGATTGCCGGATTTGGCTCGTTAAGCGCCAGAACAGAGGACTGAATATGCAGGAGAAGATGATGAGTTATCAGCCACAAACCGAAgc                                           <  2:1549053/93‑1 (MQ=255)
                                              cGGAGTGGATTGCCGGATTTGGCTCGTTAAGCGCCAGAACAGAGGACTGAATATGCAGGAGAAGATGATGAGTTATCAGCCACAAACCGAAgc                                           <  1:7602487/93‑1 (MQ=255)
                                                gAGTGGATTGCCGGATTTGGCTCGTTAAGCGCCAGAACAGAGGACTGAATATGCAGGAGAAGATGATGAGTTATCAGCCACAAACCGAAgccg                                         >  1:431739/1‑93 (MQ=255)
                                                 aGTGGATTGCCGGATTTGGCTCGTTAAGCGCCAGAACAGAGGACTGAATATGCAGGAGAAGATGATGAGTTATCAGCCACAAACCGAAgccgc                                        <  2:1334741/93‑1 (MQ=255)
                                                   tGGATTGCCGGATTTGGCTCGTTAAGCGCCAGAACAGAGGACTGAATATGCAGGAGAAGATGATGAGTTATCAGCCACAAAACGAAGACGcca                                      >  2:8863643/1‑93 (MQ=255)
                                                    ggATTGCCGGATTTGGCTCGTTAAGCGCCAGAACAGAGGACTGAATATGCAGGAGAAGATGATGAGTTATCAGCCACAAACCGAAGCCGccac                                     >  2:632553/1‑93 (MQ=255)
                                                    ggATTGCCGGATTTGGCTCGTTAAGCGCCAGAACAGAGGACTGAATATGCAGGAGAAGATGATGAGTTATCAGCCACAAACCGAAGCCGccac                                     >  2:9672869/1‑93 (MQ=255)
                                                      gttGCCGGATTTGGCTCGTTAAGCGCCAGAACAGAGGACTGAATATGCAGGAGAAGATGATGAGTTATCAGCCACAAACCGAAGCCGccacca                                   <  2:69280/92‑1 (MQ=255)
                                                        tGCCGGATTTGGCTCGTTAAGCGCCAGAACAGAGGACTGAATATGCAGGAGAAGATGATGAGTTATCAGCCACAAACCGAAGCCGCCACCAGc                                 <  1:4934782/93‑1 (MQ=255)
                                                        tGCCGGATTTGGCTCGTTAAGCGCCAGAACAGAGGACTGAATATGCAGGAGAAGATGATGAGTTATCAGCCACAAACCGAAGCCGCCACCAGc                                 <  2:217553/93‑1 (MQ=255)
                                                          ccGGATTTGGCTCGTTAAGCGCCAGAACAGAGGACTGAATATGCAGGAGAAGATGATGAGTTATCAGCCACAAACCGAAGCCGCCACCAGCCg                               <  1:885561/93‑1 (MQ=255)
                                                               tttGGCTCGTTAAGCGCCAGAACAGAGGACTGAATATGCAGGAGAAGATGATGAGTTATCAGCCACAAACCGAAGCCGCCACCAGCCGTTTTc                          <  2:2943390/93‑1 (MQ=255)
                                                               tttGGCTCGTTAAGCGCCAGAACAGAGGACTGAATATGCAGGAGAAGATGATGAGTTATCAGCCACAAACCGAAGCCGCCACCAGCCGTTTTc                          >  2:7162588/1‑93 (MQ=255)
                                                               tttGGCTCGTTAAGCGCCAGAACAGAGGACTGAATATGCAGGAGAAGATGATGAGTTATCAGCCACAAACCGAAGCCGCCACCAGCCGTTTTc                          <  2:5664712/93‑1 (MQ=255)
                                                               tttGGCTCGTTAAGCGCCAGAACAGAGGACTGAATATGCAGGAGAAGATGATGAGTTATCAGCCACAAACCGAAGCCGCCACCAGCCGTTTTc                          <  2:2066038/93‑1 (MQ=255)
                                                                ttGGCTCGTTAAGCGCCAGAACAGAGGACTGAATATGCAGGAGAAGATGATGAGTTATCAGCCACAAACCGAAGCCGCCACCAGCCGTTTTCt                         >  1:350183/1‑93 (MQ=255)
                                                                ttGGCTCGTTAAGCGCCAGAACAGAGGACTGAATATGCAGGAGAAGATGATGAGTTATCAGCCACAAACCGAAGCCGCCACCAGCCGTTTTCt                         >  2:7827482/1‑93 (MQ=255)
                                                                ttGGCTCGTTAAGCGCAAGAACAGAGGACTGAATATGCAGGAGAAGATGATGAGTTATCAGCCACAAACCGAAGCCGCCACCAGCCGTTTTCt                         <  1:10001852/93‑1 (MQ=255)
                                                                    cTCGTTAAGCGCCAGAACAGAGGACTGAATATGCAGGAGAAGATGATGAGTTATCAGCCACAAACCGAAGCCGCCACCAGCCGTTTTCTGAAt                     <  2:8747345/93‑1 (MQ=255)
                                                                    cTCGTTAAGCGCCAGAACAGAGGACTGAATATGCAGGAGAAGATGATGAGTTATCAGCCACAAACCGAAGCCGCCACCAGCCGTTTTCTGAAt                     <  2:3726538/93‑1 (MQ=255)
                                                                    cTCGTTAAGCGCCAGAACAGAGGACTGAATATGCAGGAGAAGATGATGAGTTATCAGCCACAAACCGAAGCCGCCACCAGCCGTTTTCTGAAt                     <  2:266909/93‑1 (MQ=255)
                                                                     tCGTTAAGCGCCAGAACAGAGGACTGAATATGCAGGAGAAGATGATGAGTTATCAGCCACAAACCGAAGCCGCCACCAGCCGTTTTCTGAATg                    <  1:7876244/93‑1 (MQ=255)
                                                                        ttAAGCGCCAGAACAGAGGACTGAATATGCAGGAGAAGATGATGAGTTATCAGCCACAAACCGAAGCCGCCACCAGCCGTTTTCTGAATGTag                 >  1:4813627/1‑93 (MQ=255)
                                                                        ttAAGCGCCAGAACAGAGGACTGAATATGCAGGAGAAGATGATGAGTTATCAGCCACAAACCGAAGCCGCCACCAGCCGTTTTCTGAATGTag                 >  1:9918168/1‑93 (MQ=255)
                                                                        ttAAGCGCCAGAACAGAGGACTGAATATGCAGGAGAAGATGATGAGTTATCAGCCACAAACCGAAGCCGCCACCAGCCGTTTTCTGAATGTag                 >  1:8251842/1‑93 (MQ=255)
                                                                        ttAAGCGCCAGAACAGAGGACTGAATATGCAGGAGAAGATGATGAGTTATCAGCCACAAACCGAAGCCGCCACCAGCCGTTTTCTGAATGTag                 >  1:1484046/1‑93 (MQ=255)
                                                                         tAAGCGCCAGAACAGAGGACTTAATATGCAGGAGAAGATGATGAGTTATCAGCCACAAACCGAAGCCGCCACCAGCCGTTTTCTGAATGTaga                <  1:2678477/93‑1 (MQ=255)
                                                                         tAAGCGCCAGAACAGAGGACTGAATATGCAGGAGAAGATGATGAGTTATCAGCCACAAACCGAAGCCGCCACCAGCCGTTTTCTGAATGTaga                <  2:182605/93‑1 (MQ=255)
                                                                         tAAGCGCCAGAACAGAGGACTGAATATGCAGGAGAAGATGATGAGTTATCAGCCACAAACCGAAGCCGCCACCAGCCGTTTTCTGAATGTaga                <  1:2409575/93‑1 (MQ=255)
                                                                         tAAGCGCCAGAACAGAGGACTGAATATGCAGGAGAAGATGATGAGTTATCAGCCACAAACCGAAGCCGCCACCAGCCGTTTTCTGAATGTaga                <  1:9309902/93‑1 (MQ=255)
                                                                         tAAGCGCCAGAACAGAGGACTGAATATGCAGGAGAAGATGATGAGTTATCAGCCACAAACCGAAGCCGCCACCAGCCGTTTTCTGAATGTaga                >  2:2358518/1‑93 (MQ=255)
                                                                         tAAGCGCCAGAACAGAGGACTGAATATGCAGGAGAAGATGATGAGTTATCAGCCACAAACCGAAGCCGCCACCAGCCGTTTTCTGAATGTaga                <  1:87055/93‑1 (MQ=255)
                                                                          aaGCGCCAGAACAGAGGACTGAATATGCAGGAGAAGATGATGAGTTATCAGCCACAAACCGAAGCCGCCACCAGCCGTTTTCTGAATGTagaa               >  2:6717480/1‑93 (MQ=255)
                                                                           aGCGCCAGAACAGAGGACTGAATATGCAGGAGAAGATGATGAGTTATCAGCCACAAACCGAAGCCGCCACCAGCCGTTTTCTGAATGTagaag              <  1:1254740/93‑1 (MQ=255)
                                                                             cgcCAGAACAGAGGACTGAATATGCAGGAGAAGATGATGAGTTATCAGCCACAAACCGAAGCCGCCACCAGCCGTTTTCTGAATGTAgaagaa            <  2:1924924/93‑1 (MQ=255)
                                                                              gcCAGAACAGAGGACTGAATATGCAGGAGAAGATGATGAGTTATCAGCCACAACCCGAAGCCGCCACCAGCCGTTTTCTGAATGTAGaagaag           <  2:4629673/93‑1 (MQ=255)
                                                                              gcCAGAACAGAGGACTGAATATGCAGGAGAAGATGATGAGTTATCAGCAACAAACCGAAGCCGCCACCAGCCGTTTTCTGAATGTAGaagaag           >  2:7702402/1‑93 (MQ=255)
                                                                                cAGAACAGAGGACTGAATATGCAGGAGAAGATGATGAGTTATCAGCCACAAACCGAAGCCGCCACCAGCCGTTTTCTGAATGTAGAAGAAGCg         <  2:7053952/93‑1 (MQ=255)
                                                                                cAGAACAGAGGACTGAATATGCAGGAGAAGATGATGAGTTATCAGCCACAAACCGAAGCCGCCACCAGCCGTTTTCTGAATGTAGAAGAAGCg         <  1:1196884/93‑1 (MQ=255)
                                                                                 aGAACAGAGGACTGAATATGCAGGAGAAGATGATGAGTTATCAGCCACAAACCGAAGCCGCCACCAGCCGTTTTCTGAATGTAGAAGAAGCgg        >  1:324213/1‑93 (MQ=255)
                                                                                 aGAACAGAGGACTGAATATGCAGGAGAAGATGATGAGTTATCAGCCACAAACCGAAGCCGCCACCAGCCGTTTTCTGAATGTAGAAGAAGCgg        >  1:6371051/1‑93 (MQ=255)
                                                                                    aCAGAGGACTGAATATGCAGGAGAAGATGATGAGTTATCAGCCACAAACCGAAGCCGCCACCAGCCGTTTTCTGAATGTAGAAGAAGCgggca     >  1:3307687/1‑91 (MQ=255)
                                                                                    aCAGAGGACTGAATATGCAGGAGAAGATGATGAGTTATCAGCCACAAACCGAAGCCGCCACCAGCCGTTTTCTGAATGTAGAAGAAGCGGGTa     >  2:4904910/1‑93 (MQ=255)
                                                                                     cAGAGGACTGAATATGCAGGAGAAGATGATGAGTTATCAGCCACAAACCGAAGCCGCCACCAGCCGTTTTCTGAATGTAGAAGAAGCGGGTaa    >  1:6880202/1‑93 (MQ=255)
                                                                                      agagGACTGAATATGCAGGAGAAGATGATGAGTTATCAGCCACAAACCGAAGCCGCCACCAGCCGTTTTCTGAATGTAGAAGAAGCGGGTaaa   <  1:2299917/93‑1 (MQ=255)
                                                                                      agagGACTGAATATGCAGGAGAAGATGATGAGTTATCAGCCACAAACCAAAGCCGCCACCAGCCGTTTTCTGAATGTAGAAGAAGCGGGTaaa   <  1:214307/93‑1 (MQ=255)
                                                                                       gaggACTGAATATGCAGGAAAAGATGATGAGTAATCAGCCACAAACCGAAGCCGCCCCCAGCCGTTTTCTGAATGTAGAAGAAGCGGGTaaaa  <  1:6226197/93‑1 (MQ=255)
                                                                                          |                                                                                         
GCGTTTGGCAACTGGCGTAGCGAAGGGCGTTATTACCGCCCAATCCCGGAGTGGATTGCCGGATTTGGCTCGTTAAGCGCCAGAACAGAGAACTGAATATGCAGGAGAAGATGATGAGTTATCAGCCACAAACCGAAGCCGCCACCAGCCGTTTTCTGAATGTAGAAGAAGCGGGTAAAA  >  NC_000913/370350‑370529

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