Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 508,943 T→A 6.7% M813L (ATG→TTG)  copA ← Cu(+) exporting P‑type ATPase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913508,9430TA 6.7% 235.1 / 5.8 74M813L (ATG→TTG) copACu(+) exporting P‑type ATPase
Reads supporting (aligned to +/- strand):  ref base T (33/36);  new base A (3/2);  total (36/38)
Fisher's exact test for biased strand distribution p-value = 6.70e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.53e-01

TAAGGCGTTTACGCCGCATCCGGCAGTCATGCGTCGATGCCAAATGCGCCACCCTAAAGCAGCGCATCCGCAATGATGTACTTATTCCTTCGGTTTAAACCGCAGCAACCGGTTGGCGTTACTCACTACGGTAATCGACGAGAGCGCCATTGCCGCTCCGGCAACTACCGGGTTAAGCAGTGTTCCAGTGAACGGCCACAAAATACCGGCGGCGACCGGAATACCGATACTGTTGTAGATAAACGCACCGAGCAGGTTCTGCTTCATGTTGTGCAGCGTTGCGCGGGAAATAGCGAGCG  >  NC_000913/508794‑509092
                                                                                                                                                     |                                                                                                                                                     
tAAGGCGTTTACGCCGCATCCGGCAGTCATGCGTCGATGCCAAATGCGCCACCCTAAAGCAGCGCATCCGCAATGATGTACTTATTCCTTCGGTTTAAACCGCAGCAACCGGTTGGCGTTACTCACTACGGTAATCGACGAGAGCGCCAt                                                                                                                                                       >  2:1186016/1‑150 (MQ=255)
   ggCGTTTACGCCGCATCCGGCAGTCATGCGTCGATGCCAAATGCGCCACCCTAAAGCAGCGCATCCGCAATGATGTACTTATTCCTTCGGTTTAAACCGCAGCAACCGGTTGGCGTTACTCACTACGGTAATCGACGAGAGCGCCATTGc                                                                                                                                                    <  1:63030/150‑1 (MQ=255)
      gTTTACGCCGCATCCGGCAGTCATGCGTCGATGCCAAATGCGCCACCCTAAAGCAGCGCATCCGCAATGATGTACTTATTCCTTCGGTTTAAACCGCAGCAACCGGTTGGCGTTACTCACTACGGTAATCGACGAGAGCGCCATTGCCGc                                                                                                                                                 >  1:316991/1‑150 (MQ=255)
        ttACGCCGCATCCGGCAGTCATGCGTCGATGCCAAATGCGCCACCCTAAAGCAGCGCATCCGCAATGATGTACTTATTCCTTCGGTTTAAACCGCAGCAACCGGTTGGCGTTACTCACTACGGTAATCGACGAGAGCGCCATTGCCGCTc                                                                                                                                               >  2:1161222/1‑150 (MQ=255)
         tACGCCGCATCCGGCAGTCATGCGTCGATGCCAAATGCGCCACCCTAAAGCAGCGCATCCGCAATGATGTACTTATTCCTTCGGTTTAAACCGCAGCAACCGGTTGGCGTTACTCACTACGGTAATCGACGAGAGCGGCAATGGCGCTct                                                                                                                                              <  1:1097037/150‑2 (MQ=255)
          aCGCCGCATCCGGCAGTCATGCGTCGATGCCAAATGCGCCACCCTAAAGCAGCGCATCCGCAATGATGTACTTATTCCTTCGGTTTAAACCGCAGCAACCGGTTGGCGTTACTCACTACGGTAATCGACGAGAGCGCCATTGCCGCTCCg                                                                                                                                             <  1:894376/150‑1 (MQ=255)
           cgccgcATCCGGCAGTCATGCGTCGATGCCAAATGCGCCACCCTAAAGCAGCGCATCCGCAATGATGTACTTATTCCTTCGGTTTAAACCGCAGCAACCGGTTGGCGTTACTCACTACGGTAATCGACGAGAGCGCCATTGCCGCTCCgg                                                                                                                                            >  2:1410930/1‑150 (MQ=255)
           cgccgcATCCGGCAGTCATGCGTCGATGCCAAATGCGCCACCCTAAAGCAGCGCATCCGCAATGATGTACTTATTCCTTCGGTTTAAACCGCAGCAACCGGTTGGCGTTACTCACTACGGTAATCGACGAGAGCGCCATTGCCGCTCCgg                                                                                                                                            <  2:316991/150‑1 (MQ=255)
           cgccgcATCCGGCAGTCATGCGTCGATGCCAAATGCGCCACCCTAAAGCAGCGCATCCGCAATGATGTACTTATTCCTTCGGTTTAAACCGCAGCAACCGGTTGGCGTTACTCACTACGGTAATCGACGAGAGCGCCATTGCCGCTCCgg                                                                                                                                            <  2:702753/150‑1 (MQ=255)
            gccgcATCCGGCAGTCATGCGTCGATGCCAAATGCGCCACCCTAAAGCAGCGCATCCGCAATGATGTACTTATTCCTTCGGTTTAAACCGCAGCAACCGGTTGGCGTTACTCACTACGGTAATCGACGAGAGCGCCATTGCCGCTCCGGc                                                                                                                                           <  1:887797/150‑1 (MQ=255)
             ccgcATCCGGCAGTCATGCGTCGATGCCAAATGCGCCACCCTAAAGCAGCGCATCCGCAATGATGTACTTATTCCTTCGGTTTAAACCGCAGCAACCGGTTGGCGTTACTCACTACGGTAATCGACGAGAGCGGCAATGGCGCTctcgtc                                                                                                                                          <  1:726802/150‑6 (MQ=255)
             ccgcATCCGGCAGTCATGCGTCGATGCCAAATGCGCCACCCTAAAGCAGCGCATCCGCAATGATGTACTTATTCCTTCGGTTTAAACCGCAGCAACCGGTTGGCGTTACTCACTACGGTAATCGACGAGAGCGCCATTGCCGCTCCGGCa                                                                                                                                          >  2:542527/1‑150 (MQ=255)
              cgcATCCGGCAGTCATGCGTCGATGCCAAATGCGCCACCCTAAAGCAGCGCATCCGCAATGATGTACTTATTCCTTCGGTTTAAACCGCAGCAACCGGTTGGCGTTACTCACTACGGTAATCGACGAGAGCGCCATTGCCGCTCCGGCaa                                                                                                                                         <  1:914358/150‑1 (MQ=255)
                cATCCGGCAGTCATGCGTCGATGCCAAATGCGCCACCCTAAAGCAGCGCATCCGCAATGATGTACTTATTCCTTCGGTTTAAACCGCAGCAACCGGTTGGCGTTACTCACTACGGTAATCGACGAGAGCGCCATTGCCGCTCCGGCAACt                                                                                                                                       <  2:238028/150‑1 (MQ=255)
                   ccGGCAGTCATGCGTCGATGCCAAATGCGCCACCCTAAAGCAGCGCATCCGCAATGATGTACTTATTCCTTCGGTTTAAACCGCAGCAACCGGTTGGCGTTACTCACTACGGTAATCGACGAGAGCGCCATTGCCGCTCCGGCAACTAcc                                                                                                                                    >  2:660671/1‑150 (MQ=255)
                     ggCAGTCATGCGTCGATGCCAAATGCGCCACCCTAAAGCAGCGCATCCGCAATGATGTACTTATTCCTTCGGTTTAAACCGCAGCAACCGGTTGGCGTTACTCACTACGGTAATCGACGAGAGCGCCATTGCCGCTCCGGCAACTACCgg                                                                                                                                  <  1:493202/150‑1 (MQ=255)
                           cATGCGTCGATGCCAAATGCGCCACCCTAAAGCAGCGCATCCGCAATGATGTACTTATTCCTTCGGTTTAAACCGCAGCAACCGGTTGGCGTTACTCACTACGGTAATCGACGAGAGCGCCATTGCCGCTCCGGCAACTACCGGGTTAAg                                                                                                                            >  2:26046/1‑150 (MQ=255)
                             tGCGTCGATGCCAAATGCGCCACCCTAAAGCAGCGCATCCGCAATGATGTACTTATTCCTTCGGTTTAAACCGCAGCAACCGGTTGGCGTTACTCACTACGGTAATCGACGAGAGCGCCATTGCCGCTCCGGCAACTACCGGGTTAAGCa                                                                                                                          >  1:534074/1‑150 (MQ=255)
                                    aTGCCAAATGCGCCACCCTAAAGCAGCGCATCCGCAATGATGTACTTATTCCTTCGGTTTAAACCGCAGCAACCGGTTGGCGTTACTCACTACGGTAATCGACGAGAGCGCCATTGCCGCTCCGGCAACTACCGGGTTAAGCAGTGTTcc                                                                                                                   >  2:1194911/1‑150 (MQ=255)
                                       ccAAATGCGCCACCCTAAAGCAGCGCATCCGCAATGATGTACTTATTCCTTCGGTTTAAACCGCAGCAACCGGTTGGCGTTACTCACTACGGTAATCGACGAGAGCGCCATTGCCGCTCCGGCAACTACCGGGTTAAGCAGTGTTCCAGt                                                                                                                <  1:1194911/150‑1 (MQ=255)
                                            tGCGCCACCCTAAAGCAGCGCATCCGCAATGATGTACTTATTCCTTCGGTTTAAACCGCTGCAACCGGTTGGCGTTACTCACTACGGTAATCGACGAGAGCGCCATTGCCGCTCCGGCAACTACCGGGTTAAGCAGTGTTCCAGTGAACg                                                                                                           <  1:1161222/150‑1 (MQ=255)
                                             gcgcCACCCTAAAGCAGCGCATCCGCAATGATGTACTTATTCCTTCGGTTTAAACCGCAGCAACCGGTTGGCGTTACTCACTACGGTAATCGACGAGAGCGCCATTGCCGCTCCGGCAACTACCGGGTTAAGCAGTGTTCCAGTGAACgg                                                                                                          <  2:835629/150‑1 (MQ=255)
                                              cgcCACCCTAAAGCAGCGCATCCGCAATGATGTACTTATTCCTTCGGTTTAAACCGCAGCAACCGGTTGGCGTTACTCACTACGGTAATCGACGAGAGCGCCATTGCCGCTCCGGCAACTACCGGGTTAAGCAGTGTTCCAGTGAACGGc                                                                                                         >  1:885614/1‑150 (MQ=255)
                                                   cccTAAAGCAGCGCATCCGCAATGATGTACTTATTCCTTCGGTTTAAACCGCAGCAACCGGTTGGCGTTACTCACTACGGTAATCGACGAGAGCGCCATTGCCGCTCCGGCAACTACCGGGTTAAGCAGTGTTCCAGTGAACGGCCACaa                                                                                                    >  2:345703/1‑150 (MQ=255)
                                                    ccTAAAGCAGCGCATCCGCAATGATGTACTTATTCCTTCGGTTTAAACCGCAGCAACCGGTTGGCGTTACTCACTACGGTAATCGACGAGAGCGCCATTGCCGCTCCGGCAACTACCGGGTTAAGCAGTGTTCCAGTGAACGGCCACaaa                                                                                                   >  1:637305/1‑150 (MQ=255)
                                                      tAAAGCAGCGCATCCGCAATGATGTACTTATTCCTTCGGTTTAAACCGCAGCAACCGGTTGGCGTTACTCACTACGGTAATCGACGAGAGCGCCATTGCCGCTCCGGCAACTACCGGGTTAAGCAGTGTTCCAGTGAACGGCAACAAAAt                                                                                                 >  1:1311341/1‑150 (MQ=255)
                                                      tAAAGCAGCGCATCCGCAAAGAGGTACTTATTACTTCGGTTTAAACCGCAGCAACCGGTTGGCGTTACTCACTACGGTAATCGACGAGAGCGCCATTGCCGCTCCGGCAACTACCGGGTTAAGCAGTGTTCCAGTGAACGGCCACAAAAt                                                                                                 >  2:1074219/1‑150 (MQ=255)
                                                        aaGCAGCGCATCCGCAATGATGTACTTATTCCTTCGGTTTAAACCGCAGCAACCGGTTGGCGTTACTCACTACGGTAATCGACGAGAGCGCCATTGCCGCTCCGGCAACTACCGGGTTAAGCAGTGTTCCAGTGAACGGCCACAAAATAc                                                                                               <  1:824465/150‑1 (MQ=255)
                                                          gcagcGCATCCGCAATGATGTACTTATTCCTTCGGTTTAAACCGCAGCAACCGGTTGGCGTTACTCACTACGGTAATCGACGAGAGCGCCATTGCCGCTCCGGCAACTACCGGGTTAAGCAGTGTTCCAGTGAACGGCCACAAAATACcg                                                                                             >  1:532777/1‑150 (MQ=255)
                                                              cgcATCCGCAATGATGTACTTATTCCTTCGGTTTAAACCGCAGCAACCGGTTGGCGTTACTCACTACGGTAATCGACGAGAGCGCCATTGCCGCTCCGGCAACTACCGGGTTAAGCAGTGTTCCAGTGAACGGCCACAAAATACcggcgg                                                                                         >  1:13352/1‑150 (MQ=255)
                                                               gcATCCGCAATGATGTACTTATTCCTTCGGTTTAAACCGCAGCAACCGGTTGGCGTTACTCACTACCGTAATCGACGAGAGCGCCATTGCCGCTCCGGCAACTACCGGGTTAAGCAGTGTTCCAGTGAACGGCCACAAAATACCggcggc                                                                                        <  1:256741/150‑1 (MQ=255)
                                                                cATCCGCAATGATGTACTTATTCCTTCGGTTTAAACCGCAGAAACCGGTTGTAGTTACTCACTACGGTAATCGACGAGAGCGCCATTGCCGCTCCGCAAACTACCGGGTTAAGCAGTGTTCCAGAGAACGGCCACAAAATACCGgcggcg                                                                                       >  1:379369/1‑150 (MQ=255)
                                                                   ccGCAATGATGTACTTATTCCTTCGGTTTAAACCGCAGCAACCGGTTGGCGTTACTCACTACGGTAATCGACGAGAGCGCCATTGCCGCTCCGGCAACTACCGGGTTAAGCAGTGTTCCAGTGAACGGCCACAAAATACCGGCGGCGAcc                                                                                    >  2:775127/1‑150 (MQ=255)
                                                                     gCAATGATGTACTTATTCCTTCGGTTTAAACCGCAGCAACCGGTTGGCGTTACTCACTACGGTAATCGACGAGAGCGCCATTGCCGCTCCGGCAACTACCGGGTTAAGCAGTGTTCCAGTGAACGGCCACAAAATACCGGCGGCGACCgg                                                                                  <  2:1361156/150‑1 (MQ=255)
                                                                      cAATGATGTACTTATTCCTTCGGTTTAAACCGCAGCAACCGGTTGGCGTTACTCACTACGGTAATCGACGAGAGCGCCATTGCCGCTCCGGCAACTACCGGGTTAAGCAGTGTTCCAGTGAACGGCCACAAAATACCGGCGGCGACCGGa                                                                                 <  2:1002144/150‑1 (MQ=255)
                                                                        atgatgTACTTATTCCTTCGGTTTAAACCGCAGCAACCGGTTGGCGTTACTCACTACGGTAATCGACGAGAGCGCCATTGCCGCTCCGGCAACTACCGGGTTAAGCAGTGTTCCAGTGAACGGCCACAAAATACCGGCGGCGACCGGAAt                                                                               <  2:814772/150‑1 (MQ=255)
                                                                          gatgTACTTATTCCTTCGGTTTAAACCGCAGCAACCGGTTGGCGTTACTCACTACGGTAATCGACGAGAGCGCCATTGCCGCTCCGGCAACTACCGGGTTAAGCAGTGTTCCAGTGAACGGCCACAAAATACCGGCGGCGACCGGAATAc                                                                             <  2:16428/150‑1 (MQ=255)
                                                                           atgTACTTATTCCTTCGGTTTAAACCGCAGCAACCGGTTGGCGTTACTCACTACGGTAATCGACGAGAGCGCCATTGCCGCTCCGGCAACTACCGGGTTAAGCAGTGTTCCAGTGAACGGCCACAAAATACCGGCGGCGACCGGAATAcc                                                                            <  2:616917/150‑1 (MQ=255)
                                                                            tgTACTTATTCCTTCGGTTTAAACCGCAGCAACCGGTTGGCGTTACTCACTACGGTAATCGACGAGAGCGCCATTGCCGCTCCGGCAACTACCGGGTTAAGCAGTGTTCCAGTGAACGGCCACAAAATACCGGCGGCGACCGGAATACCg                                                                           >  1:1076077/1‑150 (MQ=255)
                                                                             gTACTTATTCCTTCGGTTTAAACCGCAGCAACCGGTTGGCGTTACTCACTACGGTAATCGACGAGAGCGCCATTGCCGCTCCGGCAACTACCGGGTTAAGCAGTGTTCCAGTGAACGGCCACAAAATACCGGCGGCGACCGGAATACCGa                                                                          <  1:866086/150‑1 (MQ=255)
                                                                                cTTATTCCTTCGGTTTAAACCGCAGCAACCGGTTGGCGTTACTCACTACGGTAATCGACGAGAGCGCCATTGCCGCTCCGGCAACTACCGGGTTAAGCAGTGTTCCAGTGAACGGCCACAAAATACCGGCGGCGACCGGAATACCGATAc                                                                       >  2:660962/1‑150 (MQ=255)
                                                                                     tCCTTCGGTTTAAACCGCAGCAACCGGTTGGCGTTACTCACTACGGTAATCGACGAGAGCGCCATTGCCGCTCCGGCAACTACCGGGTTAAGCAGTGTTCCAGTGAACGGCCACAAAATACCGGCGGCGACCGGAATACCGATACtgttg                                                                  <  1:1410930/150‑1 (MQ=255)
                                                                                      ccTTCGGTTTAAACCGCAGCAACCGGTTGGCGTTACTCACTACGGTAATCGACGAGAGCGCCATTGCCGCTCCGGCAACTACCGGGTTAAGCAGTGTTCCAGTGAACGGCCACAAAATACCGGCGGCGACCGGAATACCGATACtgttgt                                                                 >  2:1378371/1‑150 (MQ=255)
                                                                                       cTTCGGTTTAAACCGCAGCAACCGGTTGGCGTTACTCACTACGGTAATCGACGAGAGCGCCATTGCCGCTCCGGCAACTACCGGGTTAAGCAGTGTTCCAGTGAACGGCCACAAAATACCGGCGGCGACCGGAATACCGATACTGTTGTa                                                                >  2:1360572/1‑150 (MQ=255)
                                                                                         tCGGTTTAAACCGCAGCAACCGGTTGGCGTTACTCACTACGGTAATCGACGAGAGCGCCATTGCCGCTCCGGCAACTACCGGGTTAAGCAGTGTTCCAGTGAACGGCCACAAAATACCGGCGGCGACCGGAATACCGATACTGTTGTAGa                                                              >  1:869683/1‑150 (MQ=255)
                                                                                          cGGTTTAAACCGCAGCAACCGGTTGGCGTTACTCACTACGGTAATCGACGAGAGCGCCATTGCCGCTCCGGCAACTACCGGGTTAAGCAGTGTTCCAGTGAACGGCCACAAAATACCGGCGGCGACCGGAATACCGATACTGTTGTAGAt                                                             >  2:492372/1‑150 (MQ=255)
                                                                                            gTTTAAACCGCAGCAACCGGTTGGCGTTACTCACTACGGTAATCGACGAGAGCGCCATTGCCGCTCCGGCAACTACCGGGTTAAGCAGTGTTCCAGTGAACGGCCACAAAATACCGGCGGCGACCGGAATACCGATACTGTTGTAGATaa                                                           <  1:660671/150‑1 (MQ=255)
                                                                                             tttAAACCGCAGCAACCGGTTGGCGTTACTCACTACGGTAATCGACGAGAGCGCCATTGCCGCTCCGGCAACTACCGGGTTAAGCAGTGTTCCAGTGAACGGCCACAAAATACCGGCGGCGACCGGAATACCGATACTGTTGTAGATaaa                                                          >  2:582485/1‑150 (MQ=255)
                                                                                              ttAAACCGCAGCAACCGGTTGGCGTTACTCACTACGGTAATCGACGAGAGCGCCATTGCCGCTCCGGCAACTACCGGGTTAAGCAGTGTTCCAGTGAACGGCCACAAAATACCGGCGGCGACCGGAATACCGATACTGTTGTAGATAAAc                                                         >  2:1380646/1‑150 (MQ=255)
                                                                                               tAAACCGCAGCAACCGGTTGGCGTTACTCACTACGGTAATCGACGAGAGCGCCATTGCCGCTCCGGCAACTACCGGGTTAAGCAGTGTTCCAGTGAACGGCCACAAAATACCGGCGGCGACCGGAATACCGATACTGTTGTAGATAAACg                                                        <  2:641260/150‑1 (MQ=255)
                                                                                                   ccGCAGCAACCGGTTGGCGTTACTCACTACGGTAATCGACGAGAGCGCCATTGCCGCTCCGGCAACTACCGGGTTAAGCAGTGTTCCAGTGAACGGCCACAAAATACCGGCGGCGACCGGAATACCGATACTGTTGTAGATAAACGCAcc                                                    >  1:687062/1‑150 (MQ=255)
                                                                                                      cagcaACCGGTTGGCGTTACTCACTACGGTAATCGACGAGAGCGCCATTGCCGCTCCGGCAACTACCGGGTTAAGCAGTGTTCCAGTGAACGGCCACAAAATACCGGCGGCGACCGGAATACCGATACTGTTGTAGATAAACGCACCGAg                                                 >  1:1292725/1‑150 (MQ=255)
                                                                                                       agcaACCGGTTGGCGTTACTCACTACGGTAATCGACGAGAGCGCCATTGCCGCTCCGGCAACTACCGGGTTAAGCAGTGTTCCAGTGAACGGCCACAAAATACCGGCGGCGACCGGAATACCGATACTGTTGTAGATAAACGCACCGAGc                                                <  1:542527/150‑1 (MQ=255)
                                                                                                        gcaACCGGTTGGCGTTACTCACTACGGTAATCGACGAGAGCGCCATTGCCGCTCCGGCAACTACCGGGTTAAGCAGTGTTCCAGTGAACGGCCACAAAATACCGGCGGCGACCGGAATACCGATACTGTTGTAGATAAACGCACCGAGCa                                               >  1:1264773/1‑150 (MQ=255)
                                                                                                                ttGGCGTTACTCACTACGGTAATCGACGAGAGCGCCATTGCCGCTCCGGCAACTACCGGGTTAAGCAGTGTTCCAGTGAACGGCCACAAAATACCGGCGGCGACCGGAATACCGATACTGTTGTAGATAAACGCACCGAGCAGGTTCTGc                                       <  1:443122/150‑1 (MQ=255)
                                                                                                                 tGGCGTTACTCACTACGGTAATCGACGAGAGCGCCATTGCCGCTCCGGCAACTACCGGGTTAAGCAGTGTTCCAGTGAACGGCCACAAAATACCGGCGGCGACCGGAATACCGATACTGTTGTAGATAAACGCACCGAGCAGGTTCTGCt                                      <  2:637305/150‑1 (MQ=255)
                                                                                                                    cGTTACTCACTACGGTAATCGACGAGAGCGCCATTGCCGCTCCGGCAACTACCGGGTTAAGCAGTGTTCCAGTGAACGGCCACAAAATACCGGCGGCGACCGGAATACCGATACTGTTGTAGATAAACGCACCGAGCAGGTTCTGCTTCa                                   <  1:335250/150‑1 (MQ=255)
                                                                                                                       tACTCACTACGGTAATCGACGAGAGCGCCATTGCCGCTCCGGCAACTACCGGGTTAAGCAGTGTTCCAGTGAACGGCCACAAAATACCGGCGGCGACCGGAATACCGATACTGTTGTAGATAAACGCACCGAGCAGGTTCTGCTTCAtgt                                <  1:660962/150‑1 (MQ=255)
                                                                                                                         cTCACTACGGTAATCGACGAGAGCGCCATTGCCGCTCCGGCAACTACCGGGTTAAGCAGTGTTCCAGTGAACGGCCACAAAATACCGGCGGCGACCGGAATACCGATACTGTTGTAGATAAACGCACCGAGCAGGTTCTGCTTCAtgttg                              <  1:1360572/150‑1 (MQ=255)
                                                                                                                          tCACTACGGTAATCGACGAGAGCGCCATTGCCGCTCCGGCAACTACCGGGTTAAGCAGTGTTCCAGTGAACGGCCACAAAATACCGGCGGCGACCGGAATACCGATACTGTTGTAGATAAACGCACCGAGCAGGTTCTGCTTCAtgttgt                             >  2:235216/1‑150 (MQ=255)
                                                                                                                            aCTACGGTAATCGACGAGAGCGCCATTGCCGCTCCGGCAACTACCGGGTTAAGCAGTGTTCCAGTGAACGGCCACAAAATACCGGCGGCGACCGGAATACCGATACTGTTGTAGATAAACGCACCGAGCAGGTTCTGCTTCATGTTGTGc                           <  1:492372/150‑1 (MQ=255)
                                                                                                                            aCTACGGTAATCGACGAGAGCGCCATTGCCGCTCCGGCAACTACCGGGTTAAGCAGTGTTCCAGTGAACGGCCACAAAATACCGGCGGCGACCGGAATACCGATACTGTTGTAGATAAACGCACCGAGCAGGTTCTGCTTCATGTTGTGc                           >  1:1382646/1‑150 (MQ=255)
                                                                                                                                cGGTAATCGACGAGAGCGCCATTGCCGCTCCGGCAACTACCGGGTTAAGCAGTGTTCCAGTTAAGGGCCACAAAATACCGGCGGCGACCGGAATACCTATACTGTTGTAGATAAACGCACCGAGCAGGTTCTGCTTCATGTTGTGCAGCg                       <  1:26046/150‑1 (MQ=255)
                                                                                                                                  gTAATCGACGAGAGCGCCATTGCCGCTCCGGCAACTACCGGGTTAAGCAGTGTTCCAGTGAACGGCCACAAAATACCGGCGGCGACCGGAATACCGATACTGTTGTAGATAAACGCACCGAGCAGGTTCTGCTTCATGTTGTGCAGCGtt                     >  2:1196493/1‑150 (MQ=255)
                                                                                                                                        tgccggagCGGCAATGGCGCTCTGGCAACTACCGGGTTAAGCAGTGTTCCAGTGAACGGCCACAAAATACCGGCGGCGACCGGAATACCGATACTGTTGTAGATAAACGCACCGAGCAGGTTCTGCTTCATGTTGTGCAGCGTTGCGCgg               >  1:565884/6‑150 (MQ=255)
                                                                                                                                        tgccggagCGGCAATGGCGCTCCGGCAACTACCGGGTTAAGCAGTGTTCCAGTGAACGGCCCCAAAATACCGGCGGCGCCCGGAATACCGATACTGTTGTAGATAAACGCCCCGAGCAGGTTCTGCTTCATGTTGTGCAGCGTTGCGCgg               >  1:394644/6‑150 (MQ=255)
                                                                                                                                        gacgaGAGCGCCATTGCCGCTCCGGCAACTACCGGGTTAAGCAGTGTTCCAGTGAACGGCCACAAAATACCGGCGGCGACCGGAATACCGATACTGTTGTAGATAAACGCACCGAGCAGGTTCTGCTTCATGTTGTGCAGCGTTGCGCgg               >  2:1110964/1‑150 (MQ=255)
                                                                                                                                          cgaGAGCGCCATTGCCGCTCCGGCAACTACCGGGTTAAGCAGTGTTCCAGTGAACGGCCACAAAATACCGGCGGCGACCGGAATACCGATACTGTTGTAGATAAACGCACCGAGCAGGTTCTGCTTCATGTTGTGCAGCGTTGCGCGGGa             <  1:690188/150‑1 (MQ=255)
                                                                                                                                          cgaGAGCGCCATTGCCGCTCCGGCAACTACCGGGTTAAGCAGTGTTCCAGTGAACGGCCACAAAATACCGGCGGCGACCGGAATACCGATACTGTTGTAGATAAACGCACCGAGCAGGTTCTGCTTCATGTTGTGCAGCGTTGCGCGGGa             <  1:1196493/150‑1 (MQ=255)
                                                                                                                                          ccggagCGGCAATGGCGCTCCGGCAACTACCGGGTTAAGCAGTGTTCCAGTGAACGGCCACAAAATACCGGCGGCGACCGGAATACCGATACTGTTGTAGATAAACGCACCGAGCAGGTTCTGCTTCATGTTGTGCAGCGTTGCGCGGGa             >  2:1371615/4‑150 (MQ=255)
                                                                                                                                              agCGCCATTGCCGCTCCGGCAACTACCGGGTTAAGCAGTGTTCCAGTGAACGGCCACAAAATACCGGCGGCGACCGGAATACCGATACTGTTGTAGATAAACGCACCGAGCAGGTTCTGCTTCATGTTGTGCAGCGTTGCGCGGGAAATa         <  1:1418498/150‑1 (MQ=255)
                                                                                                                                              agCGCCATTGCCGCTCCGGCAACTACCGGGTTAAGCAGTGTTCCAGTGAACGGCCACAAAATACCGGCGGCGACCGGAATACCGATACTGTTGTAGATAAACGCACCGAGCAGGTTCTGCTTCATGTTGTGCAGCGTTGCGCGGGAAATa         <  1:345703/150‑1 (MQ=255)
                                                                                                                                               gcgcCATTGCCGCTCCGGCAACTACCGGGTTAAGCAGTGTTCCAGTGAACGGCCACAAAATACCGGCGGCGACCGGAATACCGATACTGTTGTAGATAAACGCACCGAGCAGGTTCTGCTTCATGTTGTGCAGCGTTGCGCGGGAAATag        >  1:985821/1‑150 (MQ=255)
                                                                                                                                                 gcCATTGCCGCTCCGGCAACTACCGGGTTAAGCAGTGTTCCAGTGAACGGCCACAAAATACCGGCGGCGACCGGAATACCGATACTGTTGTAGATAAACGCACCGAGCAGGTTCTGCTTCATGTTGTGCAGCGTTGCGCGGGAAATagcg      <  1:582485/150‑1 (MQ=255)
                                                                                                                                                   cATTGCCGCTCCGGCAACTACCGGGTTAAGCAGTGTTCCAGTGAACGGCCACAAAATACCGGCGGCGACCGGAATACCGATACTGTTGTAGATAAACGCACCGAGCAGGTTCTGCTTCATGTTGTGCAGCGTTGCGCGGGAAATagcgag    <  2:885614/150‑1 (MQ=255)
                                                                                                                                                     ttGCCGCTCCGGCAACTACCGGGTTAAGCAGTGTTCCAGTGAACGGCCACAAAATACCGGCGGCGACCGGAATACCGATACTGTTGTAGATAAACGCACCGAGCAGGTTCTGCTTCATGTTGTGCAGCGTTGCGCGGGAGATagcgagcg  <  2:508787/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
TAAGGCGTTTACGCCGCATCCGGCAGTCATGCGTCGATGCCAAATGCGCCACCCTAAAGCAGCGCATCCGCAATGATGTACTTATTCCTTCGGTTTAAACCGCAGCAACCGGTTGGCGTTACTCACTACGGTAATCGACGAGAGCGCCATTGCCGCTCCGGCAACTACCGGGTTAAGCAGTGTTCCAGTGAACGGCCACAAAATACCGGCGGCGACCGGAATACCGATACTGTTGTAGATAAACGCACCGAGCAGGTTCTGCTTCATGTTGTGCAGCGTTGCGCGGGAAATAGCGAGCG  >  NC_000913/508794‑509092

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