Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,227,555 A→T 6.1% intergenic (+44/‑176) 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,5550AT 6.1% 324.7 / 8.5 99intergenic (+44/‑176)metH/yjbBhomocysteine‑N5‑methyltetrahydrofolate transmethylase, B12‑dependent/putative Na+/Pi‑cotransporter
Reads supporting (aligned to +/- strand):  ref base A (47/46);  new base T (3/3);  total (50/49)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.63e-01

Only 100 of 102 total aligned reads displayed.

CAAATTCAGCGCGATCAGGTTGAAGATTATGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGAGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGG‑AATAATCTTAGCCGGAAAAACGAATGTCAAAAGGCACAGTAAAGATGCCAGTCAGCGTGTTTCCGGTCCTTTGGTCATTCAATTGCTTCAG  >  NC_000913/4227407‑4227699
                                                                                                                                                    |                                                                                                                                                 
caAATTCAGCGCGATCAGGTTGAAGATTATGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGAGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCaa                                                                                                                                                   >  2:271561/1‑149 (MQ=255)
 aaaTTCAGCGCGATCAGGTTGAAGATAATGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGAGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAAt                                                                                                                                                  <  2:194264/149‑1 (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)
                aGGTTGAAGATTATGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGAGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGtctc                                                                                                                                    <  3:153604/148‑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)
                             agCCCGCCGTAAAGGTATGAGCGTTACCGAAGTGGAGCGCTGGCTGGCACCGAATCTGGGGGATGACGCGGACTGATTCACAAATCCGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGaa                                                                                                                      <  6:297662/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                                                                                                                                           <  7:129049/80‑1 (MQ=255)
                                                                             aCCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCACTCATTGAGCGcc                                                                                                                                           >  8:129049/1‑80 (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                                                                                                                                         <  2:229735/74‑1 (MQ=38)
                                                                                     tGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCACTCATTGAGCGCCCt                                                                                                                                         >  1:229735/1‑74 (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                         >  2:294672/1‑149 (MQ=255)
                                                                                                                         ttttCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGG‑AATAATCTTAGCCGGAAAAACGAATGTCAAAAGGCACAGTAAAGATGCCAGTCAGCGTGTTTCCGGTc                         >  6:103184/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                                                                                                                   >  3:171/1‑50 (MQ=255)
                                                                                                                                   aacaaacaGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAaca                                                                                                                   <  4:171/50‑1 (MQ=255)
                                                                                                                                   aacaaacaGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGG‑AATAATCTTAGCCGGAAAAACGAATGTCAAAAGGCACAGTAAAGATGCCAGTCAGCGTGTTTCCGGTCCTTTGGTCAt               >  6:7858/1‑149 (MQ=255)
                                                                                                                                   aacaaacaGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGG‑AATAATCTTAGCCGGAAAAACGAATGTCAAAAGGCACAGTAAAGATGCCAGTCAGCGTGTTTCCGGTCCTTTGGTCAt               >  1:131146/1‑149 (MQ=255)
                                                                                                                                   aacaaacaGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGG‑AATAATCTTAGCCGGAAAAACGAATGTCAAAAGGCACAGTAAAGATGCCAGTCAGCGTGTTTCCGGTCCTTTGGTCAt               >  3:230159/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)
                                                                                                                                        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)
                                                                                                                                                 tCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGG‑AATAATCTTAGCCGGAAAAACGAATGTCAAAAGGCACAGTAAAGATGCCAGTCAGCTTGTTTCCTGTCCTTTGGTCATTCAATTGCTTCat  >  6:376733/1‑147 (MQ=255)
                                                                                                                                                   aaTGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGG‑AATAATCTTAGCCGGAAAAACGAATGTCAAAAGGCACAGTAAAGATGCCAGTCAGCGTGTTTCCGGTCCTTTGGTCATTCAATTGc       <  8:179439/141‑1 (MQ=255)
                                                                                                                                                   aaTGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGG‑AATAATCTTAGCCGGAAAAACGAATGTCAAAAGGCACAGTAAAGATGCCAGTCAGCGTGTTTCCGGTCCTTTGGTCATTCAATTGc       >  7:179439/1‑141 (MQ=255)
                                                                                                                                                    |                                                                                                                                                 
CAAATTCAGCGCGATCAGGTTGAAGATTATGCCCGCCGTAAAGGTATGAGCGTTACCGAAGTTGAGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGG‑AATAATCTTAGCCGGAAAAACGAATGTCAAAAGGCACAGTAAAGATGCCAGTCAGCGTGTTTCCGGTCCTTTGGTCATTCAATTGCTTCAG  >  NC_000913/4227407‑4227699

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: ‑