Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,308,411 A→G G133G (GGT→GGC b3164 ← polynucleotide phosphorylase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,308,4110AG96.3% 628.3 / 8.2 215G133G (GGT→GGCb3164polynucleotide phosphorylase
Reads supporting (aligned to +/- strand):  ref base A (1/7);  new base G (90/117);  total (91/124)
Fisher's exact test for biased strand distribution p-value = 1.42e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.57e-01

Only 100 of 221 total aligned reads displayed.

CAGTACGTACTGGTCATTGATGTAACCTACGCGGGCAGCACCAATCGGGCCATTGAACGGAATACCAGACAGAGACAGCGCTGCGGAAGCACCAATCATCGCGACGATATCCGGGTTAACTTGCGGGTTAACAGAAACCACGGTGGCGATAACCTGAACTTCGTTGACGAAGCCTTCCGG  >  NC_000913/3308321‑3308500
                                                                                          |                                                                                         
cAGTACGTACTGGTCATTGATGTAACCTACGCGGGCAGCACCAATCGGGCCATTGAACGGAATACCAGACAGAGACAGCGCTGCGGAAGCGcc                                                                                         <  1:583984/93‑1 (MQ=255)
 aGTACGTACTGGTCATTGATGTAACCTACGCGGGCAGCACCAATCGGGCCATTGAACGGAATACCAGACAGAGACAGCGCTGCGGAAGCGCCa                                                                                        >  1:2279933/1‑93 (MQ=255)
 aGTACGTACTGGTCATTGATGTAACCTACGCGGGCAGCACCAATCGGGCCATTGAACGGAATACCAGACAGAGACAGCGCTGCGGAAGCGCCa                                                                                        >  2:3428212/1‑93 (MQ=255)
  gtacgtacTGGTCATTGATGTAACCTACGCGGGCAGCACCAATCGGGCCATTGAACGGAATACCAGACAGAGACAGCGCTGCGGAAGCGCCaa                                                                                       <  2:1498476/93‑1 (MQ=255)
  gtacgtacTGGTCATTGATGTAACCTACGCGGGCAGCACCAATCGGGCCATTGAACGGAATACCAGACAGAGACAGCGCTGCGGAAGCGCCaa                                                                                       <  2:4348602/93‑1 (MQ=255)
  gtacgtacTGGTCATTGATGTAACCTACGCGGGCAGCACCAATCGGGCCATTGAACGGAATACCAGACAGAGACAGCGCTGCGGAAGCGCCaa                                                                                       <  2:2921239/93‑1 (MQ=255)
   tacgtacTGGTCATTGATGTAACCTACGCGGGCAGCACCAATCGGGCCATTGAACGGAATACCAGACAGAGACAGCGCTGCGGAAGCGCCAat                                                                                      <  1:35387/93‑1 (MQ=255)
   tacgtacTGGTCATTGATGTAACCTACGCGGGCAGCACCAATCGGGCCATTGAACGGAATACCAGACAGAGACAGCGCTGCGGAAGCGCCAat                                                                                      >  1:4549031/1‑93 (MQ=255)
    acgtacTGATCATTGATGTAACCTACGCGTGCAGCACCAATCGGGCCATTGAACGGAATACCAGACAGAGACAGCGCTGCGGAAGCGCCAatc                                                                                     <  2:2433426/93‑1 (MQ=255)
      gtacTGGTCATTGATGTAACCTACGCGGGCAGTACCAATCGGGCCATTGAACGGAATACCAGACAGAGACAGCGCTGCGGAAGCGCCAatcat                                                                                   <  1:4430701/93‑1 (MQ=255)
      gtacTGGTCATTGATGTAACCTACGCGGGCAGCACCAATCGGGCCATTGAACGGAATACCAGACAGAGACAGCGCTGCGGAAGCGCCAatcat                                                                                   >  2:300298/1‑93 (MQ=255)
         cTGGTCATTGATGTAACCTACGCGGTCAGCACCAATCGGACCATTGAACGGAATCCCAGACAGAGACAGCGCTGCGGAAGCGCCAATCATcgc                                                                                <  2:943356/93‑1 (MQ=255)
         cTGGTCATTGATGTAACCTACGCGGGCAGCACCAATCGGGCCATTGAACGGAATACCAGACAGAGACAGCGCTGCGGAAGCGCCAATCATcgc                                                                                >  2:2021663/1‑93 (MQ=255)
             tCATTGATGTAACCTACGCGGGCAGCACCAATCGGGCCATTGAACGGAATACCAGACAGAGACAGCGCTGCGGAAGCGCCAATCATCGcgacg                                                                            >  1:1339200/1‑93 (MQ=255)
              cATTTATGTAACCTACGCGGGCAGCACCAATCGGGCCATTGAACGGAATACCAGACAGAGACAGCGCTGCGGAAGCGCCAATCATCGcgacga                                                                           <  2:1911296/93‑1 (MQ=255)
              cATTGATGTAACCTACGCGGGCAGCACCAATCGGGCCATTGAACGGAATACCAGACAGAGACAGCGCTGCGGACGCGCCAATCATCGcgacga                                                                           <  2:4558215/93‑1 (MQ=255)
              cATTGATGTAACCTACGCGGGCAGCACCAATCGGGCCATTGAACGGAATACCAGACAGAGACAGCGCTGCGGAAGCGCCAATCATCGcgacga                                                                           >  2:4006772/1‑93 (MQ=255)
                ttGATGTAACCTACGCGGGCAGCACCAATCGGGCCATTGAACGGAATACCAGACAGAGACAGCGCTGCGGAAGCGCCAATCATCGCGACGata                                                                         >  2:4114183/1‑93 (MQ=255)
                ttGATGTAACCTACGCGGGCAGCACCAATCGGGCCATTGAACGGAATACCAGACAGAGACAGCGCTGCGGAAGCGCCAATCATCGCGACGata                                                                         >  2:859783/1‑93 (MQ=255)
                ttGATGTAACCTACGCGGGCAGCACCAATCGGGCCATTGAACGGAATACCAGACAGAGACAGCGCTGCGGAAGCGCCAATCATCGCGACGata                                                                         >  1:1178076/1‑93 (MQ=255)
                  gATGTAACCTACGCGGGCAGCACCAATCGGGCCATTGAACGGAATACCAGACAGAGACAGCGCTGCGGAAGCGCCAATCATCGCGACGATATc                                                                       <  1:3190355/93‑1 (MQ=255)
                  gATGTAACCTACGCGGGCAGCACCAATCGGGCCATTGAACGGAATACCAGACAGAGACAGCGCTGCGGAAGCGCCAATCATCGCGACGATATc                                                                       <  1:163025/93‑1 (MQ=255)
                  gATGGAACCTACGCGGGCAGCACCAATCGGGCCATTGAACGGAATACCAGACAGAGACAGCGCTGCGGAAGCGCCAATCATCGCGACGATATc                                                                       <  2:3735789/93‑1 (MQ=255)
                   aTGTAACCTACGCGGGCAGCACCAATCGGGCCATTGAACGGAATACCAGACAGAGACAGCGCTGCGGAAGCGCCAATCATCGCGACGATATcc                                                                      <  1:3317515/93‑1 (MQ=255)
                   aTGTAACCTACGCGGGCAGCACCAATCGGGCCATTGAACGGAATACCAGACAGAGACAGCGCTGCGGAAGCGCCAATCATCGCGACGATATcc                                                                      <  1:3841110/93‑1 (MQ=255)
                    tGTAACCTACGCGGGCAGCACCAATCGGGCCATTGAACGGAATACCAGACAGAGACAGCGCTGCGGAAGCGCCAATCATCGCGACGATATCCg                                                                     <  2:2671777/93‑1 (MQ=255)
                    tGTAACCTACGCGGGCAGCACCAATCGGGCCATTGAACGGAATACCAGACAGAGACAGCGCTGCGGAAGCGCCAATCATCGCGACGATATCCg                                                                     >  2:4132074/1‑93 (MQ=255)
                    tGTAACCTACGCGGGCAGCACCAATCGGGCCATTGAACGGAATACCAGACAGAGACAGCGCTGCGGAAGCGCCAATCATCGCGACGATATCCg                                                                     <  2:3542055/93‑1 (MQ=255)
                     gTAACCTACGCGGGCAGCACCAATCGGGCCATTGAACGGAATACCAGACAGAGACAGCGCTGCGGAAGCGCCAATCATCGCGACGATATCCgg                                                                    <  1:3289138/93‑1 (MQ=255)
                      tAACCTACGCGGGCAGCACCAATCGGGCCATTGAACGGAATACCAGACAGAGACAGCGCTGCGGAAGCGCCAATCATCGCGACGATATCCggg                                                                   <  1:2784744/93‑1 (MQ=255)
                        aCCTACGCGGGCAGCACCAATCGGGCCATTGAACGGAATACCAGACAGAGACAGCGCTGCGGAAGCGCCAATCATCGCGACGATATCCGGGtt                                                                 >  1:3089629/1‑93 (MQ=255)
                        aCCTACGCGGGCAGCACCAATCGGGCCATTGAACGGAATACCAGACAGAGACAGCGCTGCGGAAGCGCCAATCATCGCGACGATATCCGGGtt                                                                 >  2:1863598/1‑93 (MQ=255)
                        aCCTACGCGGGCAGCACCAATCGGGCCATTGAACGGAATACCAGACAGAGACAGCGCTGCGGAAGCGCCAATCATCGCGACGATATCCGGGtt                                                                 <  1:2033552/93‑1 (MQ=255)
                         ccTACGCGGGCAGCACCAATCGGGCCATTGAACGGAATACCAGACAGAGACAGCGCTGCGGAAGCGCCAATCATCGCGACGATATCCGGGTTa                                                                >  2:4358295/1‑93 (MQ=255)
                          cTACGCGGGCAGCACCAATCGGGCCATTGAACGGAATACCAGACAGAGACAGCGCTGCGGAAGCGCCAATCATCGCGACGATATCCGGGTTaa                                                               >  1:1514110/1‑93 (MQ=255)
                          cTACGCGGGCAGCACCAATCGGGCCATTGAACGGAATACCAGACAGAGACAGCGCTGCGGAAGCGCCAATCATCGCGACGATATCCGGGTTaa                                                               >  1:2164676/1‑93 (MQ=255)
                          cTACGCGGGCAGCACCAATCGGGCCATTGAACGGAATACCAGACAGAGACAGCGCTGCGGAAGCGCCAATCATCGCGACGATATCCGGGTTaa                                                               >  2:4632864/1‑93 (MQ=255)
                           tACGCGGGCAGCACCAATCGGGCCATTGAACGGAATACCAGACAGAGACAGCGCTGCGGAAGCGCCAATCATCGCGACGATATCCGGGTTAAc                                                              <  2:3903653/93‑1 (MQ=255)
                           tACGCGGCCAGCACCAATCGGGCCATTGAACGGAATACCAGACAGAGACAGCGCTGCGGAAGCGCCAATCATCGCGACGATATCCGGGTTAAc                                                              <  2:3956254/93‑1 (MQ=255)
                             cgcgGGCAGCACCAATCGGGCCATTGAACGGAATACCAGACAGAGACAGCGCTGCGGAAGCGCCAATCATCGCGACGATATCCGGGTTAACtt                                                            <  1:3330069/93‑1 (MQ=255)
                               cgGGCAGCACCAATCGGGCCATTGAACGGAATACCAGACAGAGACAGCGCTGCGGAAGCGCCAATCATCGCGACGATATCCGGGTTAACTTGc                                                          <  1:2703788/93‑1 (MQ=255)
                               cgGGCAGCACCAATCGGGCCATTGAACGGAATACCAGACAGAGACAGCGCTGCGGAAGCGCCAATCATCGCGACGATATCCGGGTTAACTTGc                                                          <  1:1016533/93‑1 (MQ=255)
                                  gcagcaCCAATCGGGCCATTGAACGGAATACCAGACAGAGACAGCGCTGCGGAAGCGCCAATCATCGCGACGATATCCGGGTTAACTTGCggg                                                       <  1:2395771/93‑1 (MQ=255)
                                  gcagcaCCAATCGGGCCATTGAACGGAATACCAGACAGAGACAGCGCTGCGGAAGCACCAATCATCGCGACGATATCCGGGTTAACTTGCggg                                                       <  2:1593097/93‑1 (MQ=255)
                                   cagcaCCAATCGGGCCATTGAACGGAATACCAGACAGAGACAGCGCTGCGGAAGCGCCAATCATCGCGACGATATCCGGGTTAACTTGCGGGt                                                      <  2:1770049/93‑1 (MQ=255)
                                     gcaCCAATCGGGCCATTGAACGGAATACCAGACAGAGACAGCGCTGCGGAAGCGCCAATCATCGCGACGATATCCGGGTTAACTTGCGGGTTa                                                    >  1:2999042/1‑93 (MQ=255)
                                     gcaCCAATCGGGCCATTGAACGGAATACCAGACAGAGACAGCGCTGCGGAAGCGCCAATCATCGCGACGATATCCGGGTTAACTTGCGGGTTa                                                    >  2:3203299/1‑93 (MQ=255)
                                     gcaCCAATCGGGCCATTGAACGGAATACCAGACAGAGACAGCGCTGCGGAAGCGCCAATCATCGCGACGATATCCGGGTTAACTTGCGGGTTa                                                    >  2:546259/1‑93 (MQ=255)
                                     gcaCCAATCGGGCCATTGAACGGAATACCAGACAGAGACAGCGCTACGGAAGCGCCAATCATCGCGACGATATCCGGGTTAACTTGCGGGTTa                                                    >  1:3440955/1‑93 (MQ=255)
                                      caCCAATCGGGCCATTGAACGGAATACCAGACAGAGACAGCGCTGCGGAAGCGCCAATCATCGCGACGATATCCGGGTTAACTTGCGGGTTaa                                                   <  1:1756105/93‑1 (MQ=255)
                                       aCCAATCGGGCCATTGAACGGAATACCAGACAGAGACAGCGCTGCGGAAGCGCCAATCATCGCGACGATATCCGGGTTAACTTGCGGGTTAAc                                                  <  2:1262096/93‑1 (MQ=255)
                                        ccAAGCGGGCCATTGAACGGAATACCAGACAGAGACAGCGCTGCGGAAGCGCCAATCATCGCGACGATATCCGGGTTAACTTGCGGGTTAACa                                                 <  2:1090196/93‑1 (MQ=255)
                                          aaTCGGGCCATTGAACGGAATACCAGACAGAGACAGCGCTGCGGAAGCGCCAATCATCGCGACGATATCCGGGTTAACTTGCGGGTTAACAGa                                               <  1:398318/93‑1 (MQ=255)
                                          aaTCGGGCCATTGAACGGAATACCAGACAGAGACAGCGCTGCGGAAGCGCCAATCATCGCGACGATATCCGGGTTAACTTGCGGGTTAACAGa                                               <  1:3596616/93‑1 (MQ=255)
                                          aaTCGGGCCATTGAACGGAATACCAGACAGAGACAGCGCTGCGGAAGCGCCAATCATCGCGACGATATCCGGGTTAACTTGCGGGTTAACAGa                                               <  2:1202195/93‑1 (MQ=255)
                                          aaTCGGGCCATTGAACGGAATACCAGACAGAGACAGCGCTGCGGAAGCGCCAATCATCGCGACGATATCCGGGTTAACTTGCGGGTTAACAGa                                               <  1:4575171/93‑1 (MQ=255)
                                           aTCGCGCGATTGAACGGAATACCAGACAGAGACAGCGCTGCGGAAGCGCCAATCATCGCGACGATATCCGGGTTAACTTGCGGGTTAACAGaa                                              <  2:2801225/93‑1 (MQ=255)
                                            tCGGGCCATTGAACGGAATACCAGACAGAGACAGCGCTGCGGAAGCGCCAATCATCGCGACGATATCCGGGTTAACTTGCGGGTTAACAGaaa                                             <  1:4094620/93‑1 (MQ=255)
                                            tCGGGCCATTGAACGGAATACCAGACAGAGACAGCGCTGCGGAAGCGCCAATCATCGCGACGATATCCGGGTTAACTTGCGGGTTAACAGaaa                                             <  2:2301868/93‑1 (MQ=255)
                                               ggCCATTGAACGGAATACCAGACAGAGACAGCGCTGCGGAAGCGCCAATCATCGCGACGATATCCGGGTTAACTTGCGGGTTAACAGAAACCa                                          <  1:4260160/93‑1 (MQ=255)
                                               ggCCATTGAACGGAATACCAGACAGAGACAGCGCTGCGGAAGCGCCAATCATCGCGACGATATCCGGGTTAACTTGCGGGTTAACAGAAACCa                                          <  1:3687503/93‑1 (MQ=255)
                                                gCCATTGAACGGAATACCAGACAGAGACAGCGCTGCGGAAGCGCCAATCATCGCGACGATATCCGGGTTAACTTGCGGGTTAACAGAAACCAc                                         <  1:3251329/93‑1 (MQ=255)
                                                 ccATTGAACGGAATACCAGACAGAGACAGCGCTGCGGAAGCGCCAATCATCGCGACGATATCCGGGTTAACTTGCGGGTTAACAGAAACCACg                                        >  1:89960/1‑93 (MQ=255)
                                                  cATTGAACGGAATACCAGACAGAGACAGCGCTGCGGAAGCGCCAATCATCGCGACGATATCCGGGTTAACTTGCGGGTTAACAGAAACCACgg                                       >  2:1315069/1‑93 (MQ=255)
                                                   aTTGAACGGAATACCAGACAGAGACAGCGCTGCGGAAGCGCCAATCATCGCGACGATATCCGGGTTAACTTGCGGGTTAACAGAAACCACGGt                                      <  2:262564/93‑1 (MQ=255)
                                                    ttGAACGGAATACCAGACAGAGACAGCGCTGCGGAAGCGCCAATCATCGCGACGATATCCGGGTTAACTTGCGGGTTAACAGAACCCACGGTg                                     <  2:3632261/93‑1 (MQ=255)
                                                    ttGAACGGAATACCAGACAGAGACAGCGCTGCGGAAGCGCCAATCATCGCGACGATATCCGGGTTAACTTGCGGGTTAACAGAAACCACGGTg                                     <  2:2716313/93‑1 (MQ=255)
                                                           gAATACCAGACAGAGACAGCGCTGCGGAAGCGCCAATCATCGCGACGATATCCGGGTTAACTTGCGGGTTAACAGAAACCACGGTGGCGATaa                              <  1:1955112/93‑1 (MQ=255)
                                                           gAATACCAGACAGAGACAGCGCTGCGGAAGCGCCAATCATCGCGACGATATCCGGGTTAACTTGCGGGTTAACAGAAACCACGGTGGCGATaa                              >  2:3181576/1‑93 (MQ=255)
                                                              tACCAGACAGAGACAGCGCTGCGGAAGCGCCAATCATCGCGACGATATCCGGGTTAACTTGCGGGTTAACAGAAACCACGGTGGCGATAACCt                           >  2:2643777/1‑93 (MQ=255)
                                                               aCCAGACAGAGACAGCGCTGCGGAAGCGCCAATCATCGCGACGATATCCGGGTTAACTTGCGGGTTAACAGAAACCACGGTGGCGATAACCTg                          <  1:287793/93‑1 (MQ=255)
                                                                ccAGACAGAGACAGCGCTGCGGAAGCGCCAATCATCGCGACGATATCCGGGTTAACTTGCGGGTTAACAGAAACCACGGTGGCGATAACCTGa                         <  2:580722/93‑1 (MQ=255)
                                                                  agacagagacagCGCTGCGGAAGCGCCAATCATCGCGACGATATCCGGGTTAACTTGCGGGTTAACAGAAACCACGGTGGCGATAACCTGAAc                       <  2:2063727/93‑1 (MQ=255)
                                                                   gacagagacagCGCTGCGGAAGCGCCAATCATCGCGACGATATCCGGGTTAACTTGCGGGTTAACAGAAACCACGGTGGCGATAACCTGAACt                      <  2:142833/93‑1 (MQ=255)
                                                                    acagagacagCGCTGCGGAAGCGCCAATCATCGCGACGATATCCGGGTTAACTTGCGGGTTAACAGAAACCACGGTGGCGATAACCTGAACtt                     <  1:1966546/93‑1 (MQ=255)
                                                                    acagagacagCGCTGCGGAAGCGCCAATCATCGCGACGATATCCGGGTTAACTTGCGGGTTAACAGAAACCACGGTGGCGATAACCTGAACtt                     <  1:260422/93‑1 (MQ=255)
                                                                    acagagacagCGCTGCGGAAGCGCCAATCATCGCGACGATATCCGGGTTAACTTGCGGGTTAACAGAAACCACGGTGGCGATAACCTGAACtt                     <  2:4063268/93‑1 (MQ=255)
                                                                     cagagacagCGCTGCGGAAGCGCCAATCATCGCGACGATATCCGGGTTAACTTGCGGGTTAACAGAAACCACGGTGGCGATAACCTGAACTTc                    <  1:2897360/93‑1 (MQ=255)
                                                                       gagacagCGCTGCGGAAGCGCCAATCATCGCGACGATATCCGGGTTAACTTGCGGGTTAACAGAAACCACGGTGGCGATAACCTGAACTTCGt                  <  2:2875993/93‑1 (MQ=255)
                                                                       gagacagCGCTGCGGAAGCGCCAATCATCGCGACGATATCCGGGTTAACTTGCGGGTTAACAGAAACCACGGTGGCGATAACCTGAACTTCGt                  <  1:2745943/93‑1 (MQ=255)
                                                                       gagacagCGCTGCGGAAGCGCCAATCATCGCGACGATATCCGGGTTAACTTGCGGGTTAACAGAAACCACGGTGGCGATAACCTGAACTTCGt                  >  1:717445/1‑93 (MQ=255)
                                                                        agacagCGCTGCGGAAGCACCAATCATCGCGACGATATCCGGGTTAACTTGCGGGTTAACAGAAACCACGGTGGCGATAACCTGAACTTCGtt                 <  2:484641/93‑1 (MQ=255)
                                                                         gacagCGCTGCGGAAGCGCCAATCATCGCGACGATATCCGGGTTAACTTGCGGGTTAACAGAAACCACGGTGGCGATAACCTGAACTTCGTTg                <  2:3426022/93‑1 (MQ=255)
                                                                         gacagCGCTGCGGAAGCGCCAATCATCGCGACGATATCCGGGTTAACTTGCGGGTTAACAGAAACCACGGTGGCGATAACCTGAACTTCGTTg                <  2:1058610/93‑1 (MQ=255)
                                                                          acagCGCTGCGGAAGCGCCAATCATCGCGACGATATCCGGGTTAACTTGCGGGTTAACAGAAACCACGGTGGCGATAACCTGAACTTCGTTGa               >  1:617157/1‑93 (MQ=255)
                                                                           cagCGCTGCGGAAGCGCCAATCATCGCGACGATATCCGGGTTAACTTGCGGGTTAACAGAAACCACGGTGGCGATAACCTGAACTTCGTTGAc              <  1:830732/93‑1 (MQ=255)
                                                                             gcgcTGCGGAAGCGCCAATCATCGCGACGATATCCGGGTTAACTTGCGGGTTAACAGAAACCACGGTGGCGATAACCTGAACTTCGTTGACGa            <  2:913373/93‑1 (MQ=255)
                                                                               gcTGCGGAAGCGCCAATCATCGCGACGATATCCGGGTTAACTTGCGGGTTAACAGAAACCACGGTGGCGATAACCTGAACTTCGTTGACGAAg          <  1:3920982/93‑1 (MQ=255)
                                                                                cTGCGGAAGCGCCAATCATCGCGACGATATCCGGGTTAACTTGCGGGTTAACAGAAACCACGGTGGCGATAACCTGAACTTCGTTGACGAAGc         >  2:4261640/1‑93 (MQ=255)
                                                                                cTGCGGAAGCACCAATCATCGCGACGATATCCGGGTTAACTTGCGGGTTAACAGAAACCACGGTGGCGATAACCTGAACTTCGTTGACGAAGc         >  2:237285/1‑93 (MQ=255)
                                                                                 tGCGGAAGCGCCAATCATCGCGACGATATCCGGGTTAACTTGCGGGTTAACAGAAACCACGGTGGCGATAACCTGAACTTCGTTGACGAAGcc        >  2:180794/1‑93 (MQ=255)
                                                                                 tGCGGAAGCGCCAATCATCGCGACGATATCCGGGTTAACTTGCGGGTTAACAGAAACCACGGTGGCGATAACCTGAACTTCGTTGACGAAGcc        >  2:716703/1‑93 (MQ=255)
                                                                                 tGCGGAAGCGCCAATCATCGCGACGATATCCGGGTTAACTTGCGGGTTAACAGAAACCACGGTGGCGATAACCTGAACTTCGTTGACGAAGcc        >  1:3901240/1‑93 (MQ=255)
                                                                                   cGGAAGCGCCAATCATCGCGACGATATCCGGGTTAACTTGCGGGTTAACAGAAACCACGGTGGCGATAACCTGAACTTCGTTGACGAAGCCtt      >  2:1629833/1‑93 (MQ=255)
                                                                                    ggAAGCGCCAATCATCGCGACGATATCCGGGTTAACTTGCGGGTTAACAGAAACCACGGTGGCGATAACCTGAACTTCGTTGACGAAGCCTTc     <  1:2805826/93‑1 (MQ=255)
                                                                                     gAAGCGCCAATCATCGCGACGATATCCGGGTTAACTTGCGGGTTAACAGAAACCACGGTGGCGATCAACTGAACTTCGTTGACGAAGCCTTcc    >  1:3766433/1‑93 (MQ=255)
                                                                                     gAAGCGCCAATCATCGCGACGATATCCGGGTTAACTTGCGGGTTAACAGAAACCACGGTGGCGATAACCTGAACTTCGTTGACGAAGCCTTcc    >  1:1567395/1‑93 (MQ=255)
                                                                                     gAAGCGCCAATCATCGCGACGATATCCGGGTTAACTTGCGGGTTAACAGAAACCACGGTGGCGATAACCTGAACTTCGTTGACGAAGCCTTcc    <  1:4313168/93‑1 (MQ=255)
                                                                                     gAAGCCCCAATCATCGCGACGATATCCGGGTTAACTTGCGGGTTAACAGAAACCACGGTGGCGATAACCTGAACTTCGTTGACGAAGCCTTcc    <  2:1679485/93‑1 (MQ=255)
                                                                                      aaGCGCCAATCATCGCGACGATATCCGGGTTAACTTGCGGGTTAACAGAAACCACGGTGGCGATAACCTGAACTTCGTTGACGAAGCCTTCCg   >  1:237732/1‑93 (MQ=255)
                                                                                          |                                                                                         
CAGTACGTACTGGTCATTGATGTAACCTACGCGGGCAGCACCAATCGGGCCATTGAACGGAATACCAGACAGAGACAGCGCTGCGGAAGCACCAATCATCGCGACGATATCCGGGTTAACTTGCGGGTTAACAGAAACCACGGTGGCGATAACCTGAACTTCGTTGACGAAGCCTTCCGG  >  NC_000913/3308321‑3308500

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