Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,632,588 A→T 5.7% intergenic (‑53/+16) 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,5880AT 5.7% 312.0 / 6.0 116intergenic (‑53/+16)guaA/guaBGMP synthetase (glutamine aminotransferase)/IMP dehydrogenase
Reads supporting (aligned to +/- strand):  ref base A (50/57);  new base T (3/4);  total (53/63)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.77e-01

Only 100 of 120 total aligned reads displayed.

AAACACCCAGCTCACGCACGCGGCGCGCAACCAGTTGAGTGTACTGAGAACCGAAGTCCAGAATGAGGATGCGATGCTTATGAATGTTTTCCGTCATTGACGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACATGAAGTC‑GGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAGTAATGGTCACGTCGTGAACGTGGCTTTCCTGAATGCCCGCACCGCTGATACGTACAAACTCCGCTTTAGTACGCAGTTCGTCGATAGTACCACAGCCGG  >  NC_000913/2632439‑2632737
                                                                                                                                                     |                                                                                                                                                      
aaaCACCCAGCTCACGCACGCGGCGCGCAACCAGTTGAGTGTACTGAGAACCGAAGTCCAGAATGAGGATGCGATGCTTATGAATGTTTTCCGTCATTGCCGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACATGaa                                                                                                                                                        >  1:1038146/1‑150 (MQ=255)
    acCCAGCTCACGCACGCGGCGCGCAACCAGTTGAGTGTACTGAGAACCGAAGTCCAGAATGAGGATGCGATGCTTATGAATGTTTTCCGTCATTGACGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACATGAAGTC‑g                                                                                                                                                   >  2:1310458/1‑150 (MQ=255)
     cccAGCTCACGCACGCGGCGCGCAACCAGTTGAGTGTACTGAGAACCGAAGTCCAGAATGAGGATGCGATGCTTATGAATGTTTTCCGTCATTGACGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACATGAAGTC‑gg                                                                                                                                                  >  2:603916/1‑150 (MQ=255)
     cccAGCTCACGCACGCGGCGCGCAACCAGTTGAGTGTACTGAGAACCGAAGTCCAGAATGAGGATGCGATGCTTATGAATGTTTTCCGTCATTGACGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACATGAAGTC‑gg                                                                                                                                                  >  2:1599241/1‑150 (MQ=255)
        aGCTCACGCACGCGGCGCGCAACCAGTTGAGTGTACTGAGAACCGAAGTCCAGAATGAGGATGCGATGCTTATGAATGTTTTCCGTCATTGACGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACATGAAGTC‑GGGCg                                                                                                                                               <  1:1053797/150‑1 (MQ=255)
        aGCTCACGCACGCGGCGCGCAACCAGTTGAGTGTACTGAGAACCGAAGTCCAGAATGAGGATGCGATGCTTATGAATGTTTTCCGTCATTGACGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACATGAAGTC‑GGGCg                                                                                                                                               <  1:320921/150‑1 (MQ=255)
        aGCTCACGCACGCGGCGCGCAACCAGTTGAGTGTACTGAGAACCGAAGTCCAGAATGAGGATGCGATGCTTATGAATGTTCTCCGTCATTGACGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACATGAAGTC‑GGGCg                                                                                                                                               <  1:231735/150‑1 (MQ=255)
          cTCACGCACGCGGCGCGCAACCAGTTGAGTGTACTGAGAACCGAAGTCCAGAATGAGGATGCGATGCTTATGAATGTTTTCCGTCATTGACGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACATGAAGTC‑GGGCGaa                                                                                                                                             >  2:137103/1‑150 (MQ=255)
           tCACGCACGCGGCGCGCAACCAGTTGAGTGTACTGAGAACCGAAGTCCAGAATGAGGATGCGATGCTTATGAATGTTTTCCGTCATTGACGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACTTCATGTC‑GGGCGatt                                                                                                                                            <  1:1633064/150‑3 (MQ=255)
           tCACGCACGCGGCGCGCAACCAGTTGAGTGTACTGAGAACCGAAGTCCAGAATGAGGATGCGATGCTTATGAATGTTTTCCGTCATTGACGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACTTCATGTC‑GGGCGatt                                                                                                                                            <  1:540803/150‑3 (MQ=255)
            cacgcaccCGGCGCGCAACCAGGTGAGTGTACTGAGAACCGAAGTCCAGAATGAGGATGCGATGCTTATGAATGTTTTCCGTCATTGACGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACATGAAGTC‑GGGCGAaga                                                                                                                                           <  1:1064337/150‑1 (MQ=255)
               gcacgcGGCGCGCAACCAGTTGAGTGTACTGAGAACCGAAGTCCAGAATGAGGATGCGATGCTTATGAATGTTTTCCGTCATTGACGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACATGAAGTC‑GGGCGAAGAGaa                                                                                                                                        <  1:170436/150‑1 (MQ=255)
                cacgcGGCGCGCAACCAGTTGAGTGTACTGAGAACCGAAGTCCAGAATGAGGATGCGATGCTTATGAATGTTTTCCGTCATTGACGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACATGAAGTC‑GGGCGAAGAGAAt                                                                                                                                       >  2:1085296/1‑150 (MQ=255)
                 acgcGGCGCGCAACCATTTGAGTGTACTGAGAACCGAAGTCCAGAATGAGGATGCGATGCTTATGAATGTTTTCCGTCATTGACGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACATGAAGTC‑GGGCGAAGAGAATc                                                                                                                                      <  2:712364/150‑1 (MQ=255)
                  cgcgGCGCGCAACCAGTTGAGTGTACTGAGAACCGAAGTCCAGAATGAGGATGCGATGCTTATGAATGTTTTCCGTCATTGACGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACATGAAGTC‑GGGCGAAGAGAATCa                                                                                                                                     >  1:32188/1‑150 (MQ=255)
                    cggcgCGCAACCAGTTGAGTGTACTGAGAACCGAAGTCCAGAATGAGGATGCGATGCTTATGAATGTTTTCCGTCATTGACGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACATGAAGTC‑GGGCGAAGAGAATCAgg                                                                                                                                   <  1:862496/150‑1 (MQ=255)
                         cgcAACCAGTTGAGTGTACTGAGAACCGAAGTCCAGAATGAGGATGCGATGCTTATGAATGTTTTCCGTCATTGACGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACATGAAGTC‑GGGCGAAGAGAATCAGGAGccc                                                                                                                              <  2:1074952/150‑1 (MQ=255)
                             aCCAGTTGAGTGTACTGAGAACCGAAGTCCAGAATGAGGATGCGATGCTTATGAATGTTTTCCGTCATTGACGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACATGAAGTC‑GGGCGAAGAGAATCAGGAGCCCAGAc                                                                                                                          >  1:723607/1‑150 (MQ=255)
                                    gAGTGTACTGAGAACCGAAGTCCAGAATGAGGATGCGATGCTTATGAATGTTTTCCGTCATTGCCGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACATGAAGTC‑GGGCGAAGAGAATCAGGAGCCCAGCCGGTAGtt                                                                                                                   >  1:392810/1‑150 (MQ=255)
                                         tACTGAGAACCGAAGTCCAGAATGAGGATGCGATGCTTATGAATGTTTTCCGTCATTGACGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACATGAAGTC‑GGGCGAAGAGAATCAGGAGCCCAGACGGTAg                                                                                                                     >  1:668477/1‑143 (MQ=255)
                                         tACTGAGAACCGAAGTCCAGAATGAGGATGCGATGCTTATGAATGTTTTCCGTCATTGACGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACATGAAGTC‑GGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCgggg                                                                                                              <  2:1685259/150‑1 (MQ=255)
                                            tGAGAACCGAAGTCCATAATGAGCATGCGAGGCATATGAATGTTTTCCGTCATTTACGCTTATTCCGAGGCAAGTTAAACAGATAATATAAATCGCCCGACATGACCTC‑GGGCGAAGAGAATAAGTAGCCCAGACGGTAGTTCTGGGAct                                                                                                           <  2:183326/150‑1 (MQ=255)
                                            tGAGAACCGAAGTCCAGAATGAGGATGCGATGCTTATGAATGTTTTCCGTCATTGACGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACATGAAGTC‑GGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGAct                                                                                                           >  1:666540/1‑150 (MQ=255)
                                             gagaACCGAAGTCCAGAATGAGGATGCGATGCTTATGAATGTTTTCCGTCATTGACGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACATGAAGTC‑GGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGActc                                                                                                          <  2:18553/150‑1 (MQ=255)
                                               gaACCGAAGTCCAGAATGAGGATGCGATGCTTATGAATGTTTTCCGTCATTGACGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACATGAAGTC‑GGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCtt                                                                                                        >  2:148495/1‑150 (MQ=255)
                                                aaCCGAAGTCCAGAATGAGGATGCGATGCTTATGAATGTTTTCCGTCATTGACGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACATGAAGTC‑GGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCttt                                                                                                       <  2:511171/150‑1 (MQ=255)
                                                  ccGAAGTCCAGAATGAGGATGCGATGCTTATGAATGTTTTCCGTCATTGACGCTTATTCCGAGGCAATTGAAACAGATAATATAAATCGCCCGACATGAAGTC‑GGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAg                                                                                                     <  2:1077921/150‑1 (MQ=255)
                                                  ccGAAGTCCAGAATGAGGATGCGATGCTTATGAATGTTTTCCGTCATTGACGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACATGGAGTC‑GGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAg                                                                                                     <  1:396130/150‑1 (MQ=255)
                                                  ccGAAGTCCAGAATGAGGATGCGATGCTTATGAATGTTTTCCGTCATTGACGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACATGAAGTC‑GGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGAGGACTCTTTAg                                                                                                     <  2:156217/150‑1 (MQ=255)
                                                    gAAGTCCAGAATGAGGATGCGATGCTTATGAATGTTTTCCGTCATTGACGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACATGAAGTC‑GGGCGAAGAGAATCAGGAGCCCAGCCGGTAGTTCGGGGACTCTTTAGTa                                                                                                   >  1:92011/1‑150 (MQ=255)
                                                     aaGTCCAGAATGAGGATGCGATGCTTATGAATGTTTTCCGTCATTGACGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACATGAAGTC‑GGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAGTaa                                                                                                  <  1:1326546/150‑1 (MQ=255)
                                                       gTCCAGAATGAGGATGCGATGCTTATGAATGTTTTCCGTCATTGACGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACATGAAGTC‑GGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAGTAATg                                                                                                >  2:74356/1‑150 (MQ=255)
                                                       gTCCAGAATGAGGATGCGATGCTTATGAATGTTTTCCGTCATTGACGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACATGAAGTC‑GGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAGTAATg                                                                                                <  2:1600919/150‑1 (MQ=255)
                                                       gTCCAGAATGAGGATGCGATGCTTATGAATGTTTTCCGTCATTGACGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACATGAAGTC‑GGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAGTAATg                                                                                                <  2:1655503/150‑1 (MQ=255)
                                                        tCCAGAATGAGGATGCGATGCTTATGAATGTTTTCCGTCATTGACGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACATGAAGTC‑GGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAGTAATgg                                                                                               <  1:74356/150‑1 (MQ=255)
                                                         ccAGAATGAGGATGCGATGCTTATGAATGTTTTCCGTCATTGACGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACATGAAGTC‑GGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAGTAATGGt                                                                                              >  2:795667/1‑150 (MQ=255)
                                                          cAGAATGAGGATGCGATGCTTATGAATGTTTTCCGTCATTGACGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACATGAAGTC‑GGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAGTAATGGTc                                                                                             >  2:653272/1‑150 (MQ=255)
                                                           aGAATGAGGATGCGATGCTTATGAATGTTTTCCGTCATTGACGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACATGAAGTC‑GGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAGTAATGGTCa                                                                                            >  2:1621719/1‑150 (MQ=255)
                                                             aaTGAGGATGCGATGCTTATGAATGTTTTCCGTCATTGACGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACATGAAGTC‑GGGTGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAGTAATGGTCAcg                                                                                          <  2:896743/150‑1 (MQ=255)
                                                                gaggagGCGATGCTTATGAATGTTTTCCGTCATTGACGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACATGAAGTC‑GGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAGTAATGGTCAcgtcg                                                                                       <  1:1460185/150‑1 (MQ=255)
                                                                 aGGATGCGATGCATATGAATGTTTTCCGTCATTGACGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACATGAAGTC‑GGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAGTAATGGTCAcgtcgt                                                                                      <  2:1412782/150‑1 (MQ=255)
                                                                  ggATGCGATGCTTATGAATGTTTTCCGTCATTGACGCTTATTCCGAGGCAAGTGAAACAGATAAAATAAATCGCCCGACATGAAGTC‑GGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGTGACTCTTTAGTAATGGTCACGTCGTg                                                                                     <  2:469838/150‑1 (MQ=255)
                                                                   gatgcgatgcTTATGAATGTTTTCCGTCATTGACGCTTATTCCGAGGCATGTGAAACAGATAATATAAATCGCCCGACATGAAGTC‑GGGCGAAGAGATTCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAGTAATGGTCACGTCGTGa                                                                                    <  2:583968/150‑1 (MQ=255)
                                                                   gatgcgatgcTTATGAATGTTTTCCGACATTGACGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACATGCAGTC‑GGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAGTAATGGTCACGTCGTGa                                                                                    <  2:835224/150‑1 (MQ=255)
                                                                       cgatgcTTATGAATGTTTTCCGTCATTGACGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACATGAAGTC‑GGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTGTTTAGTAATGGTCACGTCGTGAACGt                                                                                >  2:1536967/1‑150 (MQ=255)
                                                                        gatgcTTATGAATGTTTTCCGTCATTGACGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACATGAAGTC‑GGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAGTAATGGTCACGTCGTGAACGTg                                                                               >  1:1329849/1‑150 (MQ=255)
                                                                        gatgcTTATGAATGTTTTCCGTCATTGACGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACATGAAGTC‑GGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAGTAATGGTCACGTCGTGAACGTg                                                                               <  1:1621719/150‑1 (MQ=255)
                                                                         atgcTTATGAATGTTTTCCGTCATTGACGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACATGAAGTC‑GGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAGTAATGGTCACGTCGTGAACGTgg                                                                              <  2:417941/150‑1 (MQ=255)
                                                                           gcTTATGAATGTTTTCCGTCATTGACGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACATGAAGTC‑GGGCGAAGAGAATCAGGAGCACAGACGGTAGTTCTGGGACTCTTTAGTAATGGTCACGTCGTGAACGTGGCt                                                                            >  1:148295/1‑150 (MQ=255)
                                                                           gcTTATGAATGTTTTCCGTCATTGACGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACATGAAGTC‑GGGCGAAGAGAATCAGGAGACCAGACGGTAGTTCGGGGACTCTTTAGTAATGGTCACGTCGTGATCTTGGCt                                                                            >  1:1653795/1‑150 (MQ=255)
                                                                              tATGAATGTTTTCCGTCATTGACGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACATGAAGTC‑GGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAGTAATGGTCACGTCGTGAACGTGGCTTTc                                                                         >  2:747026/1‑150 (MQ=255)
                                                                                 taaTGTTTTCCGTCATTGACGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACATGAAGTC‑GGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAGTAATGGTCACGTCGTGAACGTGGCTTTCCTg                                                                      <  2:740750/149‑1 (MQ=255)
                                                                                 taaTGTTTTCCGTCATTGACGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACATGAAGTC‑GGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAGTAATGGTCACGTCGTGAACGTGGCTTTCCTg                                                                      <  1:499527/149‑1 (MQ=255)
                                                                                   aTGTTTTCCGTCATTGACGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACATGAAGTC‑GGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAGTAATGGTCACGTCGTGAACGTGGCTTTCCTGaa                                                                    >  1:1042481/1‑150 (MQ=255)
                                                                                    ggTTTTCCGGCATTGACGCTTATTCCGAGGCCAGTGACACAGATAATATAATTCGCCCGACCTGTAGTC‑GGGCGAAGAGAATCAGGAGCCCAGTCGGTAGTTCTTGTACTCTTTAGTAATGGTCACGTCGTGAACGTGGCTTTCCTGAAt                                                                   <  1:254767/149‑1 (MQ=255)
                                                                                                gtgtcGCTTATTCCGACGCAAGTGAAACAGATAATATAAATCGCCCGACATGAAGTC‑GGGCGAAGAGAATCAGGAGCCCAGACGGTAGTGCGGGGACTCTCTAGTAATGGTCACGTCGTGAACGTGGCTTTCCTGAATGCCCGCACCGcg                                                       <  1:137178/146‑2 (MQ=255)
                                                                                                  gACGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACATGAAGTC‑TGGCGAAGAGAATCAGGAGCCAAGACGGTAGTTCGGTGACTCTTTAGTAATTGTCACGTCGTGAACGTGGCTTTCCTGAATGCCCGCACCGCTGa                                                     >  1:1043394/1‑150 (MQ=255)
                                                                                                  gACGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACATGAAGTC‑GGGCGAAGAGAAACAGGAGCCCAGCCGGTAGTTCGGGGACTCTTTAGTAATGGTCACGTCGTGAACGTGGCTTGACTGAATGACCGCACCGCTGa                                                     >  1:857536/1‑150 (MQ=255)
                                                                                                   aCGCTTATTCCGAGGCAAGTGAAACAGATAATTTAAATCGCCCGACATGCAGTC‑GGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAGTAATGGTCACGTCGTGAACGTGGCTTTCCTGAATGCCCGCACCGCTGAt                                                    <  2:1516255/150‑1 (MQ=255)
                                                                                                   aCGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACATGAAGTC‑GGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAGTAATGGTCACGTCGTGAACGTGGCTTTCCTGAATGCCCGCACCGCTGAt                                                    <  1:747026/150‑1 (MQ=255)
                                                                                                   aCGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACATGAAGTC‑GGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACGCTTTAGTAATGGTCACGTCGGGAACGTGGCTTTCCTGAATGCCCGCACCGCTGAt                                                    <  2:116193/150‑1 (MQ=255)
                                                                                                     gCTTAATCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACATGAAGTC‑GGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAGTAATGGTCACGTCGTGAACGTGGCTTTCCTGAATGCCCGCACCGCTGATAc                                                  <  1:1136945/150‑1 (MQ=255)
                                                                                                        tATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACATGAAGTC‑GGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAGTAATGGTCACGTCGTGAACGTGGCTTTCCTGAATGCCCGCACCGCTGATACGTa                                               <  2:1701688/150‑1 (MQ=255)
                                                                                                          ttCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACATGAAGTC‑GGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAGTAATGGTCACGTCGTGAACGTGGCTTTCCTGAATGCCCGCACCGCTGATACGTACa                                             >  2:297602/1‑150 (MQ=255)
                                                                                                           tCCGATGCAAGTGAATCAGATAATATAAATCGCCCGACATGAAGTC‑GGGCGAAGAGAATCAGGAGCCCAGACGGTAGTACGGGGACTCTTTAGTAATGGTCACGTCGAGAACGAGGCTTTCCTGAATGCCCGCACCTCTCATAAGTaaaa                                            >  1:59746/1‑150 (MQ=255)
                                                                                                           tCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACATGAAGTC‑GGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAGTAATGGTCACGTCGTGAACGTGGCTTTCCTGAATGCCCGCACCGCTGATACGTACaa                                            >  1:810762/1‑150 (MQ=255)
                                                                                                           tCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACATGAAGTC‑GGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAGTAATGGTCACGTCGTGAACGTGGCTTTCCTGAATGCCCGCACCGCTGATACGTACaa                                            <  1:603916/150‑1 (MQ=255)
                                                                                                              gAGGCAAGTGAAACAGATAATATAAATCGCCCGACATGAAGTC‑GGGCGAAGAGAATCAGGAGCCAAGACGGTAGTTCGGGGACTCTTTAGTAATGGTAACGTCGTGAACGAGGCTATCCTTAATGCCCGCACCGATGATACGTACAAACt                                         >  1:574423/1‑150 (MQ=255)
                                                                                                               aGGCAAGTTACACAGATAATATAAATGGCCCGACATGAAGTC‑GGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAGTAATGGTCACGTCGTGAACGTGGCTTTCCTGAATGCCCGCACCGCTGATACGTACAAACTc                                        <  1:137103/150‑1 (MQ=255)
                                                                                                               aGGCAAGTGAAACAGATAATATAAATCGCCCTACATGATGTC‑GGGCGAAGAGCAGCAGGAGCCCAGACTGTAGTTCGGGGACGCTTTAGTAATGGTCACGTCGTGAACGTGGCTTTCCTGAATGTCCGCACCGCTGATTCGTACTAACTc                                        <  2:291270/150‑1 (MQ=255)
                                                                                                                gtcAAGTGAAACAGATAATATAAATCGCCCGACATGAAGTC‑GGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAGTAATGGTCACGTCGTTAACGTGGCTTTCCTGAATGCCCGCACCGCTGATACGTACAAACTcc                                       <  2:428328/148‑1 (MQ=255)
                                                                                                                ggCAAGTGAAACAGATAATATAAATCGCCCGACATGAAGTC‑GGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAGTAATGGTCACGTCGTGAACGTGGCTTTCCTGAATGCCCGCACCGCTGATACGTACAAACTcc                                       <  2:1303693/150‑1 (MQ=255)
                                                                                                                gacAAGTGAAACAGATAATATAAATCGCCCGACATGAAGTC‑GGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAGTAATGGTCACGTCGTGAACGTGGCTTTCCTGAATGCCCGCACCGCTGATACGTACAAACTcc                                       >  2:1201312/3‑150 (MQ=255)
                                                                                                                 gCAAGTGAAACAGATAATATAAATCGCCCGACATGAAGTC‑GGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAGTAATGGTCACGTCGTGAACGTGGCTTTCCTGAATGCCCGCACCGCTGATACGTACAAACTCCg                                      >  2:1610970/1‑150 (MQ=255)
                                                                                                                     gTGAAACAGATAATATAAATCGCCCGACATGAAGTC‑GGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAGTAATGGTCACGTCGTGAACGTGGCTTTCCTGAATGCCCGCACCGCTGATACGTACAAACTCCGCttt                                  >  2:881923/1‑150 (MQ=255)
                                                                                                                     gTGAAACAGATAATATAAATCGCCCGACATGAAGTC‑GGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAGTAATGGTCACGTCGTGAACGTGGCTTTCCTGAATGCCCGCACCGCTGATACGTACAAACTCCGCttt                                  <  1:800204/150‑1 (MQ=255)
                                                                                                                     gTGAAACAGATAATATAAATCGCCCGACATGAAGTC‑GGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAGTAATGGTCACGTCGTGAACGTGGCTTTCCTGAATGCCCGCACCGCTGATACGTACAAACTCCGCttt                                  >  2:366591/1‑150 (MQ=255)
                                                                                                                      tGAAACAGATAATATAAATCGCCCGACATGAAGTCGGGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAGTAATGGTCACGTCGTGAACGTGGCTTTCCTGAATGCCCGCACCGCTGATACGTACAAACTCCGCttt                                  <  1:1486234/150‑1 (MQ=255)
                                                                                                                      tGAAACAGATAATATAAATCGCCCGACATGAAGTC‑GGGCGAAGAGAATCAGGAGCCCAGACTGTAGTTCGTTGACTCTTTAGTAATGGTCACGTCGTGAACGTGGCTTTCCTGAATGCCCGCACCGCTGATACGTACAAACTCCGCTTTa                                 <  1:1085296/150‑1 (MQ=255)
                                                                                                                         aaCAGATAATATAAATCGCCCGACATGAAGTC‑GGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAGTAATGGTCACGTCGTGAACGTGGCTTTCCTGAATGCCCGCACCGCTGATACGTACAAACTCCGCTTTAGTa                              <  1:1031180/150‑1 (MQ=255)
                                                                                                                           cAGATAATATAAATCGCCCGACATGAAGTC‑GGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAGTAATGGTCACGTCGTGAACGTGGCTTTCCTGAATGCCCGCACCGCTGATACGTACAACCTCCGCTTTAGTACg                            >  1:1025948/1‑150 (MQ=255)
                                                                                                                           cAGATAATATAAATCGCCCGACATGAAGTC‑GGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAGTAATGGTCACGTCGTGAACGTGGCTTTCCTGAATGCCCGCACCGCTGATACGTACAAACTCCGCTTTAGTACg                            >  2:1353169/1‑150 (MQ=255)
                                                                                                                           cAGATAATATAAAACGCCCGACATCACGTC‑GGGCGAAGAGACTCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAGTAATGGTCACGTCGTGAACGTGGCTTTCCTGAATGCCCGCACCGCTGATACGTACAAACTCCGCTTTAGTACg                            <  1:864704/150‑1 (MQ=255)
                                                                                                                              ataatat‑aatCGCCCGACATGAAGTC‑GGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAGTAATGGTCACGTCGTGAACGTGGCTTTCATGAATGCCCGCACCGCTGATACGTACAAACTCCGCTTTAGTACGCAGt                        >  2:655224/1‑150 (MQ=255)
                                                                                                                                 atatAAATCGCCCGACATGAAGTC‑GGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAGTAATGGTCACGTCGTGAACGTGGCTTTCCTGAATGCCCGCACCGCTGATACATACAAACTCCGCTTTAGTACGCAGttg                      >  1:987841/1‑149 (MQ=255)
                                                                                                                                   ataAATCGCCCGACATGAAGTC‑GGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAGTAATGGTCACGTCGTGAACGTGGCTTTCCTGAATGCCCGCACCGCTGATACGTACAAACTCCGCTTTAGTACGCAGTtcgt                    >  2:1352056/1‑150 (MQ=255)
                                                                                                                                    taAATCGCCCGACATGAAGTC‑GGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACGCTTTAGTAATGGTCACGTCGTGCACGTGGCTTTCCTGAATGCCCGCACCGCTGATACGTACAAACTCCGCTGTTGTAGGGAGTtcgtg                   >  1:442052/1‑149 (MQ=255)
                                                                                                                                     tcttCGCCCGACTTCATGTC‑GGGCGAAGAGAATCAGGAGCCCTGACCGGAGTTCGGGGACTCATAAGTAATGCTCAGGTCGTGGACGTCGCGTTCCAGCAAGCACCCACCGCTGTTACGTACAAACGCCGCTTTATGAGGCAGTtggtcc                  >  2:331566/4‑149 (MQ=255)
                                                                                                                                     aaaTCGCCCGACATGAAGTC‑GGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAGTAATGGTCACGTCGTGAACGTGGCTTTCCTGAATGCCCGCACCGCTGATACGTACAAACTCCGCTTTAGTACGCAGTtcgtcg                  >  2:1488988/1‑150 (MQ=255)
                                                                                                                                      cttCGCCCGACTTAATGTC‑GGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAGTAATGGTCACGTCGTGAACGTGGCTTTCCTGAATGCCCGCACCGCTGATACGTACAAACTCCGCTTTAGTACGCAGTTCGTCGa                 >  1:1070609/3‑150 (MQ=255)
                                                                                                                                      aaTCGCCCGACATGAAGTC‑GGGCGAAGAGATGCTGGAGCCCTGACGGGTGTTCGGGGACTCTTTAGTAATGGTCACGTGGGGCACGTGGGTTGCCGGAATGCCCGCACCGCTGATACGTACCCACTGCGCGGAAGTACGCCGTTCGGCGa                 >  2:350907/1‑150 (MQ=255)
                                                                                                                                       aTCGCCCGACATGAAGTC‑GGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTTGTAATGGTCACGTCGTGAACGTGGCTTTCCTGAATGCCCGCACCGCTGATACGTACAAACTCCGCTTTAGTACGCAGGTCGTCGAt                >  2:88007/1‑150 (MQ=255)
                                                                                                                                         cGCCCGACTTCATGTC‑GGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAGTAATGGTCACGTCGTGAACGTGGCTTTCCTGAATGCCCGCACCGCTGATACGTACAAACTCCGGTTTAGTACGCAGTTCGTCGATat              >  2:905463/1‑149 (MQ=255)
                                                                                                                                         cGCCCGACATGAAGTC‑GGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAGTAATGGTCACGTCGTGAACGTGGCTTTCCTGAATGCCCGCACCGCTGATACGTACAAAc                                          <  2:289713/122‑1 (MQ=255)
                                                                                                                                         cGCCCGACATGAAGTC‑GGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAGTAATGGTCACGTCGTGAACGTGGCTTTCCTGAATGCCCGCACCGCTGATACGTACAAAc                                          >  1:289713/1‑122 (MQ=255)
                                                                                                                                         cGCCCGACATGAAGTC‑GGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAGTAATGGTCACGTCGTGAACGTGGCTTTCCTGAATGCCCGCACCGCTGATACGTACAAACTCCGCTTTAGTACGCAGTTCGTCGATAg              >  1:1335380/1‑150 (MQ=255)
                                                                                                                                            ccGACATGAAGGC‑GGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAGTAATGGTCACGTCGTGAACGTGGCTTTCCTGAATGCCCGCACCGCTGATACGTACAAACTCCGCTTTAGTACGCAGTTCGTCGATAGTAc           <  1:1421240/150‑1 (MQ=255)
                                                                                                                                              gACATGAAGTC‑GGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAGTAATGGTCACGTCGTGAACGTGGCTTTCCTGAATGCCCGCACCGCTGATACGTACAAACTCCGCTTTAGTACGCAGTTCGTCGATAGTACca         >  1:1472685/1‑150 (MQ=255)
                                                                                                                                                  tGAAGTC‑GGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAGTAATGGTCACGTCGTGAACGTGGCTTTCCTGAATGCCCGCACCGCTGATACGTACAAACTCCGCTTTAGTACGCAGTTCGTCGATAGTACCACAGc     >  1:18866/1‑150 (MQ=255)
                                                                                                                                                    aaGTC‑GGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAGTAATGGTCACGTCGTGAACGTGGCTTTCCTGAATGCCCGCACCGCTGATACGTACAAACTCCGCTTTAGTACGCAGTTCGTCGATAGTACCACAGCCg   <  2:593229/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                      
AAACACCCAGCTCACGCACGCGGCGCGCAACCAGTTGAGTGTACTGAGAACCGAAGTCCAGAATGAGGATGCGATGCTTATGAATGTTTTCCGTCATTGACGCTTATTCCGAGGCAAGTGAAACAGATAATATAAATCGCCCGACATGAAGTC‑GGGCGAAGAGAATCAGGAGCCCAGACGGTAGTTCGGGGACTCTTTAGTAATGGTCACGTCGTGAACGTGGCTTTCCTGAATGCCCGCACCGCTGATACGTACAAACTCCGCTTTAGTACGCAGTTCGTCGATAGTACCACAGCCGG  >  NC_000913/2632439‑2632737

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