Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 242,344 C→G 6.6% P320P (CCG→CCC)  fadE ← acyl coenzyme A dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913242,3440CG 6.6% 208.8 / 5.9 61P320P (CCG→CCC) fadEacyl coenzyme A dehydrogenase
Reads supporting (aligned to +/- strand):  ref base C (29/28);  new base G (2/2);  total (31/30)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.57e-01

GGCACTCCACCAGCATCCGCCAGCCTTGCCCGGCCATTTTCGGCCCGCCGATGATGTAATCGATCGGCACGAAGACATCTTTACCGCGCGTCGGTCCGTTCTGGAACGGTACGTTCAGCGGGAAGTGGCGACGACCAATTTCCACGCCCGGCGTGGTGGTTGGGATCAGCGCACAGGTAATGCCTAAATCTTCTGCACCGCCGAGTAATTTTTCCGGGTCGGAGAGTTTAAACGCCAGCCCAAGCACGGTCGCAATCGGTGCCAGCGTAATGTAGCGTTTGTTCCAG  >  NC_000913/242196‑242482
                                                                                                                                                    |                                                                                                                                          
ggCACTCCACCAGCATCCGCCAGCCTTGCCCGGCCATTTTCGGCCCGCCGATGATGTAATCGATCGGCACGAAGACATCTTTACCGCGCGTCGGTCCGTTCTGGAACGGTACGTTCAGCGGGAAGTGGCGACGACCAATTTCCACGccc                                                                                                                                            >  1:162941/1‑149 (MQ=255)
     tCCACCAGCATCCGCCAGCCTTGCCCGGCCATTTTCGGCCCGCCGATGATGTAATCGATCGGCACGAAGACATCTTTACCGCGCGTCGGTCCGTTCTGGAACGGTACGTTCAGCGGGAAGTGGCGACGACCAATTTCCACGCCCGGCgt                                                                                                                                       >  2:173319/1‑149 (MQ=255)
         ccaGCATCCGCCAGCCTTGCCCGGCCATTTTCGGCCCGCCGATGATGTAATCGATCGGCACGAAGACATCTTTACCGCGCGTCGGTCCGTTCTGGAACGGTACGTTCAGCGGGAAGTGGCGACGACCAATTTCCACGCCCGGCgtggtg                                                                                                                                   <  2:182484/149‑1 (MQ=255)
          caGCATCCGCCAGCCTTGCCCGGCCATTTTCGGCCCGCCGATGATGTAATCGATCGGCACGAAGACATCTTTACCGCGCGTCGGTCCGTTCTGGAACGGTACGTTCAGCGGGAAGTGGCGACGACCAATTTCCACGCCCGGCGtggtgg                                                                                                                                  >  3:12639/1‑149 (MQ=255)
           aGCATCCGCCAGCCTTGCCCGGCCATTTTCGGCCCGCCGATGATGTAATCGATCGGCACGAAGACATCTTTACCGCGCGTCGGTCCGTTCTGGAACGGTACGTTCAGCGGGAAGTGGCGACGACCAATTTCCACGCCGGGCgtggaaat                                                                                                                                 <  4:126280/149‑5 (MQ=255)
            gCATCCGCCAGCCTTGCCCGGCCATTTTCGGCCCGCCGATGATGTAATCGATCGGCACGAAGACATCTTTACCGCGCGTCGGTCCGTTCTGGAACGGTACGTTCAGCGGGAAGTGGCGACGACCAATTTCCACGCCCGGCGTGGTGGtt                                                                                                                                <  4:267284/149‑1 (MQ=255)
            gCATCCGCCAGCCTTGCCCGGCCATTTTCGGCCCGCCGATGATGTAATCGATCGGCACGAAGACATCTTTACCGCGCGTCGGTCCGTTCTGGAACGGTACGTTCAGCGGGAAGTGGCGACGACCAATTTCCACGCCCGGCGTGGTGGtt                                                                                                                                >  2:347799/1‑149 (MQ=255)
              aTCCGCCAGCCTTGCCCGGCCATTTTCGGCCCGCCGATGATGTAATCGATCGGCACGAAGACATCTTTACCGCGCGTCGGTCCGTTCTGGAACGGTACGTTCAGCGGGAAGTGGCGACGACCAATTTCCACGCCCGGCGTGGTGGTTgg                                                                                                                              <  1:287237/149‑1 (MQ=255)
                 cGCCAGCCTTGCCCGGCCATTTTCGGCCCGCCGATGATGTAATCGATCGGCACGAAGACATCTTTACCGCGCGTCGGTCCGTTCTGGAACGGTACGTTCAGCGGGAAGTGGCGACGACCAATTTCCACGCCCGGCGTGGTGGTTGGGAt                                                                                                                           <  3:210190/149‑1 (MQ=255)
                             ccGGCCATTTTCGGCCCGCCGATGATGTAATCGATCGGCACGAAGACATCTTTACCGCGCGTCGGTCCGTTCTGGAACGGTACGTTCAGCGGGAAGTGGCGACGACCAATTTCCACGCCCGGCGTGGTGGTTGGGATCAGCGCACAGGt                                                                                                               >  3:4647/1‑149 (MQ=255)
                                gCCATTTTCGGCCCGCCGATGATGTAATCGATCGGCACGAAGACATCTTTACCGCGCGTCGGTCCGTTCTGGAACGGTACGTTCAGCGGGAAGTGGCGACGACCAATTTCCACGCCCGGCGTGGTGGTTGGGATCAGCGCACAGGTAAt                                                                                                            >  1:90280/1‑149 (MQ=255)
                                gCCATTTTCGGCCCGCCGATGATGTAATCGATCGGCACGAAGACATCTTTACCGCGCGTCGGTCCGTTCTGGAACGGTACGTTCAGCGGGAAGTGGCGACGACCAATTTCCACGCCCGGCGTGGTGGTTGGGATCAGCGCACAGGTAAt                                                                                                            <  1:422187/149‑1 (MQ=255)
                                  cATTTTCGGCCCGCCGATGATGTAATCGATCGGCACGAAGACATCTTTACCGCGCGTCGGTCCGTTCTGGAACGGTACGTTCAGCGGGAAGTGGCGACGACCAATTTCCACGCCCGGCGTGGTGGTTGGGATCAGCGCACAGGTAATGc                                                                                                          >  1:322372/1‑149 (MQ=255)
                                                 gatgatGTAATCGATCGGCACGAAGACATCTTTACCGCGCGTCGGTCCGTTCTGGAACGGTACGTTCAGCGGGAAGTGGCGACGACCAATTTCCACGCCCGGCGTGGTGGTTGGGATCAGCGCACAg                                                                                                                 >  1:403073/1‑127 (MQ=255)
                                                 gatgatGTAATCGATCGGCACGAAGACATCTTTACCGCGCGTCGGTCCGTTCTGGAACGGTACGTTCAGCGGGAAGTGGCGACGACCAATTTCCACGCCCGGCGTGGTGGTTGGGATCAGCGCACAg                                                                                                                 <  2:403073/127‑1 (MQ=255)
                                                      tgTAATCGATCGGCACGAAGACATCTTTACCGCGCGTCGGTCCGTTCTGGAACGGTACGTTCAGCGGGAAGTGGCGACGACCAATTTCCACGCCCGGCGTGGTGGTTGGGATCAGCGCACAGGTAATGCCTAAATCTTCTGCAccgccg                                                                                      >  1:266036/1‑149 (MQ=255)
                                                       gTAATCGATCGGCACGAAGACATCTTTACCGCGCGTCGGTCCGTTCTGGAACGGTACGTTCAGCGGGAAGTGGCGACGACCAATTTCCACGCCCGGCGTGGTGGtt                                                                                                                                <  4:83755/106‑1 (MQ=255)
                                                       gTAATCGATCGGCACGAAGACATCTTTACCGCGCGTCGGTCCGTTCTGGAACGGTACGTTCAGCGGGAAGTGGCGACGACCAATTTCCACGCCCGGCGTGGTGGtt                                                                                                                                >  3:83755/1‑106 (MQ=255)
                                                       gTAATCGATCGGCACGAAGACATCTTTACCGCGCGTCGGTCCGTTCTGGAACGGTACGTTCAGCGGGAAGTGGCGACGACCAATTTCCACGCCCGGCGTGGTGGTTGGGATCAGCGCACAGGTAATGCCTAAATCTTCTGCACCGCCGa                                                                                     <  2:356109/149‑1 (MQ=255)
                                                       gTAATCGATCGGCACGAAGACATCTTTACCGCGCGTCGGTCCGTTCTGGAACGGTACGTTCAGCGGGAAGTGGCGACGACCAATTTCCACGCCCGGCGTGGTGGTTGGGATCAGCGCACAGGTAATGCCTAAATCTTCTGCACCGCCGa                                                                                     >  4:246165/1‑149 (MQ=255)
                                                        tAATCGATCGGCACGAAGACATCTTTACCGCGCGTCGGTCCGTTCTGGAACGGTACGTTCAGCGGGAAGTGGCGACGACCAATTTCCACGCCCGGCGTGGTGGTTGGGATCAGCGCACAGGTAATGCCTAAATCTTCTGCACCGCCGAg                                                                                    >  1:175099/1‑149 (MQ=255)
                                                        tAATCGATCGGCACGAAGACATCTTTACCGCGCGTCGGTCCGTTCTGGAACGGTACGTTCAGCGGGAAGTGGCGACGACCAATTTCCACGCCCGGCGTGGTGGTTGGGATCAGCGCACAGGTAATGCCTAAATCTTCTGCACCGCCGAg                                                                                    <  3:246165/149‑1 (MQ=255)
                                                              atcgGCACGAAGACATCTTTACCGCGCGTCGGTCCGTTCTGAAACGGTACGTTCAGCGGGAAGTGGCGACGACCAATTTCCACGCCCGGCGTGGTGGTTGGGATCAGCGCACAGGTAATGCCTAAATCTTCTGCACCGCCGAGTAAttt                                                                              >  3:310817/1‑149 (MQ=255)
                                                               tcgGCACGAAGACATCTTTACCGCGCGTCGGTCCGTTCTGGAACGGTACGTTCAGCGGGAAGTGGCGACGACCAATTTCCACGCCCGGCGTGGTGGTTGGGATCAGCGCACAGGTAATGCCTAAATCTTCTGCACCGCCGAGTAAtttt                                                                             >  1:361802/1‑149 (MQ=255)
                                                               tcgGCACGAAGACATCTTTACCGCGCGTCGGTCCGTTCTGGAACGGTACGTTCAGCGGGAAGTGGCGACGACCAATTTCCACGCCCGGCGTGGTGGTTGGGATCAGCGCACAGGTAATGCCTAAATCTTCTGCACCGCCGAGTAAtttt                                                                             <  3:187564/149‑1 (MQ=255)
                                                                     cGAAGACATCTTTACCGCGCGTCGGTCCGTTCTGGAACGGTACGTTCAGCGGGAAGTGGCGACGACCAATTTCCACGCCCGGCGTGGTGGTTGGGATCAGCGCACAGGTAATGCCTAAATCTTCTGCACCGCCGAGTAATTTTTCCggg                                                                       >  1:125873/1‑149 (MQ=255)
                                                                      gAAGACATCTTTACCGCGCGTCGGTCCGTTCTGGAACGGTACGTTCAGCGGGAAGTGGCGACGACCAATTTCCACGCCCGGCGTGGTGGTTGGGATCAGCGCACAGGTAATGCCTAAATCTTCTGCACCGCCGAGTAATTTTTCCGGGt                                                                      <  1:173319/149‑1 (MQ=255)
                                                                            aTCTTTACCGCGCGTCGGTCCGTTCTGGAACGGTACGTTCAGCGGGAAGTGGCGACGACCAATTTCCACGCCCGGCGTGGTGGTTGGGATCAGCGCACAGGTaa                                                                                                             <  4:360422/104‑1 (MQ=255)
                                                                            aTCTTTACCGCGCGTCGGTCCGTTCTGGAACGGTACGTTCAGCGGGAAGTGGCGACGACCAATTTCCACGCCCGGCGTGGTGGTTGGGATCAGCGCACAGGTaa                                                                                                             >  3:360422/1‑104 (MQ=255)
                                                                            aTCTTTACCGCGCGTCGGTCCGTTCTGGAACGGTACGTTCAGCGGGAAGTGGCGACGACCAATTTCCACGCCCGGCGTGGTGGTTGGGATCAGCGCACAGGTAATGCCTAAATCTTCTGCACCGCCGAGTAATTTTTCCGGGTCGgaga                                                                <  2:266036/149‑1 (MQ=255)
                                                                            aTCTTTACCGCGCGTCGGTCCGTTCTGGAACGGTACGTTCAGCGGGAAGTGGCGACGACCAATTTCCACGCCCGGCGTGGTGGTTGGGATCAGCGCACAGGTAATGCCTAAATCTTCTGCACCGCCGAGTAATTTTTCCGGGTCGgaga                                                                >  2:308866/1‑149 (MQ=255)
                                                                                 tACCGCGCGTCGGTCCGTTCTGGAACGGTACGTTCAGCGGGAAGTGGCGACGACCAATTTCCACGCCCGGCGTGGTGGTTGGGATCAGCGCACAGGTAATGCCTAAATCTTCTGCACCGCCGAGTAATTTTTCCGGGTCGGAGAGTTTa                                                           <  2:233985/149‑1 (MQ=255)
                                                                                    cgcgcgTCGGTCCGTTCTGGAACGGTACGTTCAGCGGGAAGTGGCGACGACCAATTTCCACGCCCGGCGTGGTGGTTGGGATCAGCGCACAGGTAATGCCTAAATCTTCTGCACCGCCGAGTAATTTTTCCGGGTCGGAGAGTTTAAAc                                                        <  2:175099/149‑1 (MQ=255)
                                                                                      cgcgTCGGTCCGTTCTGGAACGGTACGTTCAGCGGGAAGTGGCGACGACCAATTTCCACGCCCGGCGTGGTGGTTGGGATCAGCGCACAGGTAATGCCTAAATCTTCTGCACCGCCGAGTAATTTTTCCGGGTCGGAGAGTTTAAACGc                                                      <  2:322372/149‑1 (MQ=255)
                                                                                       gcgTCGGTCCGTTCTGGAACGGTACGTTCAGCGGGAAGTGGCGACGACCAATTTCCACGCCCGGCGTGGTGGTTGGGATCAGCGCACAGGTAATGCCTAAATCTTCTGCACCGCCGAGTAATTTTTCCGGGTCGGAGAGTTTAAACGcc                                                     <  4:427118/149‑1 (MQ=255)
                                                                                              tCCGTTCTGGAACGGTACGTTCAGCGGGAAGTGGCGACGACCAATTTCCACGCCCGGCGTGGTGGTTGGGATCAGCGCACAGGTAATGCCTAAATCTTCTGCACCGCCGAGTAATTTTTCCGGGTCGGAGAGTTTAAACGCCAGCCCaa                                              <  2:181755/149‑1 (MQ=255)
                                                                                              tCCGTTCTGGAACGGTACGTTCAGCGGGAAGTGGCGACGACCAATTTCCACGCCCGGCGTGGTGGTTGGGATCAGCGCACAGGGAATGCCTAAATCTTCTGCACCGCCGAGTAATTTTTCCGGGTCGGAGAGTTTAAACGCCAGCCCaa                                              >  2:107619/1‑149 (MQ=255)
                                                                                                      ggAACGGTACGTTCAGCGGGAAGTGGCGACGACCAATTTCCACGCCCGGCGTGGTGGTTGGGATCAGCGCACAGGTAATGCCTAAATCTTCTGCACCGCCGAGTAATTTTTCCGGGTCGGAGAGTTTAAACGCCAGCCCAAGCACGGTc                                      <  2:162941/149‑1 (MQ=255)
                                                                                                      ggAACGGTACGTTCAGCGGGAAGTGGCGACGACCAATTTCCACGCCCGGCGTGGTGGTTGGGATCAGCGCACAGGTAATGCCTAAATCTTCTGCACCGCCGAGTAATTTTTCCGGGTCGGAGAGTTTAAACGCCAGCCCAAGCACGGTc                                      >  4:332761/1‑149 (MQ=255)
                                                                                                      ggAACGGTACGTTCAGCGGGAAGTGGCGACGACCAATTTCCACGCCCGGCGTGGTGGTTGGGATCAGCGCACAGGTAATGCCTAAATCTTCTGCACCGCCGAGTAATTTTTCCGGGTCGGAGAGTTTAAACGCCAGCCCAAGCACGGTc                                      <  4:4647/149‑1 (MQ=255)
                                                                                                            gTACGTTCAGCGGGAAGTGGCGACGACCATTTTCCACGCCCGGCGTGGTGGTTGGGATCAGCGCACAGGTAATGCCTAAATCTTCTGCACCGCCGAGTAATCTTTCCGGGTCGGAGGGTTTAAACGCCAGCCCAAGCACTGTCGCAATc                                >  3:100235/1‑149 (MQ=255)
                                                                                                               cGTTCAGCGGGAAGTGGCGACGACCAATTTCCACGCCCGGCGTGGTGGTTGGGATCAGCGCACAGGTAATGCCTAAATCTTCTGCACCGCCGAGTAATTTTTCCGGGTCGGAGAGTTTAAACGCCAGCCCAAGCACGGTCGCAATCGGt                             <  2:90280/149‑1 (MQ=255)
                                                                                                                  tCAGCGGGAAGTGGCGACGACCAATTTCCACGCCCGGCGTGGTGGTTGGGATCAGCGCACAGGTAATGCCTAAATCTTCTGCACCGCCGAGTAATTTTTCCGGGTCGGAGAGTTTAAACGCCAGCCCAAGCACGGTCGCAATCGGTGcc                          >  4:88664/1‑149 (MQ=255)
                                                                                                                  tCAGCGGGAAGTGGCGACGACCAATTTCCACGCCCGGCGTGGTGGTTGGGATCAGCGCACAGGTAATGCCTAAATCTTCTGCACCGCCGAGTAATTTTTCCGGGTCGGAGAGTTTAAACGCCAGCCCAAGCACGGTCGCAATCGGTGcc                          >  1:240370/1‑149 (MQ=255)
                                                                                                                             tGGCGACGACCAATTTCCACGCCCGGCGTGGTGGTTGGGATCAGCGCACAGGTAATGCCTAAATCTTCTGCACCGCCGAGTAATTTTTCCGGGTCGGAGAGTTTAAACGCCAGCCCAAGCACGGTCGCAATCGGTGCCAGCGTAATGTa               <  3:200286/149‑1 (MQ=255)
                                                                                                                                 gacgacCAATTTCCACGCCCGGCGTGGTGGTTGGGATCAGCGCACAGGTAATGCCTAAATCTTCTGCACCGCCGAGTAATTTTTCCGGGTCGGAGAGTTTAAACGCCAGCCCAAGCACGGTCGCAATCGGTGCCAGCGTAATGTAGCGt           >  2:141612/1‑149 (MQ=255)
                                                                                                                                    gacCAATTTCCACGCCCGGCGTGGTGGTTGGGATCAGCGCACAGGTAATGCCTAAATCTTCTGCACCGCCGAGTAATTTTTCCGGGTCGGAGAGTTTAAACGCCAGCCCAAGCACGGTCGCAATCGGt                             <  3:102306/128‑1 (MQ=255)
                                                                                                                                    gacCAATTTCCACGCCCGGCGTGGTGGTTGGGATCAGCGCACAGGTAATGCCTAAATCTTCTGCACCGCCGAGTAATTTTTCCGGGTCGGAGAGTTTAAACGCCAGCCCAAGCACGGTCGCAATCGGt                             >  4:102306/1‑128 (MQ=255)
                                                                                                                                    gacCAATTTCCACGCCCGGCGTGGTGGTTGGGATCAGCGCACAGGTAATGCCTAAATCTTCTGCACCGCCGAGTAATTTTTCCGGGTCGGAGAGTTTAAACGCCAGCCCAAGCACGGTCGCAATCGGTGCCAGCGTAATGTAGCGTTTg        >  2:189849/1‑149 (MQ=255)
                                                                                                                                     acCAATTTCCACGCCCGGCGTGGTGGTTGGAATCAGCGCACAGGTAATGCCTAAATCTTCTGCACCGCCGAGTAATTTTTCCGGGTCGGAGAGTTTAAACGCCAGCCCAAGCACGGTCGCAATCGGTGCCAGCGTAATGTAGCGTTTGt       <  1:308866/149‑1 (MQ=255)
                                                                                                                                       gtcggtCCACGCCGGGCGTGGTGGTTGGGATCAGCGCACAGGTAATGCCTAAATCTTCTGCACCGCCGAGTAATTTTTCCGGGTCGGAGAGTTTAAACGCCAGCCCAAGCACGGTCGCAATCGGTGCCAGCGTAATGTAGCGTTTGTTc     <  4:317893/144‑1 (MQ=255)
                                                                                                                                       cAATTTCCACGCCCGGCGTGGTGGTTGGGATCAGCGCACAGGTAATGCCTAAATCTTCTGCACCGCCGAGTAATTTTTCCGGGTCGGAGAGTTTAAACGCCAGCCCAAGCACGGTCGCAATCGGTGCCAGCGTAATGTAGCGTTTGTTc     >  2:213162/1‑149 (MQ=255)
                                                                                                                                       cAATTTCCACGCCCGGCGTGGTGGTTGGGATCAGCGCACAGGTAATGCCTAAATCTTCTGCACCGCCGAGTAATTTTTCCGGGTCGGAGAGTTTAAACGCCAGCCCAAGCACGGTCGCAATCGGTGCCAGCGTAATGTAGCGTTTGTTc     >  4:285755/1‑149 (MQ=255)
                                                                                                                                       cAATTTCCACGCCCGGCGTGGTGGTTGGGATCAGCGCACAGGTAATGCCTAAATCTTCTGCACCGCCGAGTAATTTTTCCGGGTCGGAGAGTTTAAACGCCAGCCCAAGCACGGTCGCAATCGGTGCCAGCGTAATGTAGCGTTTGTTc     >  2:343011/1‑149 (MQ=255)
                                                                                                                                       cAATTTCCACGCCCGGCGTGGTGGTTGGGATCAGCGCACAGGTAATGCCTAAATCTTCTGCACCGCCGAGTAATTTTTCCGGGTCGGAGAGTTTAAACGCCAGCCCAAGCACGGTCGCAATCGGTGCCAGCGTAATGTAGCGTTTGTTc     >  1:15750/1‑149 (MQ=255)
                                                                                                                                        aaccaccACGCCGGGCGTGGTGGTTGGGATCAGCGCACAGGTAATGCCTAAATCTTCTGCACCGCCGAGTAATTTTTCCGGGTCGGAGAGTTTAAACGCCAGCCCAAGCACGGTCGCAATCGGTGCCAGCGTAATGTAGCGTTTGTTcc    >  1:240849/6‑149 (MQ=255)
                                                                                                                                         accaccACGCCGGGCGTGGTGGTTGGGATCAGCGCACAGGTAATGCCTAAATCTTCTGCACCGCCGAGTAATTTTTCCGGGTCGGAGAGTTTAAACGCCAGCCCAAGCACGGTCGCAATCGGTGCCAGCGTAATGTAGCGTTTGTTCCa   >  2:248987/5‑149 (MQ=255)
                                                                                                                                         aTTTCCACGCCCGGCGTGGTGGTTGGGATCAGCGCACAGGTAATGCCTAAATCTTCTGCACCGCCGAGTAATTTTTCCGGGTCGGAGAGTTTAAACGCCAGCCCAAGCACGGTCGCAATCGGTGCCAGCGTAATGTAGCGTTTGTTCCa   <  1:347799/149‑1 (MQ=255)
                                                                                                                                          tttCCACGCCCGGCGTGGTGGTTGGGATCAGCGCACAGGTAATGCCTAAATCTTCTGCACCGCCGAGTAATTTTTCCGGGTCGGAGAGTTTAAACGCCAGCCCAAGCACGGTCGCAATCGGTGCCAGCGTAATGTAGCGTTTGTTCCAg  <  4:12639/149‑1 (MQ=255)
                                                                                                                                          tttACACGCCCGGCGTGGTGGTTGGGATTAGCGCACCGGTAATGCCTAAATCTTCTGCACCGCCGAGTAATTTTTCCGGGTCGGAGAGTTTAAACGCCAGCCCAAGCGCGGTCGCAATCGGTGCCAGCGTAATGTAGCGTTTGTTCCAg  <  4:100235/149‑1 (MQ=255)
                                                                                                                                             ccACGCCCGGCGTGGTGGTT‑GGATCAGCGCACAGGTAATGCCTAAATCTTCTGCACCGCCGAGTAATTTTTCCGGGTCGGAGAGTTTa                                                           <  3:342138/88‑1 (MQ=255)
                                                                                                                                             ccACGCCCGGCGTGGTGGTT‑GGATCAGCGCACAGGTAATGCCTAAATCTTCTGCACCGCCGAGTAATTTTTCCGGGTCGGAGAGTTTa                                                           >  4:342138/1‑88 (MQ=255)
                                                                                                                                                    |                                                                                                                                          
GGCACTCCACCAGCATCCGCCAGCCTTGCCCGGCCATTTTCGGCCCGCCGATGATGTAATCGATCGGCACGAAGACATCTTTACCGCGCGTCGGTCCGTTCTGGAACGGTACGTTCAGCGGGAAGTGGCGACGACCAATTTCCACGCCCGGCGTGGTGGTTGGGATCAGCGCACAGGTAATGCCTAAATCTTCTGCACCGCCGAGTAATTTTTCCGGGTCGGAGAGTTTAAACGCCAGCCCAAGCACGGTCGCAATCGGTGCCAGCGTAATGTAGCGTTTGTTCCAG  >  NC_000913/242196‑242482

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑