Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,227,550 G→A 17.3% 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,5500GA17.3% 454.2 / 74.6 190intergenic (+39/‑181)metH/yjbBhomocysteine‑N5‑methyltetrahydrofolate transmethylase, B12‑dependent/putative Na+/Pi‑cotransporter
Reads supporting (aligned to +/- strand):  ref base G (83/74);  new base A (15/18);  total (98/92)
Fisher's exact test for biased strand distribution p-value = 4.51e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.78e-01

Only 100 of 197 total aligned reads displayed.

ATGAGCGTTACCGAAGTTGAGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATGTCAAAAGGCACAGTA  >  NC_000913/4227452‑4227649
                                                                                                  |                                                                                                   
 tGAGCGTTACCGAAGTTGAGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCt                                                                                                   >  2:2450526/1‑100 (MQ=255)
   aGCGTTACCGAAGTTGAGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCCCTTTTCCTTACAACAAACAGGGCGCTc                                                                                                  >  2:602549/1‑99 (MQ=255)
   aGCGTTACCGAAGTTGAGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCa                                                                                                 <  2:949924/100‑1 (MQ=255)
   aGCGTTACCGAAGTTGAGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCa                                                                                                 >  2:182840/1‑100 (MQ=255)
     cGTTACCGAAGTTGAGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAAt                                                                                               <  1:1068855/100‑1 (MQ=255)
       ttACCGAAGTTGAGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGa                                                                                             >  1:271607/1‑100 (MQ=255)
          ccGAAGTTGAGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGCTTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTg                                                                                          >  2:3524988/1‑100 (MQ=255)
          ccGAAGTTGAGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTg                                                                                          <  1:3117477/100‑1 (MQ=255)
          ccGAAGTTGAGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTg                                                                                          <  1:535146/100‑1 (MQ=255)
          ccGAAGTTGAGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCACTCATTGAgcg                                                                                          >  2:3711431/1‑98 (MQ=255)
             aaGTTGAGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCACTCATTGAGCGccc                                                                                       <  2:2102849/100‑1 (MQ=255)
             aaGTTGAGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTCCTTACAACAAACAGGGCACTCATTGAGCGCCCt                                                                                      <  2:685543/100‑1 (MQ=255)
               gTTGAGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTg                                                                                     <  2:2842071/100‑1 (MQ=255)
               gTTGAGCGCTGGCTGGAACCGAATCTGGGGTATGACGTGGACTGATTCCCAACACTATCACTTTTCCTTACAACAAACAGAGCGCTCAATGAGTGCCCt                                                                                      >  2:370441/1‑99 (MQ=255)
                ttGAGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCACTCATTGAGCGCCCTGt                                                                                    <  2:2640306/100‑1 (MQ=255)
                 tGAGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGtc                                                                                   >  1:86086/1‑100 (MQ=255)
                 tGAGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGtc                                                                                   >  2:597179/1‑100 (MQ=255)
                 tGAGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCACTCATTGAGCGCCCTGtt                                                                                   <  1:3711452/100‑2 (MQ=39)
                  gAGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGtct                                                                                  <  1:618662/100‑1 (MQ=255)
                  gAGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCACTCATTGAGCGCCCTGttt                                                                                  <  1:3209454/100‑3 (MQ=39)
                   aGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGtct                                                                                  <  2:1204923/99‑1 (MQ=255)
                   aGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGAGCTCAATGAGTGCCCTGtctc                                                                                 >  1:3071712/1‑100 (MQ=255)
                   aGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTACTTACAACAAACAGGGCGCTCAATGAGTGCCCTGtctc                                                                                 >  1:2081673/1‑100 (MQ=255)
                   aGCGCTGCCGCGCACCGAATCTGGGATATGACCCGGACTCATTCACAACTCTTTCACTTTTCCTTACAACACACAGGGCACTCATTGAGCGCCCTGtt                                                                                   <  2:2542223/98‑2 (MQ=19)
                    gcgcTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCACTCATTGAGCGCCCTGtttgt                                                                                <  1:3509421/100‑5 (MQ=38)
                     cgcTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCACTCATTGAGCGCCCTGtttgtt                                                                               <  1:225161/100‑6 (MQ=38)
                      gctggctgGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCACTCATTGAGCGCCCTGtttgttg                                                                              <  2:2810192/100‑7 (MQ=37)
                      gctggctgGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCACTCATTGAGCGCCCTGtttgttg                                                                              <  1:1334850/100‑7 (MQ=37)
                       ctggctggCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTtta                                                                             <  1:4004505/100‑1 (MQ=255)
                        tggctggcACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTttat                                                                            <  1:138307/100‑1 (MQ=255)
                         ggctggcACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTttatt                                                                           <  1:7457/100‑1 (MQ=255)
                         ggctggcACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTctatt                                                                           <  2:1657282/100‑1 (MQ=255)
                          gctggcACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTttatta                                                                          >  1:2147580/1‑100 (MQ=255)
                           ctggcACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTaa                                                                         >  2:1692826/1‑100 (MQ=255)
                           ctggcACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTaa                                                                         >  1:3138599/1‑100 (MQ=255)
                           ctggcACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTaa                                                                         >  1:1797794/1‑100 (MQ=255)
                           ctggcACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTaa                                                                         >  1:525865/1‑100 (MQ=255)
                            tggcACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAat                                                                        <  1:1175790/100‑1 (MQ=255)
                               cACCGAATCTTGGGTATTACGCGGACTGATTCACAAATCTGTCACTTGTCCTTACATCAAACAGGGCGCTCAATGCGTGCCCTGTCTCTctatta                                                                          <  2:2576090/95‑1 (MQ=255)
                               cACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATg                                                                     >  2:4098570/1‑100 (MQ=255)
                                 ccGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGaa                                                                   <  2:1573933/100‑1 (MQ=255)
                                 ccGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGaa                                                                   >  1:3385900/1‑100 (MQ=255)
                                 ccGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGaa                                                                   <  1:3031083/100‑1 (MQ=255)
                                  cGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGaaa                                                                  >  1:158737/1‑100 (MQ=255)
                                  cGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGaaa                                                                  <  2:3116534/100‑1 (MQ=255)
                                       cTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAacacct                                                             <  1:2824065/100‑3 (MQ=255)
                                       cTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACtt                                                             >  1:2361947/1‑100 (MQ=255)
                                        tGGTGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTta                                                            <  2:1991903/100‑1 (MQ=255)
                                        tGGGGTATGACGCGGACTGATTCACCACTCTGTCACTTTCCCTTCCAACAAACAGGCCGCTCAATGAGTGCCCTGTCTCTTCATTAATATGAAACACtt                                                             <  2:3344148/99‑1 (MQ=255)
                                         ggggTATTACGCGGACTGATTCACAAATCTGTCATTTTTCCTTACAATAAACAGGGGGCTCAATGAGTGCAGTGTCTCTTTATTAATATGAAACACTtag                                                           <  2:174497/100‑2 (MQ=255)
                                           ggTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATAc                                                         <  2:2134227/100‑1 (MQ=255)
                                            gTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACt                                                        <  2:525864/100‑1 (MQ=255)
                                             tATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTg                                                       >  1:1005822/1‑100 (MQ=255)
                                             tATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTg                                                       >  2:3041008/1‑100 (MQ=255)
                                              aTGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTgg                                                      <  2:111661/100‑1 (MQ=255)
                                              aTGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTgg                                                      <  1:4014765/100‑1 (MQ=255)
                                               tGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACCACAAACCGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGa                                                     >  2:1025305/1‑100 (MQ=255)
                                               tGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGa                                                     <  1:2190309/100‑1 (MQ=255)
                                                gACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGaa                                                    >  2:4046287/1‑100 (MQ=255)
                                                gACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGaa                                                    <  1:3181943/100‑1 (MQ=255)
                                                  cgcgGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAAc                                                  <  1:1966253/100‑1 (MQ=255)
                                                     ggACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAgg                                               <  1:1606643/100‑1 (MQ=255)
                                                        cTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTg                                            <  2:3470303/100‑1 (MQ=255)
                                                         tGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTgg                                           <  2:184779/100‑1 (MQ=255)
                                                         tGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTgg                                           <  2:592108/100‑1 (MQ=255)
                                                         atcttCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTgg                                           <  2:847199/97‑1 (MQ=255)
                                                          gATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGa                                          >  1:3949874/1‑100 (MQ=255)
                                                          gATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGa                                          >  2:2988836/1‑100 (MQ=255)
                                                           aTTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGaa                                         <  2:641816/100‑1 (MQ=255)
                                                           aTTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGaa                                         <  2:393664/100‑1 (MQ=255)
                                                           aTTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGaa                                         <  1:838429/100‑1 (MQ=255)
                                                            ttCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATCTGAAACACTTATACTGGAAACAGGCTGGaat                                        >  1:444165/1‑100 (MQ=255)
                                                            ttCACAAATCTGTAACTTTTCCTTACAACAAACAGGGCGATAAATGAGTGCACTGTCTCTTTATTAATCTGAAACACTTATACTGGAAACAAGCTAGa                                          >  2:2933249/1‑98 (MQ=255)
                                                              cacaAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGaataa                                      >  1:2356062/1‑100 (MQ=255)
                                                              cacaAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGaataa                                      >  2:2348092/1‑100 (MQ=255)
                                                                caAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATc                                    >  1:1702016/1‑100 (MQ=255)
                                                                      tGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGcc                              <  2:3004508/100‑1 (MQ=255)
                                                                       gTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCg                             <  2:2463535/100‑1 (MQ=255)
                                                                       gTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCg                             >  1:1253615/1‑100 (MQ=255)
                                                                        tCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGAGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCgg                            >  1:2963075/1‑100 (MQ=255)
                                                                           cTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGGAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGaaa                         >  2:2675685/1‑100 (MQ=255)
                                                                             tttCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGaaaaa                       >  2:228448/1‑100 (MQ=255)
                                                                               tCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACg                     <  2:3245704/100‑1 (MQ=255)
                                                                                ccTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGa                    >  1:761026/1‑100 (MQ=255)
                                                                                 cTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGaa                   <  1:3194286/100‑1 (MQ=255)
                                                                                 cTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGaa                   >  1:3432343/1‑100 (MQ=255)
                                                                                   ataaagagacaGGGCACTCATTGAGCGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAAt                  >  2:3498689/9‑99 (MQ=37)
                                                                                     caacaaACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATGTc               >  2:1383535/1‑100 (MQ=255)
                                                                                        caaacaGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATGTCaaa            >  1:3643596/1‑100 (MQ=255)
                                                                                        caaacaGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATGTCaaa            >  2:3892342/1‑100 (MQ=255)
                                                                                          gacaGGGCACTCATTGAGCGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATGTCAAAAg          >  1:1370122/2‑100 (MQ=39)
                                                                                          gacaGGGCACTCATTGAGCGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATGTCAAAAg          >  1:2257763/2‑100 (MQ=39)
                                                                                          aacaGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATGTCAAAAg          <  1:937183/100‑1 (MQ=255)
                                                                                           acaGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAACATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATGTCAAAAgg         <  1:2939470/100‑1 (MQ=255)
                                                                                           acaGGGCACTCATTGAGCGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATGTCAAAAgg         >  2:692986/1‑100 (MQ=255)
                                                                                            caGGGCGCTCAATGCGTGCCCTGTCTCTTTATTAATATGCAACACTTATACTGGAACCAGGCTGGAATAATCTTATCCGGAAAAACGAATGTCTAAAGGc        >  2:1921837/1‑100 (MQ=255)
                                                                                            caGGGCACTCATTGAGCGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATGTCAAAAGGc        >  2:1772776/1‑100 (MQ=255)
                                                                                             aGGGCACTCATTGAGCGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATGTCAAAAGGca       >  2:3153188/1‑100 (MQ=255)
                                                                                              gggCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATGTCAAAAGGcac      <  1:3981662/100‑1 (MQ=255)
                                                                                                 agcTCAATGAGTGCCCTGTCTCTGTATTAATATGAAACACTTATAATGGAAACAGGCTGGAATAATCTTAGCCGGAAAAAGGAATGTCAAAAGGcaca     >  2:1331838/2‑98 (MQ=255)
                                                                                                  |                                                                                                   
ATGAGCGTTACCGAAGTTGAGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATGTCAAAAGGCACAGTA  >  NC_000913/4227452‑4227649

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