Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,632,586 G→C 15.9% intergenic (‑51/+18) 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,5860GC15.9% 528.1 / 84.9 197intergenic (‑51/+18)guaA/guaBGMP synthetase (glutamine aminotransferase)/IMP dehydrogenase
Reads supporting (aligned to +/- strand):  ref base G (83/81);  new base C (13/18);  total (96/101)
Fisher's exact test for biased strand distribution p-value = 4.36e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.83e-01

Only 100 of 200 total aligned reads displayed.

AACCGAAGTCCAGAATGAGGATGCGATGCTTATGAATGTTTTCCGTCATTGACGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACATGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAGTAATGGTCACGTCGTGAACGTGGCTTTCCTGAATGCCCGCACCGCTG  >  NC_000913/2632487‑2632685
                                                                                                   |                                                                                                   
  ccGAAGTCCAGAATGAGGATGCGATGCTTATGAATGTTTTCCGTCATTGACGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACATGaa                                                                                                   >  1:2117459/1‑100 (MQ=255)
    gAAGTCCAGAATGAGGATGCGATGCTTATGAATGTTTTCCGTCATTGACGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACATGAAGt                                                                                                 >  1:3205472/1‑100 (MQ=255)
     aaGTCCAGAATGAGGATGCGATGCTTATGAATGTTTTCCGTCATTGACGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACATGAAGTc                                                                                                <  1:1215744/100‑1 (MQ=255)
       gTCCAGAATGAGGATGCGATGCTTATGAATGTTTTCCGTCATTGACGCTTATTCCGAGGCAAGTGAAACAGATAATATGAATCGCCCGACATGAAGTCgg                                                                                              <  2:452961/100‑1 (MQ=255)
       gTCCAGAATGAGGATGCGATGCTTATGAATGTTTTCCGTCATTGACGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACATGAAGTCgg                                                                                              >  2:976921/1‑100 (MQ=255)
       gTCCAGAATGAGGATGCGATGCTTATGAATGTTTTCCGTCATTGACGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACATGAAGTCgg                                                                                              <  1:1915947/100‑1 (MQ=255)
        tCCAGAATGAGGATGCGATGCTTATGAATGTTTTCCGTCATTGACGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACATGAAGTCggg                                                                                             <  1:3313460/100‑1 (MQ=255)
         ccAGAATGAGGATGCGATGCTTATGAATGTTTTCCGTCATTGACGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACTTCATGTCGGGc                                                                                            <  2:2817884/100‑1 (MQ=255)
         ccAGAATGAGGATGCGATGCTTATGAATGTTTTCCGTCATTGACGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACTTCATGTCGGGc                                                                                            <  2:3011426/100‑1 (MQ=255)
         ccAGAATGAGGATGCGATGCTTATGAATGTTTTCCGTCATTGACGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACTTCATGTCGGGc                                                                                            <  2:3870962/100‑1 (MQ=255)
         ccAGAATGAGGATGCGATGCTTATGAATGTTTTCCGTCATTGAAGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACTTCATGTCGGGc                                                                                            <  2:186277/100‑1 (MQ=255)
         ccAGAATGAGGATGCGATGCTTATGAATGTGTTCCGTCATTGACGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACTTCATGTCGGGc                                                                                            <  1:826531/100‑1 (MQ=255)
         ccAGAATGAGGATGCGATGATTATGAATGTTTTCCGTCATTGACGCTTATTCCGAGGCAAGTGAAACAGATGATATAAATCGCCAGACATGAAGTCggg                                                                                             >  1:2105845/1‑99 (MQ=255)
          taGAATGAGGATGCGATGCTTATGAATGTTTTCCGTCATTGACGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACTTCATGTCGGGCg                                                                                           <  1:2854387/99‑1 (MQ=255)
          cAGAATGAGGATGCGATGCTTATGAATGTTTTCCGTCATTGACGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACATGAAGTCGGGCg                                                                                           <  1:4057619/100‑1 (MQ=255)
          cAGAATGAGGATGCGATACTTATGAATGTTTTCCGTCATTGACGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACATGAAGTCGGGCg                                                                                           >  1:3190694/1‑100 (MQ=255)
           aGAATGAGGATGCGATGCTTATGAATGTTTTCCGTCATTGACGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACATGAAGTCGGGCGa                                                                                          >  2:1932907/1‑100 (MQ=255)
            gAATGAGGATGCGATGCTTATGAATGTTTTCCGTCATTGACGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACATGAAGTCGGGCGaa                                                                                         >  2:3200865/1‑100 (MQ=255)
               tGAGGCTGCGTTGCTTATGAATGTTGTCCGTCATTGACGCTTATTCCGCCGCAAGTGAAACAGATAATATACATCGCCCGACATTAAGCCGGGCGAacag                                                                                      <  1:1375752/100‑1 (MQ=255)
               tGAGGATGGGATGCTTATGAATGTTTTACGTCATTGACGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACTTCATGTCGGGCGa                                                                                          <  1:2805910/96‑1 (MQ=255)
                gAGGATGCGATGCTTATGAATGTTTTCCGTCATTGACGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACATGAAGTCGGGCGAAgaga                                                                                     >  2:3653855/1‑100 (MQ=255)
                 aGGATGCGATGCTTATGAATGTTTTCCGTCATTGACGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACATGAAGTCGGGCGAAGAGaa                                                                                    >  2:2205339/1‑100 (MQ=255)
                      gcgatgcTTATGAATGTTTTCCGTCATTGACGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACATGAAGTCGGGCGAAGAGAATCAgg                                                                               >  1:2134260/1‑100 (MQ=255)
                       cgatgcTTATGAATGTTTTCCGTCATTGACGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACATGAAGTCGGGCGAAGAGAATCAGGa                                                                              >  2:3580348/1‑100 (MQ=255)
                       cgatgcTTATGAATGTTTTCCGTCATTGACGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACATGAAGTCGGGCGAAGAGAATCAGGa                                                                              <  2:606315/100‑1 (MQ=255)
                        gatgcTTATGAATGTTTTCCGTCATTGACGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACATGAAGTCGGGCGAAGAGAATCAGGAg                                                                             <  2:2910153/100‑1 (MQ=255)
                         atgcTTATGAATGTTTTCCGTCATTGACGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACATGAAGTCGGGCGAAGAGAATCAGGAGc                                                                            <  1:4146103/100‑1 (MQ=255)
                            cTTATGAATGTTTTCCGTCATTGACGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACATGAAGTCGGGCGAAGAGAATCAGGAGCCCa                                                                         <  2:3079567/100‑1 (MQ=255)
                               aTGAATGTTTTCCGTCATTGACGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACATGAAGTCGGGCGAAGAGAATCAGGAGCCCAGAc                                                                      <  1:4293966/100‑1 (MQ=255)
                                tGAATGTTTTCCGTCATTGACGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACATGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACg                                                                     >  2:3512550/1‑100 (MQ=255)
                                tGAATGTTTTCCGTCATTGACGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACATGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACg                                                                     <  2:2656271/100‑1 (MQ=255)
                                       tttCCGTCATTGACGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACATGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTc                                                              >  2:2850509/1‑100 (MQ=255)
                                         tCCGTCATTGACGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACATGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCgg                                                            <  2:1470950/100‑1 (MQ=255)
                                          ccGTCATTGACGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACATGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCggg                                                           <  1:3978411/100‑1 (MQ=255)
                                          ccGTCATTGACGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACATGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCggg                                                           >  2:4226766/1‑100 (MQ=255)
                                          ccGTCATTGACGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACATGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCggg                                                           >  1:644675/1‑100 (MQ=255)
                                           cGTCATTGACGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCGGACATGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCgggg                                                          <  1:4402243/100‑1 (MQ=255)
                                           cGTCATTGACGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACATGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCgggg                                                          >  2:3359138/1‑100 (MQ=255)
                                            gTCATTGACGTTTATCCCGAGGAAAGTAAAACAGATAATATAAAGCGCCCGACATGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCgggg                                                          <  1:4380509/99‑1 (MQ=255)
                                            gTCATTGACGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACATGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGa                                                         >  1:1784127/1‑100 (MQ=255)
                                              cATTGACGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACATGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGAct                                                       >  2:2112198/1‑100 (MQ=255)
                                              cATTGACGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACATGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGAct                                                       >  2:52309/1‑100 (MQ=255)
                                               aTTGACGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACATGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGActc                                                      <  1:1932900/100‑1 (MQ=255)
                                                  gACGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACATGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCttt                                                   <  1:376372/100‑1 (MQ=255)
                                                    cGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCTACATGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAg                                                 <  2:2433967/100‑1 (MQ=255)
                                                    cGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACATGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAg                                                 <  1:976919/100‑1 (MQ=255)
                                                     gCTTATTGCGTGGCAAGTGGAACAGATAATATAAATCGCCCGGCTTGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAg                                                 <  1:3592147/99‑1 (MQ=255)
                                                       ttATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACATGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAGTaa                                              >  1:1236830/1‑100 (MQ=255)
                                                        tATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACATGAAGTCGGGCGAAGAGAATGAGGAGACCAGACGGTAGTTCGGGGACTCTTTAGTAAt                                             >  1:2971832/1‑100 (MQ=255)
                                                          ttCCGAGGCAATGGAAACAGATAATATAAATCGCCCGACATAAAGCGGGGCGAAGATCTTCAGGAGCCCAGACTGTAGTTCGGGGACTCTTTAGTAATgg                                           <  1:739991/100‑1 (MQ=255)
                                                          ttCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACATGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAGTAATgg                                           <  2:4116451/100‑1 (MQ=255)
                                                            ccGAGGCAAGTGAAACAGATAATATAAATCGCCCGACATGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAGTAATGGTc                                         >  2:556881/1‑100 (MQ=255)
                                                             cGAGGCAAGTGAAACAGATAATATAAATCGCCCGACCTGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAGTAATGGTCa                                        >  1:2000822/1‑100 (MQ=255)
                                                              gAGGGAAGTGAAACAGATAATATAAATCGCCCGACATGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAGTAATGGTCa                                        <  1:3828668/99‑1 (MQ=255)
                                                              gAGGCAAGTGAAACAGATAATATAAATCGCCCGACATGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAGTAATGGTCAc                                       >  1:3164058/1‑100 (MQ=255)
                                                                ggCAAGTGAAACAGATAATATAAATCGCCCGACATGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAGTAATGGTCAcgt                                     >  1:3842388/1‑100 (MQ=255)
                                                                 gCAAGTGAAACAGATAATATAAATCGCCCGACATGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAGTAATGGTCAcg                                      >  2:1761069/1‑98 (MQ=255)
                                                                   aaGTGAAACAGATAATATAAATCGCCCGACATGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAGTAATGGTCAcgtcgt                                  <  2:141530/100‑1 (MQ=255)
                                                                       gAAACAGATAATATAAATCGCGCGACATGACGTCGGGCGAAGAGAATCAGGAGCCCAGCCGGTACTTCGGGGACTCTTTAGTAATGGTCACGTCGTGa                                <  1:3055040/98‑1 (MQ=255)
                                                                       gAAACAGATAATATAAATCGCCCGACATGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAGTAATGGTCACGTCGTGAAc                              >  2:2290462/1‑100 (MQ=255)
                                                                          aCAGATAATATAAATCGCCCGACATGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAGTAATGGTCACGTCGTGAACGTg                           <  1:419476/100‑1 (MQ=255)
                                                                          aCAGATAATATAAATCGCCCGACATGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAGTAATGGTCACGTCGTGAACGTg                           <  1:2727159/100‑1 (MQ=255)
                                                                             gATAATATAAATCGCCCGACATGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAGTAATGGTCACGTCGTGAACGTGGCt                        >  2:1994957/1‑100 (MQ=255)
                                                                               taataTAAATCGCCCGACATGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAGTAATGGTCACGTCGTGAACGTGGCttt                      >  2:1185391/1‑100 (MQ=255)
                                                                               taataTAAATCGCCCGACATGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAGTAATGGTCACGTCGTGAACGTGGCttt                      <  1:1599876/100‑1 (MQ=255)
                                                                                aataTAAATCGCCCGACATGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAGTAATGGTCACGTCGTGAACGTGGCTTTc                     <  1:2217834/100‑1 (MQ=255)
                                                                                aataTAAATCGCCCGACATGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAGTAATGGTCACGTCGTGAACGTGGCTTTc                     >  1:3512803/1‑100 (MQ=255)
                                                                                 atatAAATCGCCCGACATGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAGTAATGGTCACGTCGTGAACGTGGCTTTcc                    >  2:892702/1‑100 (MQ=255)
                                                                                  tataAATCGCCCGACATGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAGTAATGGTCACGTCGTGAACGTGGCTTTCCt                   >  2:3805509/1‑100 (MQ=255)
                                                                                   ataAATCGCCCGACATGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAGTAATGGTCACGTCGTGAACGTGGCTTTCCTg                  <  1:7272/100‑1 (MQ=255)
                                                                                    taAATCGCCCGACATGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAGTAATGGTCACGTCGTGAACGTGGCTTTCCTGa                 >  2:1146632/1‑100 (MQ=255)
                                                                                     tcttCGCCCGACTTCATGTCGGGCGAGGAGAATCAGGAGCCCAGACGGTAGTGCGGGGACTCTTTAGTAATGGTCACGTCGTGAACGTGGCTTTCCTGaa                >  2:2401201/4‑100 (MQ=38)
                                                                                     tcttCGCCCGACTTCATGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTCTAGTAATGGTCACGTCGGGAACGTGGCTTTCCTgca                >  1:4237221/4‑98 (MQ=38)
                                                                                     aaaTCGCCCGACATGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAGTAATGGTCACGTCGTGAACGTGGCTTTCCTGaa                <  2:3507563/100‑1 (MQ=255)
                                                                                     aaaTCGCCCGACATGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAGTAATGGTCACGTCGTGAACGTGGCTTTCCTGa                 <  2:2872685/99‑1 (MQ=255)
                                                                                      cttCGCCCGACTTCATGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAGTAATGGTCACGTCGTGAACGTGGCTTTCCTGAAt               >  2:1684736/3‑100 (MQ=255)
                                                                                      aaTCGCCCGACATGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAGTAATGGTCACGTCGTGAACGTGGCTTTCCTGAAt               >  1:1091809/1‑100 (MQ=255)
                                                                                      aaTCGCCCGACATGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAGTAATGGTCACGTCGTGAACGTGGCTTTCCTGAAt               >  1:4270828/1‑100 (MQ=255)
                                                                                       aTCGCCCGACATGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAGTAATGGTCACGTCGTGAACGTGGCTTTCCTGAATg              <  1:1994950/100‑1 (MQ=255)
                                                                                        tCGCCCGACTTCATGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAGTAATGGTCACGTCGTGAACGTGGCTTTCCTGAATGc             >  2:2149636/1‑100 (MQ=255)
                                                                                         cGCCCGACTTCATGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAGTAATGGTCACGTCGTGAACGTGGCTTTCCTGAATGcc            >  2:3696466/1‑100 (MQ=255)
                                                                                         cGCCCGACATGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAGTAATGGTCACGTCGTGAACGTGGCTTTCCTGAATGcc            <  2:4374208/100‑1 (MQ=255)
                                                                                         cGCCCGACATGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAGTAATGGTCACGTCGTGAACGTGGCTTTCCTGAATGcc            >  1:350358/1‑100 (MQ=255)
                                                                                          gCCCGACTTCATGTCGGGTGAAGAGAATCAGGAGCCCAGACGGTAGTTCGTGGACTCTTTAGTAATGGTCACGTCGTGAACGTGGCTTTCCTGAATg              >  1:1036324/1‑97 (MQ=255)
                                                                                          gCCCGACATGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAGTAATGGTCACGTCGTGAACGTGGCTTTCCTGAATGccc           <  2:1277867/100‑1 (MQ=255)
                                                                                           cccGACATGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAGTAATGGTCACGTCGTGAACGTGGCTTTCCTGAATGCCCg          <  2:1503116/100‑1 (MQ=255)
                                                                                           cccGACATGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAGTAATGGTCACGTCGTGAACGTGGCTTTCCTGAATGCCCg          <  2:1914806/100‑1 (MQ=255)
                                                                                            ccGACATGTTGTCGTGCGGATAGGGTCAGGAGCGCAGTCGATAGTTCGGGGACACTTCAGTAATGGGCACGTCGTGAACGTGGCATCCAGGAATGccac          <  1:1882962/99‑3 (MQ=255)
                                                                                            ccGACATGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAGTAATGGTCACGTCGTGAACGTGGCTTTCCTGAATGCCCGc         >  2:2323927/1‑100 (MQ=255)
                                                                                              gACATGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAGTAATGGTCACGTCGTGAACGTGGCTTTCCTGAATGCCCGCAc       <  2:2030507/100‑1 (MQ=255)
                                                                                              gACATGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAGTAATGGTCACGTCGTGAACGTGGCTTTCCTGAATGCCCGCAc       >  2:2492235/1‑100 (MQ=255)
                                                                                              gACATGAAGTCGGGCGAAGAGAACCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAGTAATGGTCACGTCGTGGACGTGGCTTTCCTGAATGCCCGCAc       >  2:1567311/1‑100 (MQ=255)
                                                                                               aCATGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAGTAATGGTCACGTCGTGAACGTGGCTTTCCTGAATGCCCGCAcc      <  2:1109377/100‑1 (MQ=255)
                                                                                               aCATGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAGTAATGGTCACGTCGTGAACGTGGCTTTCCTGAATGCCCGCAcc      >  2:3173490/1‑100 (MQ=255)
                                                                                                  tGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAGTAATGGTCACGTCGTGAACGTGGCTTTCCTGAATGCCCGCACCGCt   >  2:1429746/1‑100 (MQ=255)
                                                                                                  tGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAGTAATGGTCACGTCGTGAACGTGGCTTTCCTGAATGCCCGCACCGCt   <  1:892699/100‑1 (MQ=255)
                                                                                                  tGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAGTAATGGTCACGTCGTGAACGTGGCTTTCCTGAATGCCCGCACCGCt   >  1:3370146/1‑100 (MQ=255)
                                                                                                  tGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAGTAATGGTCACGTCGTGAACGTGGCTTTCCTGAATGCCCGCACCGCt   <  2:676948/100‑1 (MQ=255)
                                                                                                  tGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAGTAATGGTCACGTCGTGAACGTGGCTTTCCTGAATGCCCGCACCGCt   >  1:1133734/1‑100 (MQ=255)
                                                                                                   gAAGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAGTAATGGTCACGTCGTGAACGTGGCTTTCCTGAATGCCCTCACCGCTg  <  1:2436338/100‑1 (MQ=255)
                                                                                                   |                                                                                                   
AACCGAAGTCCAGAATGAGGATGCGATGCTTATGAATGTTTTCCGTCATTGACGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACATGAAGTCGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAGTAATGGTCACGTCGTGAACGTGGCTTTCCTGAATGCCCGCACCGCTG  >  NC_000913/2632487‑2632685

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: