Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,198,437 G→A 5.2% R241H (CGC→CAC)  hemE → uroporphyrinogen decarboxylase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,198,4370GA 5.2% 244.3 / 2.2 76R241H (CGC→CAC) hemEuroporphyrinogen decarboxylase
Reads supporting (aligned to +/- strand):  ref base G (42/30);  new base A (2/2);  total (44/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.11e-01

TGAATGCGCAGATTAAAGCCGGTGCTCAGGCAGTGATGATTTTCGACACCTGGGGCGGCGTGCTTACCGGGCGCGATTATCAACAGTTCTCGCTCTATTACATGCATAAAATTGTTGATGGTTTACTGCGTGAAAACGACGGTCGCCGCGTACCGGTCACGCTGTTTACCAAAGGCGGCGGACAGTGGCTGGAAGCCATGGCAGAAACCGGTTGCGATGCGCTGGGCCTCGACTGGACAACGGACATCGCCGATGCGCGCCGCCGTGTGGGCAATAAAGTCGCGTTGCAGGG  >  NC_000913/4198293‑4198584
                                                                                                                                                |                                                                                                                                                   
tGAATGCGCAGATTAAAGCCGGTGCTCAGGCAGTGATGATTTTCGACACCTGGGGCGGCGTGCTTACCGGGCGCGATTATCAACAGTTCTCGCTCTATTACATGCATAAAATTGTTGATGGTTTACTGCGTGAAAACGACGGTcgccgc                                                                                                                                                 >  1:217719/1‑149 (MQ=255)
  aaTGCGCAGATTAAAGCCGGTGCTCAGGCAGTGATGATTTTCGACACCTGGGGCGGCGTGCTTACCGGGCGCGATTATCAACAGTTCTCGCTCTATTACATGCATAAAATTGTTGATGGTTTACTGCGTGAAAACGACGGTCGCCGCGt                                                                                                                                               >  2:1563847/1‑149 (MQ=255)
    tGCGCAGATTAAAGCCGGTGCTCAGGCAGTGATGATTTTCGACACCTGGGGCGGCGTGCTTACCGGGCGCGATTATCAACAGTTCTCGCTCTATTACATGCATAAAATTGTTGATGGTTTACTGCGTGAAAACGACGGTCGCCGCGTAc                                                                                                                                             <  1:1443688/149‑1 (MQ=255)
       gcAGATTAAAGCCGGTGCTCAGGCAGTGATGATTTTCGACACCTGGGGCGGCGTGCTTACCGGGCGCGATTATCAACAGTTCTCGCTCTATTACATGCATAAAATTGTTGATGGTTTACTGCGTGAAAACGACGGTCGCCGCGTACCgg                                                                                                                                          >  1:511678/1‑149 (MQ=255)
         aGATTAAAGCCGGTGCTCAGGCAGTGATGATTTTCGACACCTGGGGCGGCGTGCTTACCGGGCGCGATTATCAACAGTTCTCGCTCTATTACATGCATAAAATTGTTGATGGTTTACTGCGTGAAAACGACGGTCGCCGCGTACCGGTc                                                                                                                                        >  1:693929/1‑149 (MQ=255)
               aaGCCGGTGCTCAGGCAGTGATGATTTTCGACACCTGGGGCGGCGTGCTTACCGGGCGCGATTATCAACAGTTCTCGCTCTATTACATGCATAAAATTGTTGATGGTTTACTGCGTGAAAACGACGGTCGCCGCGTACCGGTCACGCTg                                                                                                                                  >  2:351204/1‑149 (MQ=255)
               aaGCCGGTGCTCAGGCAGTGATGATTTTCGACACCTGGGGCGGCGTGCTTACCGGGCGCGATTATCAACAGTTCTCGCTCTATTACATGCATAAAATTGTTGATGGTTTACTGCGTGAAAACGACGGTCGCCGCGTACCGGTCACGCTg                                                                                                                                  >  2:1378093/1‑149 (MQ=255)
                aGCCGGTGCTCAGGCAGTGATGATTTTCGACACCTGGGGCGGCGTGCTTACCGGGCGCGATTATCAACAGTTCTCGCTCTATTACATGCATAAAATTGTTGATGGTTTACTGCGTGAAAACGACGGTCGCCGCGTACCGGTCACGCTGt                                                                                                                                 <  2:524082/149‑1 (MQ=255)
                aGCCGGTGCTCAGGCAGTGATGATTTTCGACACCTGGGGCGGCGTGCTTACCGGGCGCGATTATCAACAGTTCTCGCTCTATTACATGCATAAAATTGTTGATGGTTTACTGCGTGAAAACGACGGTCGCCGCGTACCGGTCACGCTGt                                                                                                                                 <  2:946035/149‑1 (MQ=255)
                aGCCGGTGCTCAGGCAGTGATGATTTTCGACACCTGGGGCGGCGTGCTTACCGGGCGCGATTATCAACAGTTCTCGCTCTATTACATGCATAAAATTGTTGATGGTTTACTGCGTGAAAACGACGGTCGCCGCGTACCGGTCACGCTGt                                                                                                                                 <  2:1090183/149‑1 (MQ=255)
                 gCCGGTGCTCAGGCAGTGATGATTTTCGACACCTGGGGCGGCGTGCTTACCGGGCGCGATTATCAACAGTTCTCGCTCTATTACATGCATAAAATTGTTGATGGTTTACTGCGTGAAAACGACGGTCGCCGCGTACCGGTCACGCTGtt                                                                                                                                >  1:1186440/1‑149 (MQ=255)
                   cGGTGCTCAGGCAGTGATGATTTTCGACACCTGGGGCGGCGTGCTTACCGGGCGCGATTATCAACAGTTCTCGCTCTATTACATGCATAAAATTGTTGATGGTTTACTGCGTGAAAACGACGGTCGCCGCGTACCGGTCACGCTGTTTa                                                                                                                              <  2:1220009/149‑1 (MQ=255)
                        cTCAGGCAGTGATGATTTTCGACACCTGGGGCGGCGTGCTTACCGGGCGCGATTATCAACAGTTCTCGCTCTATTACATGCATAAACTTGTTGATGGTTTACTGCGTGAAAACGACGGTCGCCGCGTACCGGTCACGCTGTTTACCaaa                                                                                                                         <  2:249887/149‑1 (MQ=255)
                         tCAGGCAGTGATGATTTTCGACACCTGGGGCGGCGTGCTTACCGGGCGCGATTATCAACAGTTCTCGCTCTATTACATGCATAAAATTGTTGATGGTTTACTGCGTGAAAACGACGGTCGCCGCGTACCGGTCACGCTGTTTACCAAAg                                                                                                                        >  2:884486/1‑149 (MQ=255)
                          cAGGCAGTGATGATTTTCGACACCTGGGGCGGCGTGCTTACCGGGCGCGATTATCAACAGTTCTCGCTCTATTACATGCATAAAATTGTTGAAGGTTTACTGCGTGAAAACGACGGTCGCCGCGTACCGGTCACGCTGTTTACCAAAgg                                                                                                                       <  2:135774/149‑1 (MQ=255)
                                  gatgatTTTCGACACCTGGGGCGGCGTGCTTACCGGGCGCGATTATCAACAGTTCTCGCTCTATTACATGCATAAAATTGTTGATGGTTTACTGCGTGAAAACGACGGTCGCCGCGTACCGGTCACGCTGTTTACCAAAGGCGGCGGAc                                                                                                               >  2:1506386/1‑149 (MQ=255)
                                   atgatTTTCGACACCTGGGGCGGCGTGCTTACCGGGCGCGATTATCAACAGTTCTCGCTCTATTACATGCATAAAATTGTTGATGGTTTACTGCGTGAAAACGACGGTCGCCGCGTACCGGTCACGCTGTTTACCAAAGGCGGCGGACa                                                                                                              >  2:545806/1‑149 (MQ=255)
                                      atTTTCGACACCTGGGGCGGCGTGCTTACCGGGCGCGATTATCAACAGTTCTCGCTCTATTACATGCATAAAATTGTTGATGGTTTACTGCGTGAAAACGACGGTCGCCGCGTACCGGTCACGCTGTTTACCAAAGGCGGCGGACAGTg                                                                                                           >  1:1273121/1‑149 (MQ=255)
                                      atTTTCGACACCTGGGGCGGCGTGCTTACCGGGCGCGATTATCAACAGTTCTCGCTCTATTACATGCATAAAATTGTTGATGGTTTACTGCGTGAAAACGACGGTCGCCGCGTACCGGTCACGCTGTTTACCAAAGGCGGCGGACAGTg                                                                                                           >  2:191957/1‑149 (MQ=255)
                                       ttttCGACACCTGGGGCGGCGTGCTTACCGGGCGCGATTATCAACAGTTCTCGCTCTATTACATGCATAAAATTGTTGATGGTTTACTGCGTGAAAACGACGGTCGCCGCGTACCGGTCACGCTGTTTACCAAAGGCGGCGGACAGTgg                                                                                                          <  2:1186440/149‑1 (MQ=255)
                                          tCGACACCTGGGGCGGCGTGCTTACCGGGCGCGATTATCAACAGTTCTCGCTCTATTACATGCATAAAATTGTTGATGGTTTACTGCGTGAAAACGACGGTCGCCGCGTACCGGTCACGCTGTTTACCAAAGGCGGCGGACAGTGGCTg                                                                                                       <  1:1470457/149‑1 (MQ=255)
                                             acacCTGGGGCTGCGTGCTTAACGTGCTCGATTATCAACAGTTCTCGCTCTATTACATGCCTAAAATTGTTGATGGTTTACTGCGTGAAAACGACGGTCGCCGCGTACCGGTCACGCTGTTTACCAAAGGCGGCTGACAGTGGCTGGaa                                                                                                    <  2:411795/149‑1 (MQ=255)
                                              cacCTGGGGCGGCGTGCTTACCGGGCGCGATTATCAACAGTTCTCGCTCTATTACATGCATAAAATTGTTGATGGTTTACTGCGTGAAAACGACGGTCGCCGCGTACCGGTCACGCTGTTTACCAAAGGCGGCGGACAGTGGCTGGAAg                                                                                                   <  2:181915/149‑1 (MQ=255)
                                                 cTGGGGCGGCGTGCTTACCGGGCGCGATTATCAACAGTTCTCGCTCTATTACATGCATAAAATTGTTGATGGTTTACTGCGTGAAAACGACGGTCGCCGCGTACCGGTCACGCTGTTTACCAAAGGCGGCGGACAGTGGCTGGAAGCCa                                                                                                >  1:49928/1‑149 (MQ=255)
                                                   ggggCGGCGTGCTTACCGGGCGCGATTATCAACAGTTCTCGCTCTATTACATGCATAAAATTGTTGATGGTTTACTGCGTGAAAACGACGGTCGCCGCGTACCGGTCACGCTGTTTACCAAAGGCGGCGGACAGTGGCTGGAAGCCATg                                                                                              >  1:486478/1‑149 (MQ=255)
                                                   ggggCGGCGTGCTTACCGGGCGCGATTATCAACAGTTCTCGCTCTATTACATGCATAAAATTGTTGATGGTTTACTGCGTGAAAACGACGGTCGCCGCGTACCGGTCACGCTGTTTACCAAAGGCGGCGGACAGTGGCTGGAAGCCATg                                                                                              >  1:133628/1‑149 (MQ=255)
                                                    gggCGGCGTGCTTACCGGGCGCGATTATCAACAGTTCTCGCTCTATTACATGCATAAAATTGTTGATGGTTTACTGCGTGAAAACGACGGTCGCCGCGTACCGGTCACGCTGTTTACCAAAGGCGGCGGACAGTGGCTGGAAGCCATgg                                                                                             >  1:64767/1‑149 (MQ=255)
                                                    gggCGGCGTGCTTACCGGGCGCGATTATCAACAGTTCTCGCTCTATTACATGCATAAAATTGTTGATGGTTTACTGCGTGAAAACGACGGTCGCCGCGTACCGGTCACGCTGTTTACCAAAGGCGGCGGACAGTGGCTGGAAGCCATgg                                                                                             >  1:747151/1‑149 (MQ=255)
                                                     ggcggcGTGCTTACCGGGCGCGATTATCAACAGTTCTCGCTCTATTACATGCATAAAATTGTTGATGGTTTACTGCGTGAAAACGACGGTCGCCGCGTACCGGTCACGCTGTTTACCAAAGGCGGCGGACAGTGGCTGGAAGCCATGGc                                                                                            <  2:964153/149‑1 (MQ=255)
                                                     ggcggcGTGCTTACCGGGCGCGATTATCAACAGTTCTCGCTCTATTACATGCATAAAATTGTTGATGGTTTACTGCGTGAAAACGACGGTCGCCGCGTACCGGTCACGCTGTTTACCAAAGGCGGCGGACAGTGGCTGGAAGCCATGGc                                                                                            <  2:1453296/149‑1 (MQ=255)
                                                           gTGCTTACCGGGCGCGATTATCAACAGTTCTCGCTCTATTACATGCATAAAATTGTTGATGGTTTACTGCGTGAAAACGACGGTCGCCGCGTACCGGTCACGCTGTTTACCAAAGGCGGCGGACAGTGGCTGGAAGCCATGGCAGAAAc                                                                                      >  2:629570/1‑149 (MQ=255)
                                                             gCTTACCGGGCGCGATTATCAACAGTTCTCGCTCTATTACATGCATAAAATTGTTGATGGTTTACTGCGTGAAAACGACGGTCGCCGCGTACCGGTCACGCTGTTTACCAAAGGCGGCGGACAGTGGCTGGAAGCCATGGCAGAAACCg                                                                                    >  2:1370115/1‑149 (MQ=255)
                                                               ttACCGGGCGCGATTATCAACAGTTCTCGCTCTATTACATGCATAAAATTGTTGATGGTTTACTGCGTGAAAACGACGGTCGCCGCGTACCGGTCACGCTGTTTACCAAAGGCGGCGGACAGTGGCTGTAAGCAATGGCAGAACCCGGt                                                                                  >  1:27292/1‑149 (MQ=255)
                                                               ttACCGGGCGCGATTATCAACAGTTCTCGCTCTATTACATGCATAAAATTGTTGATGGTTTACTGCGTGAAAACGACGGTCGCCGCGTACCGGTCACGCTGTTTACCAAAGGCGGCGGACAGTGGCTGGAAGCCATGGCAGAAACCGGt                                                                                  >  2:667083/1‑149 (MQ=255)
                                                                      gcgcGATTATCAACAGATCTCGCTCTATTACATGCATAACATTGTTGATGGTTTACTGCGTGAAACCGACGGTCGCCGCGTACCGGTCACGCTGTTTACCAAAGGCGGCGGACAGTGGCTGGAAGCCATGGCAGAAACCGGTTGCGATg                                                                           <  1:1388904/149‑1 (MQ=255)
                                                                         cgATTATCAACAGTTCTCGCTCTATTACATGCATAAAATTGTTGATGGTTTACTGCGTGAAAACGACGGTCGCCGCGTACCGGTCACGCTGTTTACCAAAGGCGGCGGACAGTGGCTGGAAGCCATGGCAGAAACCGGTTGCGATgcgc                                                                        <  1:1563847/149‑1 (MQ=255)
                                                                          gATTATCAACAGTTCTCGCTCTATTACATGCATAAAATTGTTGATGGTTTACTGCGTGAAAACGACGGTCGCCGCGTACCGGTCACGCTGTTTACCAAAGGCGGCGGACAGTgg                                                                                                          >  1:1103822/1‑114 (MQ=255)
                                                                          gATTATCAACAGTTCTCGCTCTATTACATGCATAAAATTGTTGATGGTTTACTGCGTGAAAACGACGGTCGCCGCGTACCGGTCACGCTGTTTACCAAAGGCGGCGGACAGTgg                                                                                                          <  2:1103822/114‑1 (MQ=255)
                                                                          gATTATCAACAGTTCTCGATCTATTACATGCATAAAATTGTTGATGGTTTACTGCGTGAAAACGACGGTCGCCGAGTACCGGTCACGCTGTTTACCAAAGGCGGCGGACAGTGGATGGAAGCCATGGCAGAAACCTGTTGAGATGCGCt                                                                       >  1:667331/1‑149 (MQ=255)
                                                                           aTTATCAACAGTTCTCGCTCTATTACATGCATAAAATTGTTGATGGTTTACTGCGTGAAAACGACGGTCGCCGCGTACCGGTCACGCTGTTTACCAAAGGCGGCGGACAGTGGCTGGAAGCCATGGCAGAAACCGGTTGCGATGCGCTg                                                                      <  2:1586147/149‑1 (MQ=255)
                                                                            ttATCAACAGTTCTCGCTCTATTACATGCATAAAATTGTTGATGGTTTACTGCGTGAAAACGACGGTCGCCGCGTACCGGTCACGCTGTTTACCAAAGGCGGCGGACAGTGGCTGGAAGCCATGGCAGAAACCGGTTGCGATGCGCTgg                                                                     <  1:763626/149‑1 (MQ=255)
                                                                            ttATCAACAGTGCTCGCTCTATTACATGCATAAAATTGTTGATGGTTTACTGTGTGAAAACGACGGGCGCCGCGTACCGGTCACGCTGTTTACCAAAGGCGGCGGACAGTGGCTGGAAGCCATGGCAGAAACCGGTTGCGATGCGCTgg                                                                     <  2:580043/149‑1 (MQ=255)
                                                                                  aCAGTTCTCGCTCTATTACATGCATAAAATTGTTGATGGTTTACTGCGTGAAAACGACGGTCGCCGCGTACCGGTCACGCTGTTTACCAAAGGCGGCGGACAGTGGCTGGAAGCCATGGCAGAAACCGGTTGCGATGCGCTGGGCCTCg                                                               >  1:1257702/1‑149 (MQ=255)
                                                                                   cAGTTCGCGCTCTATTACATGCATAAAATTGTTGATGGTTTACTGCGTGAAAACGACGGTCGCCGCGTACAGGTCACGCTGTTTACCAAAGGCGGCGGACAGTGGCTGGAAGCCATGGCAGAAACCGGTTGCGATGCGCTGGGCCTCGa                                                              <  2:175347/149‑1 (MQ=255)
                                                                                    aGTTCTCGCTCTATTACATGCATAAAATTGTTGATGGTTTACTGCGTGAAAACGACGGTCACCGCGTACCGGTCACGCTGTTTACCAAAGGCGGCGGACAGTGGCTGGAAGCCATGGCAGAAAc                                                                                      >  1:1173376/1‑124 (MQ=255)
                                                                                    aGTTCTCGCTCTATTACATGCATAAAATTGTTGATGGTTTACTGCGTGAAAACGACGGTCACCGCGTACCGGTCACGCTGTTTACCAAAGGCGGCGGACAGTGGCTGGAAGCCATGGCAGAAAc                                                                                      >  1:1167193/1‑124 (MQ=255)
                                                                                    aGTTCTCGCTCTATTACATGCATAAAATTGTTGATGGTTTACTGCGTGAAAACGACGGTCACCGCGTACCGGTCACGCTGTTTACCAAAGGCGGCGGACAGTGGCTGGAAGCCATGGCAGAAAc                                                                                      <  2:1167193/124‑1 (MQ=255)
                                                                                    aGTTCTCGCTCTATTACATGCATAAAATTGTTGATGGTTTACTGCGTGAAAACGACGGTCACCGCGTACCGGTCACGCTGTTTACCAAAGGCGGCGGACAGTGGCTGGAAGCCATGGCAGAAAc                                                                                      <  2:1173376/124‑1 (MQ=255)
                                                                                      ttCTCGCTCTATTACATGCATAAAATTGTTGATGGTTTACTGCGTGAAAACGACGGTCGCCGCGTACCGGTCACGCTGTTTACCAAAGGCGGCGGACAGTGGCTGGAAGCCATGGCAGAAACCGGTTGCGATGCGCTGGGCCTCGACTg                                                           >  2:758640/1‑149 (MQ=255)
                                                                                      ttCTCGCTCTATTACATGCATAAAATTGTTGATGGTTTACTGCGTGAAAACGACGGGCGCCGCGTACCGGTCACGCTGTTTACCAAAGGCGGCGGACAGTGGCTGGAAGCCATGGCAGAAACCGGTTGCGATGCGCTGGGCCTCGACTg                                                           <  2:1071828/149‑1 (MQ=255)
                                                                                        ctcGCTCTATTACATGCATAAACTTGTAGATGGTTTACTGCGTGAAAACGACGGTCGCCGCGTACCGGTCACGCGGTTTACCAAAGGCGGCGGACAGTGGCTGGAAGCCATGGCAGAAACCGGTTGCGATGCGCTGGGCCTCGACTGGa                                                         <  1:184933/149‑1 (MQ=255)
                                                                                        ctcGCTCTATTACATGCATAAAATTGTTGATGGTTTACTGCGTGAAAACGACGGTCGCCGCGTACCGGTCACGCTGTTTACCAAAGGCGGCGGACAGTGGCTGGAAGCCATGGCAGAACCCGGTTGCGATGCGCAGGGCCTCGACTGGa                                                         >  2:263511/1‑149 (MQ=255)
                                                                                             tctATTACATGCATAAAATTGTTGATGGTTTACTGCGTGAAAACGACGGTCGCCGCGTACCGGTCACGCTGTTTACCAAAGGCGGCGGACAGTGGCTGGAAGCCATGGCAGAAACCGGTTGCGATGCGCTGGGCCTCGACTGGACAACg                                                    >  1:1476074/1‑149 (MQ=255)
                                                                                                 ttACATGCATAAAATTGTTGATGGTTTACTGCGTGAAAACGACGGTCGCCGCGTACCGGTCACGCTGTTTACCAAAGGCGGCGGACAGTGGCTGGAAGCCATGGCAGAAACCGGTTGCGATGCGCTGGGCCTCGACTGGACAACGGACa                                                <  2:511678/149‑1 (MQ=255)
                                                                                                 ttACATGCATAAAATTGTTGATGGTTTACTGCGTGAAAACGACGGTCGCCGCGTACCGGTCACGCTGTTTACCAAAGGCGGCGGACAGTGGCTGGAAGCCATGGCAGAAACCGGTTGCGATGCGCTGGGCCTCGACTGGACAACGGACa                                                >  2:1256725/1‑149 (MQ=255)
                                                                                                     aTGCATAAAATTGTTGATGGTTTACTGCGTGAAAACGACGGTCGCCGCGTACCGGTCACGCTGTTTACCAAAGGCGGCGGACAGTGGCTGGAAGCCATGGCAGAAACCGGTTGCGATGCGCTGGGCCTCGACTGGACAACGGACATCGc                                            <  2:217719/149‑1 (MQ=255)
                                                                                                          tAAAATTGTTGATGGTTTACTGCGTGAAAACGACGGTCGCCGCGTACCGGTCACGCTGTTTACCAAAGGCGGCGGACAGTGGCTGGAAGCCATGGCAGAACCCGGTTGCGATGCGCTGGGCCTCGACTGGACAACGGACATCGCCGATg                                       >  2:1378956/1‑149 (MQ=255)
                                                                                                          tAAAATTGTTGATGGTTTACTGCGTGAAAACGACGGTCGCCGCGTACCGGTCACGCTGTTTACCAAAGGCGGCGGACAGTGGCTGGAAGCCATGGCAGAAACCGGTTGCGATGCGCTGGGCCTCGACTGGACAACGGACATCGCCGATg                                       >  1:1377884/1‑149 (MQ=255)
                                                                                                                         tttACTGCGTGAAAACGACGGTCGCCGCGTACCGGTCACGCTGTTTACCAAAGGCGGCGGACAGTGGCTGGAAGCCATGGCAGAAACCGGTTGCGATGCGCTGGGCCTCGACTGGACAACGGACATCGCCGATGCGCGCCGCCGTGTgg                        >  2:1184353/1‑149 (MQ=255)
                                                                                                                               gCGTGAAAACGACGGTCGCCGCGTACCGGTCACGCTGTTTACCAAAGGCGGCGGACAGTGGCTGGAAGCCATGGCAGAAACCGGTTGCGATGCGCTGGGCCTCGACTGGACAACGGACATCGCCGATGCGCGCCGCCGTGTGGGCAATa                  <  1:351204/149‑1 (MQ=255)
                                                                                                                                 gTGAAAACGACGGTCGCCGCGTACCGGTCACGCTGTTTACCAAAGGCGGCGGACAGTGGCTGGAAGCCATGGCAGAAACCGGTTGCGATGCGCTGGGCCTCGACTGGACAACGGACATCGCCGATGCGCGCCGCCGTGTGGGCAATaaa                <  1:1152143/149‑1 (MQ=255)
                                                                                                                                 gTGAAAACGACGGTCGCCGCGTACCGGTCACGCTGTTTACCAAAGGCGGCGGACAGTGGCTGGAAGCCATGGCAGAAACCGGTTGCGATGCGCTGGGCCTCGACTGGACAACGGACATCGCCGATGCGCGCCGCCGTGTGGGCAATaaa                >  2:1224260/1‑149 (MQ=255)
                                                                                                                                  tGAAAACGACGGTCGCCGCGTACCGGTCACGCTGTTTACCAAAGGCGGCGGACAGTGGCTGGAAGCCATGGCAGAAACCGGTTGCGATGCGCTGGGCCTCGACTGGACAACGGACATCGCCGATGCGCGCCGCCGTGTGGGCAATAAAg               >  1:993674/1‑149 (MQ=255)
                                                                                                                                  tGAAAACGACGGTCGCCGCGTACCGGTCACGCTGTTTACCAAAGGCGGCGGACAGTGGCTGGAAGCCATGGCAGAAACCGGTTGCGATGCGCTGGGCCTCGACTGGACAACGGACATCGCCGATGCGCGCCGCCGTGTGGGCAATAAAg               >  1:1102978/1‑149 (MQ=255)
                                                                                                                                  tGAAAACGACGGTCGCCGCGTACCGGTCACGCTGTTTACCAAAGGCGGCGGACAGTGGCTGGAAGCCATGGCAGAAACCGGTTGCGATGCGCTGGGCCTCGACTGGACAACGGACATCGCCGATGCGCGCCGCCGTGTGGGCAATAAAg               >  1:1183010/1‑149 (MQ=255)
                                                                                                                                    aaaaCGACGGTCGCCGCGTACCGGTCACGCTGTTTACCAAAGGCGGCGGACAGTGGCTGGAAGCCATGGCAGAAACCGGTTGCGATGCGCTGGGCCTCGCCTGGACAACGGACATCGCCGATGCGCGCCGCCGTGTGGGCAATAAAGTc             >  1:1178247/1‑149 (MQ=255)
                                                                                                                                    aaaaCGACGGTCGCCGCGTACCGGTCACGCTGTTTACCAAAGGCGGCGGACAGTGGCTGGAAGCCATGGCAGAAACCGGTTGCGATGCGCTGGGCCTCGACTGGACAACGGACATCGCCGATGCGCGCCGCCGTGTGGGCAATAAAGTc             <  1:629570/149‑1 (MQ=255)
                                                                                                                                    aaaaCGACGGTCGCCGCGTACCGGTCACGCTGTTTACCAAAGGCGGCGGACAGTGGCTGGAAGCCATGGCAGAAACCGGTTGCGATGCGCTGGGCCTCGACTGGACAACGGACATCGCCGATGCGCGCCGCCGTGTGGGCAATAAAGTc             >  2:1571139/1‑149 (MQ=255)
                                                                                                                                     aaaCGACGGTCGCCGCGTACTGGTCACGCTGTTTACCAAAGGCGGCGGACAGTGGCTGGAAGCCATGGCAGAAACCGGTTGCGATGCGCTGGGCCTCGACTGGACAACGGACATCGCCGATGCGCGCCGCCGTGTGGGCAATAAAGTcg            >  1:1255005/1‑149 (MQ=255)
                                                                                                                                     aaaCGACGGTCGCCGCGTACCGGTCACGCTGTTTACCAAAGGCGGCGGACAGTGGCTGGAAGCCATGGCAGAAACCGGTTGCGATGCGCTGGGCCTCGACTGGACAACGGACATCGCCGATGCGCGCCGCCGTGTGGGCAATAAAGTcg            <  2:64767/149‑1 (MQ=255)
                                                                                                                                       acgacgGTCGCCGCGTACCTGTCACGCTGTTTACCAAAGGCGGCGGACAGTGGCTGGAAGCCATGGCAGAAACCGGTTGCGATGCGCTTGGCCTCGACTGTACAACGGACATCGCAGATGCGCGCCGCCGTGTCGGCAATAAAGTcgag          >  2:239556/1‑147 (MQ=255)
                                                                                                                                       acgacgGGCGCCGCGTACCGGTCACGCTGTTTACCAAAGGCGGCGGACAGTGGCTGGAAGCCATGGCAGAAACCGGTTGCGATGCGCTGGGCCTCGACTGGACAACGGACATCGCCGATGCGCGCCGCCGTGTGGGCAATAAAGTcgcg          <  1:1378093/149‑1 (MQ=255)
                                                                                                                                         gacgGTCGCCGCGTCCCGGTCACTCTGTTTACCAAAGGCGGCGGACAGTGGCTGGAAGCCATGGCAGACACCGGTTGCGATGCGCTGGGCCTCGACTGGACAACGGACATCGCCGATTCGCGCCGCCGGGTGGGCAATAAAGTCGCGtt        >  2:466030/1‑149 (MQ=255)
                                                                                                                                         gacgGTCGCCGCGTACCGGTCACGCTGTTTACCAAAGGCGGCGGACAGTGGCTGGAAGCCATGGCAGAAACCGGTTGCGATGCGCTGGGCCTCGACTGGACAACGGCCATCGCCGATGCGCGCCGCCGTGTGGGCAATAAAGTCGCGtt        >  2:12781/1‑149 (MQ=255)
                                                                                                                                         gacgGTCGCCGCGTACCGGTCACGCTGTTTACCAAAGGCGGCGGACAGTGGCTGGAAGCCATGGCAGAAACCGGTTGCGAAGCGCTGGGCCTCGACTGGACAACGGACATCGCCGATGCGCGCCGCCGTGTGGGCAATAAAGTCGCGtt        >  2:460084/1‑149 (MQ=255)
                                                                                                                                          acgGTCGCCGCGTACCGGTCACGCTGTTTACCAAAGGCGGCGGACAGTGGCTGGAAGCCATGGCAGAAACCGGTTGCGATGCGCTGGGCCTCGACTGGACAACGGACATCGCCGATGCGCGCCGCCGTGTGGGCAATAAAGTCGCGTTg       >  2:767075/1‑149 (MQ=255)
                                                                                                                                              tCGCCGCGTACCGGTCACGCTGTTTACCAAAGGCGGCGGACAGTGGCTGGAAGCCATGGCAGAAACCGGTTGCGATGCGCTGGGCCTCGACTGGACAACGGACATCGCCGATGCGCGCCGCCGTGTGGGCAATAAAGTCGCGTTGCAgg   <  2:486478/149‑1 (MQ=255)
                                                                                                                                               cgccgcGTACCGGTCACGCTGTTTACCAAAGGCGGCGGACAGTGGCTGGAAGCCATGGCAGAAACCGGTTGCGATGCGCTGGGCCTCGACTGGACAACGGACATCGCCGATGCGCGCCGCCGTGTGGGCAATAAAGTCGCGTTGCAggg  <  2:1476074/149‑1 (MQ=255)
                                                                                                                                               cgccgcGTACCGGTCACGCTGTTTACCAAAGGCGGCGGACAGTGGCTGGAAGCCATGGCAGAAACCGGTTGCGATGCGCTGGGCCTCGACTGGACAACGGACATCGCCGATGCGCGCCGCCGTGTGGGCAATAAAGTCGCGTTGCAggg  <  2:1273121/149‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                                   
TGAATGCGCAGATTAAAGCCGGTGCTCAGGCAGTGATGATTTTCGACACCTGGGGCGGCGTGCTTACCGGGCGCGATTATCAACAGTTCTCGCTCTATTACATGCATAAAATTGTTGATGGTTTACTGCGTGAAAACGACGGTCGCCGCGTACCGGTCACGCTGTTTACCAAAGGCGGCGGACAGTGGCTGGAAGCCATGGCAGAAACCGGTTGCGATGCGCTGGGCCTCGACTGGACAACGGACATCGCCGATGCGCGCCGCCGTGTGGGCAATAAAGTCGCGTTGCAGGG  >  NC_000913/4198293‑4198584

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