Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,632,584 A→T 15.1% intergenic (‑49/+20) guaA ← / ← guaB GMP synthetase (glutamine aminotransferase)/IMP dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,632,5840AT15.1% 515.4 / 76.7 199intergenic (‑49/+20)guaA/guaBGMP synthetase (glutamine aminotransferase)/IMP dehydrogenase
Reads supporting (aligned to +/- strand):  ref base A (78/90);  new base T (12/18);  total (91/108)
Fisher's exact test for biased strand distribution p-value = 5.56e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.86e-01

Only 100 of 202 total aligned reads displayed.

AGAACCGAAGTCCAGAATGAGGATGCGATGCTTATGAATGTTTTCCGTCATTGACGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACATGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAGTAATGGTCACGTCGTGAACGTGGCTTTCCTGAATGCCCGCACCGC  >  NC_000913/2632485‑2632683
                                                                                                   |                                                                                                   
 gaACCAAAGTCCAGAATGAGGATGCGATGCTTATGAATGTTTTCCGTCATTGACGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACAt                                                                                                    >  1:3460841/1‑100 (MQ=255)
   aCCGAAGTCCAGAATGAGGATGCGATGCTTATGAATGTTTTCCGTCATTGACGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACATGa                                                                                                  <  2:3532817/100‑1 (MQ=255)
    ccGAAGTCCAGAATGAGGATGCGATGCTTATGAATGTTTTCCGTCATTGACGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACATGaa                                                                                                 >  1:2614922/1‑100 (MQ=255)
    ccGAAGTCCAGAATGAGGATGCGATGCTTATGAATGTTTTCCGTCATTGACGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACATGaa                                                                                                 <  2:2065574/100‑1 (MQ=255)
    ccGAAGTCCAGAATGAGGATGCGATGCTTATGAATGTTTTCCGTCATTGACGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACATGaa                                                                                                 >  1:1240306/1‑100 (MQ=255)
    ccGAAGTCCAGAATGAGGATGCGATGCTTATGAATGTGTTCCGTCATTGACGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACATGa                                                                                                  >  1:2983569/1‑99 (MQ=255)
          tCCAGAATGAGGATGCGATGCTTATGAATGTTTTCCGTCATTGACGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACATGAAGTCggg                                                                                           <  2:2714462/100‑1 (MQ=255)
          tCCAGAATGAGGATGCGATGCTTATGAATGTTTTCCGTCATTGACGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACATGAAGTCggc                                                                                           <  2:795307/100‑2 (MQ=255)
           ccAGAATGAGGATGCGATGCTTATGAATGTTTTCCGTCATTGACGCTTATTCCGAGGCAAGTGAAACAGATAATATTAATCGCCCGACATGAAGTCGGGc                                                                                          <  2:3727957/100‑1 (MQ=255)
           ccAGAATGAGGATGCGATGCTTATGAATGTTTTCCGTCATTGACGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACTTCATGTCGGGc                                                                                          <  1:4277518/100‑1 (MQ=255)
           ccAGAATGAGGATGCGATGCTTATGAATGTTTTCCGTCATTGACGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACTTCATGTCGGGc                                                                                          <  2:3993601/100‑1 (MQ=255)
           ccAGAATGAGGATGCGATGCTTATGAATGTTTTCCGTCATTGACGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACATGAAGTCGGGc                                                                                          >  2:4045200/1‑100 (MQ=255)
            cAGAATGAGGATGCGATGCTTATGAATGTTTTCCGTCATTGACGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACTTCATGTCGGGCg                                                                                         <  2:84627/100‑1 (MQ=255)
            cAGAATGAGGATGCGATGCTTATGAATGTTTTCCGTCATTGACGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACTTCATGTCGGGCg                                                                                         <  2:84293/100‑1 (MQ=255)
            cAGAATGAGGATGCGATGCTTATGAATGTTTTCCGTCATTGACGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACATGAAGTCGGGCg                                                                                         <  2:4334833/100‑1 (MQ=255)
             aGAATGAGGATGCGATGCTTATGAATGTTTTCCGTCATTGACGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACTTCATGTCGGGCGa                                                                                        <  1:1412075/100‑1 (MQ=255)
              gAATGAGGATGCGATGCTTATGAATGTTTTCCGTCATTGACGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACTTCATGTCGGGCGat                                                                                       <  2:2663475/100‑2 (MQ=255)
               aaTGAGGATGCGATGCTTATGAATGTTTTCCGTCATTGACGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACTTCATGTCGGGCGatt                                                                                      <  2:449080/100‑3 (MQ=255)
               aaTGAGGATGCGATGCTTATGAATGTTTTCCGTCATTGACGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACATGAAGTCGGGCGAag                                                                                      <  2:2339544/100‑1 (MQ=255)
                  gAGGATGCGATGCTTATGAATGTTTTCCGTCATTGACGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACATGAAGTCGGGCGAagag                                                                                    <  1:4192368/99‑1 (MQ=255)
                  gAGGATGCGATGCTTATGAATGTTTTCCGTCATTGACGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACATGAAGTCGGGCGAAgaga                                                                                   <  2:1972769/100‑1 (MQ=255)
                   aGGATGCGATGCTTATGAATGTTTTCCGTCATTGACGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACATGAAGTCGGGCGAAGAGaa                                                                                  >  2:1179012/1‑100 (MQ=255)
                   aGGATGCGATGCTTATGAATGTTTACCGTCATTGACGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACTTCATGTCGGGCGatttata                                                                                  <  1:1671265/100‑7 (MQ=37)
                     gatgcgatgcTTATGAATGTTTTCCGTCATTGACGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACTTCATGTCGGGCGatttatatt                                                                                <  2:3911404/100‑9 (MQ=37)
                     gatgcgatgcTTATGAATGTTTTCCGTCATTGACGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACATGAAGTCGGGCGAAGAGAATc                                                                                >  1:902718/1‑100 (MQ=255)
                      atgcgatgcTTATGAATGTTTTCCGTCATTGACGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACATGAAGTCGGGCGAAGAGAATCa                                                                               <  2:3814425/100‑1 (MQ=255)
                      atgcgatgcTTATGAAAGTTTGCCGTCATTGACGCGTATTCCGAGGCAAGAGAAAAAGATAATATAAATCGCCTGACATGAAGTCGGGCGAag                                                                                      <  1:4206658/93‑1 (MQ=255)
                       tgcgatgcTTATGAATGTTTTCCGTCATTGACGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACATGAAGTCGGGCGAAGAGAATCAg                                                                              <  1:1420170/100‑1 (MQ=255)
                       tgcgatgcTTATGAATGTTTTCCGTCATTGACGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACATGAAGTCGGGCGAAGAGAATCAg                                                                              <  2:97312/100‑1 (MQ=255)
                         cgatgcTTATGAATGTTTTCCGTCATTGACGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACATGAAGTCGGGCGAAGAGAATCAgg                                                                             <  1:1072038/99‑1 (MQ=255)
                         cgatgcTTATGAATGTTTTCCGTCATTGACGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACATGAAGTCGGGCGAAGAGAATCAGGa                                                                            >  1:3657834/1‑100 (MQ=255)
                             gcTTATGAATTTTTTCCGTCATTGACGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCTGACATGAAGTCGGGCGAAGAGAATCAGGAGccc                                                                        <  1:2890581/100‑1 (MQ=255)
                             gcTTATGAATGTTTTCCGTCATTGACGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACATGAAGTCGGGCGAAGAGAATCAGGAGccc                                                                        <  2:3061775/100‑1 (MQ=255)
                              cTTATGAATGTTTTCCGTCATTGACGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACATGAAGTCGGGCGAAGAGAATCAGGAGCCCa                                                                       <  2:4070907/100‑1 (MQ=255)
                                   naaTGTTTTCCGTCATTGACGCTTATTCCGAGGCAAGTGAAAAAGATTATAAAAATCGCCCGACATGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACgg                                                                  >  2:4365986/2‑100 (MQ=255)
                                    aaTGTTTTCCGTCATTGACGCTTATTCCGAGGCAAGTGAAACAGGTAATATAAATCGCCCGACATGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACgct                                                                 <  2:2976464/100‑3 (MQ=255)
                                     ctGTTTTCCGTCATTGACGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACATGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTa                                                                <  1:3994119/99‑1 (MQ=255)
                                     aTGTTTTCCGTCATTGACGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACATGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTa                                                                <  1:1179007/100‑1 (MQ=255)
                                      tGTTTTCCGTCATTGACGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACATGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAg                                                               >  2:2040639/1‑100 (MQ=255)
                                      tGTTTTCCGTCATTGACGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACATGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAg                                                               >  2:576918/1‑100 (MQ=255)
                                            ccGTCATTGACGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACATGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCggg                                                         <  2:2924387/100‑1 (MQ=255)
                                            ccGTCATTGACGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACATGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCggg                                                         >  2:22014/1‑100 (MQ=255)
                                            ccGTCATTGACGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACATGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCggg                                                         >  1:482716/1‑100 (MQ=255)
                                            ccGTCATTGACGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACATGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCggg                                                         >  1:551884/1‑100 (MQ=255)
                                             cGTCATTGACGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACATGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCgggg                                                        <  2:272179/100‑1 (MQ=255)
                                             cGTCATTGACGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACATGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCgggg                                                        >  1:3382240/1‑100 (MQ=255)
                                             cGTCATTGACGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACATGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCgggg                                                        >  2:536626/1‑100 (MQ=255)
                                              gTCATTGACGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACATGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGa                                                       >  1:2915284/1‑100 (MQ=255)
                                               tCATTGACGCTTATTCCGAGGCAAGTGAAACGGATAATATAAATCGCCCGACATGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGAc                                                      <  2:2643044/100‑1 (MQ=255)
                                               tCATTGACGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACATGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGAc                                                      <  2:1557579/100‑1 (MQ=255)
                                               tCATTGAAGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACATGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGAc                                                      <  2:688214/100‑1 (MQ=255)
                                                    gACGCTTATTCCGAGGCAAGTGCAACAGATACTATAACTCGCCCGACATGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCttt                                                 <  1:2968869/100‑1 (MQ=255)
                                                    gACGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACATGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGActctat                                                 <  2:2988377/100‑3 (MQ=255)
                                                    gACGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACATGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCttt                                                 <  1:3130025/100‑1 (MQ=255)
                                                      cGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACATGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAg                                               >  2:3684897/1‑100 (MQ=255)
                                                      cGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACATGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAg                                               <  2:3098072/100‑1 (MQ=255)
                                                      cGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACATGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAg                                               >  2:3166989/1‑100 (MQ=255)
                                                      cGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACATGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAg                                               <  2:2625275/100‑1 (MQ=255)
                                                      cGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACATGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAg                                               <  2:1278745/100‑1 (MQ=255)
                                                        cTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACATGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAGTa                                             >  1:4250681/1‑100 (MQ=255)
                                                          tATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACATGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAGTAAt                                           >  2:129289/1‑100 (MQ=255)
                                                          tATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACATGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAGTAAt                                           >  2:4100627/1‑100 (MQ=255)
                                                             tCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACATGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAGTAATGGt                                        >  1:2213745/1‑100 (MQ=255)
                                                              ccGAGGCAAGTGAAACAGATAATATAAATCGCCCGACATGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAGTAATGGTc                                       >  1:757051/1‑100 (MQ=255)
                                                              ccGAGGCAAGTGAAACAGATAATATAAATCGCCCGACATGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAGTAATGGTc                                       >  2:2853522/1‑100 (MQ=255)
                                                                gAGGCAAGTGAAACAGATAATATAAATCGCCCGACATGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAGTAATGGTCAc                                     >  1:297411/1‑100 (MQ=255)
                                                                  ggCAAGTGAAACAGATAATATAAATCGCCCGACATGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAGTAATGGTCAcgt                                   >  1:133713/1‑100 (MQ=255)
                                                                  ggCAAGTGAAACAGATAATATAAATCGCCCGACATGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAGTAATGGTCAcgt                                   <  2:2875660/100‑1 (MQ=255)
                                                                   gCAAGTGAAACAGATAATATAAATCGCCCGACATGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAGTAATGGTCAcgtc                                  <  2:775111/100‑1 (MQ=255)
                                                                    cAAGTGAAACAGATAATATAAATCGCCCGACATGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAGTAATGGTCAcgtcg                                 >  1:223661/1‑100 (MQ=255)
                                                                     aaGTGAAACAGATAATATAAATCGCCCGACATGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAGTAATGGTCAcgtcgt                                <  1:971686/100‑1 (MQ=255)
                                                                      aGTGAAACAGATAATATAAATCGCCCGACATGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAGTAATGGTCACGTCGTg                               >  2:3413080/1‑100 (MQ=255)
                                                                      aGTGAAACAGATAATATAAATCGCCCGACATGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAGTAATGGTCACGTCGTg                               <  2:202168/100‑1 (MQ=255)
                                                                       gTGAAACAGATAATATAATTCGCCCTACATGAAGTCGGGCGAATAGAATCAGGAGCCCAGACTGTAGTTCGGGCACTCTTTAGTAATGGTCAcgtcgt                                <  1:1293095/98‑1 (MQ=255)
                                                                           aaCAGATAATATAAATCGCCCGACATGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAGTAATGGTCACGTCGTGAACGt                          >  2:759823/1‑100 (MQ=255)
                                                                            aCAGATAATATAAATCGCCCGACATGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAGTAATGGTCACGTCGTGAACGTg                         <  2:1012914/100‑1 (MQ=255)
                                                                             cAGATAATATAAATCGCCCGACATGAAGGCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAGTAATGGTCACGTCGTGAACGTgg                        <  1:2025686/100‑1 (MQ=255)
                                                                              aGATAATATAAATCGCCCGACATGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAGTAATGGTCACGTCGTGAACGTGGc                       >  1:2092797/1‑100 (MQ=255)
                                                                              aGATAATATAAATCGCCCGACATGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAGTAATGGTCACGTCGTGAACGTGGc                       <  1:3935249/100‑1 (MQ=255)
                                                                               gATAATATAAATCGCCCGACATGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAGTAATGGTCACGTCGTGAACGTGGCt                      >  2:3258968/1‑100 (MQ=255)
                                                                                 taataTAAATCGCCCGACATGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAGTAATGGTCACGTCGTGAACGTGGCttt                    >  1:2326143/1‑100 (MQ=255)
                                                                                  aataTAAAGCGCCCGACATGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAGTAATGGTCACGTCGTGAACGTGGCTTTc                   <  1:4355757/100‑1 (MQ=255)
                                                                                   atatAAATCGCCCGACATGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAGTAATGGTCACGTCGTGAACATGGCTTTcc                  >  1:2803575/1‑100 (MQ=255)
                                                                                    tataAATCGCCCGACATGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAGTAATGGTCACGTCGTGAACGTGGCTTTCCt                 >  1:2582988/1‑100 (MQ=255)
                                                                                     tctcttCGCCCGACTTCATGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAGTAATGGTCACGTCGAAAACGTGGCTTTcat                 >  1:685701/6‑97 (MQ=37)
                                                                                     atgaaTCGCCCGACATGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGTCTCTTTAGTAATGGTCACGGCGTGAACGTGGCttt                    <  1:2180510/93‑1 (MQ=255)
                                                                                       tcttCGCCCGACTTCATGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAGTAATGGTCACGTCGTGAACGTGGCTTTGCTGa               >  1:2383732/4‑99 (MQ=38)
                                                                                       tcttCGCCCGACTTCATGGCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAGTAATGGTCACGTCGTGAACGTGGCTTTCCTGaa              <  1:205380/97‑1 (MQ=255)
                                                                                       aaaTCGCCCGACATGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAGTAATGGTCACGTCGTGAACGTGGCTTTCCTGaa              <  1:1953807/100‑1 (MQ=255)
                                                                                        cttCGCCCGACTTCATGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAGTAATGGTCACGTCGTGAACGTGGCTTTCCTGAAt             >  2:2505715/3‑100 (MQ=255)
                                                                                        aaTCGCCCGACATGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAGTAATGGTCACGTCGTGAACGTGGCTTTCCTGAAt             <  2:1735208/100‑1 (MQ=255)
                                                                                          tCGCCCGACTTCATGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAGTAATGGTCACGTCGTGAACGTGGCTTTCCTGAATGc           >  2:3548956/1‑100 (MQ=255)
                                                                                           cGCCCGACATGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAGTAATGGTCACGTCGTGAACGTGGCTTTCCTGAATGcc          >  2:2550927/1‑100 (MQ=255)
                                                                                             cccGACATGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAGTAATGGTCACGTCGTGAACGTGGCTTTCCTGAATGcc          >  1:1802858/1‑98 (MQ=255)
                                                                                             cccGACATGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAGTAATGGTCACGTCGTGAACGTGGCTTTCCTGAATGCCCg        <  1:4028167/100‑1 (MQ=255)
                                                                                             cccGACATGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAGTAATGGTCACGTCGTGAACGTGGCTTTCCTGAATGCCCg        <  1:439216/100‑1 (MQ=255)
                                                                                             cccGACATGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAGTAATGGTCACGTCGTGAACGTGGCTTTCCTGAATGCCCg        <  2:2935412/100‑1 (MQ=255)
                                                                                              ccGACATGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAGTAATGGTCACGTCGTGAACGTGGCTTTCCTGAATGCCCGc       >  1:3634332/1‑100 (MQ=255)
                                                                                              ccGACATGAAGTCCGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAGTAATGGTCACGTCGTGAACGTGGCTTTCCTGAATGCCCg        <  1:2999207/99‑1 (MQ=255)
                                                                                               cGACATGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAGTAATGGTCACGTCGTGTACGTGGCTTTCCTGAATGCCCGc       <  1:4146873/99‑1 (MQ=255)
                                                                                                   |                                                                                                   
AGAACCGAAGTCCAGAATGAGGATGCGATGCTTATGAATGTTTTCCGTCATTGACGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACATGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAGTAATGGTCACGTCGTGAACGTGGCTTTCCTGAATGCCCGCACCGC  >  NC_000913/2632485‑2632683

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: