Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,448,298 T→C D359G (GAT→GGT)  tynA ← copper amine oxidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,448,2980TC99.5% 601.0 / ‑4.9 193D359G (GAT→GGT) tynAcopper amine oxidase
Reads supporting (aligned to +/- strand):  ref base T (0/1);  new base C (95/96);  total (96/97)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.05e-01

Only 100 of 205 total aligned reads displayed.

TCGTACATGACTTTGCGTTTGGTGCCATTGTCGTTATAAGTCACGGTGGAGATCATCGGCCCGACGCGAGAGTTCATGCTGAGGTGAAAATCCCAGTTCCGCCAGTGAATCATATCGCCAGTAATGGTGTAATTTTTACCTTCAGGCTCAATGATTTGCATAGGCTTAACTGCCGGAGCA  >  NC_000913/1448208‑1448387
                                                                                          |                                                                                         
tCGTACATGACTTTGCGTTTGGTGCCATTGTCGTTATAAGTCACGGTGGAGATCATCGGCCCGACGCGAGAGTTCATGCTGAGGTGAAAAccc                                                                                         <  1:3847600/93‑1 (MQ=255)
 cGTACATGACTTTGCGTTTGGTGCCATTGTCGTTATAAGTCACGGTGGAGATCATCGGCCCGACGCGAGAGTTCATGCTGAGGTGAAAAcccc                                                                                        <  2:5513624/93‑1 (MQ=255)
   tACATGACTTTGCGTTTGGTGCCATTGTCGTTATAAGTCACGGTGGAGATCATCGGCCCGACGCGAGAGTTCATGCTGAGGTGAAAACCCCAg                                                                                      <  2:5352230/93‑1 (MQ=255)
   tACATGACTTTGCGTTTGGTGCCATTGTCGTTATAAGTCACGGTGGAGATCATCGGCCCGACGCGAGAGTTCATGCTGAGGTGAAAACCCCAg                                                                                      <  1:534613/93‑1 (MQ=255)
    aCATGACTTTGCGTTTGGTGCCATTGTCGTTATAAGTCACGGTGGAGATCATCGGCCCGACGCGAGAGTTCATGCTGAGGTGAAAACCCCAGt                                                                                     >  1:5133255/1‑93 (MQ=255)
     cATGACTTTGCGTTTGGTGCCATTGTCGTTATAAGTCACGGTGGAGATCATCGGCCCGACGCGAGAGTTCATGCTGAGGTGAAAACCCCAGtt                                                                                    >  2:4701993/1‑93 (MQ=255)
         aCTTTGCGTTTGGTGCCATTGTCGTTATAAGTCACGGTGGAGATCATCGGCCCGACGCGAGAGTTCATGCTGAGGTGAAAACCCCAGTTCCGc                                                                                <  2:371787/93‑1 (MQ=255)
           tttGCGTTTGGTGCCATTGTCGTTATAAGTCACGGTGGAGATCATCGGCCCGACGCGAGAGTTCATGCTGAGGTGAAAACCCCAGTTCCGCCa                                                                              <  1:5456526/93‑1 (MQ=255)
            ttGCGTTTGGTGCCATTGTCGTTATAAGTCACGGTGGAGATCATCGGCCCGACGCGAGAGTTCATGCTGAGGTGAAAACCCCAGTTCCGCCAg                                                                             <  2:4092752/93‑1 (MQ=255)
              gCGTTTGGTGCCATTGTCGTTATAAGTCACGGTGGAGATCATCGGCCCGACGCGAGAGTTCATGCTGAGGTGAAAACCCCAGTTCCGCCAGTg                                                                           <  2:2443985/93‑1 (MQ=255)
               ggTTTGGTGCCATTGTCGGTATAAGGCACGGTGGAGATCATCGGCCCGACGCGAGAGTTCATGCTGAGGTGAAAACCCCAGTTCCGCCAGTGa                                                                          <  2:1417383/92‑1 (MQ=255)
               cGTTTGGTGCCATTGTCGTTATAAGTCACGGTGGAGATCATCGGCCCGACGCGAGAGTTCATGCTGAGGTGAAAACCCCAGTTCCGCCAGTGa                                                                          >  2:55751/1‑93 (MQ=255)
                gTTTGGTGCCATTGTCGTTATAAGTCACGGTGGAGATCATCGGCCCGACGCGAGAGTTCATGCTGAGGTGAAAACCCCAGTTCCGCCAGTGaa                                                                         <  1:7591771/93‑1 (MQ=255)
                gTTTGGTGCCATTGTCGTTATAAGTCACGGTGGAGATCATCGGCCCGACGCGAGAGTTCATGCTGAGGTGAAAACCCCAGTTCCGCCAGTGaa                                                                         >  1:3586867/1‑93 (MQ=255)
                 tttGGTGCCATTGTCGTTATAAGTCACGGTGGAGATCATCGGCCCGACGCGAGAGTTCATGCTGAGGTGAAAACCCCAGTTCCGCCAGTGAAt                                                                        <  1:4579630/93‑1 (MQ=255)
                 tttGGTGCCATTGTCGTTATAAGTCACGGTGGAGATCATCGGCCCGACGCGAGAGTTCATGCTGAGGTGAAAACCCCAGTTCCGCCAGTGAAt                                                                        >  2:1912960/1‑93 (MQ=255)
                 tttGGTGCCATTGTCGTTATAAGTCACGGTGGAGATCATCGGCCCGACGCGAGAGTTCATGCTGAGGTGAAAACCCCAGTTCCGCCAGTGAAt                                                                        >  1:3733925/1‑93 (MQ=255)
                  ttGGTGCCATTGTCGTTATAAGTCACGGTGGAGATCATCGGCCCGACGCGAGAGTTCATGCTGAGGTGAAAACCCCAGTTCCGCCAGTGAAta                                                                       >  1:1530927/1‑92 (MQ=255)
                   tGGTGCCATTGTCGTTATAAGTCACGGTGGAGATCATCGGCCCGACGCGAGAGTTCATGCTGAGGTGAAAACCCCAGTTCCGCCAGTGAATCa                                                                      <  2:5766635/93‑1 (MQ=255)
                    ggTGCCATTGTCGTTATAAGTCACGGTGGAGATCATCGGCCCGACGCGAGAGTTCATGCTGAGGTGAAAACCCCAGTTCCGCCAGTGAATCat                                                                     >  2:7857003/1‑93 (MQ=255)
                     gTGCCATTGTCGTTATAAGTCACGGTGGAGATCATCGGCCCGACGCGAGAGTTCATGCTGAGGTGAAAACCCCAGTTCCGCCAGTGAATCata                                                                    <  2:423313/93‑1 (MQ=255)
                     gTGCCATTGTCGTTATAAGTCACGGTGGAGATCATCGGCCCGACGCGAGAGTTCATGCTGAGGTGAAAACCCCAGTTCCGCCAGTGAATCata                                                                    <  1:704376/93‑1 (MQ=255)
                      tGCCATTGTCGTTATAAGTCACGGTGGAGATCATCGGCCCGACGCGAGAGTTCATGCTGAGGTGAAAACCCCAGTTCCGCCAGTGAATCatat                                                                   <  1:7137280/93‑1 (MQ=255)
                      tGCCATTGTCGTTATAAGTCACGGTGGAGATCATCGGCCCGACGCGAGAGTTCATGCTGAGGTGAAAACCCCAGTTCCGCCAGTGAATCatat                                                                   >  2:3969109/1‑93 (MQ=255)
                       gCCATTGTCGTTATAAGTCACGGTGGAGATCATCGGCCCGACGCGAGAGTTCATGCTGAGGTGAAAACCCCAGTTCCGCCAGTGAATCATATc                                                                  <  1:5841049/93‑1 (MQ=255)
                       gCCATTGTCGTTATAAGTCACGGTGGAGATCATCGGCCCGACGCGAGAGTTCATGCTGAGGTGAAAACCCCAGTTCCGCCAGTGAATCATATc                                                                  <  1:4500487/93‑1 (MQ=255)
                        ccATTGTCGTTATAAGTCACGGTGGAGATCATCGGCCCGACGCGAGAGTTCATGCTGAGGTGAAAACCCCAGTTCCGCCAGTGAATCATATCg                                                                 <  2:4678905/93‑1 (MQ=255)
                         cATTGTCGTTATAAGTCACGGTGGAGATCATCGGCCCGACGCGAGAGTTCATGCTGAGGTGAAAACCCCAGTTCCGCCAGTGAATCATATCGc                                                                <  1:6756155/93‑1 (MQ=255)
                         cATTGTCGTTATAAGTCACGGTGGAGATCATCGGCCCGACGCGAGAGTTCATGCTGAGGTGAAAACCCCAGTTCCGCCAGTGAATCATATCGc                                                                >  1:297715/1‑93 (MQ=255)
                          aTTGTCGTTATAAGTCACGGTGGAGATCATCGGCCCGACGCGAGAGTTCATGCTGAGGTGAAAACCCCAGTTCCGCCAGTGAATCATATCGcc                                                               >  2:2608248/1‑93 (MQ=255)
                           ttGTCGTTATAAGTCACGGTGGAGATCATCGGCCCGACGCGAGAGTTCATGCTGAGGTGAAAACCCCAGTTCCGCCAGTGAATCATATCGCCa                                                              <  1:2009138/93‑1 (MQ=255)
                           ttGTCGTTATAAGTCACGGTGGAGATCATCGGCCCGACGCGAGAGTTCATGCTGAGGTGAAAACCCCAGTTCCGCCAGTGAATCATATCGCCa                                                              >  1:3466165/1‑93 (MQ=255)
                              tCGTTATAAGTCACGGTGGAGATCATCGGCCCGACGCGAGAGTTCATGCTGAGGTGAAAACCCCAGTTCCGCCAGTGAATCATATCGCCAGTa                                                           <  2:2040670/93‑1 (MQ=255)
                              tCGTTATAAGTCACGGTGGAGATCATCGGCCCGACGCGAGAGTTCATGCTGAGGTGAAAACCCCAGTTCCGCCAGTGAATCATATCGCCAGTa                                                           >  2:3983144/1‑93 (MQ=255)
                               cGTTATAAGTCACGGTGGAGATCATCGGCCCGACGCGAGAGTTCATGCTGAGGTGAAAACCCCAGTTCCGCCAGTGAATCATATCGCCAGTaa                                                          <  2:3484284/93‑1 (MQ=255)
                                gTTATAAGTCACGGTGGTGATCATCGGCCCGACGCGTGAGTTCATGCTGAGGCGTAAACCTCAGTTTCGCCAGTGAATCATATCGCCAGTGAt                                                         >  2:1030647/1‑93 (MQ=255)
                                gTTATAAGTCACGGTGGAGATCATCGGCCCGACGCGAGAGTTCATGCTGAGGTGAAAACCCCAGTTCCGCCAGTGAATCATATCGCCAGTAAt                                                         >  2:2453792/1‑93 (MQ=255)
                                gTTATAAGTCACGGTGGAGATCATCGGCCCGACGCGAGAGTTCATGCTGAGGTGAAAACCCCAGTTCCGCCAGTGAATCATATCGCCAGTAAt                                                         <  1:1032401/93‑1 (MQ=255)
                                gTTATAAGTCACGGTGGAGATCATCGGCCCGACGCGAGAGTTCATGCTGAGGTGAAAACCCCAGTTCCGCCAGTGAATCATATCGCCAGTAAt                                                         >  2:6721737/1‑93 (MQ=255)
                                 ttATAAGTCACGGTGGAGATCATCGGCCCGACGCGAGAGTTCATGCTGAGGTGAAAACCCCAGTTCCGCCAGTGAATCATATCGCCAGTAATg                                                        <  1:2096923/93‑1 (MQ=255)
                                 ttATAAGTCACGGTGGAGATCATCGGCCCGACGCGAGAGTTCATGCTGAGGTGAAAACCCCAGTTCCGCCAGTGAATCATATCGCCACTAATg                                                        >  1:4225309/1‑93 (MQ=255)
                                 ttATAAGTCACGGTGGAGATCATCGGCCCGACGAGAGAGTTCATGCTGAGGTGAAAACCCCAGTTCCGCCAGTGAATCATATCGCCAGTAATg                                                        <  1:2165462/93‑1 (MQ=255)
                                  tataAGTCACGGTGGAGATCATCGGCCCGACGCGAGAGTTCATGCTGAGGTGAAAACCCCAGTTCCGCCAGTGAATCATATCGCCAGTAATgg                                                       <  2:1713650/93‑1 (MQ=255)
                                   ataAGTCACGGTGGAGATCATCGGCCCGACGCGAGAGTTCATGCTGAGGTGAAAACCCCAGTTCCGCCAGTGAATCATATCGCCAGTAATGgt                                                      <  1:3293956/93‑1 (MQ=255)
                                    taAGTCACGGTGGAGATCATCGGCCCGACGCGAGCGTTCATGCTGAGGGGAAAAATCCCGTTCCACCAGTGAATAATATAGCCAGAAATGgta                                                     >  2:5967465/1‑92 (MQ=255)
                                    taAGTCACGGTGGAGATCATCGGCCCGACGCGAGAGTTCATGCTGAGGTGAAAACCCCAGTTCCGCCAGTGAATCATATCGCCAGTAATGgtg                                                     >  2:4859457/1‑93 (MQ=255)
                                    taAGTCACGGTGGAGATCATCGGCCCGACGCGAGAGTTCATGCTGAGGTGAAAACCCCAGTTCCGCCAGTGAATCATATCGCCAGTAATGgtg                                                     >  2:4326091/1‑93 (MQ=255)
                                    taAGTCACGGTGGAGATCATCGGCCCGACGCGAGAGTTCATGCTGAGGTGAAAACCCCAGTTCCGCCAGTGAATCATATCGCCAGTAATGgtg                                                     <  1:2508396/93‑1 (MQ=255)
                                       gTCACGGTGGAGATCATCGGCCCGACGCGAGAGTTCATGCTGAGGTGAAAACCCCAGTTCCGCCAGTGAATCATATCGCCAGTAATGGTGTaa                                                  >  1:2570412/1‑93 (MQ=255)
                                            ggTGGAGATCATCGGCCCGACGCGAGAGTTCATGCTGAGGTGAAAACCCCAGTTCCGCCAGTGAATCATATCGCCAGTAATGGTGTAAttttt                                             >  1:5256958/1‑93 (MQ=255)
                                             gTGGAGATCATCGGCCCGACGCGAGAGTTCATGCTGAGGTGAAAACCCCAGTTCCGCCAGTGAATCATATCGCCAGTAATGGTGTAATTTTTa                                            <  1:6232844/93‑1 (MQ=255)
                                             gTGGAGATCATCGGCCCGACGCGAGAGTTCATGCTGAGGTGAAAACCCCAGTTCCGCCAGTGAATCATATCGCCAGTAATGGTGTAATTTTTa                                            <  2:2739950/93‑1 (MQ=255)
                                             gTGGAGATCATCGGCCCGACGCGAGAGTTCATGCTGAGGTGAAAACCCCAGTTCCGCCAGTGAATCATATCGCCAGTAATGGTGTAATTTTTa                                            >  1:2605446/1‑93 (MQ=255)
                                                gagaTCATCGGCCCGACGCGAGAGTTCATGCTGAGGTGAAAACCCCAGTTCCGCCAGTGAATCATATCGCCAGTAATGGTGTAATTTTTACCt                                         >  2:3174997/1‑93 (MQ=255)
                                                gagaTCATCGGCCCGACGCGAGAGTTCATGCTGAGGTGAAAACCCCAGTTCCGCCAGTGAATCATATCGCCAGTAATGGTGTAATTTTTACCt                                         <  1:1758499/93‑1 (MQ=255)
                                                gagaTCATCGGCCCGACGCGAGAGTTCATGCTGAGGAGAAAACCCCAGTTCCGCCAGTGAATCATATCGCCAGTAATGGTGTAATTTTTACCt                                         >  2:1188959/1‑93 (MQ=255)
                                                  gaTCATCGGCCCGACGCGAGAGTTCATGCTGAGGTGAAAACCCCAGTTCCGCCAGTGAATCATATCGCCAGTAATGGTGTAATTTTTACCTTc                                       <  2:1739003/93‑1 (MQ=255)
                                                  gaTCATCGGCCCGACGCGAGAGTTCATGCTGAGGTGAAAACCCCAGTTCCGCCAGTGAATCATATCGCCAGTAATGGTGTAATTTTTACCTTc                                       >  1:5900075/1‑93 (MQ=255)
                                                   atcatcGGCCCGACGCGAGAGTTCATGCTGAGGTGAAAACCCCAGTTCCGCCAGTGAATCATATCGCCAGTAATGTTGGAATTTTTTCCTTcc                                      >  1:827379/1‑92 (MQ=255)
                                                   atcatcGGCCCGACGCGAGAGTTCATGCTGAGGTGAAAACCCCAGTTCCGCCAGTGAATCATATCGCCAGTAATGGTGTAATTTTTACCTTCa                                      <  1:2783741/93‑1 (MQ=255)
                                                   atcatcGGCCCGACGCGAGAGTTCATGCTGAGGTGAAAACCCCAGTTCCGCCAGTGAATCATATCGCCAGTAATGGTGTAATTTTTACCTTCa                                      <  1:2936966/93‑1 (MQ=255)
                                                   atcatcGGCCCGACGCGAGAGTTCATGCTGAGGTGAAAACCCCAGTTCCGCCAGTGAATCATATCGCCAGTAATGGTGTAATTTTTACCTTCa                                      >  2:874048/1‑93 (MQ=255)
                                                    tcatcGGCCCGACGCGAGAGTTCATGCTGAGGTGAAAACCCCAGTTCCGCCAGTGAATCATATCGCCAGTAATGGTGTAATTTTTACCTTCAg                                     <  1:3376454/93‑1 (MQ=255)
                                                      agcGGCCCGACGCGAGAGTTCATGCTGAGGTGAAAACCCCAGTTCCGCCAGTGAATCATATCGCCAGTAATGGTGTAATTTTTACCTTCAGGc                                   <  1:992341/91‑1 (MQ=255)
                                                       tcGGCCCGACGCGAGAGTTCATGCTGAGGTGAAAACCCCAGTTCCGCCAGTGAATCATATCGCCAGTAATGGTGTAATTTTTACCTTCAGGCt                                  >  2:2996994/1‑93 (MQ=255)
                                                        cGGCCCGACGCGAGAGTTCATGCTGAGGTGAAAACCCCAGTTCCGCCAGTGAATCATATCGCCAGTAATGGTGTAATTTTTACCTTCAGGCTc                                 <  2:6536077/93‑1 (MQ=255)
                                                         ggCCCGACGCGAGAGTTCATGCTGAGGTGAAAACCCCAGTTCCGCCAGTGAATCATATCGCCAGTAATGGTGTAATTTTTACCTTCAGGCTCa                                >  1:4911600/1‑93 (MQ=255)
                                                           cccGACGCGAGAGTTCATGCTGAGGTGAAAACCCCAGTTCCGCCAGTGAATCATATCGCCAGTAATGGTGTAATTTTTACCTTCAGGCTCAAt                              <  1:2219256/93‑1 (MQ=255)
                                                            ccGACGCGCGAGTTCATGCTGAGGTGAAAACCCCAGTTCCGCCAGTGAATCATATCGCCAGTAATGGTGTAATTTTTACCTTCAGGCTCAATg                             >  2:3261596/1‑93 (MQ=255)
                                                            ccGACGCGAGAGTTCATGCTGAGGTGAAAACCCCAGTTCCGCCAGTGAATCATATCGCCAGTAATGGTGTAATTTTTACCTTCAGGCTCAAta                             >  1:6403752/1‑92 (MQ=255)
                                                            ccGACGCGAGAGTTCATGCTGAGGTGAAAACCCCAGTTCCGCCAGTGAATCATATCGCCAGTAATGGTGTAATTTTTACCTTCAGGCTCAATg                             >  1:737467/1‑93 (MQ=255)
                                                             cGACGCGAGAGTTCATGCTGAGGTGAAAACCCCAGTTCCGCCAGTGAATCATATCGCCAGTAATGGTGTAATTTTTACCTTCAGGCTCAATGa                            <  1:7100412/93‑1 (MQ=255)
                                                               aCGCGAGAGTTCATGCTGAGGTGAAAACCCCAGTTCCGCCAGTGAATCATATCGCCAGTAATGGTGTAATTTTTACCTTCAGGCTCAATGAtt                          <  2:6809893/93‑1 (MQ=255)
                                                                cgcgAGAGTTCATGCTGAGGTGAAATCCCCAGTTCCGCCAGTGAATCATATCGCCAGTAATGGTGTAATTTTTACCTTCAGGCTCAATGAttt                         <  1:7240478/93‑1 (MQ=255)
                                                                cgcgAGAGTTCATGCTGAGGTGAAAGCCCCAGTTCCGCCAGTGAATCATATCGCCAGTAATGGTGTAATTTTTACCTTCAGGCTCAATGAttt                         <  1:51363/93‑1 (MQ=255)
                                                                cgcgAGAGTTCATGCTGAGGTGAAAACCCCAGTTCCGCCAGTGAATCATATCGCCAGTAATGGTGTAATTTTTACCTTCAGGCTCAATGAttt                         <  2:344130/93‑1 (MQ=255)
                                                                 gcgAGAGTTCATGCTGAGGTGAAAACCCCAGTTCCGCCAGTGAATCATATCGCCAGTAATGGTGTAATTTTTACCTTCAGGCTCAATGATTTg                        >  1:6511164/1‑93 (MQ=255)
                                                                  cgAGAGTTCATGCTGAGGTGAAAACCCCAGTTCCGCCAGTGAATCATATCGCCAGTAATGGTGTAATTTTTACCTTCAGGCTCAATGATTTGc                       >  1:6978257/1‑93 (MQ=255)
                                                                    agagTTCATGCTGAGGTGAAAACCCCAGTTCCGCCAGTGAATCATATCGCCAGTAATGGTGTAATTTTTACCTTCAGGCTCAATGATTTGCAt                     >  2:2195321/1‑93 (MQ=255)
                                                                    agagTTCATGCTGAGGTGAAAACCCCAGTTCCGCCAGTGAATCATATCGCCAGTAATGGTGTAATTTTTACCTTCAGGCTCAATGATTTGCAt                     >  1:1933468/1‑93 (MQ=255)
                                                                     gagTTCATGCTGAGGTGAAAACCCCAGTTCCGCCAGTGAATCATATCGCCAGTAATGGTGTAATTTTTACCTTCAGGCTCAATGATTTGCATa                    >  1:888123/1‑93 (MQ=255)
                                                                       gTTCATGCTGAGGTGAAAACCCCAGTTCCGCCAGTGAATCATATCGCCAGTAATGGTGTAATTTTTACCTTCAGGCTCAATGATTTGCATAgg                  >  2:3116334/1‑93 (MQ=255)
                                                                        ttCATGCTGAGGTGAAAACCCCAGTTCCGCCAGTGAATCATATCGCCAGTAATGGTGTAATTTTTACCTTCAGGCTCAATGATTTGCATAGGc                 >  2:1287832/1‑93 (MQ=255)
                                                                         tCATGCTGAGGTGAAAACCCCAGTTCCGCCAGTGAATCATATCGCCAGTAATGGTGTAATTTTTACCTTCAGGCTCAATGATTTGCATAGGCt                >  2:6356016/1‑93 (MQ=255)
                                                                         tCATGCTGAGGTGAAAACCCCAGTTCCGCCAGTGAATCATATCGCCAGTAATGGTGTAATTTTTACCTTCAGGCTCAATGATTTGCATAGGCt                >  1:5559941/1‑93 (MQ=255)
                                                                          cATGCTGAGGTGAAAACCCCAGTTCCGCCAGTGAATCATATCGCCAGTAATGGTGTAATTTTTACCTTCAGGCTCAATGATTTGCATAGGCtt               <  1:1217114/93‑1 (MQ=255)
                                                                           aTGCTGAGGTGAAAACCCCAGTTCCGCCAGTGAATCATATCGCCAGTAATGGTGTAATTTTTACCTTCAGGCTCAATGATTTGCATAGGCTTa              >  2:7094959/1‑93 (MQ=255)
                                                                           aTGCTGAGGTGAAAACCCCAGTTCCGCCAGTGAATCATATCGCCAGTAATGGTGTAATTGTTACCTTCAGGCTCAATGATTCGCATAGGCTTa              >  2:7013089/1‑93 (MQ=255)
                                                                            tGCTGAGGTGAAAACCCCAGTTCCGCCAGTGAATCATATCGCCAGTAATGGTGTAATTTTTACCTTCAGGCTCAATGATTTGCATAGGCTTaa             >  1:4478218/1‑93 (MQ=255)
                                                                            tGCTGAGGTGAAAACCCCAGTTCCGCCAGTGAATCATATCGCCAGTAATGGTGTAATTTTTACCTTCAGGCTCAATGATTTGCATAGGCTTaa             <  2:2337582/93‑1 (MQ=255)
                                                                             gCTGAGGTGAAAACCCCAGTTCCGCCAGTGAATCATATCGCCAGTAATGGTGTAATTTTTACCTTCAGGCTCAATGATTTGCATAGGCTTAAc            >  1:6872508/1‑93 (MQ=255)
                                                                             gCTGAGGTGAAAACCCCAGTTCCGCCAGTGAATCATATCGCCAGTAATGGTGTAATTTTTACCTTCAGGCTCAATGATTTGCATAGGCTTAAc            <  2:6949688/93‑1 (MQ=255)
                                                                              cTGAGGTGAAAACCCCAGTTCCGCCAGTGAATCATATCGCCAGTAATGGTGTAATTTTTACCTTCAGGCTCAATGATTTGCATAGGCTTAACt           >  2:1118650/1‑93 (MQ=255)
                                                                              cTGAGGTGAAAACCCCAGTTCCGCCAGTGAATCATATCGCCAGTAATGGTGTAATTTTTACCTTCAGGCTCAATGATTTGCATAGGCTTAACt           <  2:7275697/93‑1 (MQ=255)
                                                                                 aGGTGAAAACCCCAGTTCCGCCAGTGAATCATATCGCCAGTAATGGTGTAATTTTTACCTTCAGGCTCAATGATTTGCATAGGCTTAACTGcc        <  2:5040287/93‑1 (MQ=255)
                                                                                  ggTAAAAACCCCAGTTCCGCCAGTGAATCATATCGCCAGTAATGGTGTAATTTTTACCTTCAGGCTCAATGATTTGCATAGGCTTAACTGCCg       <  2:3299625/93‑1 (MQ=255)
                                                                                    tGAAAACCCCAGTTCCGCCAGTGAATCATATCGCCAGTAATGGTGTAATTTTTACCTTCAGGCTCAATGATTTGCATAGGCTTAACTGCCGGa     <  1:3977042/93‑1 (MQ=255)
                                                                                    tGAAAACCCCAGTTCCGCCAGTGAATCATATCGCCAGTAATGGTGTAATTTTTACCTTCAGGCTCAATGATTTGCATAGGCTTAACTGCCGGa     <  1:7632499/93‑1 (MQ=255)
                                                                                      aaaaCCCCAGTTCCGCCAGTGAATCATATCGCCAGTAATGGTGTAATTTTTACCTTCAGGCTCAATGATTTGCATAGGCTTAACTGCCGGAGc   >  2:2587987/1‑93 (MQ=255)
                                                                                       aaaCCCCAGTTCCGCCAGTGAATCATATCGCCAGTAATGGTGTAATTTTTACCTTCAGGCTCAATGATTTGCATAGGCTTAACTGCCGGAGCa  <  2:7705420/93‑1 (MQ=255)
                                                                                          |                                                                                         
TCGTACATGACTTTGCGTTTGGTGCCATTGTCGTTATAAGTCACGGTGGAGATCATCGGCCCGACGCGAGAGTTCATGCTGAGGTGAAAATCCCAGTTCCGCCAGTGAATCATATCGCCAGTAATGGTGTAATTTTTACCTTCAGGCTCAATGATTTGCATAGGCTTAACTGCCGGAGCA  >  NC_000913/1448208‑1448387

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