Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,227,550 G→A 7.0% intergenic (+39/‑181) metH → / → yjbB homocysteine‑N5‑methyltetrahydrofolate transmethylase, B12‑dependent/putative Na+/Pi‑cotransporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,227,5500GA 7.0% 310.0 / 9.2 100intergenic (+39/‑181)metH/yjbBhomocysteine‑N5‑methyltetrahydrofolate transmethylase, B12‑dependent/putative Na+/Pi‑cotransporter
Reads supporting (aligned to +/- strand):  ref base G (48/45);  new base A (4/3);  total (52/48)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.99e-01

Only 100 of 103 total aligned reads displayed.

CACAAATTCAGCGCGATCAGGTTGAAGATTATGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGAGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGG‑AATAATCTTAGCCGGAAAAACGAATGTCAAAAGGCACAGTAAAGATGCCAGTCAGCGTGTTTCCGGTCCTTTGGTCATTCAATTGCTT  >  NC_000913/4227405‑4227696
                                                                                                                                                 |                                                                                                                                                   
cacaAATTCAGCGCGATCAGGTTGAAGATTATGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGAGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCt                                                                                                                                                   >  1:203108/1‑148 (MQ=255)
 acaAATTCAGCGCGATCAGGTTGAAGATTATGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGAGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTc                                                                                                                                                  >  8:395349/1‑148 (MQ=255)
  caAATTCAGCGCGATCAGGTTGAAGATTATGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGAGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCaa                                                                                                                                                >  2:271561/1‑149 (MQ=255)
   aaaTTCAGCGCGATCAGGTTGAAGATTATGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGAGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCACTCatt                                                                                                                                               >  2:223724/1‑147 (MQ=255)
   aaaTTCAGCGCGATCAGGTTGAAGATAATGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGAGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAAt                                                                                                                                               <  2:194264/149‑1 (MQ=255)
   aaaTTCAGCGCGATCAGGTTGAACATTATGCCCGCCGTAAAGGTATGAGCGTTACCGAACTTGAGCGCCGGCTGCCACCGAATCTGGGGTATGACGCTGACTGATTCACCAATCTGCCACTCTTCCTTCCAACAAACAGGgcgcgccat                                                                                                                                               <  4:220109/149‑6 (MQ=255)
     aTTCAGCGCGATCAGGTTGAAGATTATGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGAGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGa                                                                                                                                             <  4:300591/149‑1 (MQ=255)
     aTTCAGCGCGATCAGGTTGAAGATTATGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGAGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTGCCTTACAACAAACAGGGCGCTCAATGa                                                                                                                                             <  5:187121/149‑1 (MQ=255)
         aGCGCGATCAGGTTGAAGATTATGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGAGCGCTGGCTGGCACCGCATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGCGTGc                                                                                                                                         >  2:370814/1‑149 (MQ=255)
           cgcgATCAGGTTGAAGATTATGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGAGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGccc                                                                                                                                       <  1:319308/149‑1 (MQ=255)
            gcgATCAGGTTGAAGATTATGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGAGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCt                                                                                                                                      >  8:48185/1‑149 (MQ=255)
              gATCAGGTTGAAGATTATGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGAGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGTGCGCTCAATGAGTGCCCTGt                                                                                                                                    <  5:23239/149‑1 (MQ=255)
                  aGGTTGAAGATTATGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGAGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCACTCATTGAGCGCCCTGtttg                                                                                                                                 <  1:223724/148‑4 (MQ=255)
                   ggTTGAAGATTATGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGAGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCtt                                                                                                                               <  2:75649/149‑1 (MQ=255)
                    gTTGAAGATTATGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGAGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCtt                                                                                                                               >  4:63540/1‑148 (MQ=255)
                     ttGTAGATTATGCGCGCAGTTAAGGTATGAGCGTTACCGAAGTTGAGCGGTGGCTGGCAGCGTATCTGGGGTTTGACGCCGACTGATTCAATAAGCTGTAACTTTTCCTTTCATCTAACAGGTCGCTCAATGAGTGCCCTGTCTCTtta                                                                                                                             <  1:137081/149‑1 (MQ=255)
                      tGAAGATTATGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGAGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTGTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTttat                                                                                                                            >  6:42771/1‑149 (MQ=255)
                         aGATTATGCCCGCCGTAATGGTATGAGCGTTACCGAAGTTGAGCGCTGGCTGGAAACGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTaa                                                                                                                         <  1:118954/149‑1 (MQ=255)
                         aGATTATGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGAGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTaa                                                                                                                         <  7:196021/149‑1 (MQ=255)
                          gATTATGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGAGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAat                                                                                                                        <  2:203108/149‑1 (MQ=255)
                            ttATGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGAGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAatat                                                                                                                      <  4:72586/149‑1 (MQ=255)
                               tGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGAGCGCTGGCTGGCACCGAATCTGGGGTATGAAGAGGACTGATTCACAAATCTCTCACTTTTCCTTACACCAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGa                                                                                                                    <  1:265502/148‑1 (MQ=255)
                                gCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGAGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGAATGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGaaa                                                                                                                  >  6:120214/1‑149 (MQ=255)
                                 cccGCCGTAAAGGTATGAGCGTTACCGAAGTTGAGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGaaa                                                                                                                  <  8:66223/148‑1 (MQ=255)
                                  ccgccgTAAAGGTATGAGCGTTACCGAAGTTGAGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAaca                                                                                                                >  5:306658/1‑149 (MQ=255)
                                  ccgccgTAAAGGTATGAGCGTTACCGAAGTTGAGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAaca                                                                                                                >  4:165748/1‑149 (MQ=255)
                                   cgccgTAAAGGTATGAGCGTTACCGAAGTTGAGCGCTGGCTGGCTCCGACTCTGGGGTATGACCCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATCTGAAacac                                                                                                               <  8:389173/149‑1 (MQ=255)
                                    gccgTAAAGGTATGAGCGTTACCGAAGTTGAGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACt                                                                                                              <  2:212410/149‑1 (MQ=255)
                                      cgTAAAGGTATGAGCGTTACCGAAGTTGCGCGCTGGCTGCCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACCAACAGGGCGCTCAATGAGTGACCTGTCTCTTTATTAATATGAAACACTta                                                                                                            <  2:118995/149‑1 (MQ=255)
                                      cgTAAAGGTATGAGCGTTACCGAAGTTGAGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTta                                                                                                            >  3:245197/1‑149 (MQ=255)
                                        tAAAGGTATGAGCGTTACCGAAGTTGAGCGCTGGCGGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTTTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTtata                                                                                                          <  5:209623/149‑1 (MQ=255)
                                          aaGGTATGAGCGTTACCGAAGTTGAGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATAc                                                                                                         <  4:245197/148‑1 (MQ=255)
                                            ggTATGAGCGTTACCGAAGTTGAGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTgg                                                                                                      <  8:73832/149‑1 (MQ=255)
                                              taggAGCGTTACCGAAGGTGAGCGCTGGCTGGCACCGAAGCTGGGGTATGACGCGGACTGATTCACAAATCTGTAACTTTTCCTTACAACAAACAGCGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGaa                                                                                                    <  8:366564/146‑1 (MQ=255)
                                                tGAGCGTTACCGAAGTTGAGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGaaa                                                                                                   >  8:44836/1‑148 (MQ=255)
                                                tGAGCGTTACCGAAGTTGAGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAAc                                                                                                  <  3:156785/149‑1 (MQ=255)
                                                     gTTACCGAAGTTGAGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCGTACAACAAACAGGGCGCTCAATGAGTGCCCGGTCTCTTTATTAATATGAAACACTGATACTGGAAACAGGCt                                                                                             >  6:219115/1‑149 (MQ=255)
                                                         ccGAAGTTGAGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGG‑aa                                                                                        <  8:47834/149‑1 (MQ=255)
                                                         ccGAAGTTGAGCGC‑GGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGG‑aa                                                                                        >  1:222830/1‑148 (MQ=255)
                                                          cGAAGTTGAGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGG‑aa                                                                                        >  7:288819/1‑148 (MQ=255)
                                                          cGAAGTTGAGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGG‑aa                                                                                        >  1:277283/1‑148 (MQ=255)
                                                          cGAAGTTGAGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGAGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGG‑aa                                                                                        >  2:96455/1‑148 (MQ=255)
                                                                       tggctggcACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGG‑AATAATCTTAGCCgg                                                                           >  7:238356/1‑148 (MQ=255)
                                                                         gctggcACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGG‑AATAATCTTAGCCGGaa                                                                         >  8:351606/1‑148 (MQ=255)
                                                                         gctggcACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGG‑AATAATCTTAGCCGGa                                                                          <  7:395349/147‑1 (MQ=255)
                                                                           tggcACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGACACAGGCTGG‑AATAATCTTAGCCGGaaaaa                                                                      >  1:274462/1‑149 (MQ=255)
                                                                              cACCGAATCTTGGGTATGACGCGGACTGATTCACAAATCTGGCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACATATACTGGAAACAGGCTGG‑AATAATCTTAGCCGGAAAAACGa                                                                   <  2:224537/149‑1 (MQ=255)
                                                                               aCCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCACTCATTGAGCGcc                                                                                                                                        >  8:129049/1‑80 (MQ=255)
                                                                               aCCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCACTCATTGAGCGcc                                                                                                                                        <  7:129049/80‑1 (MQ=255)
                                                                                ccGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGT‑AATAATCTTAGCCGGAAAAACGAAt                                                                 <  4:263885/149‑1 (MQ=255)
                                                                                 cGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGG‑AATAATCTTAGCCGGAAAAACGAATg                                                                <  7:48185/149‑1 (MQ=255)
                                                                                  gAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGG‑AATAATCTTAGCCGGAAAAACGAATGt                                                               >  5:328166/1‑149 (MQ=255)
                                                                                   aaTCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGG‑AATAATCTTAGCCGGAAAAACGAATGTc                                                              <  6:306658/149‑1 (MQ=255)
                                                                                       tGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCACTCATTGAGCGCCCt                                                                                                                                      >  1:229735/1‑74 (MQ=38)
                                                                                       tGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCACTCATTGAGCGCCCt                                                                                                                                      <  2:229735/74‑1 (MQ=38)
                                                                                           gTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGG‑AATAATCTTAGCCGGAAAAACGAATGTCAAAAGGca                                                      >  6:355637/1‑149 (MQ=255)
                                                                                              tGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGG‑AATAATCTTAGCCGGAAAAACGAATGTCAAAAGGcaca                                                    <  1:271561/148‑1 (MQ=255)
                                                                                              tGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGG‑AATAATCTTAGCCGGAAAAACGAATGTCAAAAGGcaca                                                    <  8:43028/148‑1 (MQ=255)
                                                                                              tGACGCGGAATGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGG‑AATAATCTTAGCCGGAAAAACGAATGTCAAAAGGcaca                                                    <  5:120214/148‑1 (MQ=255)
                                                                                                       cTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGG‑AATAATCTTAGCCGGAAAAACGAATGTCAAAAGGCACAGTAAAGATGc                                          <  8:238356/149‑1 (MQ=255)
                                                                                                          aTTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGG‑AATAATCTTAGCCGGAAAAACGAATGTCAAAAGGCACAGGAAAGATGCCa                                        <  7:54675/148‑1 (MQ=255)
                                                                                                           ttCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGG‑AATAATCTTAGCCGGAAAAACGAATGTCAAAAGGCACAGTAAAGATGCCAGt                                      <  8:276190/149‑1 (MQ=255)
                                                                                                             cacaAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGG‑AATAATCTTAGCCGGAAAAACGAATGTCAAAAGGCACAGTAAAGATGCCAGt                                      >  2:279790/1‑147 (MQ=255)
                                                                                                             cacaAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGG‑AATAATCTTAGCCGGAAAAACGAATGTCAAAAGGCACAGTAAAGATGCCAGt                                      >  3:318529/1‑147 (MQ=255)
                                                                                                              acaAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGG‑AATAATCTTAGCCGGAAAAACGAATGTCAAAAg                                                         >  5:48750/1‑127 (MQ=255)
                                                                                                              acaAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGG‑AATAATCTTAGCCGGAAAAACGAATGTCAAAAg                                                         <  6:48750/127‑1 (MQ=255)
                                                                                                              acaAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGG‑AATAATCTTAGCCGGAAAAACGAATGTCAAAAGGCACAGTAAAGATGCCAGTCAg                                   >  2:153649/1‑149 (MQ=255)
                                                                                                              acaAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTAG‑AATAATCTTAGCCGGAAAAACGAATGTCAAAAGGCACAGTAAAGATGCCAGTCAg                                   >  1:320883/1‑149 (MQ=255)
                                                                                                                     tGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGG‑ATTAATCTTAGCCGGAAAAACGAATGTCAAAAGGCACAGTAAAGATGCCAGTCAGCGTGttt                            >  1:130417/1‑149 (MQ=255)
                                                                                                                      gTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGG‑AATAATCTTAGCCGGAAAAACGAATGTCAAAAGGCACAGTAAAGATGCCAGTCAGCGTGTTTc                           >  7:73943/1‑149 (MQ=255)
                                                                                                                       tCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGG‑AATAATCTTAGCCGGAAAAACGAATGTCAAAAGGCACAGTAAAGATGCCAGTCAGCGTGTTTcc                          >  4:272029/1‑149 (MQ=255)
                                                                                                                       tCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGG‑AATAATCTTAGCCGGAAAAACGAATGTCAAAAGGCACAGTAAAGATGCCAGTCAGCGTGTTTcc                          >  7:11318/1‑149 (MQ=255)
                                                                                                                       tCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGG‑AATAATCTTAGCCGGAAAAACGAATGTCAAAAGGCACAGTAAAGATGCCAGTCAGCGTGTTTcc                          >  6:156567/1‑149 (MQ=255)
                                                                                                                           ttttCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGG‑AATAATCTTAGCCGGAAAAACGAATGTCAAAAGGCACAGTAAAGATGCCAGTCAGCGTGTTTCCGGTc                      >  6:103184/1‑149 (MQ=255)
                                                                                                                           ttttCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGG‑AATAATCTTAGCCGGAAAAACGAATGTCAAAAGGCACAGTAAAGATGCCAGTCAGCGTGTTTCCGGTc                      >  2:294672/1‑149 (MQ=255)
                                                                                                                               ccTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGG‑AATAATCTTAGCCGGAAAAACGAATGTCAAAAGGCACAGTAAAGATGCCAGTCa                                    >  1:344875/1‑131 (MQ=255)
                                                                                                                               ccTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGG‑AATAATCTTAGCCGGAAAAACGAATGTCAAAAGGCACAGTAAAGATGCCAGTCa                                    <  2:344875/131‑1 (MQ=255)
                                                                                                                                cTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGG‑AATAATCTTAGCCGGAAAAACGAATGTCAAAAGGCACAGTAAAGATGCCAGTCAGCGTGTTTCCGGTCCttt                  <  5:42771/148‑1 (MQ=255)
                                                                                                                                cTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGG‑AATAATCTTAGCCGGAAAAACGAATGTCAAAAGGCACAGTAAAGATGCCAGTCAGCGTGTTTCCGGTCCTTTg                 >  1:57942/1‑149 (MQ=255)
                                                                                                                                cTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGG‑AATAATCTTAGCCGGAAAAACGAATGTCAAAAGGCACAGTAAAGATGCCAGTCAGCGTGTTTCCGGTCCTTTg                 <  3:63540/149‑1 (MQ=255)
                                                                                                                                 ttACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGG‑AATAATCTTAGCCGGAAAAACGAATGTCAAAAGGCACAGTAAAGATGCCAGTCAGCGTGTTTCCGGTCCTTTgg                >  7:303673/1‑149 (MQ=255)
                                                                                                                                 ttACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGG‑AATAATCTTAGCCGGAAAAACGAATGTCAAAAGGCACAGTAAAGATGCCAGTCAGCGTGTTTCCGGTCCTTTgg                >  3:248674/1‑149 (MQ=255)
                                                                                                                                  tACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGG‑AATAATCTTAGCCGGATAAACGAATGTCAAAAGGCACAGTAAAGATGCCAGTCAGCGTGTATCCGGTCCTTTGGt               >  7:336085/1‑149 (MQ=255)
                                                                                                                                   acaacaAACAGGGCGATCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACATATACTGGAAACAGGCTGG‑AATAATCTTAGCCGGAAAAACGAATGTCAAAAGGCACAGTAAAGTTGCCAGTCAGCGTGTTTCCGGTCCTTTGGt               <  4:46045/148‑1 (MQ=255)
                                                                                                                                    caacaaACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGG‑AATAATCTTAGCCGGAAAAACGAATGTCAAAAGGCACAGTAAAGATGCCAGTCAGCGTGTTTCCGGTCCTTTGGTCa             >  6:252782/1‑149 (MQ=255)
                                                                                                                                     aacaaacaGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAaca                                                                                                                <  4:171/50‑1 (MQ=255)
                                                                                                                                     aacaaacaGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAaca                                                                                                                >  3:171/1‑50 (MQ=255)
                                                                                                                                     aacaaacaGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGG‑AATAATCTTAGCCGGAAAAACGAATGTCAAAAGGCACAGTAAAGATGCCAGTCAGCGTGTTTCCGGTCCTTTGGTCAt            >  3:230159/1‑149 (MQ=255)
                                                                                                                                     aacaaacaGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGG‑AATAATCTTAGCCGGAAAAACGAATGTCAAAAGGCACAGTAAAGATGCCAGTCAGCGTGTTTCCGGTCCTTTGGTCAt            >  6:7858/1‑149 (MQ=255)
                                                                                                                                     aacaaacaGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGG‑AATAATCTTAGCCGGAAAAACGAATGTCAAAAGGCACAGTAAAGATGCCAGTCAGCGTGTTTCCGGTCCTTTGGTCAt            >  1:131146/1‑149 (MQ=255)
                                                                                                                                      acaaacaGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGG‑AATAATCTTAGCCGGAAAAACGAATGTCAAAAGGCACAGTAAAGATGCCAGTCAGCGTGTTTCCGGTCCTTTGGTCAtt           >  6:341736/1‑149 (MQ=255)
                                                                                                                                       caaacaGGGCGCTCAATGAGTGCCCTGTCTCTTTTTAAATTTAAAAC‑CTTATTATGGAAACAGGGAGGAAATATATTTACCCGAAAAAACTAAT‑TAAAAAAGCAAATGAAAAATTGCAGGAAACGTTGTTCCGGTTCGTTGGTGATTc          >  8:130014/1‑148 (MQ=255)
                                                                                                                                       caaacaGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGG‑AATAATCTTAGCCGGAAAAACGAATGTCAAAAGGCACAGTAAAGATGCCAGTCAGCGTGTTTCCGGTCCTTTGGTCATTc          <  1:279790/149‑1 (MQ=255)
                                                                                                                                        aaacaGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGG‑AATAATCTTAGCCGGAAAAACGAATGTCAAAAGGCACAGTAAAGATGCCAGTCAGCGTGTTTCCGGTCCTTTGGTCATTc          >  5:28613/1‑148 (MQ=255)
                                                                                                                                         gacaGGGCACTCATTGAGCGCCCTGTCTATTTATTAATATGAAACACTTATACTGGAAACAGGCTGG‑AATAATCTTAGCCGGAAAAACGAATGTCAAAAGGCACAGTAAAGATGCCAGTCAGCGTGTTTCCGGACGTTTGTTCATTaaa        >  5:16798/2‑149 (MQ=255)
                                                                                                                                          acaGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGG‑AATAATCTTAGCCGGAAAAACGAATGTCAAAAGGCACAGTAAAGATGCCAGTCAGCGTGTTTCCGGTCCTTTGGTCATTCAAt       <  7:351606/149‑1 (MQ=255)
                                                                                                                                            aGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGG‑AATAATCTTAGCCGGAAAAACGAATGTCAAAAGGCACAGTAAAGATGCCAGTCAGCGTGTTTCCGGTCCTTTGGTCATTCAATTg     >  1:20048/1‑149 (MQ=255)
                                                                                                                                              ggCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGG‑AATAATCTTAGCCGGAAAAACGAATGTCAAAAGGCACAGTAAAGATGCCAGTCAGCGTGTTTCCGGTCCTTTGGTCATTCAATTGCt   <  6:230151/149‑1 (MQ=255)
                                                                                                                                               gcgcTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGG‑AATAATCTTAGCCGGAAAAACGAATGTCAAAAGGCACAGTAAAGATGCCAGTCAGCGTGTTTCCGGTCCTTTGGTCATTCAATTGCtt  >  1:208402/1‑149 (MQ=255)
                                                                                                                                               gcgcTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGG‑AATAATCTTAGCCGGAAAAACGAATGTCAAAAGGCACAGTAAAGATGCCAGTCAGCGTGTTTCCGGTCCTTTGGTCATTCAATTGCtt  <  6:328166/149‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                   
CACAAATTCAGCGCGATCAGGTTGAAGATTATGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGAGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGG‑AATAATCTTAGCCGGAAAAACGAATGTCAAAAGGCACAGTAAAGATGCCAGTCAGCGTGTTTCCGGTCCTTTGGTCATTCAATTGCTT  >  NC_000913/4227405‑4227696

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