Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,340,460 A→T 11.1% intergenic (+114/+28) PP_3812 → / ← PP_5632 lipase/allo‑isoleucine metabolism protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,340,4600AT11.1% 321.3 / 28.7 117intergenic (+114/+28)PP_3812/PP_5632lipase/allo‑isoleucine metabolism protein
Reads supporting (aligned to +/- strand):  ref base A (51/53);  new base T (7/6);  total (58/59)
Fisher's exact test for biased strand distribution p-value = 7.77e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.62e-01

Only 100 of 118 total aligned reads displayed.

GGCGGGCAACGAACTTAAAGTGCATCTGCAGTGAGGCGATGAGCAATAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTTCGCGGGTAAACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCCACAATGATCACCTCAGCCTTCGCCACCCCTTGGCGATGGCTGCATGCCTCTTGGTACAGCGCGGAGCGGTAACAGGCGACGGCCTGCGCGTACGAATCGAAC  >  NC_002947/4340313‑4340608
                                                                                                                                                    |                                                                                                                                                    
ggCGGGCAACGAACTTAAAGTGCATCTGCAGTGAGGCGATGAGCAATAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTTCGCGGGTAAAc                                                                                                                                                    <  1:1761380/150‑1 (MQ=255)
   gggCAACGAACTTAAAGTGCATCTGCAGTGAGGCGATGAGCAATAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTTCGCGGGTAAACCcg                                                                                                                                                 <  1:2892146/150‑1 (MQ=255)
    ggCAACGAACTTAAAGTGCATCTGCAGTGAGGCGATGAGCAATAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTTCGCGGGTAAACCcgc                                                                                                                                                <  2:4743054/150‑1 (MQ=255)
    ggCAACGAACTTAAAGTGCATCTGCAGTGAGGCGATGAGCAATAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTTCGCGGGTAAACCcgc                                                                                                                                                <  2:2625284/150‑1 (MQ=255)
    ggCAACGAACTTAAAGTGCATCTGCAGTGAGGCGATGAGCAATAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTTCGCGGGTAAACCcgc                                                                                                                                                <  2:3378032/150‑1 (MQ=255)
       aaCGAACTTAAAGTGCATCTGCAGTGAGGCGATGAGCAATAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTTCGCGGGTTTACCCGCGaa                                                                                                                                             >  1:2371606/1‑150 (MQ=255)
       aaCGAACTTAAAGTGCATCTGCAGTGAGGCGATGAGCAATAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTTCGCGGGTTTAACCGCGaa                                                                                                                                             >  1:2371611/1‑150 (MQ=255)
       aaCGAACTTAAAGTGCATCTGCAGTGAGGCGATGAGCAATAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTTCGCGGGTAAACCCGCGaa                                                                                                                                             <  1:1129053/150‑1 (MQ=255)
        aCGAACTTAAAGTGCATCTGCAGTGAGGCGATGAGCAATAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTTCGCGGGTTTACCCGCGaaa                                                                                                                                            <  1:4575024/150‑1 (MQ=255)
        aCGAACTTAAAGTGCATCTGCAGTGAGGCGATGAGCAATAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTTCGCGGGTAAACCCGCGaaa                                                                                                                                            <  2:314853/150‑1 (MQ=255)
         cGAACTTAAAGTGCATCTGCAGTGAGGCGATGAGCAATCGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCGGCCAGCGGTGCGGCAGGTCCGGCCCTTTCGCGGGTTTACCCGCGAAAg                                                                                                                                           <  2:444048/150‑1 (MQ=255)
          gAACTTAAAGTGCATCTGCAGTGAGGCGATGAGCAATAGACCCGCGTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTTCGCGGGTTTACCCGCGAAAgg                                                                                                                                          <  2:1488017/150‑1 (MQ=255)
          gAACTTAAAGTGCATCTGCAGTGAGGCGATGAGCAATAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTTCGCGGGTTTACCCGCGAAAgg                                                                                                                                          <  2:3537193/150‑1 (MQ=255)
          gAACTTAAAGTGCATCTGCAGTGAGGCGATGAGCAATAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTTCGCGGGGAAACCCCCGAAAgg                                                                                                                                          >  2:2364454/1‑150 (MQ=255)
          gAACTTAAAGTGCATCTGCAGTGAGGCGATGAGCAAGAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCGTTCGCGGGTTTACCCGCGAAAgg                                                                                                                                          <  2:4188439/150‑1 (MQ=255)
            aCTTAAAGTGCATCTGCAGTGAGGCGATGAGCAATAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGGCCGTTGGCCGGGGGAACCCCCGAAgggc                                                                                                                                         >  1:487585/1‑149 (MQ=255)
            aCTTAAAGTGCATCTGCAGTGAGGCGATGAGCAATAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTTCGCGGGTTTACCCGCGAAAGGGc                                                                                                                                        <  1:2955807/150‑1 (MQ=255)
                aaaGTGCATCTGCAGTGAGGCGATGAGCAATAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCGTTCGCGGGGGAATCCCCCAAAGGGCCCGc                                                                                                                                    >  1:2216686/1‑150 (MQ=255)
                       aTCTGCAGTGAGGCGATGAGCAATAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTTCGCGGGTAAACCCGCGAAAGGGCCGGCAAATCCt                                                                                                                             >  2:3100504/1‑150 (MQ=255)
                         ctgctgTGTGGCTATGAGCATTAGACCCGCGTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCTGGCTTCAATTAATCTGCCAGCGGTGCGGCTGGTCGGGCCCTTTGGGGGTTAAACCCGCGAAAGGGCCGGCAAATCCTGc                                                                                                                           <  1:1571862/150‑1 (MQ=255)
                          tGCAGTGAGGCGATGAGCAATAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCGTTCGCGGGGGAAACCCCGCAAGGGCCGGCCAAACCAGCt                                                                                                                          >  1:2792956/1‑150 (MQ=255)
                             cgTGAGGCGATGAGCAATAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTAAATTAATCTGCCAGCGGTGCGGCTGGTCCGGCCCTTTCGCGGGTAAACCCGCGAAAGGGCCGGCAAATCCTGCTCAg                                                                                                                       <  1:421035/149‑1 (MQ=255)
                                   gCGATGAGCAATAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTTCGCGGGTAAACCCGCGAAAGGGCCGGCAAATCCAGCTCAGGGCTCa                                                                                                                 >  1:2528390/1‑150 (MQ=255)
                                     gATGAGCACTAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACCGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTTCGCGGGTAAACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCaaa                                                                                                               <  2:4144905/150‑1 (MQ=255)
                                      aTGAGCAATAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTTCGCGGGTAAACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAAc                                                                                                              <  2:1993115/150‑1 (MQ=255)
                                           cAATAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTTCGCGGGTAAACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTc                                                                                                         <  2:2184160/150‑1 (MQ=255)
                                             aTAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCGTTCGCGGGTAAACCCGCGCAAGGGCCGGCAAAACCAGCTCAGGGCTCAAACCCTTcag                                                                                                       >  2:4816636/1‑148 (MQ=255)
                                                gACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTTCGCGGGTAAACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCATCCACaa                                                                                                    >  2:4885315/1‑150 (MQ=255)
                                                  cccGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACCGCCCGCAGGGTTCAATGATTCTGCCAGCGGTGCGGCTGGTCCGGCCCTTTCGCGGGTAAACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCCACAATg                                                                                                  <  2:4892241/150‑1 (MQ=255)
                                                   ccGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTTCGCGGGTAAACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCCACAATGa                                                                                                 >  2:2821761/1‑150 (MQ=255)
                                                    cgcgTGGGGGGGTCTTTTGCTTCTTTTTGCGGAAACAGCCCGCGGGCTTCAATGAATCTGCCAGCGGTGCGGCTGGTCCGGCCCTTTCGCGGGTAAACCCGCGAATGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCCACAATGAt                                                                                                <  1:1222847/150‑1 (MQ=255)
                                                           gggtttattttGCTTCTTTTTGCGTAAACCGCCCGCAGGATTCATTTATTCTGCCCGGGGTGCGGCTGGTCCGGCCCTTTCGGGGGTAAACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCCACAATGATCACCTCa                                                                                         <  1:915710/143‑1 (MQ=255)
                                                                 cTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTTCGCGGGTAAACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCCACAATGATCACCTCAGCCTTc                                                                                   <  2:688978/150‑1 (MQ=255)
                                                                  ttttGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTTCGCGGGGAAACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCGCAAACCCGTCCACAAAGATCACCTCAGCCATCg                                                                                  >  2:710820/1‑150 (MQ=255)
                                                                  ttttGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTGTCGCGGGTAAACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAAACCTTCCACAATGATCACCTCAGCCTTCg                                                                                  >  2:2750099/1‑150 (MQ=255)
                                                                   tttGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTTCGCGGGTAAACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCCACAATGATCACCTCAGCCTTCGc                                                                                 <  2:4929373/150‑1 (MQ=255)
                                                                      gCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTTCGCGGGTAAACCCGCGAAAGGGCCCGCAAAAACAGCTCACGGCTCAAAACCGTCCACAAAGATCACCACAGCCATAGCCAc                                                                              >  2:1906107/1‑150 (MQ=255)
                                                                      gCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTTCGCGGGTAAACCCCCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCCACAATGATAACCTCAGCCTTAGCCAc                                                                              >  2:587875/1‑150 (MQ=255)
                                                                        ttCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTTCGCGGGTAAACCCGCGAAAGGGGCGGCAAAACCTGCTCCGGGCTCAAACCCTTCCACAAATTTCACCTCAGCCGTCGCCAccc                                                                            >  1:3812943/1‑150 (MQ=255)
                                                                        ttCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTTCGCGGGTAAACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCCACAATGATCACCTCAGCCTTCGCCAccc                                                                            <  2:317381/150‑1 (MQ=255)
                                                                        ttCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCATTCGCGGGTTAACCCGCGGAAGGGCCGGCAAATACAGCTCAGGGCTCAAAACCGAACACAAAGATCACCTCAGCCAACGCCAccc                                                                            >  1:423265/1‑150 (MQ=255)
                                                                        tccTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTTCGCGGGTAAACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCCACAATGATCACCTCAGCCTTCGCCAccc                                                                            <  2:4524647/148‑1 (MQ=255)
                                                                           tATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTTCGCGGGTAAACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCCACAATGATCACCTCAGCCTTCGCCAccccgt                                                                         >  2:4133997/1‑148 (MQ=255)
                                                                           tATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTTCGCGGGTAAACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCCACAATGATCACCTCAGCCTTCGCCACCCCtt                                                                         <  1:1665735/150‑1 (MQ=255)
                                                                             tttGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTTCGCGGGTAAACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCCACAATGATCACCTCAGCCGTCGCCACCCCTTgg                                                                       >  2:1547442/1‑150 (MQ=255)
                                                                             tttGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCGTTCGCGGGGTAACCCGCGAAAGGGCCGGCAAATCCAGCTCAGGGCGCAAACCCATCCACAAAGATCACCACACCCGTCGCCACCCCATgg                                                                       >  2:172882/1‑150 (MQ=255)
                                                                             tttGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCATTCGCGGGTAAACCCGCGAAAGGGCCGGCAAAACCTGCTCAGGGCTCAAAACCATCCACCAAGATCAACACAGCCATCGCCACCCCATgg                                                                       >  1:1577632/1‑150 (MQ=255)
                                                                              tggCGTAAACAGCCGGCAGGCTTCAATGAATCTGCCCGGGGTGCGGCCGGCCCGGCCCTTGGGCGGGTAAACCCGCGAACGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCCACAATGATCACCTCAGCCTTCGCCACCCCTTGGc                                                                      <  1:3934994/148‑1 (MQ=255)
                                                                                  ggAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTTCGCGGGTAAACCCGCGAAAGGGCCGGCCAATCCTGCTCAGGGCTCAAACCCTTCCACAATGATCACCGCAGCCTTCGCCACCCCTTGGCGATg                                                                  >  2:781937/1‑150 (MQ=255)
                                                                                  ggAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTTCGCGGGTAAACCCCCGAAAGGGCCGGCAAATACTGCTCAGGGCTCAAAACCTTACACAATGATAACCTCAGCCTTCGCCACCCCATGGCgcag                                                                  >  2:4420756/1‑147 (MQ=255)
                                                                                  ggAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCATTCGCGGGGAAAACCCCGAAAGGGCCGGCCAAAACGGCGCAGGGCTCAAAACCCTCCACAAAGTTAACCACAGCCATCGCCACCCCCTGGCGGTg                                                                  >  1:3681271/1‑150 (MQ=255)
                                                                                   gAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTTCGCGGGTAAACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCCACAATGATCACCTCAGCCTTCGCCACCCCTTGGCGATgg                                                                 <  2:5006402/150‑1 (MQ=255)
                                                                                        aGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCGTTCGCGGGTAAACCCGCGAAAGGGCCGGCCAATCCTGCTCAGGGCGCAAACCCTTCCACAAAGATCACCGCAGCCGGCGCCACCCCATGGCGAGGGCTGCa                                                            >  2:2342314/1‑150 (MQ=255)
                                                                                             gCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCGTTCGCGGGTAAACCCGCGAAAGGGCCCGCAAATCCTGCTCAGGGCGCAAAACCGTCCACCAGGGTCACCGCAGCCTTCGCCCCCCCGTGGCGATtgctgccggccg                                                       >  2:4500210/1‑144 (MQ=255)
                                                                                                ggCTTCAATGAATGTGCCAGGGGTGCGGCAGGTCCGGCCCTTTCGCGGGAAAACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCCACAATGATCACCTCAGCCTTCGCCACCCCTTGGCGATGGCTGCATGCCTCtt                                                    <  1:897423/150‑1 (MQ=255)
                                                                                                ggCTTCAATGAAGCGGCCAGCGGTGCGGCCGGTCCGGCCCTGTCGCGGGTAAACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCCACAATGATCACCTCAGCCTTCGCCACCCCTTGGCGATGGCTGCATGCCTCtt                                                    <  1:3934671/150‑1 (MQ=255)
                                                                                                 gCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTTCGCGGGTAAACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCCACAATGATCACCTCAGCCTTCGCCACCCCTTGGCGATGGCTGCATGCCTCTTg                                                   <  2:819121/150‑1 (MQ=255)
                                                                                                 gCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTTCGCGGGTAAACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCCACAATGATCACCTCAGCCTTCGCCACCCCTTGGCGATGGCTGCATGCCTCTTg                                                   <  2:265889/150‑1 (MQ=255)
                                                                                                 gCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTTCGCGGGTAAACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCCACAATGATCACCTCAGCCTTCGCCACCCCTTGGCGATGGCTGCATGCCTCTTg                                                   <  2:265879/150‑1 (MQ=255)
                                                                                                 gCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTTCGCGGGTAAACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCCACAATGATCACCTCAGCCTTCGCCACCCCTTGGCGATGGCTGCATGCCTCTTg                                                   <  1:5003245/150‑1 (MQ=255)
                                                                                                 gCTTAAATGAATCTGCCAGCGGTGCGGCTGGTCCGGCCATTTCGCGGGTAAACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCCACAATGATCACCTCAGCCTTCGCCACCCCTTGGCGATGGCTGCATGCCTCTTg                                                   <  1:4307592/150‑1 (MQ=255)
                                                                                                  cTTCAATGAAGCGGCCCGCGGTGCGGCAGGTCCGGCCCTTGCGCGGGTAAACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCCACAATGATCACCTCAGCCTTCGCCACCCCTTGGCGATGGCTGCATGCCTCTTgg                                                  <  1:3043573/150‑1 (MQ=255)
                                                                                                  cTTAAATGAATCTGCCCGGGGTGCGGCAGGTCCGGCCCTTTCGCGGTTAAACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCCACAATGATCACCTCAGCCTTCGCCACCCCTTGGCGATGGCTGCATGCCTCTTgg                                                  <  1:3043601/150‑1 (MQ=255)
                                                                                                   ttCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTTCGCGGGTAAACCCGCGaaa                                                                                                                                            >  1:3974891/1‑59 (MQ=255)
                                                                                                   ttCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTTCGCGGGTAAACCCGCGaaa                                                                                                                                            <  2:3974891/59‑1 (MQ=255)
                                                                                                   ttCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTTCGCGGGTAAACCCGCGGAAGGGCCGGCCAATCCTGCTCAGGGCGCAAACCCATCCACAATGATCACCGCAGCCGTCGCCACCCCGTGGCGATGGCTGCATGCCTCTTGGt                                                 >  2:3043601/1‑150 (MQ=255)
                                                                                                   ttCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTTCGCGGGGAAACCCGCGAAAGGGCCGGCAAATACTGCTCAGGGCTCAAACCCGTCCACAATGATAACCTCAGCCTTCGCCCCCCCTTGGCGATGGCTGCATGCctcggggt                                                 >  2:3043573/1‑145 (MQ=255)
                                                                                                   ttCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCGTTCGCGGGTAAACCCGCGAAAGGGCCGGCAAATCCTGCGCAGGGCTCAAACCCGTCCACAAAGATAACCGCAGCCATCGCCACCCCGTGGCGATGGCTGCATGCCGCGTGGt                                                 >  2:1824682/1‑150 (MQ=255)
                                                                                                   ttCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCATTCGCCGGTAAAACCGCGAAAGGGCCGGCAAATCCAGCTCAGGGCTCAAAACCGTCCACAATGATCACCTCAGCCGTCGCCACCCCATGGCGATGGCTGCATGCCGCGTGGt                                                 >  2:1824641/1‑150 (MQ=255)
                                                                                                    tCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTTCGCGGGTAAACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCCACAATGATCACCTCAGCCTTCGCCCCCCCGTGGCGATGGCTGCATGCCGCGTGGTa                                                >  2:5057683/1‑150 (MQ=255)
                                                                                                    tCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTTCGCGGGTAAACCCGCGAAAGGGCCGGCAAATCCAGCTCAGGGCTCAAAACCTTCCACAATGATCACCTCAGCCTTCGCCACCCCGTGGCGATGGCTGCATGCCTCTTGGTa                                                >  2:3680694/1‑150 (MQ=255)
                                                                                                        tcaagcTGCCCGCGGTGCGGCAGGTCCGGCCCTTGCGCGGGTAAACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCCACAATGATCACCTCAGCCTTCGCCACCCCTTGGCGATGGCTGCATGCCTCTTGGTACAgc                                            <  1:710820/145‑1 (MQ=255)
                                                                                                        tcaaTCTGCCAGCGGTGCGGCAGGTCCGGCCCTTTCGCGGGAAAACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCCACAATGATCACCTCAGCCTTCGCCACCCCTTGGCGATGGCTGCATGCCTCTTGGTACAgc                                            <  2:4084588/148‑1 (MQ=255)
                                                                                                           aTCTGCCAGCGGTGCGGCAGGTCCGGCCCTTTCGCGGGTAAACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCCACAATGATCACCTCAGCCTTCGCCACCCCTTGGCGATGGCTGCATGCCTCTTGGTACAGcgcg                                         <  2:1410874/150‑1 (MQ=255)
                                                                                                             cTGCCCGCGGTGGGGCTGGTCCGGCCCTTGCGGGGTTAACCCCGCGAACGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCCACAATGATCACCTCAGCCTTCGCCACCCCTTGGCGATGGCTGCATGCCTCTTGGTACAGCGCGGa                                       <  1:468141/150‑1 (MQ=255)
                                                                                                             cTGCCAGCGGTGCGGCCGGTCCGGCCCTTTCGCGGTTAAACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCCACAATGATCACCTCAGCCTTCGCCACCCCTTGGCGATGGCTGCATGCCTCTTGGTACAGCGCGGa                                       <  1:2364454/150‑1 (MQ=255)
                                                                                                              tGCCAGCGGTGCGGCAGGTCCGGCCCTTTGGCGGGTTAACCCGCGAAAGGGGCGGGAAATACTTCTGAGGGGTGAAACCCGTCCACAAAGATAACCTCAGCCTTAGCCACCCCGTGGCGATGGGTGGATGGCGGTTGGTTCAGCGCgggg                                      >  2:1374711/1‑148 (MQ=255)
                                                                                                              tGCCAGCGGTGCGGCAGGTCCGGCCCTTTGGCGGGTAAACCCGCGGAAGGGCCGGGAAATCCTTGTGAGGGCTCAAAACCGTACACAATGATCAACTAAGCCTTGGCCACCCCTTGGCGTTGGGTGCATGGCTGTTGGTACAGCGCgggg                                      >  2:1374706/1‑148 (MQ=255)
                                                                                                                ccAGCGGTGCGGCAGGTCCGGCCCTTTCGCGGGTAAATCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCCACAATGATCACCTCAGCCTTCGCCACCCCTTGGCGATGGCTGCATGCCTCTTGGTACAGCGCGGAGCg                                    <  2:2216686/150‑1 (MQ=255)
                                                                                                                      gTGCGGCCGGTCCGGCCCTTTCGCGGGTAAACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCCACAATGATCACCTCAGCCTTCGCCACCCCTTGGCGATGGCTGCATGCCTCTTGGTACAGCGCGGAGCGGTAACa                              <  2:1694863/150‑1 (MQ=255)
                                                                                                                      gTGCGGCAGGTCCGGCCCTTTCGCGGGTAAACCCCCGAAAGGGCCCGCAAATCCAGCTCAGGGCTAAAACCCATCCACAAAGAACACCACAGCCGTCGCCACCCCAAGGCGAAGGCTGCATGCCGCTTGGTACAGCCCGGCGCGGTAACa                              >  2:3517166/1‑150 (MQ=255)
                                                                                                                      gTGCGGCAGGGCCGGCCCTTTCGCGGGTAAACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCCACAATGATCACCTCAGCCTTCGCCACCCCTTGGCGATGGCTGCATGCCTCTTGGTACAGCGCGGAGCGGTAACa                              <  2:4765513/150‑1 (MQ=255)
                                                                                                                       tGCGGCAGGTCCGGCCCTTTCGCGGGTAAACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCCACAATGATCACCTCAGCCTTCGCCACCCCTTGGCGATGGCTGCATGCCTCTTGGTACAGCGCGGAGCGGTAACAg                             <  2:2792956/150‑1 (MQ=255)
                                                                                                                                    gCCCTTTCGCGGGTAAACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCCACAATGATCACCTCAGCCTTCGCCACCCCTTGGCGCTGGCTGCATGCCTCTTGGTACAGCGCGGAGCGGTAACAGGCGACGGCCTgcg                >  2:539461/1‑150 (MQ=255)
                                                                                                                                       cTTTCGCGGGTTTACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCCACAATGATCACCTCAGCCTTCGCCACCCCTTGGCGATGGCTGCATGCCTCTTGGTACAGCGCGGAGCGGTAACAGGCGACGGCCTGCGCGt             >  1:112285/1‑150 (MQ=255)
                                                                                                                                        tttCGCGGGTTTACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCCACAATGATCACCTCAGCCTTCGCCACCCCTTGGCGATGGCTGCATGCCTCTTGGTACAGCGCGGAGCGGTAACAGGCGACGGCCTGCGCGTa            >  1:1212946/1‑150 (MQ=255)
                                                                                                                                        tttCGCGGGTAAACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCCACAATGATCACCTCAGCCTTCGCCACCCCTTGGCGATGGCTGCATGCCTCTTGGTACAGCGCGGAGCGGTAACAGGCGACGGCCTGCGCGTa            <  2:2354073/150‑1 (MQ=255)
                                                                                                                                        tttCGCGGGTAAACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCCACAATGATCACCTCAGCCTTCGCCACCCCTTGGCGATGGCTGCATGCCTCTTGGTACAGCGCGGAGCGGTAACAGGCGACGGCCTGCGCGTa            >  2:3786065/1‑150 (MQ=255)
                                                                                                                                        tttCGCGGGTAAACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCCACAATGATCACCTCAGCCTTCGCCACCCCTTGGCGATGGCTGCATGCCTCTTGGTACAGCGCGGAGCGGTAACAGGCGACGGCCTGCGCGTa            >  2:4198579/1‑150 (MQ=255)
                                                                                                                                        tttCGCGGGTAAACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCCACAATGATCACCTCAGCCTTCGCCACCCCTTGGCGATGGCTGCATGCCGCTTGGTACAGCGCGGAGCGGTAACAGGCGACGGCCTGCGCGTa            >  2:4668520/1‑150 (MQ=255)
                                                                                                                                        tttCGCGGGTAAACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCCACAATGACCACCTCAGCCTTCGCCACCCCTTGGCGATGGCTGCATGCCTCTTGGTACAGCGCGGAGCGGTAACAGGCGACGGCCTGCGCGTa            >  2:643664/1‑150 (MQ=255)
                                                                                                                                        tttCGCGGGTAAACACGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACACTTCCACAATGATCACCTCAGCCTTCGCCACCCCTTGGAGATGGCTGCATGCCTCTTGGTCCAGCGCGGAGCGGTAACAGGCGACGGCCTGCGCGTa            >  2:1914479/1‑150 (MQ=255)
                                                                                                                                         tgcgcgGGTAAACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCCACAATGATCACCTCAGCCTTCGCCACCCCTTGGCGATGGCTGCATGCCTCTTGGTACAGCGCGGAGCGGTAACAGGCGACGGCCTGCGCGTAc           <  1:1582575/148‑1 (MQ=255)
                                                                                                                                            gcgGGTAAACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCCACAATGATCACCTCAGCCTTCGCCACCCCTTGGCGATGGCTGCATGCCTCTTGGTACAGCGCGGAGCGGTAACAGGCGACGGCCTGCGCGTACGaa        >  1:1466457/1‑150 (MQ=255)
                                                                                                                                              ggggTAACCCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCCACAATGATCACCTCAGCCTTCGCCACCCCTTGGCGATGGCTGCATGCCTCTTGGTACAGCGCGGAGCGGTAACAGGCGACGGCCTGCGCGTACGAAt       <  1:139804/149‑1 (MQ=255)
                                                                                                                                               gggTAAACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCCACAATGATCACCTCAGCCTTCGCCACCCCTTGGCGATGGCTGCATGCCTCTTGGTACAGCGCGGAGCGGTAACAGGCGACGGCCTGCGCGTACGAATc      >  1:325626/1‑150 (MQ=255)
                                                                                                                                               gggTAAACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCCACAATGATCACCTCAGCCTTCGCCACCCCTTGGCGATGGCTGCATGCCTCTTGGTACAGCGCGGAGCGGTAACAGGCGACGGCCTGCGCGTACGAATc      <  2:487585/150‑1 (MQ=255)
                                                                                                                                               gggTAAACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCCACAATGATCACCTCAGCCTTCGCCACCCCTTGGCGATGGCTGCATGCCTCTTGGTACAGCGCGGAGCGGTAACAGGCGACGGCCTGCGCGTACGAATc      >  1:1186092/1‑150 (MQ=255)
                                                                                                                                                 gTAAACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCCACAATGATCACCTCAGCCTTCGCCACCCCTTGGCGATGGCTGCATGCCTCTTGGTACAGCGCGGAGCGGTAACAGGCGACGGCCTGCGCGTACGAATCGa    <  1:4287501/150‑1 (MQ=255)
                                                                                                                                                   aaaCCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCCACAATGATCACCTCAGCCTTCGCCACCCCTTGGCGATGGCTGCATGCCTCTTGGTACAGCGCGGAGCGGTAACAGGCGACGGCCTGCGCGTACGAATCGAAc  >  2:5013476/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
GGCGGGCAACGAACTTAAAGTGCATCTGCAGTGAGGCGATGAGCAATAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTTCGCGGGTAAACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCCACAATGATCACCTCAGCCTTCGCCACCCCTTGGCGATGGCTGCATGCCTCTTGGTACAGCGCGGAGCGGTAACAGGCGACGGCCTGCGCGTACGAATCGAAC  >  NC_002947/4340313‑4340608

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