Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913298,9600AC50.7% 76.2 / 113.7 70Y659D (TAT→GAT) paoCPaoABC aldehyde oxidoreductase, Moco‑containing subunit
Reads supporting (aligned to +/- strand):  ref base A (28/6);  new base C (0/35);  total (29/41)
Fisher's exact test for biased strand distribution p-value = 8.02e-14
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.00e-11

ATACACCGCGTTGGCGATAGCCGCGCTCACGCCGCACAGGCCCAGCTCACCGACACCTTTGGCCTTCATCGGGGAGGATATGGGGTCGGTATCATCCAGGAAAATCACCTCCTGTTTTGGGATATCCGCATGAACCGGCACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTGCGCGCAGTTTTCGGATTCAGGATGCGTCCTGCAGCAC  >  NC_000913/298814‑299103
                                                                                                                                                       |                                                                                                                                               
ataCACCGCGTTGGCGATAGCCGCGCTCACGCCGCACAGGCCCAGCTCACCGACACCTTTGGCCTTCATCGGGGAGGATATGGGGTCGGTATCATCCAGGAAAATCACCTCCTGTTTTGGGATATCCGCATGAACCGGCACCTCATAcccc                                                                                                                                             <  2:644099/151‑1 (MQ=255)
  acacCGCGTTGGCGATAGCCGCGCTCACGCCGCACAGGCCCAGCTCACCGACACCTTTGGCCTTCATCGGGGAGGATATGGGGTCGGTATCATCCAGGAAAATCACCTCCTGTTTTGGGATATCCGCATGAACCGGCACCTCATACCCCg                                                                                                                                            <  1:1345842/150‑1 (MQ=255)
  acacCGCGTTGGCGATAGCCGCGCTCACGCCGCACAGGCCCAGCTCACCGACACCTTTGGCCTTCATCGGGGAGGATATGGGGTCGGTATCATCCAGGAAAATCACCTCCTGTTTTGGGATATCCGCATGAACCGGCACCTCATACCCCGc                                                                                                                                           <  1:377873/151‑1 (MQ=255)
    acCGCGTTGGCGATAGCCGCGCTCACGCCGCACAGGCCCAGCTCACCGACACCTTTGGCCTTCATCGGGGAGGATATGGGGTCGGTATCATCCAGGAAAATCACCTCCTGTTTTGGGATATCCGCATGAACCGGCACCTCATACCCCGcc                                                                                                                                          >  1:1294220/1‑150 (MQ=255)
        cgTTGGCGATAGCCGCGCTCACGCCGCACAGGCCCAGCTCACCGACACCTTTGGCCTTCATCGGGGAGGATATGGGGTCGGTATCATCCAGGAAAATCACCTCCTGTTTTGGGATATCCGCATGAACCGGCACCTCATACCCCGCCATATc                                                                                                                                     >  1:43259/1‑151 (MQ=255)
         gTTGGCGATAGCCGCGCTCACGCCGCACAGGCCCAGCTCACCGACACCTTTGGCCTTCATCGGGGAGGATATGGGGTCGGTATCATCCAGGAAAATCACCTCCTGTTTTGGGATATCCGCATGAACCGGCACCTCATACCCCGCCATATCg                                                                                                                                    >  1:402257/1‑151 (MQ=255)
          ttGGCGATAGCCGCGCTCACGCCGCACAGGCCCAGCTCACCGACACCTTTGGCCTTCATCGGGGAGGATATGGGGTCGGTATCATCCAGGAAAATCACCTCCTGTTTTGGGATATCCGCATGAACCGGCACCTCATACCCCGCCa                                                                                                                                         >  2:1223203/1‑145 (MQ=255)
          ttGGCGATAGCCGCGCTCACGCCGCACAGGCCCAGCTCACCGACACCTCTGGCCTTCATCGGGGAGGATATGGGGTCGGTATCATCCAGGAAAATCACCTCCTGTTTTGGGATATCCGCATGAACCGGCACCTCATACCCCGCCa                                                                                                                                         <  1:1222361/145‑1 (MQ=255)
                atggCCGCGCTCACGCCGCACAGGCCCAGCCCACCGACCCCTTTGGCCTTCATCGGGGAGGATAGGGGTTGGGTATCACCCAGGAAAATCCCCTCCTGTTTGGGGATATCCGCATGAACCGGCACCTCATCCCCCGCCATATCGTGAtt                                                                                                                               <  2:79895/146‑1 (MQ=255)
                   gCCGCGCTCCCGCCGCACAGGCCCAGCTCACCGACACCTTTGGCCTTCATCGGGGAGGATATGGGGTCGGTATCATCCAGGAAAATCCCCTCCTGTTTTGGGATATCCGCATGAACCGGCCCCTCATCCCCCGCCATATCGTGATTAAc                                                                                                                            <  2:988219/149‑1 (MQ=255)
                   gCCGCGCTCACGCCGCACAGGCCCAGCTCACCGACACCTTTGGCCTTCATCGGGGCGGATATGGGGTCGGTATCATCCAGGAAAATCACCTCCTGTTTTGGGATATCCGCATGAACCGGCACCTCATACCCCGCCATATCGTGATTAACg                                                                                                                           >  1:1069645/1‑150 (MQ=255)
                    ccGCCCTCCCGCCGCACAGCCCCAGCTCCCCGCCCCCTTGGGCCTTCATCGGGGAGGATAGGGGGTCGGTATCACCCAGAAAAACCCCCTCCTTTTTTGGAATATCCGCAGGAACCGCCCCCTCATCCCCCGCCATATCGTGATTAACg                                                                                                                           <  1:882702/149‑1 (MQ=255)
                     cccgcaggcccgccacgGGCCCACCTGCCCGACACCTTTGCCCTTCATCGGGAAGGTAAGGGGTTCGGTACAAACCGGAAAACTCCCCTCCTTTGGGGGGTTGCCCGTATGCCCCGCCACCTAATACCCCGCCATATCGTGATTAACGAAg                                                                                                                        <  1:1016620/138‑1 (MQ=255)
                         cTCCCGCCGCACGGGCCCACCTCCCGGACCCTTTTGGCCTTCACCGGGGAGGATAGGGGGTCGGTTTCATCCAGAAAAATCCCCTCCTGTTTTGGGATATCCGAAAGAACCGGCCCCTAATCCCCCGCCATATCGTGATTAACGAAGTAg                                                                                                                     <  2:585918/150‑1 (MQ=255)
                           cACGCCGCACAGGCCCAGCTCACCGACCCCTTTGGCCTTCATCGGGGAGGATATGGGGTCGGTATCACCCAGAAAAACCCCCCCCTGTTTTGGGATATCCGCAGGAACCGGCACCTCATCCCCCGCCATATCGTGATTAACGAAGTAGccc                                                                                                                  <  2:984854/151‑1 (MQ=255)
                           cACGCCGCACAGGCCCAGCTCACCGACACCTTTGGCCTTCATCGGGGAGGATATGGGGTCGGTATCATCCAGGAAAATCACCTCCTGTTTTGGGATATCCGCATGAACCGGCACCTCATACCCCGCCATATCGTGATTAACGAAGTAGccc                                                                                                                  >  1:636096/1‑151 (MQ=255)
                              cccgcATAAGCCCCGGCCACCGACACTTTCGGCCTTCGTGGGGGAGGCTGAGGTTGTGGAAACCCCCAGGAAAATCCCCTCCTTTTGTGAGATTCCGCTAGGTACGGCCACAGAATCCCCCGCCATATCGTGATTAACGAAGTAGCCCa                                                                                                                 <  1:19106/148‑1 (MQ=255)
                                cgcACAGGCCCAGCCCCCCGACCCTTTGGGCCTTCATCGGGGAGGATAGGGGGTGGTTATCATCCAGAAAAATCCCCTCCTTTTTTGGGATATCCGCATGACCCGCCCCCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAAc                                                                                                              <  1:1306045/150‑1 (MQ=255)
                                   acaGGCCCAGCTCACCGACACCTTTGGCCTTCATCGGGGAGGATATGGGGTCGGTATCATCCAGGAAAATCACCTCCTGTTTTGGGATATCCGCATGAACCGGCACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTc                                                                                                          >  1:928543/1‑151 (MQ=255)
                                      ggCCCAGCTCCCCGACCCCTTTGGCCTTCATCGGGGAGGATAGGGGGTCGGTATCACCCAGAAAAACCCCCTCCTGTTTTGGGATACCCGCATGAACCGGCCCCTCATCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACgg                                                                                                            <  2:1086778/146‑1 (MQ=255)
                                      ggCCCAGCTCACCGACACCTTTGGCCTTCATCGGGGAGGATAGGGGGTCGGTATCATCCAGAAAAATCCCCTCCTTTTTTGGGATATCCGCATGAACCGGCCCCTCATCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTc                                                                                                          <  2:420866/148‑1 (MQ=255)
                                        cccAGCCCCCCGCCCCCTTTGGCCTTCATCGGGGAGGATAGGGGGTCGGTATCATCCAGGAAAATCCCCTCCTTTTTTGGGATATCCGCATGAACCGGCACCTCATCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATcc                                                                                                      <  2:916002/150‑1 (MQ=255)
                                        cccAGCCCCCCGCCACCTTTGGCCTTAAGGGGGGGGGATAGGGGGTCGGTATCACCCAGAAAAACCCCCTCCTGTTTTGGAAAACCCGCATGAACCGGCCCCTCATCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCa                                                                                                     <  1:27753/151‑1 (MQ=255)
                                            gCTCACCGACACCTTTGGCCTTCATCGGGGAGGATATGGGGTCGGTATCATCCAGGAAAATCACCTCCTGTTTTGGGATATCCGCATGAACCGGCACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGc                                                                                                 >  1:751242/1‑151 (MQ=255)
                                               cACCGCCCCTTTTGCCCTTCATGGGGGAGGATAGGGGGTCGTTACCACCAAGAAAACTCCCCTCTTTTTTGGGTATATCCGCATGAACCGGCCCCTCATCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCa                                                                                               <  2:265051/150‑1 (MQ=255)
                                                 ccgcccccTTTGCCCTCCATGGGGGAGGATAGGGGGCATGTACCACCCAGGAAAATCCCCTCTTGTTTTGGGTTATCCGCATGACCGGCCCCCCCATCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGct                                                                                            <  1:651930/145‑1 (MQ=255)
                                                  cGACACCTTTGGCCTTCATCGGGGAGGATATGGGGTCGGTATCATCCAGGAAAATCACCTCCTGTTTTGGGATATCCGCATGAACCGGCACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGctc                                                                                           >  1:121353/1‑151 (MQ=255)
                                                      cccTTTGGCCCTTCAGGGGGGGGGATGGGGGGGCGTTATCACCCAGAAAAACCCCCCCTTTTTTGGGTATACCCGAAGGACCGGGCCCCTCATCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGctcc                                                                                          <  1:73765/148‑1 (MQ=255)
                                                            ggCCTTCATGGGGGAGGAATGGGGGTAGATAGCACCCAGAAAAACCCCCTCCTGTTTGGGGATACCCGAATGACCCGGCCCCTCATCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAg                                                                                  <  1:1269810/150‑1 (MQ=255)
                                                             gcccTTCATGGGGGAGGATAGGGGGTCGTTATCACCCAGAAAAATCCCCCCCTGTTTGGGGATACCCGCATGACCCGCCCCCTCATCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGt                                                                                 <  2:838740/149‑1 (MQ=255)
                                                             gcccTTCAGCGGGGAGGATAGGGGTTCGTTATCCCCCAGAAAAATCCCCTCCTTTTTTGGGATATCCGCATGACCCCCCCCCTCATCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCa                                                                                   <  2:477981/147‑1 (MQ=255)
                                                                 ttCATCGGGGAGGATATGGGGTCGGTATCATCCAGGAAAATCACCTCCTGTTTTGGGATATCCGCATGAACCGGCACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCcg                                                                             >  1:1273360/1‑151 (MQ=255)
                                                                       ggggAGGATAGGGGGTGGTTATCACCAAGGAAAATCCCCTCCTGTTTTGGGATACCCGCATAACCCGCCCCCTCATCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCcgc                                                                            <  2:409509/146‑1 (MQ=255)
                                                                       ggggAGAATTGGGGGTGGGTACCACCAAGAAAACCCCCCCCTTTTTGTGGTATATCCGGTTAACCCGGCCCCTAACCCCCCGCCATATCGTGATTAGCGAAGTAGCCCAAGCGGTCATCCTCCGCCAGCTCCTCCATCAGTGc                                                                               <  1:213657/143‑1 (MQ=255)
                                                                       ggggAGAATTGGGGGCCGGTTTCACCCGGAAAAACCCCCCCCTTTTGAGGGATACCGCAAGAACCGGCCCCCTCATCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCcgc                                                                            <  2:121372/146‑1 (MQ=255)
                                                                        gggAGGATAGGGGGTGGGTATCACCCAGAAAAATCCCCTCTTGTTTTGGGATACCCGCATGACCCGGCCCCTCATCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCgcccat                                                                      <  2:468769/151‑1 (MQ=255)
                                                                        gggAGGATAGGGGGTCGGTATCACCCAGAAAAATCCCCCCCTGTTTTGGGATATCCGCATGACCCGCCCCCTCATCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCgcccat                                                                      <  2:1233507/151‑1 (MQ=255)
                                                                                  ggggTGGGTATCACCCAGAAAAATCCCCTCCTTTTTTGGGATACCCGCATAACCCGGCCCCTCATCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGt                                                             <  2:1320131/150‑1 (MQ=255)
                                                                                  ggggTCGTTATCACCGAGGAAAACCCCCCCTTTTTTTGGGATATCCACAGGACCGGCCCCCTCAACCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCgcccatgcccat                                                                <  2:94029/147‑1 (MQ=255)
                                                                                  ggggTCGGTATCATCCAGGAAAATCACCTCCTGTTTTGGGATATCCGCATGAACCGGCACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTc                                                            >  1:1368278/1‑151 (MQ=255)
                                                                                         gTACCACCCAGGAAACTCCCCTCCTGTTTGGGTATACCCGCATGACCCGGCCCCTCATCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTgcgc                                                     <  1:696477/151‑1 (MQ=255)
                                                                                          tATCATCCAGGAAAATCACCTCCTGTTTTGGGATATCCGCATGAACCGGCACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGcc                                                    >  2:1229410/1‑151 (MQ=255)
                                                                                                ccAGAAAAATCCCCTCTTTTTTTGGGATACCCGCATGACCCGGCCCCTCATCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTgcgc                                                     <  2:1065382/144‑1 (MQ=255)
                                                                                                ccAGAAAAATCCCCTCCTTTTTGGGGATCTCCGCATGACCCGCCACCTCATCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGAc                                              <  2:87355/151‑1 (MQ=255)
                                                                                                  aGGAAAATCACCTCCTGTTTTGGGATATCCGCATGAACCGGCACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCt                                            >  1:1178795/1‑151 (MQ=255)
                                                                                                     aaaaTCCCCTCCTGTTTGGGGATCCCCCCATGACCCGCCCCCTCACCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATg                                                <  1:102854/144‑1 (MQ=255)
                                                                                                       aaccccctcctGTTTTGGGATATCCGAATGACCCGCCCCCTCATCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTg                                       <  2:528718/146‑1 (MQ=255)
                                                                                                        accccctcttGTTTTGGGATACCCGCAGAACCCGGCCCCTCATCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTg                                       <  2:585478/146‑1 (MQ=255)
                                                                                                         tCACCTCCTGTTTTGGGATATCCGCATGAACCGGCACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTgcg                                     >  1:1065351/1‑151 (MQ=255)
                                                                                                         tCACCTCCTGTTTTGGGATATCCGCATGAACCGGCACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTgcg                                     >  2:658472/1‑151 (MQ=255)
                                                                                                          cACCTCCTGTTTTGGGATATCCGCATGAACCGGCACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTgcg                                     >  1:581665/1‑150 (MQ=255)
                                                                                                           aCCTCCTGTTTTGGGATATCCGCATGAACCGGCACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTgcgcg                                   >  1:395882/1‑151 (MQ=255)
                                                                                                             ctcctGTTTTGGGATATCCGCATGAACCGGCACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTgcgcgc                                  >  2:726682/1‑150 (MQ=255)
                                                                                                             ctcctGTTTTGGGATATCCGCATGAACCGGCACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTgcgcgc                                  >  1:1240883/1‑150 (MQ=255)
                                                                                                             ctcctGTTTTGGGATATCCGCATGAACCGGCACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAATCATTGCGCCAATGACCTGGCTgcgcgc                                  >  2:690327/1‑150 (MQ=255)
                                                                                                                 tGTTTTGGGATATCCGCATGAACCGGCACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTGCGCGCAGttt                             >  1:1381215/1‑151 (MQ=255)
                                                                                                                   ttttGTTATCTCCGCTTGTACCGGCACCTTTCTCCCCGCCATATCGTGCGTGACGAGGTAACCCAGCCGGTCATCTACCGGCAGCTCCTGCAGCAGGGCGGCGGCCTTGCCCATAGTCATTGCGCCAAGGTCCTGGCGGCGCGCGGttt                             >  1:56818/1‑149 (MQ=255)
                                                                                                                    tttGGGATATCCGCATGAACCGGCACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTGCGCGCAGTTTTCg                          >  1:159232/1‑151 (MQ=255)
                                                                                                                       tGGGATATCCGCATAAACCGGCACCTCATACCCCGCCATATCGTGCTTAACGAAGTAGACCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCAGGCCCATAGTCATTGCGCCAATGACCTGGCTGCGCGCAGGTTTCgg                         >  1:492882/1‑150 (MQ=255)
                                                                                                                             tatCCGCATGAACCGGCAACTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCACTGACCTGGCCGCGCGCAGTTTTCGGATTc                     >  1:376573/1‑149 (MQ=255)
                                                                                                                                ccGCATGAACCGGCCCCTCATCCCCCGCCATATCGTGATCAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTGCGCGCAGTTTTCGGATTCa                    <  2:392565/147‑1 (MQ=255)
                                                                                                                                ccGCATGAACCGCCCCCTCATCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTGCGCGCAGTTTTCGGATTCAgg                  <  2:395803/149‑1 (MQ=255)
                                                                                                                                  gCATGAACCGGCACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTGCGCGCAGTTTTCGGATTCAGGATGc              >  2:1044251/1‑151 (MQ=255)
                                                                                                                                     catcaaccccccccTCATCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTGCGCGCAGTTTTCGGATTCAGGATg               <  1:351050/137‑1 (MQ=255)
                                                                                                                                         aCCGGCACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTGCGCGCAGTTTTCGGATTCAGGATGCGTCCTg        >  1:450130/1‑151 (MQ=255)
                                                                                                                                          ccGGCCCGTCAACCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTGCGCGCAGTTTTCGGATTCAGGATGCGTCCt         <  1:1132032/149‑1 (MQ=255)
                                                                                                                                            cgcccccTCATCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTGCGCGCAGTTTTCGGATTCAGGATGCGTcc          <  1:397130/142‑1 (MQ=255)
                                                                                                                                            cgcccccTCATCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTGCGCGCAGTTTTCGGATTCAGGATGCGTCCt         <  1:719108/143‑1 (MQ=255)
                                                                                                                                                aCCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTGCGCGCAGTTTTCGGATTCAGGATGCGTCCTGCAGcac  >  2:757136/1‑151 (MQ=255)
                                                                                                                                                 ccTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTGCGCGCAGTTTTCGGATTCAGGATGCGTCCTGCAGcac  >  1:871811/1‑150 (MQ=255)
                                                                                                                                                   tcccaCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTGCGCGCAGTTTTCGGATTCAGGATGCGTCCTGCAGcac  <  1:914293/144‑1 (MQ=255)
                                                                                                                                                       |                                                                                                                                               
ATACACCGCGTTGGCGATAGCCGCGCTCACGCCGCACAGGCCCAGCTCACCGACACCTTTGGCCTTCATCGGGGAGGATATGGGGTCGGTATCATCCAGGAAAATCACCTCCTGTTTTGGGATATCCGCATGAACCGGCACCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTGCGCGCAGTTTTCGGATTCAGGATGCGTCCTGCAGCAC  >  NC_000913/298814‑299103

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: