Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 796,681 C→A 63.8% A43D (GCC→GAC)  modC → molybdate ABC transporter ATPase; chlorate resistance protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913796,6810CA63.8% 96.1 / 130.7 120A43D (GCC→GAC) modCmolybdate ABC transporter ATPase; chlorate resistance protein
Reads supporting (aligned to +/- strand):  ref base C (21/22);  new base A (40/36);  total (61/59)
Fisher's exact test for biased strand distribution p-value = 7.07e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.14e-01

Only 100 of 122 total aligned reads displayed.

TGGGCAACCATTGCCTGACTATTAATGAAACGCTGCCCGCCAATGGCATCACTGCTATCTTTGGCGTCTCCGGTGCCGGAAAAACTTCGCTGATTAACGCCATCAGTGGACTGACGCGCCCGCAAAAAGGGCGGATTGTCCTCAATGGGCGGGTACTAAATGATGCCGAAAAAGGTATCTGCCTGACGCCGGAAAAGCG  >  NC_000913/796582‑796780
                                                                                                   |                                                                                                   
tGGGCAACCATTGCCTGACTATTAATGAAACGCTGCCCGCCAATGGCATCACTGCTATCTTTGGCGTCTCCGGTGCCGGAAAAACTTCGCTGATTAACGc                                                                                                     >  2:2832342/1‑100 (MQ=255)
  ggCAACCATTGCCTGACTATTAATGAAACGCTGCCCGCCAATGGCATCACTGCTATCTTTTGCGTCTCCGGTGCCGGAAAAACTTCGCTGATTAACGACa                                                                                                   >  2:316219/1‑100 (MQ=255)
  ggCAACCATTGCCTGACTATTAATGAAACGCTGCCCGCCAATGGCATCACTGCTATCTTTGGCGTCTCCGGTGCCGGAAAAACTTCGCTGATTAACGACa                                                                                                   >  2:1890021/1‑100 (MQ=255)
  ggCAACCATTGCCTGACTATTAATGAAACGCTGCCCGCCAATGGCATCACTGCTATCTTTGGCGTCTCCGGTGCCGGAAAAACTTCGCTGATTAACGACa                                                                                                   >  1:588104/1‑100 (MQ=255)
    cAACCATTGCCTGACTATTAATGAAACGCTGCCCGCCAATGGCATCACTGCTATCTTTGGCGTCTCCGGTGCCGGAAAAACTTCGCTGATTAACGACATc                                                                                                 <  1:2922945/100‑1 (MQ=255)
    cAACCATTGCCTGACTATTAATGAAACGCTGCCCGCCAATGGCATCACTGCTATCTTTGGCGTCTCCGGTGCCGGAAAAACTTCGCTGATTAACGACATc                                                                                                 <  2:3798512/100‑1 (MQ=255)
       ccATTGCCTGACTATTAATGAAACGCTGCCCGCCAATGGCATCACTGCTATCTTTGGCGTCTCCGGTGCCGGAAAAACTTCGCTGATTAACGACATCAGt                                                                                              <  1:4113097/100‑1 (MQ=255)
         aTTGCCTGACTATTAATGAAACGCTGCCCGCCAATGGCATGACTGCTATCTTTGGCGTCTCCGGTGCCGGAAAAACTTCGCTGATTAACGACATCAGTg                                                                                             <  2:2336877/99‑1 (MQ=255)
         aTTGCCTGACTATTAATGAAACGCTGCCCGCCAATGGCATCACTGCTATCTTTGGCGTCTCCGGTGCCGGAAAAACTTCGCTGATTAACGACATCAGTgg                                                                                            <  2:671491/100‑1 (MQ=255)
          ttGCTTGATTATTAATGAAACGCTGCCCGTCAATGGCATTACTTCTATCTTTGGCGTCTCCGGTGCCGGAAAAACTTCGCTGATTAACGACATCAGTgg                                                                                            <  1:3454129/99‑1 (MQ=255)
           tGCCTGACTATTAATGAAACGCTGCCCGCCAATGGCATCACTGCTATCTTTGGCGTCTCCGGTGCCGGAAAAACTTCGCTGATTAACGCCATCAGTGGAc                                                                                          <  1:309981/100‑1 (MQ=255)
           tGCCTGACTATTAATGAAACGCTGCCCGCCAATGGCATCACTGCTATCTTTGGCGTCTCCGGTGCCGGAAAAACTTCGCTGATTAACGCCATCAGTGGAc                                                                                          >  2:3272617/1‑100 (MQ=255)
           tGCCTGACTATTAATGAAACGCTGCCCGCCAATGGCATCACTGCTATCTTTGGCGTCTCCGGTGCCGGAAAAACTTCGCTGATTAACGACATCAGTGGAc                                                                                          <  2:4133830/100‑1 (MQ=255)
           tGCCTGACTATTAATGAAACGCTGCCCGCCAATGGCATCACTGCTATCTTTGGCGTCTCCGGTGCCGGAAAAACTTCGCTGATTAACGACATCAGTGGAc                                                                                          >  2:666897/1‑100 (MQ=255)
            gCCTGACTATTAATGAAACGCTGCCCGCCAATGGCATCACTGCTATCTTTGGCGTCTCCGGTGCCGGAAAAACTTCGCTGATTAACGCCATCAGTGGACt                                                                                         >  2:2057754/1‑100 (MQ=255)
             ccTGACTATTAATGAAACGCTGCCCGCCAATGGCATCACTGCTATCTTTGGCGTCTCCGGTGCCGGAAAAACTTCGCTGATTAACGACATCAGTGGACTg                                                                                        >  2:3072361/1‑100 (MQ=255)
             ccTGACTATTAATGAAACGCTGCCCGCCAATGGCATCACTGCTATCTTTGGCGTCTCCGGTGCCGGAAAAACTTCGCTGATTAACGACATCAGTGGACTg                                                                                        >  1:778185/1‑100 (MQ=255)
             ccTGACTATTAATGAAACGCTGCCCGCCAATGGCATCACTGCTATCTTTGCCGTCTCCGGTGCCGGAAAAACTTCGCTGATTAACGACATCAGTGGACTg                                                                                        >  1:1600903/1‑100 (MQ=255)
              cTGACTATTAATGAAACGCTGCCCGCCAATGGCATCACTGCTATCTTTGGCGTCTCCGGTGCCGGAAAAACTTCGCTGATTAACGCCATCAGTGGACTGa                                                                                       >  1:9108/1‑100 (MQ=255)
              cTGACTATTAATGAAACGCTGCCCGCCAATGGCATCACTGCTATCTTTGGCGTCTCCGGTGCCGGAAAAACTTCGCTGATTAACGACATCAGTGGACTGa                                                                                       <  2:1043100/100‑1 (MQ=255)
                    aTTAATGAAACGCTGCCCGCCAGTGGCATCACTGCTATCTTTGGCGTCTCCGCTGCCGGAAAAACCTCGCTGACTAACGCCATCAGTGGGCTGACGCGcc                                                                                 >  2:899489/1‑100 (MQ=255)
                    aTTAATGAAACGCTGCCCGCCAATGGCATCACTGCTATCTTTGGCGTCTCCGGTGCCGGAAAAACTTCGCTGATTAACGGCATCAGTGGACTGACGCGcc                                                                                 <  2:2514999/100‑1 (MQ=255)
                    aTTAATGAAACGCTGCCCGCCAATGGCATCACTGCTATCTTTGGCGTCTCCGGTGCCGGAAAAACTTCGCTGATTAACGACATCAGTGGACTGACGCGcc                                                                                 >  1:626029/1‑100 (MQ=255)
                    aTTAATGAAACGCTGCCCGCCAATGGCATCACTGCTATCTTTGGCGTCTCCGGTGCCGGAAAAACTTCGCTGATTAACGACATCAGTGGACTGACGCGcc                                                                                 >  1:2526765/1‑100 (MQ=255)
                    aTTAATGAAACGCTGCCCGCCAATGGCATCACTGCTATCTTTGGCGTCTCCGGTGCCGGAAAAACTTCGCTGATTAACGACATCAGTGGACTGACGCGcc                                                                                 >  1:167889/1‑100 (MQ=255)
                     ttAATGAAACGCTGCCCGCCAATGGCATCACTGCTATCTTTGGCGTCTCCGGTGCCGGAAAAACTTCGCTGATTAACGACATCAGTGGACTGACGCGccc                                                                                >  1:2465628/1‑100 (MQ=255)
                      tAATGAAACGCTGCCCGCCAATGGCATCACTGCTATCTTTGGCGTCTCCGGTGCCGGAAAAACTTCGCTGATTAACGACATCAGTGGACTGACGCGCCCg                                                                               <  1:3816474/100‑1 (MQ=255)
                       aaTGAAACGCTGCCCGCCAATGGCATCACTGCTATCTTTGGCGTCTCCGGTGCCGGAAAAACTTCGCTGATTAACGACATCAGTGGACTGACGCGCCCGc                                                                              >  2:709895/1‑100 (MQ=255)
                       aaTGAAACGCTGCCCGCCAATGGCATCACTGCTATCTTTGGCGTCTCCGGTGCCAGAAAAACTTCGCTGATTAACGACATCAGTGGACTGACGCGCCCGc                                                                              <  1:1813960/100‑1 (MQ=255)
                       aaTGAAACGCTGCCCGCCAATGGCATCACTGCTATCTTTGGCGTCTCCGGTGCCAGAAAAACTTCGCTGATTAACGACATCAGTGGACTGACGCGCCCGc                                                                              <  1:1814023/100‑1 (MQ=255)
                       aaTGAAACGATGCCCGCCAATGGCATCACTGCTATCGTTGGCGTCTCCGGTGCCGGAAAAACTTCGCTGATTAACGACATCAGTGGACTGACGCGCCCGc                                                                              >  1:2910625/1‑100 (MQ=255)
                        aTGAAACGCTGCCCGCCAATGGCATCACTGCTATCTTTGGCGTCTCCGGTGCCGGAAAAACTTCGCTGATTAACGCCATCAGTGGACTGACGCGCCCGCa                                                                             >  1:1277330/1‑100 (MQ=255)
                        aTGAAACGCTGCCCGCCAATGGCATCACTGCTATCTTTGGCGTCTCCGGTGCCGGAAAAACTTCGCTGATTAACGACATCAGTGGACTGACGcgcacgca                                                                             >  2:1286498/1‑100 (MQ=255)
                         tGAAACGCTGCCCGCCAATGGCATCACTGCTATCTTTGGCGTCTCCGGTGCCGGAAAAACTTCGCTGATTAACGCCATCAGTGACTGACGCGCCCGCaaa                                                                           >  2:2182241/1‑100 (MQ=255)
                         tGAAACGCTGCCCGCCAATGGCATCACTGCTATCTTTGGCGTCTCCGGTGCCGGAAAAACTTCGCTGATTAACGACATCAGTGGACTGACGCGCCCGCaa                                                                            >  2:2795459/1‑100 (MQ=255)
                          gAAACGCTGCCCGCCAATGGCATCACTGCTATCTTTGGCGTCTCCGGTGCCGGAAAAACTTCGCTGATTAACGCCATCAGTGGACTGACGCGCCCGCaaa                                                                           >  2:1615754/1‑100 (MQ=255)
                           aaaCGCTGCCCGCCAATGGCATCACTGCTATCTTTGGCGTCTCCGGTGCCGGAAAAACTTCGCTGATTAACGCCATCAGTGGACTGACGCGCCCGCaaaa                                                                          >  1:2042587/1‑100 (MQ=255)
                              cGCTGCCCGCCAATGGCATCACTGCTATCTTTGGCGTCTCCGGTGCCGGAAAAACTTCGCTGATTAACGACATCAGTGGACTGACGCGCCCGCAAAAAgg                                                                       >  2:2531247/1‑100 (MQ=255)
                              cGCTGCCCGCCAATGGCATCACTGCTATCTTTGGCGTCTCCGGTGCCGGAAAAACTTCGCTGATTAACGACATCAGTGGACTGACGCGCCCGCAAAAAg                                                                        >  1:580783/1‑99 (MQ=255)
                               gCTGCCCGCCAATGGCATCACTGCTATCTTTGGCGTCTCCGGTGCCGGAAAAACTTCGCTGATTAACGCCATCAGTGGACTGACGCGCCCGCAAAAAggg                                                                      <  1:3498883/100‑1 (MQ=255)
                               gCTGCCCGCCAATGGCATCACTGCTATCTTTGGCGTCTCCGGTGCCGGAAAAACTTCGCTGATTAACGCCATCAGTGGACTGACGCGCCCGCAAAAAggg                                                                      <  1:2237036/100‑1 (MQ=255)
                                cTGCCCGCCAATGGCATCACTGCTATCTTTGGCGTCTCCGGTGCCGGAAAAACTTCGCTGATTAACGCCATCAGTGGACTGACGCGCCCGCAAAAAGGGc                                                                     >  1:2732285/1‑100 (MQ=255)
                                 tGCCCGCCAATGGCATCACTGCTATCTTTGGCGTCTCCGGTGCCGGAAAAACTTCGCTGATTAACGCCATCAGTGGACTGACGCGCCCGCAAAAAGGGCg                                                                    <  1:1136065/100‑1 (MQ=255)
                                  gCCGGCCAATGGCATCACTGCTATCTTTGGCTTCTCCGGTGCCGGAAAAACTTCGCTGATTAACGACATCAGTGGACTGACGCGCCCGCAAAAAGGGCg                                                                    <  1:3833819/99‑1 (MQ=255)
                                  gCCCGCCAATGGCATCACTGCTATCTTTGGCGTCTCCGGTGCCGGAAAAACTTCGCTGATTAACGCCATCAGTGGACTGACGCGCCCGCAAAAAGGGCgg                                                                   >  1:1562081/1‑100 (MQ=255)
                                   cccGCAAATGGCATCACTGCTACCTATTGCGTCACCGCTGTGAAAACACCTTCATTGATTCACGCCCTCCGTGGACTGACGCACACGTAAAAAcgg                                                                      >  2:921453/1‑96 (MQ=255)
                                     cGCCAATGGCATCGCTGCTATCTTTGGCGTCTCCGGTGCCGGAAAAACTTCGCTGATTAACGCCATCAGTGGACTGACGCGCCCGCAAAAAGGGCgg                                                                   <  2:2222453/97‑1 (MQ=255)
                                     cGCCAATGGCATCACTGCTATCTTTGGCGTCTCCGGTGCCGGAAAAACTTCGCTGATTAACGACATCCGTGGACTGACGCGCCCGCAAAAAGGGCGGa                                                                  >  2:3375825/1‑98 (MQ=255)
                                      gCCAATGGCATCACTGCTATCTTTGGCGTCTCCGGTGCCGGAAAAACTTCGCTGATTAACGCCATCAGTGGACTGACGCGCCCGCAAAAAGGGCGGATTg                                                               <  1:2756996/100‑1 (MQ=255)
                                      gCCAATGGCATCACTGCTATCTTTGGCGTCTCCGGTGCCGGAAAAACTTCGCTGATTAACGCCATCAGTGGACTGACGCGCCCGCAAAAAGGGCGGATTg                                                               <  2:668402/100‑1 (MQ=255)
                                      gCCAATGGCATCACTGCTATCTTTGGCGTCTCCGGTGCCGGAAAAACTTCGCTGATTAACGACATCAGTGGACTGACGCGCCCGCAAAAAGGGCGGATTg                                                               >  1:2918420/1‑100 (MQ=255)
                                      gCCAATGGCATCACTGCTATCTTTGGCGTCTCCGGTGCCGGAAAAACTTCGCTGATTAACGACATCAGTGGACTGACGCGCCCGCAAAAAGGGCGGATTg                                                               <  1:4057507/100‑1 (MQ=255)
                                       ccAATGGCATCACTGCTATCTTTGGCGTCTCCGGTGCCGGAAAAACTTCGCTGATTAACGACATCAGTGGACTGACGCGCCCGCAAAAAGGGCGGATTGt                                                              >  1:2237088/1‑100 (MQ=255)
                                       ccAATGGCATCACTGCTATCTTTGGCGTCTCCGGTGCCGGAAAAACTTCGCTGATTAACGACATCAGTGGACTGACGCGCCCGCAAAAAGGGCGGATTGt                                                              >  2:2019179/1‑100 (MQ=255)
                                       ccAATGGCATCACTGCTAACTTTGGCGTATCCGGTGCCGGAAAAACTCCGCTGATTAACGCCATCAGTGGACTGACGCGCCCGCAAAAAGGGCGGa                                                                  >  2:4121285/1‑96 (MQ=255)
                                         aaTGGCATCACTGCTATCTTTGGCGTCTCCGGTGCCGGAAAAACTTCGCTGATTAACGACATCAGTGGACTGACGCGCCCGCAAAAAGGGCGGATTGTcc                                                            <  1:110913/100‑1 (MQ=255)
                                         aaTGGCATCACTGCTATCTTTGGCGTCTCCGGTGCCGGAAAAACTTCGCTGATTAACGACATCAGTGGACTGACGCGCCCGCAAAAAGGGCGGATTGTcc                                                            >  2:1048687/1‑100 (MQ=255)
                                         aaTGGCATCACTGCTATCTTTGGCGTCTCCGGTGCCCGAAAAACTTCGCTGATTAACGCCATCAGTGGACTGACGCGCCCGCAAAAAGGGCGGATTGTcc                                                            <  1:2889953/100‑1 (MQ=255)
                                             gCATCACTGCTATCTTTGGCGTCTCCGGTGCCGGAACAACTTCGCTGATTAACGCCATCAGTGGACTGGCGCGCCCGCAAAAAGGGCGGATTGTCCTCaa                                                        >  2:2783232/1‑100 (MQ=255)
                                              cATCACTGCTATCTTTGGCGTCTCCGGTGCCGGAAAAACTTCGCTGATTAACGCCATCAGTGGACTGACGCGCCCGCAAAAAGGGCGGATTGTCCTCAAt                                                       >  2:2143685/1‑100 (MQ=255)
                                                   cTGCTATCTTTGGCGTCTCCGGTGCCGGAAAAACTTCGCTGATTAACGCCATCAGTGGACTGACGCGCCCGCAACAAGGGCGGATTGACCTCAAt                                                       >  2:3366486/1‑95 (MQ=255)
                                                   cTGCTATCTTTGGCGTCTCCGGTGCCGGAAAAACTTCGCTGATTAACGACATCAGTGGACTGACGCGCCCGCAAAAAGGGCGGATTGTCCTCAATGGGCg                                                  <  1:139289/100‑1 (MQ=255)
                                                     gCTATCTTTGGCGTCTCCGGTGCCGGAAAAACTTCGCTGATTAACGACATCAGTGGACTGACGCGCCCGCAAAAAGGGCGGATTGTCCTCAATGGGCggg                                                <  1:3539544/100‑1 (MQ=255)
                                                     gCTATCTTTGGCGTCTCCGGTGCCGGAAAAACTTCGCTGATTAACGACATCAGTGGACTGACGCGCCCGCAAAAAGGGCGGATTGTCCTCAATGGGCggg                                                <  2:2628432/100‑1 (MQ=255)
                                                     gCTATCTTTGGCGTCTCCGGTGCCGGAAAAACTTCGCTGATTAACGACATCAGTGGACTGACGCGCCCGCAAAAAAGGCGGATTGTCCTCAATGGGCggg                                                >  2:3776588/1‑100 (MQ=255)
                                                      cTATCTTTGGCGTCTCCGGTGCCGGAAAAACTTCGCTGATTAACGCCATCAGTGGACTGACGCGCCCGCAAAAAGGGCGGATTGTCCTCAATGGGCGGGt                                               <  1:486032/100‑1 (MQ=255)
                                                      cTATCTTTGGCGTCTCCGGTGCCGGAAAAACGTCGCTGATTAACGACATCAGTGGACTGACGCGCCCGCAAAAAGGGCGGATTGTCCTCAATGGGCGGGt                                               <  2:2014503/100‑1 (MQ=255)
                                                       tATCTTTGGCGTCTCCGGTGCCGGAAAAACTTCGCTGATTAACGACATCAGTGGACTGACGCGCCCGCAAAAAGGGCGGATTGTCCTCAATGGGCGGGTa                                              >  1:406349/1‑100 (MQ=255)
                                                       tATCTTTGGCGTCTCCGGTGCCGGAAAAACTTCGCTGATTAACGACATCAGTGGACTGACGCGCCCGCAAAAAGGGCGGATTGTCCTCAATGGGCGGGTa                                              <  1:1743706/100‑1 (MQ=255)
                                                       tATCTTTGGCGTCTCCGGTGCCGGAAAAACTTCGCTGATTAACGACATCAGTGGACTGACGCGCCCGCAAAAAGGGCGGATTGTCCTCAATGGGCGGGTa                                              >  2:3372568/1‑100 (MQ=255)
                                                        aTCTTTGGCGTCTCCGGTGCCGGAAAAACTTCGCTGATTAACGACATCAGTGGACTGACGCGCCCGCAAAAAGGGCGGATTGTCCTCAATGGGCGGGTAc                                             <  2:3641112/100‑1 (MQ=255)
                                                        aTCTTAGGCGTCTCCGGTGCCGGAAAAACTTCGCTGATTAACGACATCAGTGTTCTAACGAAACAAAAAAAAGGGCGGATTGTCCTCAATGGGCGGGTAc                                             <  2:363778/100‑1 (MQ=255)
                                                           tttGGCGTCTCCGGTGCCGGAAAAACTTCGCTGATTAACGCCATCAGTGGACTGACGCGCCCGCAAAAAGGGCGGATTGTCCTCAATGGGCGGGTACTaa                                          >  1:417260/1‑100 (MQ=255)
                                                            ttGGCGTCTCCGGTGCCGGAAAAACTTCGCTGATTAACGCCATCAGTGGACTGACGCGCCCGCAAAAAGGGCGGATTGTCCTCAATGGGCGGGTACTaaa                                         <  2:3298494/100‑1 (MQ=255)
                                                             tGGCGTCTCCGGTGCCGGAAAAACTTCGCTGATTAACGACATCAGTGGACTGACGCGCCCGCAAAAAGGGCGGATTGTCTTCAATGGGCGGGTACTAAat                                        <  2:2526751/100‑1 (MQ=255)
                                                                 gTCTCCGGTGCCGGAAAAACTTCGCTGATTAACGACATCAGTGGACTGACGCGCCCGCAAAAAGGGCGGATTGTCCTCAATGGGCGGGTACTAAatgatg                                    <  2:2774449/100‑1 (MQ=255)
                                                                 gTCTCCGGTGCCGGAAAAACTTCGCTGATTAACGACATCAGTGGACTGACGCGCCCGCAAAAAGGGCGGATTGTCCTCAATGGGCGGGTACTAAatgatg                                    <  1:1127601/100‑1 (MQ=255)
                                                                  tctcCGGTGCCGGAAAAACTTCGCTGATTAACGCCATCAGTGGACTGACGCGCCCGCAAAAAGGGCGGATTGTCCTCAATGGGCGGGTACTAAATGATGc                                   <  1:815121/100‑1 (MQ=255)
                                                                     ccGGTGCCGGAAAAACTTCGCTGATTAACGCCATCAGTGGACTGACGCGCCCGCAAAAAGGGCGGATTGTCCTCAATGGGCGGGTACTAAATGATGCCGa                                <  1:541487/100‑1 (MQ=255)
                                                                      cGGTGCCGGAAAAACTTCGCTGATTAACGCCATCAGTGGACTGACGCGCCCGCAAAATGGGCGGATTGACCTCAATGGGCGGGTACCAAATGATGCCgta                               <  1:2792579/100‑3 (MQ=255)
                                                                         tGCCGGAAAAACTTCGCTGATTAACGACATCAGTGGACTGACGCGCCCGCAAAAAGGGCGGATTGTCCTCAATGGGCGGGTACTAAATGATGGCGaaaaa                            <  1:2169994/100‑1 (MQ=255)
                                                                         tGCCGGAAAAACTTCGCTGATTAACGACATCAGTGGACTGACGCGCCCGCAAAAAGGGCGGATTGTCCTCAATGGGCGGGTACTAAATGATGCCGaaaaa                            >  1:3337765/1‑100 (MQ=255)
                                                                          gCCGGAAAAACTTCGCTGATTAACGCCATCAGTGGACTGACGCGCCCGCAAAAAGGGCGGATTGTCCTCAATGGGCGGGTACTAAATGATGCCGAAAAAg                           <  2:2728918/100‑1 (MQ=255)
                                                                           ccGGAAAAACTTCGCTGATTAACGACATCAGTGGACTGACGCGCCCGCAAAAAGGGCGGATTGTCCTCAATGGGCGGGTACTAAATGATGCCGAAAAAgg                          >  1:1162791/1‑100 (MQ=255)
                                                                            cGGAAAAACTTCGCTGATTAACGCCATCAGTGGACTGACGCGCCCGCAAAAAGGGCGGATTGTCCTCAATGGGCGGGTACTAAATGATGCCGAAAAAGGt                         >  1:3109455/1‑100 (MQ=255)
                                                                             ggAAAAACTTCGCTGATTAACGCCATCAGTGGACTGACGCGCCCGCAAAAAGGGCGGATTGTCCTCAATGGGCGGGTACTAAATGATGCCGAAAAAGGTa                        <  2:3109436/100‑1 (MQ=255)
                                                                             ggAAAAAAGTCGCTGATTAACGACATCAGTGGACTGAAAAAAAAAAAAAAAGGGCGGATTGTCCTCAATGGGCGGGTACTAAATGATGCCGAAAAAGGTa                        <  2:3571081/100‑1 (MQ=255)
                                                                               aaaaaCTTCGCTGATTAACGACATCAGTGGACTGACGCGCCCGCAAAAAGGGCGGATTGTCCTCAATGGGCGGGTACTAAATGATGCCGAAAAAGGTATc                      >  1:821610/1‑100 (MQ=255)
                                                                                 aaaCTTCGCTGATTAACGACATCAGTGGACTGACGCGCCCGCAAAAAGGGCGGATTGTCCTCAATGGGCGGGTACTAAATGATGCCGAAAAAGGTATCTg                    >  2:1074661/1‑100 (MQ=255)
                                                                                      tCGCTGATTAACGCCATCAGTGGACTGACGCGCCCGCAACAAGGGCGGATTGTCCTCAATGGGCGGGTAGTAAATGATGCCGAAAAAGGTATCTGGCTGa               >  1:2053633/1‑100 (MQ=255)
                                                                                        gCTGATTAACGACATCAGTGGACTGACGCGCCCGCAAAAAGGGCGGATTGTCCTCAATGGGCGGGTACTAAATGATGCCGAAAAAGGTATCTGCCTGACg             >  1:1810511/1‑100 (MQ=255)
                                                                                        gCTGATTAACGACATCAGTGGACTGACGCGCCCGCAAAAAGGGCGGATTGTCCTCAATGGGCGGGTACTAAATGATGCCGAAAAAGGTATCTGCCTGACg             >  2:3612749/1‑100 (MQ=255)
                                                                                        gCTGATTAACGACATCAGTGGACTGACGCGCCCGCAAAAAGGGCGGATTGTCCTCAATGGGCGGGTACTAAATGATGCCGAAAAAGGTATCTGCCTGACg             >  2:764607/1‑100 (MQ=255)
                                                                                        gCTGATTAACGACATCAGTGGACTGACGCGCCCGCAAAAAGGGCGGATTGTCCTCAATGGGCGGGTACTAAATGATGCCGAAAAAGGTATCTGCCTGACg             >  1:107219/1‑100 (MQ=255)
                                                                                          tGATTAACGCCATCAGTGGACTGACGCGCCCGCAAAAAGGGCGGATTGTCCTCAATGGGCGGGTACTAAATGATGCCGAAAAAGGTATCTGCCTGACGcc           <  1:353404/100‑1 (MQ=255)
                                                                                           gATTAACGACATCAGTGGACTGACGCGCCCGCAAAAAGGGCGGATTGTCCTCAATGGGCGGGTACTAAATGATGCCGAAAAAGGTATCTGCCTGACGCCg          >  2:2234485/1‑100 (MQ=255)
                                                                                           gATTAACGACATCAGTGGACTGACGCGCCCGCAAAAAGGGCGGATTGTCCTCAATGGGCGGGTACTAAATGATGCCGAAAAAGGTATCTGCCTGACGCCg          >  2:951519/1‑100 (MQ=255)
                                                                                             ttAACGCCATCAGTGGACTGCCGCGCCCGCAAAAAGCGCGGATTGTCCTCAATGGGCGGGTCCTAAATGATGCCGAAAAAGGTATCTGCCTGACGCCg          >  2:3614656/1‑98 (MQ=255)
                                                                                              tAACGCCATCAGTGGACTGACGCGCCCGCAAAAAGGGCGGATTGTCCTCAATGGGCGGGTACTAAATGATGCCGAAAAAGGTATCTGCCTGACGCCGGaa       >  2:3913273/1‑100 (MQ=255)
                                                                                                   ccATCAGTGGACTGACGCGCCCGCAAAAAGGGCGGATTGTCCTCAATGGGCGGGTACTAAATGATGCCGAAAAAGGTATCTGCCTGACGCCGGAGAAGCg  >  2:857802/1‑100 (MQ=255)
                                                                                                   |                                                                                                   
TGGGCAACCATTGCCTGACTATTAATGAAACGCTGCCCGCCAATGGCATCACTGCTATCTTTGGCGTCTCCGGTGCCGGAAAAACTTCGCTGATTAACGCCATCAGTGGACTGACGCGCCCGCAAAAAGGGCGGATTGTCCTCAATGGGCGGGTACTAAATGATGCCGAAAAAGGTATCTGCCTGACGCCGGAAAAGCG  >  NC_000913/796582‑796780

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