Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,477,885 C→A 6.0% I13I (ATC→ATA dsdX → D‑serine transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,477,8850CA 6.0% 224.2 / 2.9 67I13I (ATC→ATAdsdXD‑serine transporter
Reads supporting (aligned to +/- strand):  ref base C (32/31);  new base A (2/2);  total (34/33)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.54e-01

TCACTCCATTGATGCAATTTACCTCATGTGAAAGGCAAATTTTATCGTTTGTCAGCCTGCGTTGTTTTTTTGTCCAATATCATCAGGTTAATCACAGGGGAAGGTGAGATATGCACTCTCAAATCTGGGTTGTGAGCACGCTGCTTATCAGCATCGTGTTAATTGTACTGACCATCGTGAAGTTCAAATTCCACCCGTTTCTGGCGCTGTTGCTGGCCAGCTTCTTCGTGGGAACGATGATGGGCATGGGGCCACTGGATATGGTAAATGCTATTGAAAGTGGAATTGGCGGAACGCT  >  NC_000913/2477737‑2478034
                                                                                                                                                    |                                                                                                                                                     
tCACTCCATTGATGCAATTTACCTCATGTTAAAGGCAAATTTTATCGTTTGTCAGCCTGCGTTGTTTTTTTGTCCAATATCATCAGGTTAATCACAGGGGAAGGTGAGATATGCACTCTCAAATCTGGGTTGTGAGCACGCTGCTTATCa                                                                                                                                                      <  1:1263395/150‑1 (MQ=255)
     ccATTGATGCAATTTACCTCATGTGAAAGGCAAATTTTATCGTTTGTCAGCCTGCGTTGTTTTTTTGTCCAATATCATCAGGTTAATCACAGGGGAAGGTGAGATATGCACTCTCAAATCTGGGTTGTGAGCACGCTGCTTATCAGCATc                                                                                                                                                 >  2:1022817/1‑150 (MQ=255)
       aTTGATGCAATTTACCTCATGTGAAAGGCAAATTTTATCGTTTGTCAGCCTGCGTTGTTTTTTTGTCCAATATCATCAGGTTAATCACAGGGGAAGGTGAGATATGCACTCTCAAATCTGGGTTGTGAGCACGCTGCTTATCAGCATCgt                                                                                                                                               <  2:355162/150‑1 (MQ=255)
        ttGATGCAATTTACCTCATGTGAAAGGCAAATTTTATCGTTTGTCAGCCTGCGTTGTTTTTTTGTCCAATATCATCAGGTTAATCACAGGGGAAGGTGAGATATGCACTCTCAAATCTGGGTTGTGAGCACGCTGCTTATCAGCATCgtg                                                                                                                                              <  1:822625/150‑1 (MQ=255)
          gATGCAATTTACCTCATGTGAAAGGCAAATTTTATCGTTTGTCAGCCTGCGTTGTTTTTTTGTCCAATATCATCAGGTTAATCACAGGGGAAGGTGAGATATGCACTCTCAAATCTGGGTTGTGAGCACGCTGCTTATCAGCATCGTGtt                                                                                                                                            <  2:26166/150‑1 (MQ=255)
           aTGCAATTTACCTCATGTGAAAGGCAAATTTTATCGTTTGTCAGCCTGCGTTGTTTTTTTGTCCAATATCATCAGGTTAATCACAGGGGAAGGTGAGATATGCACTCTCAAATCTGGGTTGTGAGCACGCTGCTTATCAGCATCGTGTTa                                                                                                                                           >  2:490772/1‑150 (MQ=255)
            tGCAATTTACCTCATGTGAAAGGCAAATTTTATCGTTTGTCAGCCTGCGTTGTTTTTTTGTCCAATATCATCAGGTTAATCACAGGGGAAGGTGAGATATGCACTCTCAAATCTGGGTTGTGAGCACGCTGCTTATCAGCATCGTGTTaa                                                                                                                                          >  1:606981/1‑150 (MQ=255)
             gCAATTTACCTCATGTGAAAGGCAAATTTTATCGTTTGTCAGCCTGCGTTGTTTTTTTGTCCAATATCATCAGGTTAATCACAGGGGAAGGTGAGATATGCACTCTCAAATCTGGGTTGTGAGCACGCTGCTTATCAGCATCGTGTTAAt                                                                                                                                         <  2:2665/150‑1 (MQ=255)
                aTTTACCTCATGTGAAAGGCAAATTTTATCGTTTGTCAGCCTGCGTTGTTTTTTTGTCCAATATCATCAGGTTAATCACAGGGGAAGGTGAGATATGCACTCTCAAATCTGGGTTGTGAGCACGCTGCTTATCAGCATCGTGTTAATTGt                                                                                                                                      <  1:490907/150‑1 (MQ=255)
                  ttACCTCATGTGAAAGGCAAATTTTATCGTTTGTCAGCCTGCGTTGTTTTTTTGTCCAATATCATCAGGTTAATCACAGGGGAAGGTGAGATATGCACTCTCAAATCTGGGTTGTGAGCACGCTGCTTATCAGCATCGTGTTAATTGTAc                                                                                                                                    <  2:262478/150‑1 (MQ=255)
                   tACCTCATGTGAAAGGCAAATTTTATCGTTTGTCAGCCTGCGTTGTTTTTTTGTCCAATATCATCAGGTTAATCACAGGGGAAGGTGAGATATGCACTCTCAAATCTGGGTTGTGAGCACGCTGCTTATCAGCATCGTGTTAATTGTACt                                                                                                                                   >  2:1271873/1‑150 (MQ=255)
                     ccTCATGTGAAAGGCAAATTTTATCGTTTGTCAGCCTGCGTTGTTTTTTTGTCCAATATCATCAGGTTAATCACAGGGGAAGGTGAGATATGCACTCTCAAATCTGGGTTGTGAGCACGCTGCTTATCAGCATCGTGTTAATTGTACTGa                                                                                                                                 >  1:1130181/1‑150 (MQ=255)
                       tCATGTGAAAGGCAAATTTTATCGTTTGTCAGCCTGCGTTGTTTTTTTGTCCAATATCATCAGGTTAATCACAGGGGAAGGTGAGATATGCACTCTCAAATCTGGGTTGTGAGCACGCTGCTTATCAGCATCGTGTTAATTGTACTGAcc                                                                                                                               >  1:161609/1‑150 (MQ=255)
                           gtgAAAGGCAAATTTTATCGTTTGTCAGCCTGCGTTGTTTTTTTGTCCAATATCATCAGGTTAATCACAGGGGAAGGTGAGATATGCACTCTCAAATCTGGGTTGTGAGCACGCTGCTTATCAGCATCgtg                                                                                                                                              >  2:603140/1‑131 (MQ=255)
                           gtgAAAGGCAAATTTTATCGTTTGTCAGCCTGCGTTGTTTTTTTGTCCAATATCATCAGGTTAATCACAGGGGAAGGTGAGATATGCACTCTCAAATCTGGGTTGTGAGCACGCTGCTTATCAGCATCgtg                                                                                                                                              <  1:603140/131‑1 (MQ=255)
                             gAAAGGCAAATTTTATCGTTTGTCAGCCTGCGTTGTTTTTTTGTCCAATATCATCAGGTTAATCACAGGGGAAGGTGAGATATGCACTCTCAAATCTGGGTTGTGAGCACGCTGCTTATCAGCATCGTGTTAATTGTACTGACCATCGTg                                                                                                                         >  2:525219/1‑150 (MQ=255)
                              aaaGGCAAATTTTATCGTTTGTCAGCCTGCGTTGTTTTTTTGTCCAATATCATCAGGTTAATCACAGGGGAAGGTGAGATATGCACTCTCAAATCTGGGTTGTGAGCACGCTGCTTATCAGCATCGTGTTAATTGTACTGACCATCGTGa                                                                                                                        <  2:149965/150‑1 (MQ=255)
                                       ttttATCGTTTGTCAGCCTGCGTTGTTTTTTTGTCCAATATCATCAGGTTAATCACAGGGGAAGGTGAGATATGCACTCTCAAATCTGGGTTGTGAGCACGCTGCTTATCAGCATCGTGTTAATTGTACTGACCATCGTGAAGTTCAAAt                                                                                                               <  2:659920/150‑1 (MQ=255)
                                             cGTTTGTCAGCCTGCGTTGTTTTTTTGTCCAATATCATCAGGTTAATCACAGGGGAAGGTGAGATATGCACTCTCAAATCTGGGTTGTGAGCACGCTGCTTATCAGCATCGTGTTAATTGTACTGACCATCGTGAAGTTCAAATTCCAcc                                                                                                         >  2:700702/1‑150 (MQ=255)
                                                  gTCAGCCTGCGTTGTTTTTTTGTCCAATATCATCAGGTTAATCACAGGGGAAGGTGAGATATGCACTCTCAAATCTGGGTTGTGAGCACGCTGCTTATCAGCATCGTGTTAATTGTACTGACCATCGTGAAGTTCAAATTCCACCCGttt                                                                                                    >  2:130566/1‑150 (MQ=255)
                                                  gTCAGCCTGCGTTGTTTTTTTGTCCAATATCATCAGGTTAATCACAGGGGAAGGTGAGATATGCACTCTCAAATCTGGGTTGTGAGCACGCTGCTTATCAGCATCGTGTTAATTGTACTGACCATCGTGAAGTTCAAATTCCACCCGttt                                                                                                    >  1:1067771/1‑150 (MQ=255)
                                                  gTCAGCCTGCGTTGTTTTTTTGTCCAATATCATCAGGTGAATCACAGGGGAAGGTGAGATATGCACTCTCAAATCTGGGTTGTGAGCACGCTGCTTATCAGCATCGTGTTAATACTACTGACCATCGTGAAGTTCAAATTCCACCCGttt                                                                                                    >  2:1060323/1‑150 (MQ=255)
                                                    cAGCCTGCGTTGTTTTTTTGTCCAATATCATCAGGTTAATCACAGGGGAAGGTGAGATATGCACTCTCAAATCTGGGTTGTGAGCACGCTGCTTATCAGCATCGTGTTAATTGTACTGACCATCGTGAAGTTCAAATTCCACCCGTTTCt                                                                                                  <  2:1130181/150‑1 (MQ=255)
                                                     cgCCTGCGTTGTTTTTTTGTCCAATATCATCAGGTTAATCACAGGGGAAGGTGAGATATTCACTCTCAAATCTGGGTTGTGAGCACGCTGCTTATCAGCATCGTGTTAATTGTACTGACCATCGTGAAGTTCAAATTCCACCCGTTTCTg                                                                                                 <  2:751761/149‑1 (MQ=255)
                                                     aGCCTGCGTTGTTTTTTTGTGCAATATCATCAGGTTAATCACAGGGGAAGGTGAGATATGCACTATAAAATCTGGGTTGTGAGCAAGCTGCTTAACAGCATCGTGTTAAATGTACTGACCATCGTGAAGTTCAAATTCCACCCGATTCTg                                                                                                 >  1:413533/1‑150 (MQ=255)
                                                        cTGCGTTGTTTTTTTGTCCAATATCATCAGGTTAATCACAGGGGAAGGTGAGATATGCACTCTCAAATCTGGGTTGTGAGCACGCTGCTTATCAGCATCGTGTTAATTGTACTGACCATCGTGAAGTTCAAATTCCACCCGTTTCTGgcg                                                                                              >  2:266204/1‑150 (MQ=255)
                                                         tGCGTTGTTTTTTTGTCCAATATCATCAGGTTAATCACAGGGGAAGGTGAGATATGCACTCTCAAATCTGGGTTGTGAGCACGCTGCTTATCAGCATCGTGTTAATTGTACTGACCATCGTGAAGTTCAAATTCCAc                                                                                                          <  1:482855/137‑1 (MQ=255)
                                                         tGCGTTGTTTTTTTGTCCAATATCATCAGGTTAATCACAGGGGAAGGTGAGATATGCACTCTCAAATCTGGGTTGTGAGCACGCTGCTTATCAGCATCGTGTTAATTGTACTGACCATCGTGAAGTTCAAATTCCAc                                                                                                          >  2:482855/1‑137 (MQ=255)
                                                         tGCGTTGTTTTTTTGTCCAATATCATCAGGTTAATCACAGGGGAAGGTGAGATATGCACTCTCAAATCTGGGTTGTGAGCACGCTGCTTATCAGCATCGTGTTAATTGTACTGACCATCGTGAAGTTCAAATTCCACCCGTTTCTGgcgc                                                                                             <  2:331317/150‑1 (MQ=255)
                                                               gttTTTTTGTCCAATATCATCAGGTTAATCACAGGGGAAGGTGAGATATGCACTCTCAAATCTGGGTTGTGAGCACGCTGCTTATCAGCATCGTGTTAATTGTACTGACAATCGTGAAGTTCAAATTCCACCCGTTTCTGGCGCTGTTGc                                                                                       >  2:870946/1‑150 (MQ=255)
                                                                   ttttGTCCAATATCATCAGGTTAATCACAGGGGAAGGTGAGATATGCACTCTCAAATCTGGGTTGTGAGCACGCTGCTTATCAGCATCGTGTTAATTGTACTGACCATCGTGAAGTTCAAATTCCACCCGTTTCTGGCGCTGTTGCTGGc                                                                                   <  1:1022817/150‑1 (MQ=255)
                                                                       gTCCAATATCATCAGGTTAATCACAGGGGAAGGTGAGATATGCACTCTCAAATCTGGGTTGTGAGCACGCTGCTTATCAGCATCGTGTTAATTGTACTGACCATCGTGAAGTTCAAATTCCACCCGTTTCTGGCGCTGTTGCTGGCCAGc                                                                               >  2:404656/1‑150 (MQ=255)
                                                                              atcatcAGGTTAATCACAGGGGAAGGTGAGATATGCACTCTCAAATCTGGGTTGTGAGCACGCTGCTTATCAGCATCGTGTTAATTGTACTGACCATCGTGAAGTTCAAATTCCACCCGTTTCTGGCGCTGTTGCTGGCCAGCTTCTTCg                                                                        <  1:892168/150‑1 (MQ=255)
                                                                                 atcaGGTTAATCACAGGGGAAGGTGAGATATGCACTCTCAAATCTGGGTTGTGAGCACGCTGCTTATCAGCATCGTGTTAATTGTACTGACCATCGTGAAGTTCAAATTCCACCCGTTTCTGGCGCTGTTGCTGGCCAGCTTCTTCGTgg                                                                     <  1:293570/150‑1 (MQ=255)
                                                                                   caGGTTAATCACAGGGGAAGGTGAGATATGCACTCTCAAATCTGGGTTGTGAGCACGCTGCTTATCAGCATCGTGTTAATTGTACTGACCATCGTGAAGTTCAAATTCCACCCGTTTCTGGCGCTGTTGCTGGCCAGCTTCTTCGTGGGa                                                                   >  1:290772/1‑150 (MQ=255)
                                                                                     ggTTAATCACAGGGGAAGGTGAGATATGCACTCTCAAATCTGGGTTGTGAGCACGCTGCTTATCAGCATCGTGTTAATTGTACTGACCATCGTGAAGTTCAAATTCCACCCGTTTCTGGCGCTGTTGCTGGCCAGCTTCTTCGTGGGAAc                                                                 <  2:1367332/150‑1 (MQ=255)
                                                                                        tAATCACAGGGGAAGGTGAGATATGCACTCTCAAATCTGGGTTGTGAGCACGCTGCTTATCAGCATCGTGTTAATTGTACTGACCATCGTGAAGTTCAAATTCCACCCGTTTCTGGCGCTGTTGCTGGCCAGCTTCTTCGTGGGAACgat                                                              >  2:1003580/1‑150 (MQ=255)
                                                                                        tAATCACAGGGGAAGGTGAGATATGCACTCTCAAATCTGGGTTGTGAGCACGCTGCTTATCAGCATCGTGGTAATTGTACTGACCATCGTGAAGTTCAAATTCCCCCCGTTTCTGGCGCTGTTGCTGGGCAGCTTCTTCGTGGGAACgat                                                              >  2:1342287/1‑150 (MQ=255)
                                                                                            cacaGGGGAAGGTGAGATATGCACTCTCAAATCTGGGTTGTGAGCACGCTGCTTATCAGCATCGTGTTAATTGTACTGACCATCGTGAAGTTCAAATTCCACCCGTTTCTGGCGCTGTTGCTGGCCAGCTTCTTCGTGGGAACGatgatg                                                          >  2:1315340/1‑150 (MQ=255)
                                                                                            cacaGGGGAAGGTGAGATATGCACTCTCAAATCTGGGTTGTGAGCACGCTGCTTATCAGCATCGTGTTAATTGTACTGACCATCGTGAAGTTCAAATTCCACCCGTTTCTGGCGCTGTTGCTGGCCAGCTTCTTCGTGGGAACGatgatg                                                          >  2:94950/1‑150 (MQ=255)
                                                                                              caGGGGAAGGTGAGATATGCACTCTCAAATCTGGGTTGTGAGCACGCTGCTTATCAGCATCGTGTTAATTGTACTGACCATCGTGAAGTTCAAATTCCACCCGTTTCTGGCGCTGTTGCTGGCCAGCTTCTTCGTGGGAACGATGATggg                                                        >  1:1104288/1‑150 (MQ=255)
                                                                                                  ggAAGGTGAGATATGCACTCTCAAATCTGGGTTGTGAGCACGCTGCTTATCAGCATCGTGTTAATTGTACTGACCATCGTGAAGTTCAAATTCCACCCGTTTCTGGCGCTGTTGCTGGCCAGCTTCTTCGTGGGAACGATGATGGGCATg                                                    >  1:633789/1‑150 (MQ=255)
                                                                                                   gAAGGTGAGATATGCACTCTCAAATCTGGGTTGTGAGCACGCTGCTTATCAGCATCGTGTTAATTGTACTGACCATCGTGAAGTTCAAATTCCACCCGTTTCTGGCGCTGTTGCTGGCCAGCTTCTTCGTGGGAACGATGATGGGCATgg                                                   >  1:1192502/1‑150 (MQ=255)
                                                                                                    aaGGTGAGATATGCACTCTCAAATCTGGGTTGTGAGCACGCTGCTTATCAGCATCGTGTTAATTGTACTGACCATCGTGAAGTTCAAATTCCACCCGTTTCTGGCGCTGTTGCTGGCCAGCTTCTTCGTGGGAACGATGATGGGCATggg                                                  >  1:1091895/1‑150 (MQ=255)
                                                                                                     aGGTGAGATATGCACTCTCAAATCTGGGTTGTGAGCACGCTGCTTATCAGCATAGTGTTAATTGTACTGACCATCGTGAAGTTCAAATTCCACCCGTTTCTGGCGCTGTTGCTGGACAGCTTCTTCGTGGGAACGATGATGGGCATgggg                                                 >  1:1280333/1‑150 (MQ=255)
                                                                                                               tGCACTCTCAAATCTGGGTTGTGAGCACGCTGCTTATCAGCATCGTGTTAATTGTACTGACCATCGTGAAGTTCAAATTCCACCCGTTTCTGGCGCTGTTGCTGGCCAGCTTCTTCGTGGGAACGATGATGGGCATGGGGCCACTGGata                                       <  1:525219/150‑1 (MQ=255)
                                                                                                                gCACTCTCAAATCTGGGTTGTGAGCACGCTGCTTATCAGCATCGTGTTAATTGTACTGACCATCGTGAAGTTCAAATTCCACCCGTTTCTGGCGCTGTTGCTGGCCAGCTTCTTCGTGGTAACGATGATGGGCATGGGGCCACTGGatat                                      <  1:1408949/150‑1 (MQ=255)
                                                                                                                gCACTCTCAAATCTGGGTTGTGAGCACGCTGCTTATCAGCATCGTGTTAATTGTACTGACCATCGTGAAGTTCAAATTCCACCCGTTTCTGGCGCTGTTGCTGGCCAGCTTCTTCGTGGGAACGATGATGGGCATGGGGCCACTGGatat                                      >  1:568096/1‑150 (MQ=255)
                                                                                                                 cACTCTCAAATCTGGGTTGTGAGCACGCTGCTTATCAGCATCGTGTTAATTGTACTGACCATCGTGAAGTTCAAATTCCACCCGTTTCTGGCGCTGTTGCTGGCCAGCTTCTTCGTGGGAACGATGATGGGCATGGGGCCACTGGATATg                                     >  2:240638/1‑150 (MQ=255)
                                                                                                                   ctttaaaaTTTTGGTTTTGATGACGCTGGTTATCATCATCGTGTTTATTTTTATGCACATCGTGTAGTTTAAATTCCACCCGTTTTTTGCGCTTTTGGTGGGCATCTTCTTCGTTGGAACGATGATGGGCATGGGGCCACTGGATATGGt                                   <  2:413533/145‑1 (MQ=255)
                                                                                                                     ctcAAATCTGGGTTGTGAGCACGCTGCTTATCAGCATCGTGTTAATTGTACTGACCATCGTGAAGTTCAAATTCCACCCGTTTCTGGCGCTGTTGCTGGCCAGCTTCTTCGTGGGAACGATGATGGGCATGGGGCCACTGGATATGGTaa                                 >  2:183284/1‑150 (MQ=255)
                                                                                                                       cAAATCTGGGTTGTGAGCACGCTGCTTATCAGCATCGTGTTAATTGTACTGACCATCGTGAAGTTCAAATTCCACCCGTTTCTGGCGCTGTTGCTGGCCAGCTTCTTCGTGGGAACGATGATGGGCATGGGGCCACTGGATATGGTAAAt                               >  1:757634/1‑150 (MQ=255)
                                                                                                                       cAAATCTGGGTTGTGAGCACGCTGCTTATCAGCATCGTGTTAATTGTACTGACCATCGTGAAGTTCAAATTCCACCCGTTTCTGGCGCTGTTGCTGGCCAGCTTCTTCGTGGGAACGATGATGGGCATGGGGCCACTGGATATGGTAAAt                               <  2:1192502/150‑1 (MQ=255)
                                                                                                                          aTCTGGGTTGTGAGCACGCTGCTTATCAGCATCGTGTTAATTGTACTGACCATCGTGAAGTTCAAATTCCACCCGTTTCTGGCGCTGTTGCTGGCCAGCTTCTTCGTGGGAACGATGATGGGCATGGGGCCACTGGATATGGTAAATGCt                            <  1:299332/150‑1 (MQ=255)
                                                                                                                             tGGGTTGTGAGCACGCTGCTTATCAGCATCGTGTTAATTGTACTGACCATCGTGAAGTTCAAATTCCACCCGTTTCTGGCGCTGTTGCTGGCCAGCTTCTTCGTGGGAACGATGATGGGCATGGGGCCACTGGATATGGTAAATGCTAtt                         <  1:490772/150‑1 (MQ=255)
                                                                                                                              gggTTGTGAGCACGCTGCTTATCAGCATCGTGTTAATTGTACTGACCATCGTGAAGTTCAAATTCCACCCGTTTCTGGCGCTGTTGCTGGCCAGCTTCTTCGTGGGAACGATGATGGGCATGGGGCCACTGGATATGGTAAATGCTATTg                        <  2:568096/150‑1 (MQ=255)
                                                                                                                              gggTTGTGAGCACGCTGCTTATCAGCATCGTGTTAATTGTACTGACCATCGTGAAGTTCAAATTCCACCCGTTTCTGGCGCTGTTGCTGGCCAGCTTCTTCGTGGGAACGATGATGGGCATGGGGCCACTGGATATGGTAAATGCTATTg                        <  2:757634/150‑1 (MQ=255)
                                                                                                                                  aattaacACGATGCTGATAAGCAGCGTGTTAATTGTACTGACCATCGTGAAGTTCAAATTCCACCCGTTTCTGGCGCTGTTGCTGGCCAACTTCTTCGTGGGAACGATGATGGGCATGGGGCCACTGGATATGGTAAATGCTATTGAAAg                    >  2:1336897/7‑150 (MQ=255)
                                                                                                                                   attaacACGATGCTGATAAGCAGCGTGTTAATTGTACTGACCATCGTGAAGTTCAAATTCCACCCGTTTCTGGCGCTGTTGCTGGCCAGCTTCTTCGTGGGAACGATGATGGGCATGGGGCCACTGGATATGGTAAATGCTATTGAAAGt                   >  1:1198304/6‑150 (MQ=255)
                                                                                                                                    tgAGCACGCTGCTTATCAGCATCGTGTTAATTGTACTGACCATCGTGAAGTTCAAATTCCACCCGTTTCTGGCGCTGTTGCTGGCCAGCTTCTTCGTGGGAACGATGATGGGCATGGGGCCACTGGATATGGTAAATGCTATTGAAAGTg                  <  2:161609/150‑1 (MQ=255)
                                                                                                                                    tgAGCACGCTGCTTATCAGCATCGTGTTAATTGTACTGACCATCGTGAAGTTCAAATTCCACCCGTTTCTGGCGCTGTTGCTGGCCAGCTTCTTCGTGGGAACGATGATGGGCATGGGGCCACTGGATATGGTAAATGCTATTGAAAGTg                  <  2:1104288/150‑1 (MQ=255)
                                                                                                                                       acACGATGCTGATAAGCAGCGTGTTAATTGTACTGACCATCGTGAAGTTCAAATTCCACCCGTTTCTGGCGCTGTTGCTGGCCAACTTCTTCGTGGGAACGATGATGGGCATGGGGCCACTGGATATGGTAAATGCTATTGAAAGTGGaa               <  1:1336897/149‑1 (MQ=255)
                                                                                                                                          cGCTGCTTATCAGCATCGTGTTAATTGTACTGACCATCGTGAAGTTCAAATTCCACCCGTTTCTGGCGCTGTTGCTGGCCAGCTTCTTCGTGGGAACGATGATGGGCATGGGGCCACTGGATATGGTAAATGCTATTGAAAGTGGAATTg            <  2:1067771/150‑1 (MQ=255)
                                                                                                                                            atgctGATAAGCAGCGTGTTAATTGTACTGACCATCGTGAAGTTCAAATTCCACCCGTTTCTGGCGCTGTTGCTGGCCAGCTTCTTCGTGGGAACGATGATGGGCATGGGGCCACTGGATATGGTAAATGCTATTGAAAGTGGAATTGGc          <  2:1198304/149‑1 (MQ=255)
                                                                                                                                             tgctTATCAGCATCGTGTTAATTGTACTGACCATCGTGAAGTTCAAATTCCACCCGTTTCTGGCGCTGTTGCTGGCCAGCTTCTTCGTGGGAACGATGATGGGCATGGGGCCACTGGATATGGTAAATGCTATTGAAAGTGGAATTGGCg         <  1:266204/150‑1 (MQ=255)
                                                                                                                                             tgctTATCAGCATCGTGTTAATTGTACTGACCATCGTGAAGTTCAAATTCCACCCGTTTCTGGCGCTGTTGCTGGCCAGCTTCTTCGTGGGAACGATGATGGGCATGGGGCCACTGGATATGGTAAATGCTATTGAAAGTGGAATTGGCg         <  2:648954/150‑1 (MQ=255)
                                                                                                                                                 tATCAGCATCGTGTTAATTGTACTGACCATCGTGAAGTTCAAATTCCACCCGTTTCTGGCGCTGTTGCTGGCCAGCTTCTTCGTGGGAACGATGATGGGCATGGGGCCACTGGATATGGTAAATGCTATTGAAAGTGGAATTGGCGGAAc     >  1:1032909/1‑150 (MQ=255)
                                                                                                                                                    cAGCATCGTGTTAATTGTACTGACCATCGTGAAGTTCAAATTCCACCCGTTTCTGGCGCTGTTGCTGGCCAGCTTCTTCGTGGGAACGATGATGGGCATGGGGCCACTGGATATGGTAAATGCTATTGAAAGTGGAATTGGCGGAACGCt  >  2:180889/1‑150 (MQ=255)
                                                                                                                                                    cAGCATCGTGTTAATTGTACTGACCATCGTGAAGTTCAAATTCCACCCGTTTCTGGCGCTGTTGCTGGCCAGCTTCTTCGTGGGAACGATGATGGGCATGGGGCCACTGGATATGGTAAATGCTATTGAAAGTGGAATTGGCGGAACGCt  >  2:1417381/1‑150 (MQ=255)
                                                                                                                                                    cAGCATCGTGTTAATTGTACTGACCATCGTGAAGTTCAAATTCCACCCGTTTCTGGAGCTGTTGCTGGCCAGCTTCTTCGTGGGAACGATGATGGGCATGGGGCCACTGGATATGGTAAATGCTATTGAAAGTGGAATTGGCGGAACGCt  >  1:565710/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
TCACTCCATTGATGCAATTTACCTCATGTGAAAGGCAAATTTTATCGTTTGTCAGCCTGCGTTGTTTTTTTGTCCAATATCATCAGGTTAATCACAGGGGAAGGTGAGATATGCACTCTCAAATCTGGGTTGTGAGCACGCTGCTTATCAGCATCGTGTTAATTGTACTGACCATCGTGAAGTTCAAATTCCACCCGTTTCTGGCGCTGTTGCTGGCCAGCTTCTTCGTGGGAACGATGATGGGCATGGGGCCACTGGATATGGTAAATGCTATTGAAAGTGGAATTGGCGGAACGCT  >  NC_000913/2477737‑2478034

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