Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 384,546 A→G T50A (ACC→GCC)  tauA → TauA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913384,5460AG98.8% 546.1 / ‑2.8 171T50A (ACC→GCC) tauATauA
Reads supporting (aligned to +/- strand):  ref base A (1/1);  new base G (64/105);  total (65/106)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.06e-01

Only 100 of 181 total aligned reads displayed.

ATGGCAATTTCATCGCGTAACACACTTCTTGCCGCACTGGCATTCATCGCTTTTCAGGCACAGGCGGTGAACGTCACCGTGGCGTATCAAACCTCAGCCGAACCGGCGAAAGTGGCTCAGGCCGACAACACCTTTGCTAAAGAAAGCGGAGCAACCGTGGACTGGCGTAAGTTTGACAGCGGA  >  NC_000913/384456‑384638
                                                                                          |                                                                                            
aTGGCAATTTCATCGCGTAACACACTTCTTGCCGCACTGGCATTCATCGCTTTTCAGGCACAGGCGGTGAACGTCACCGTGGCGTATCAAGcc                                                                                            <  2:2452534/93‑1 (MQ=255)
aTGGCAATTTCATCGCGTAACACACTTCTTGCCGCACTGGCATTCATCGCTTTTCAGGCACAGGCGGTGAACGTCACCGTGGCGTATCAAGcc                                                                                            <  1:3590222/93‑1 (MQ=255)
aTGGCAATTTCATCGCGTAACACACTTCTTGCCGCACTGGCATTCATCGCTTTTCAGGCACAGGCGGTGAACGTCACCGTGGCGTATCAAGcc                                                                                            <  2:1538777/93‑1 (MQ=255)
aTGGCAATTTCATCGCGTAACACACTTCTTGCCGCACTGGCATTCATCGCTTTTCAGGCACAGGCGGTGAACGTCACCGTGGCGTATCAAGcc                                                                                            >  1:1632274/1‑93 (MQ=255)
aTGGCAATTTCATCGCGTAACACACTTCTTGCCGCACTGGCATTCATCGCTTTTCAGGCACAGGCGGTGAACGTCACCGTGGCGTATCAAGcc                                                                                            <  1:1890822/93‑1 (MQ=255)
 tGGCAATTTCATCGCGTAACACACTTCTTGCCGCACTGGCATTCATCGCTTTTCAGGCACAGGCGGTGAACGTCACCGTGGCGTATCAAGCCt                                                                                           >  2:1935984/1‑93 (MQ=255)
 tGGCAATTTCATCGCGTAACACACTTCTTGCCGCACTGGCATTCATCGCTTTTCAGGCACAGGCGGTGAACGTCACCGTGGCGTATCAAGCCt                                                                                           >  2:3120918/1‑93 (MQ=255)
  ggCAATTTCATCGCGTAACACACTTCTTGCCGCACTGGCATTCATCGCTTTTCAGGCACAGGCGGTGAACGTCACCGTGGCGTATCAAGCCTc                                                                                          >  1:1987449/1‑93 (MQ=255)
  ggCAATTTCATCGCGTAACACACTTCTTGCCGCACTGGCATTCATCGCTTTTCAGGCACAGGCGGTGAACGTCACCGTGGCGTATCAAGCCTc                                                                                          >  2:3959342/1‑93 (MQ=255)
   gCAATTTCATCGCGTAACACCCTTCTTGCCGCACTGGCATTCATCGCTTTTCAGGCACAGGCGGTGAACGTCACCCTGGCGTATCCAGCCTCa                                                                                         >  2:3579288/1‑93 (MQ=255)
   gCAATTTCATCGCGTAACACACTTCTTGCCGCACTGGCATTCATCGCTTTTCAGGCACAGGCGGTGAACGTCACCGTGGCGTATCAAGCCTCa                                                                                         >  2:4254175/1‑93 (MQ=255)
   gCAATTTCATCGCGTAACACACTTCTTGCCGCACTGGCATTCATCGCTTTTCAGGCACAGGCGGTGAACGTCACCGTGGCGTATCAAGCCTCa                                                                                         >  2:3579289/1‑93 (MQ=255)
   gCAATTTCATCGCGTAACACACTTCTTGCCGCACTGGCATTCATCGCTTTTCAGGCACAGGCGGTGAACGTCACCGTGGCGTATCAAGCCTCa                                                                                         <  2:249301/93‑1 (MQ=255)
     aaTTTCATCGCGTAACACACTTCTTGCCGCACTGGCATTCATCGCTTTTCAGGCACATGCGGCGAACGTCACCGTGGCGTATCAAGCCTCAGc                                                                                       >  2:2823427/1‑93 (MQ=255)
      aTTTCATCGCGTAACACACTTCTTGCCGCACTGGCATTCATCGCTTTTCAGGCACAGGCGGTGAACGTCACCGTGGCGTATCAAGCCTCAGcc                                                                                      <  1:307042/93‑1 (MQ=255)
       tttCATCGCGTAACACACTTCTTGCCGCACTGGCATTCATCGCTTTTCAGGCACAGGCGGTGAACGTCGCCGTGGCGTATCAAGCCTCAGCCg                                                                                     <  1:3014064/93‑1 (MQ=255)
       tttCATCGCGTAACACACTTCTTGCCGCACTGGCATTCATCGCTTTTCAGGCACAGGCGGTGAACGTCACCGTGGCGTATCAAGCCTCAGCCg                                                                                     <  2:701856/93‑1 (MQ=255)
         tCATCGCGTAACACACTTCTTGCCGCACTGGCATTCATCGCTTTTCAGGCACAGGCGGTGAACGTCACCGTGGCGGATCAAGCCTCAGCCGaa                                                                                   >  1:2237773/1‑93 (MQ=255)
          cATCGCGTAACACACTTCTTGCCGCACTGGCATTCATCGCTTTTCAGGCACAGGCGGTGAACGTCACCGTGGCGTATCAAGCCTCAGCCGAAc                                                                                  <  1:3654883/93‑1 (MQ=255)
           aTCGCGTAACACACTTCTTGCCGCACTGGCATTCATCGCTTTTCAGGCACAGGCGGTGAACGTCACCGTGGCGTATCAAGCCTCAGCCGAAcc                                                                                 <  2:686462/93‑1 (MQ=255)
              gcgTAACACACTTCTTGCCGCACTGGCATTCATCGCTTTTCAGGCACAGGCGGTGAACGTCACCGTGGCGTATCAAGCCTCAGCCGAACCGGc                                                                              <  1:92296/93‑1 (MQ=255)
              gcgTAACACACTTCTTGCCGCACTGGCATTCATCGCTTTTCAGGCACAGGCGGTGAACGTCACCGTGGCGTATCAAGCCTCAGCCGAACCGGc                                                                              >  1:30542/1‑93 (MQ=255)
                gTAACACACTTCTTGCCGCACTGGCATTCATCGCTTTTCAGGCACAGGCGGTGAACGTCACCGTGGCGTATCAAGCCTCAGCCGAACCGGCGa                                                                            >  2:2745485/1‑93 (MQ=255)
                 tAACACACTTCTTGCCGCACTGGCATTCATCGCTTTTCAGGCACAGGCGGTGAACGTCACCGTGGCGTATCAAGCCTCAGCCGAACCGGCGaa                                                                           >  2:1657841/1‑93 (MQ=255)
                  aaCACACTTCTTGCCGCACTGGCATTCATCGCTTTTCAGGCACAGGCGGTGAACGTCACCGTGGCGTATCAAGCCTCAGCCGAACCGGCGaaa                                                                          <  2:1167742/93‑1 (MQ=255)
                   acacacTTCTTGCCGCACTGGCATTCATCGCTTTTCAGGCACAGGCGGTGAACGTCACCGTGGCGTATCAAGCCTCAGCCGAACCGGCGAAAg                                                                         <  2:3335750/93‑1 (MQ=255)
                   acacacTTCTTGCCGCACTGGCATTCATCGCTTTTCAGGCACAGGCGGTGAACGTCACCGTGGCGTATCAAGCCTCAGCCGAACCGGCGAAAg                                                                         <  2:3696596/93‑1 (MQ=255)
                   acacacTTCTTGCCGCACTGGCATTCATCGCTTTTCAGGCACAGGCGGTGAACGTCACCGTGGCGTATCAAGCCTCAGCCGAACCGGCGAAAg                                                                         <  1:3457564/93‑1 (MQ=255)
                   acacacTTCTTGCCGCACTGGCATTCATCGCTTTTCAGGCACAGGCGGTGAACGTCACCGTGGCGTATCAAGCCTCAGCCGAACCGGCGAAAg                                                                         >  2:1028552/1‑93 (MQ=255)
                    cacacTTCTTGCCGCACTGGCATTCATCGCTTTTCAGGCACAGGCGGTGAACGTCACCGTGGCGTATCAAGCCTCAGCCGAACCGGCGAAAGt                                                                        <  1:1013034/93‑1 (MQ=255)
                    cacacTTCTTGCCGCACTGGCATTCATCGCTTTTCAGGCACAGGCGGTGAACGTCACCGTGGCGTATCAAGCCTCAGCCGAACCGGCGAAAGt                                                                        <  1:2276198/93‑1 (MQ=255)
                     ccacTTCTTGCCGCACTGGCATTCATCGCTTTTCAGGCACAGGCGGTGAACGTCACCGTGGCGTATCAAGCCTCAGCCGAACCGGCGAAAGTg                                                                       <  2:3666765/92‑1 (MQ=255)
                      cacTTCTTGCCGCACTGGCATTCATCGCTTTTCAGGCACAGGCGGTGAACGTCACCGTGGCGTATCAAGCCTCAGCCGAACCGGCGAAAGTgg                                                                      <  2:338154/93‑1 (MQ=255)
                       ccttcttGCCGCACTGGCATTCATCGCTTTTCAGGCACAGGCGGTGAACGTCACCGTGGCGTATCAAGCCTCAGCCGAACCGGCGAAAGTGGc                                                                     <  2:3377806/92‑1 (MQ=255)
                       acTTCTTGCCGCACTGGCATTCATCGCTTTTCAGGCACAGGCGGTGAACGTCACCGTGGCGTATCAAGCCTCAGCCGAACCGGCGAAAGTGGc                                                                     <  2:1363685/93‑1 (MQ=255)
                        cttcttGCCGCACTGGCATTCATCGCTTTTCAGGCACAGGCGGTGAACGTCACCGTGGCGTATCAAGCCTCAGCCGAACCGGCGAAAGTGGCt                                                                    <  2:761597/93‑1 (MQ=255)
                         ttcttGCCGCACTGGCATTCATCGCTTTTCAGGCACAGGCGGTGAACGTCACCGTGGCGTATCAAGCCTCAGCCGAACCGGCGAAAGTGGCTc                                                                   >  2:43749/1‑93 (MQ=255)
                         ttcttGCCGCACTGGCATTCATCGCTTTTCAGGCACAGGCGGTGAACGTCACCGTGGCGTATCAAGCCTCAGCCGAACCGGCGAAAGTGGCTc                                                                   >  1:1292104/1‑93 (MQ=255)
                         ttcttGCCGCACGGGTATTCATCGCCTTTCAGGCACAGGCGTTGAACTTCACGGTGCCGTATCAAGCCTCAGCCGAACCGGCGAAAGTGGCTc                                                                   <  2:923279/93‑1 (MQ=255)
                          tcttGCCGCACTGGCATTCATCGCTTTTCAGGCACAGGCGGTGAACGTCACCGTGGCGTATCAAGCCTCAGCCGAACCGGCGAAAGTGGCTCa                                                                  >  1:2891828/1‑93 (MQ=255)
                          tcttGCCGCACTGGCATTCATCGCTTTTCAGGCACAGGCGGTGAACGTCACCGTGGCGTATCAAACCTCAGCCGAACCGGCGAAAGTGGCTCa                                                                  <  1:4267121/93‑1 (MQ=255)
                          tcttGCCGCAATGGCATTCATCGCTTTTCAGGCACAGGCGGTGAACGTCACCGTGGCGTATCAAGCCTCAGCCGAACCGGCGAAAGTGGCTCa                                                                  <  2:686813/93‑1 (MQ=255)
                           cttGCCGCACTGGCATTCATCGCTTTTCAGGCACAGGCGGTGAACGTCACCGTGGCGTATCAAGCCTCAGCCGAACCGGCGAAAGTGGCTCAg                                                                 <  2:2983351/93‑1 (MQ=255)
                           cttGCCGCACTGGCATTCATCGCTTTTCAGGCACAGGCGGTGAACGTCACCGTGGCGTATCAAGCCTCAGCCGAACCGGCGAAAGTGGCTCAg                                                                 <  1:2640278/93‑1 (MQ=255)
                             tGCCGCACTGGCATTCATCGCTTTTCAGGCACAGGCGGTGAACGTCACCGTGGCGTATCAAGCCTCAGCCGAACCGGCGAAAGTGGCTCAGGc                                                               <  2:3067408/93‑1 (MQ=255)
                             tGCCGCACTGGCATTCATCGCTTTTCAGGCACAGGCGGTGAACGTCACCGTGGCGTATCAAGCCTCAGCCGAACCGGCGAAAGTGGCTCAGGc                                                               <  2:3718630/93‑1 (MQ=255)
                             tGCCGCACTGGCATTCATCGCTTTTCAGGCACAGGCGGTGAACGTCACCGTGGCGTATCAAGCCTCAGCCGAACCGGCGAAAGTGGCTCAGGc                                                               <  1:3342712/93‑1 (MQ=255)
                              gCCGCACTGGCATTCATCGCTTTTCAGGCACAGGCGGTGAACGTCACCGTGGCGTATCAAGCCTCAGCCGAACCGGCGAAAGTGGCTCAGGcc                                                              <  2:2535147/93‑1 (MQ=255)
                               ccGCACTGGCATTCATCGCTTTTCAGGCACAGGCGGTGAACGTCACCGTGGCGTATCAAGCCTCAGCCGAACCGGCGAAAGTAGCTCACGCCg                                                             >  2:4545726/1‑93 (MQ=255)
                                cGCACTGGCATTCATCGCTTTTCAGGCACAGGCGGTGAACGTCACCGTGGCGTATCAAGCCTCAGCCGAACCGGCGAAAGTGGCTCAGGCCGa                                                            >  1:2330070/1‑93 (MQ=255)
                                 gCACTGGCATTCATCGCTTTTCAGGCACAGGCGGTGAACGTCACCGTGGCGTATCAAGCCTCAGCCGAACCGGCGAAAGTGGCTCAGGCCGac                                                           <  1:1254829/93‑1 (MQ=255)
                                 gCACTGGCATTCATCGCTTTTCAGGCACAGGCGGTGAACGTCACCGTGGCGTATCAAGCCTCAGCCGAACCGGCGAAAGTGGCTCAGGCCGac                                                           >  1:1433080/1‑93 (MQ=255)
                                 gCACTGGCATTCATCGCTTTTCAGGCACAGGCGGTGAACGTCACCGTGGCGTATCAAGCCTCAGCCGAACCGGCGAAAGTGGCTCAGGCCGac                                                           <  2:124408/93‑1 (MQ=255)
                                  cACTGGCATTCATCGCTTTTCAGGCACAGGCGGTGAACGTCACCGTGGCGTATCAAGCCTCAGCCGAACCGGCGAAAGTGGCTCAGGCCGaca                                                          <  1:116782/93‑1 (MQ=255)
                                  cACTGGCATTCATCGCTTTTCAGGCACAGGCGGTGAACGTCACCGTGGCGTATCAAGCCTCAGCCGAACCGGCGAAAGTGGCTCAGGCCGaca                                                          <  1:3600747/93‑1 (MQ=255)
                                   aCTGGCATTCATCGCTTTTCAGGCACAGGCGGTGAACGTCACCGTGGCGTATCAAGCCTCATCCGAACCGGCGAAAGTGGCTCAGGCCGacaa                                                         >  1:2339576/1‑93 (MQ=255)
                                   aCTGGCATTCATCGCTTTTCAGGCACAGGCGGTGAACGTCACCGTGGCGTATCAAGCATCAGCCGAACCGGCGACAGTGGCTCAGGCCGaacc                                                         >  2:8756/1‑90 (MQ=255)
                                    cTGGCATTCATCGCTTTTCAGGCACAGGCGGTGAACGTCACCGTGGCGTATCAAGCCTCAGCCGAACCGGCGAAAGTGGCTCAGGCCGacaac                                                        >  1:672683/1‑93 (MQ=255)
                                     tGGCATTCATCGCTTTTCAGGCACAGGCGGTGAACGTCACCGTGGCGTATCAAGCCTCAGCCGAACCGGCGAAAGTGGCTCAGGCCGacaaca                                                       >  2:4293064/1‑93 (MQ=255)
                                     tGGCATTCATCGCTTTTCAGGCACAGGCGGTGAACGTCACCGTGGCGTATCAAGCCTCAGCCGAACCGGCGAAAGTGGCTCAGGCCGacaaca                                                       <  1:4637183/93‑1 (MQ=255)
                                     tGGCATTCATCGCTTTTCAGGCACAGGCGGTGAACGTCACCGTGGCGTATCAAGCCTCAGCCGAACCGGCGAAAGTGGCTCAGGCCGacaaca                                                       >  2:1708141/1‑93 (MQ=255)
                                      ggCATTCATCGCTTTTCAGGCACAGGCGGTGAACGTCACCGTGGCGTATCAAGCCTCAGCCGAACCGGCGAAAGTGGCTCAGGCCGACAacac                                                      <  1:1841242/93‑1 (MQ=255)
                                         aTTCATCGCTTTTCAGGCACAGGCGGTGAACGTCACCGTGGCGTATCAAGCCTCAGCCGAACCGGCGAAAGTGGCTCAGGCCGACAACACCtt                                                   >  1:3548060/1‑93 (MQ=255)
                                          ttCATCGCTTTTCAGGCACAGGCGGTGAACGTCACCGTGGCGTATCAAGCCTCAGCCGAACCGGCGAAAGTGGCTCAGGCCGACAACACCttt                                                  <  1:671230/93‑1 (MQ=255)
                                          ttCATCGCTTTTCAGGCACAGGCGGTGAACGTCACCGTGGCGTATCAAGCCTCAGCCGAACCGGCGAAAGTGGCTCAGGCCGACAACACCttt                                                  >  1:4008265/1‑93 (MQ=255)
                                             ctCGCTTTTCAGGCGCAGGCGGTGAACGTCACCGTGGCGTATCAAGCCTCAGCCGAACCGGCGAAAGTGGCTCAGGCCGACAACCCCTTTGCt                                               <  1:3608901/92‑1 (MQ=255)
                                             aTCGCTTTTCAGGCACAGGCGGTGAACGTCACCGTGGCGTATCAAGCCTCAGCCGAACCGGCGAAAGTGGCTCAGGCCGACAACACCTTTGCt                                               <  2:1439943/93‑1 (MQ=255)
                                             aTCGCTTTTCAGGCACAGGCGGTGAACGTCACCGTGGCGTATCAAGCCTCAGCCGAACCGGCGAAAGTGGCTCAGGCCGACAACACCTTTGCt                                               <  2:412673/93‑1 (MQ=255)
                                             aTCGCTTTTCAGGCACAGGCGGTGAACGTCACCGTGGCGTATCAAGCCTCAGCCGAACCGGCGAAAGTGGCTCAGGCCGACAACACCTTTGCt                                               <  1:1457935/93‑1 (MQ=255)
                                              tCGCTTTTCAGGCACAGGCGGTGAACGTCACCGTGGCATATCAAGCCTCAGCCGAACCGGCGAAAGTGGCTCAGGCCGACAACACCTTTGCTa                                              <  1:4077092/93‑1 (MQ=255)
                                                gCTTTTCAGGCACAGGCGGTGAACGTCACCGTGGCGTATCAAGCCTCAGCCGAACCGGCGAAAGTGGCTCAGGCCGACAACACCTTTGCTaaa                                            <  2:2125174/93‑1 (MQ=255)
                                                gCTTTTCAGGCACAGGCGGTGAACGTCACCGTGGCGTATCAAGCCTCAGCCGAACCGGCGAAAGTGGCTCAGGCCGACAACACCTTTGCTaaa                                            <  1:3410274/93‑1 (MQ=255)
                                                  ttttCAGGCACAGGCGGTGAACGTCACCGTGGCGTATCAAGCCTCAGCCGAACCGGCGAAAGTGGCTCAGGCCGACAACACCTTTGCTaaaga                                          <  2:533541/93‑1 (MQ=255)
                                                  ttttCAGGCACAGGCGGTGAACGTCACCGTGGCGTATCAAGCCTCAGCCGAACCGGCGAAAGTGGCTCAGGCCGACAACACCTTTGCTaaaga                                          <  2:2133849/93‑1 (MQ=255)
                                                  ttttCAGGCACAGGCGGTGAACGTCACCGTGGCGTATCAAGCCTCAGCCGAACCGGCGAAAGTGGCTCAGGCCGACAACACCTTTGCTaaaga                                          <  2:2784329/93‑1 (MQ=255)
                                                  ttttCAGGCACAGGCGGTGAACGTCACCGTGGCGTATCAAGCCTCAGCCGAACCGGCGAAAGTGGCTCAGGCCGACAACACCTTTGCTaaaga                                          >  1:3522582/1‑93 (MQ=255)
                                                  ttttCAGGCACAGGCGGTGAACGTCACCGTGGCGTATCAAGCCTCAGCCGAACCGGCGAAAGTGGCTCAGGCCGACAACACCTTTGCTaaaga                                          <  1:202361/93‑1 (MQ=255)
                                                  ttttCAGGCACAGGCGGTGAACGTCACCGTGGCGTATCAAGCCTCAGCCGAACCGGCGAAAGTGGCTCAGGCCGACAACACCTTTGCTaaaga                                          <  1:2504167/93‑1 (MQ=255)
                                                     tCAGGCACAGGCGGTGAACGTCACCGTGGCGTATCAAGCCTCAGCCGAACCGGCGAAAGTGGCTCAGGCCGACAACACCTTTGCTaaagaaag                                       <  2:2273570/93‑1 (MQ=255)
                                                      cAGGCACAGGCGGTGAACGTCACCGTGGCGTATCAAGCCTCAGCCGAACCGGCGAAAGTGGCTCAGGCCGACAACACCTTTGCTAAAGAAAGc                                      <  2:2414628/93‑1 (MQ=255)
                                                        ggCACAGGCGGTGAACGTCACCGTGGCGTATCAAGCCTCAGCTGAACCGGCGAAAGTGGCTCAGGCCGACAACACCTTTGCTAAAGAAAGCgg                                    >  1:129680/1‑93 (MQ=255)
                                                          cacaGGCGGTGAACGTCACCGTGGCGTATCAAGCCTCAGCCGAACCGGCGAAAGTGGCTCAGGCCGACAACACCTTTGCTAAAGAAAGCGGAg                                  <  2:4400903/93‑1 (MQ=255)
                                                              ggCGGTGAACGTCACCGTGGCGTATCAAGCCTCAGCCGAACCGGCGAAAGTGGCTCAGGCCGACAACACCTTTGCTAAAGAAAGCGGAGCAAc                              <  1:677538/93‑1 (MQ=255)
                                                              cgCGGTGCACGTCACCGTGGCGTATCAAGCCTCAGCCGAACCGGCGAAAGTGGCTCAGGCCGACAACACCTTTGCTAAAGAAAGCGGAGCAAc                              >  2:2881953/2‑93 (MQ=255)
                                                               gCGGTGAACGTCACCGTGGCGTATCAAGCCTCAGCCGAACCGGCGAAAGTGGCTCAGGCCGACAACACCTTTGCTAAAGAAAGCGGAGCAAcc                             <  2:2613926/93‑1 (MQ=255)
                                                                  gTGAACGTCACCGTGGCGTATCAAGCCTCAGCCGAACCGGCGAAAGTGGCTCAGGCCGACAACACCTTTGCTAAAGAAAGCGGAGCAACCGTg                          >  1:3383876/1‑93 (MQ=255)
                                                                   tGAACGTCACCGTGGCGTATCAAGCCTCAGCCGAACCGGCGAAAGTGGCTCAGGCCGACAACACCTTTGCTAAAGAAAGCGGAGCAACCGTgg                         <  2:1535872/93‑1 (MQ=255)
                                                                      aCGTCACCGTGGCGTATCAAGCCTCAGCCGAACCGACGAAAGTGGCTCAGGCCGACAACACCTTTGCTAAAGAAAGCGGAGCAACCGTGGACt                      >  2:2372545/1‑93 (MQ=255)
                                                                        gTCACCGTGGCGTATCAAGCCTCAGCCGAACCGGCGAAAGTGGCTCAGGCCGACAACACCTTTGCTAAAGAAAGCGGAGCAACCGTGGACTgg                    <  1:3679118/93‑1 (MQ=255)
                                                                        gTCACCGTGGCGTATCAAGCCTCAGCCGAACCGGCGAAAGTGGCTCAGGCCGACAACACCTTTGCTAAAGAAAGCGGAGCAACCGTGGACTgg                    >  2:4012742/1‑93 (MQ=255)
                                                                        gTCACCGTGGCGTATCAAGCCTCAGCCGAACCGGCGAAAGTGGCTCAGGCCGACAACACCTTTGCTAAAGAAAGCGGAGCAACCGTGGACTgg                    <  2:4418527/93‑1 (MQ=255)
                                                                         tCACCGTGGCGTATCAAGCCTCAGCCGAACCGGCGAAAGTGGCTCAGGCCGACAACACCTTTGCTAAAGAAAGCGGAGCAACCGTGGACTGGc                   >  2:3227/1‑93 (MQ=255)
                                                                           cccGTGGCGTATCAAGCCTCAGCCGAACCGGCGAAAGTGGCTCAGGCCGACAACACCTTTGCTAAAGAAAGCGGAGCAACCGTGGACTGGCGt                 <  1:2110084/92‑1 (MQ=255)
                                                                            ccGTGGCGTATCAAGCCTCAGCCGAACCGGCGAAAGTGGCTCAGGCCGACAACACCTTTGCTAAAGAAAGCGGAGCAACCGTGGACTGGCGTa                <  2:3257863/93‑1 (MQ=255)
                                                                            ccGTGGCGTATCAAGCCTCAGCCGAACCGGCGAAAGTGGCTCAGGCCGACAACACCTTTGCTAAAGAAAGCGGAGCAACCGTGGACTGGCGTa                <  2:3833100/93‑1 (MQ=255)
                                                                             cGTGGCGTATCAAGCCTCAGCCGAACCGGCGAAAGTGGCTCAGGCCGACAACAACTTTGCTAAAGAAAGCGCAGCAACCGTAGACCGGCGTaa               >  2:1264094/1‑93 (MQ=255)
                                                                                      tCAAGCCTCAGCCGAACCGGCGAAAGTGGCTCAGGCCGACAACACCTTTGCTAAAGAAAGCGGAGCAACCGTGGACTGGCGTAAGTTTGACAg      <  1:3289301/93‑1 (MQ=255)
                                                                                       cAAGCCTCAGCCGAACCGGCGAAAGTGGCTCAGGCCGACAACACCTTTGCTAAAGAAAGCGGAGCAACCGTGGACTGGCGTAAGTTTGACAGc     <  1:465941/93‑1 (MQ=255)
                                                                                       cAAGCCTCAGCCGAACCGGCGAAAGTGGCTCAGGCCGACAACACCTTTGCTAAAGAAAGCGGAGCAACCGTGGACTGGCGTAAGTTTGACAGc     >  1:4324384/1‑93 (MQ=255)
                                                                                          aCCTCAGCCGAACCGGCGAAAGTGGCTCAGGCCGACAACACCTTTGCTAAAGAAAGCGGAGCAACCGTGGACTGGCGTAAGTTTGACAGCGGa  <  1:3873174/93‑1 (MQ=255)
                                                                                          |                                                                                            
ATGGCAATTTCATCGCGTAACACACTTCTTGCCGCACTGGCATTCATCGCTTTTCAGGCACAGGCGGTGAACGTCACCGTGGCGTATCAAACCTCAGCCGAACCGGCGAAAGTGGCTCAGGCCGACAACACCTTTGCTAAAGAAAGCGGAGCAACCGTGGACTGGCGTAAGTTTGACAGCGGA  >  NC_000913/384456‑384638

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: ‑