Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,692,292 A→G 7.0% D1084D (GAT→GAC) ‡ purL ← phosphoribosylformylglycinamide synthetase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,692,2920AG 7.0% 250.1 / 7.6 85D1084D (GAT→GAC) ‡purLphosphoribosylformylglycinamide synthetase
Reads supporting (aligned to +/- strand):  ref base A (40/39);  new base G (3/3);  total (43/42)
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.36e-01

ACCCCCAGCGCCAGCGTTTGCGGACGGTGGAAGAAGGTTGCAAACTCATCGCGTACACGGTCATTGAACAGGATTGACTTCGCCCAACCTTCACCGGCACCCAGCACATCACCGTAGGAGAAACCACCGCACGCGACCAGGGCGTGGAAATCTTCCAGGCCCGTGCGTCCGGTCAGCAGGTCACTCATATGCACGTCGATAGCATCAAAGCCTGCACGGTGGAAAGCTGCCGCCATTTCAACATGCGAGTTCACGCCCTGCTCACGCAGTACGGCAACTTTCGGACGTGCGCCAGT  >  NC_000913/2692143‑2692438
                                                                                                                                                     |                                                                                                                                                  
aCCCCCAGCGCCAGCGTTTGCGGACGGTGGAAGAAGGTTGCAAACTCATCGCGTACACGGTCATTGAACAGGATTGACTTCGCCCAACCTTCACCGGCACCCAGCACATCACCGTAGGAGAAACCACCGCACGCGACCAGGGCGTGGaaa                                                                                                                                                    >  2:1160709/1‑150 (MQ=255)
   cccAGCGCCAGCGTTTGCGGACGGTGGAAGAAGGTTGCAAACTCATCGCGTACACGGTCATTGAACAGGATTGACTTCGCCCAACCTTCACCGGCACCCAGCACATCACCGTAGGAGAAACCACCGCACGCGACCAGGGCGTGGAAATCt                                                                                                                                                 >  2:842485/1‑150 (MQ=255)
   cccAGCGCCAGCGTTTGCGGACGGTGGAAGAAGGTTGCAAACTCATCGCGTACACGGTCATTGAACAGGATTGACTTCGCCCAACCTTCACCGGCACCCAGCACATCACCGTAGGAGAAACCACCGCACGCGACCAGGGCGTGGAAATCt                                                                                                                                                 >  1:1523613/1‑150 (MQ=255)
   cccAGCGCCAGCGTTTGCGGACGGTGGAAGAAGGTTGCAAACTCATCGCGTACACGGTCATTGAACAGGATTGACTTCGCCCAACCTTCACCGGCACCCAGCACATCACCGTAGGAGAAACCACCGCACGCGACCAGGGCGTGGAAATCt                                                                                                                                                 >  1:1529270/1‑150 (MQ=255)
   cccAGCGCCAGCGTTTGCGGACGGTGGAAGAAGGTTGCAAACTCATCGCGTACACGGTCATTGAACAGGATTGACTTCGCCCAACCTTCACCGGCACCCAGCACATCACCGTAGGAGAAACCACCGCACGCGACCAGGGCGTGGAAATCt                                                                                                                                                 <  2:1382/150‑1 (MQ=255)
      agcgccagcgTTTGCGGACGGTGGAAGAAGGTTGCAAACTCATCGCGTACACGGTCATTGAACAGGATTGACTTCGCCCAACCTTCACCGGCACCCAGCACATCACCGTAGGAGAAACCACCGCACGCGACCAGGGCGTGGAAATCTTcc                                                                                                                                              <  2:809441/150‑1 (MQ=255)
      agcgccagcgTTTGCGGACGGTGGAAGAAGGTTGCAAACTCATCGCGTACACGGTCATTGAACAGGATTGACTTCGCCCAACCTTCACCGGCACCCAGCACATCACCGTAGGAGAAACCACCGCACGCGACCAGGGCGTGGAAATCTTcc                                                                                                                                              <  1:152957/150‑1 (MQ=255)
       gcgccagcgTTTGCGGACGGTGGAAGAAGGTTGCAAACTCATCGCGTACACGGTCATTGAACAGGATTGACTTCGCCCAACCTTCACCGGCACCCAGCACATCACCGTAGGAGAAACCACCGCACGCGACCAGGGCGTGGAAATCTTCCa                                                                                                                                             <  1:808812/150‑1 (MQ=255)
        cgccagcgTTTGCGGACGGTGGAAGAAGGTTGCAAACTCATCGCGTACACGGTCATTGAACAGGATTGACTTCGCCCAACCTTCACCGGCACCCAGCACATCACCGTAGGAGAAACCACCGCACGCGACCAGGGCGTGGAAATCTTCCAg                                                                                                                                            <  1:1282787/150‑1 (MQ=255)
           cagcgTTTGCGGACGGTGGAAGAAGGTTGCAAACTCATCGCGTACACGGTCATTGAACAGGATTGACTTCGCCCAACCTTCACCGGCACCCAGCACATCACCGTAGGAGAAACCACCGCACGCGACCAGGGCCTGGAAGATTTCCACGcc                                                                                                                                         <  1:1542174/150‑1 (MQ=255)
              cgTTTGCGGACGGTGGAAGAAGGTTGCAAACTCGTCGCGTACACGGTCATTGAACAGGATTGACTTCGCCCAACCTTCACCGGCACCCAGCACATCACCGTAGGAGAAACCACCGCACGCGACCAGGGCGTGGAAATCTTCCAGGCCCGt                                                                                                                                      >  2:827871/1‑150 (MQ=255)
              cgTTTGCGGACGGTGGAAGAAGGGTGCAAACTCATCGCGTACACGGTCATTGAACAGGATTGACTTCGCCCAACCTTCACCGGCACCCAGCACATCACCGTAGGAGAAACCACCGCACGCGACCAGGGCCTGGAAGATTTCCACGccctg                                                                                                                                      <  2:155517/150‑3 (MQ=255)
               gTTTGCGGACGGTGGAAGAAGGTTGCAAACTCATCGCGTACACGGTCATTGAACAGGATTGACTTCGCCCAACCTTCACCGGCACCCAGCACATCACCGTAGGAGAAACCACCGCACGCGACCAGGGCCTGGAAGATTTCCACGccctgg                                                                                                                                     <  1:1592519/150‑4 (MQ=255)
                  tGCGGACGGTGGAAGAAGGTTGCAAACTCATCGCGTACACGGTCATTGAACAGGATTGACTTCGCCCAACCTTCACCGGCACCCAGCACATCACCGTAGGAGAAACCACCGCACGCGACCAGGGCGTGGAAATCTTCCAGGCCCGTGCGt                                                                                                                                  >  2:347118/1‑150 (MQ=255)
                    cGGACGGTGGAAGAAGGTTGCAAACTCATCGCGTACACGGTCATTGAACAGGATTGACTTCGCCCAACCTTCACCGGCACCCAGCACATCACCGTAGGAGAAACCACCGCACGCGACCAGGGCGTGGAAATCTTCCAGGCCCGTGCGTcc                                                                                                                                <  2:717539/150‑1 (MQ=255)
                        cGGTGGAAGAAGGTTGCAAACTCATCGCGTACACGGTCATTGAACAGGATTGACTTCGCCCAACCTTCACCGGCACCCAGCACATCACCGTAGGAGAAACCACCGCACGCGACCAGGGCGTGGAAATCTTCCAGGCCCGTGCGTCCGGTc                                                                                                                            <  1:31567/150‑1 (MQ=255)
                           tGGAAGAAGGTTGCAAACTCATCGCGTACACTGTCATTGAACAGGATTGACTTCGACCAACCTTCACCGGCACCCAGCACATCACCGTAGGAGAAACCACCGCACGCGACCAGGGCGTGGAGATCTTCCAGGCCCGTGCGTCCGGTcagc                                                                                                                         <  2:629080/150‑1 (MQ=255)
                           tGGAAGAAGGTTGCAAACTCATCGCGTACACGGTCATTGAACAGGATTGACTTCGCCCAACCTTCACCGGCACCCAGCACATCACCGTAGGAGAAACCACCGCACGCGACCAGGGCGTGGAAATCTTCCAGGCCCGTGCGTCCGGTcagc                                                                                                                         <  1:1082121/150‑1 (MQ=255)
                            ggAAGAAGGTTGCAAACTCATCGCGTACACGGTCATTGAACAGGATTGACTTCGCCCAACCTTCACCGGCACCCAGCACATCACCGTAGGAGAAACCACCGCACGCGACCAGGGCGTGGAAATCTTCCAGGCCCGTGCGTCCGGTcagca                                                                                                                        <  1:842485/150‑1 (MQ=255)
                                 aagGTTGCAAACTCATCGCGTACACGGTCATTGAACAGGATTGACTTCGCCCAACCTTCACCGGCACCCAGCACATCACCGTAGGAGAAACCACCGCACGCGACCAGGGCGTGGAAATCTTCCAGGCCCGTGCGTCCGGTCAGCAGGTCa                                                                                                                   >  2:803493/1‑150 (MQ=255)
                                   ggTTGCAAACTCATCGCGTACACGGTCATTGAACAGGATTGACTTCGCCCAACCTTCACCGGCACCCAGCACATCACCGTAGGAGAAACCACCGCACGCGACCAGGGCGTGGAAATCTTCCAGGCCCGTGCGTCCGGTCAGCAGGTCACt                                                                                                                 >  2:1468933/1‑150 (MQ=255)
                                    ggtccAAAATCATCGCCTACACGGTAATTGAAGATGATTCACTTCGACCAAACTTCAACGGACACCAGCACATACCCGGAGGAGCACCCACCGAACGCGTCCAGGGCGTGGAAATCTTCCAGGCCCGTGAGGGCGGACAGCACGTCAcac                                                                                                                >  2:727680/5‑148 (MQ=255)
                                    gTTGCAAACTCATCGCGTACACGGTCATTGAACAGGATTGACTTCGCCCAACCTTCACCGGCACCCAGCACATCACCGTAGGAGAAACCACCGCACGCGACCAGGGCGTGGAAATCTTCCAGGCCCGTGCGTCCGGTCAGCAGGTCACTc                                                                                                                >  2:64148/1‑150 (MQ=255)
                                      tGCAAACTCATCGCGTACACGGTCATTGAACAGGATTGACTTCGCCCAACCTTCACCGGCACCCAGCACATCACCGTAGGAGAAACCACCGCACGCGACCAGGGCGTGGAAATCTTCCAGGCCCGTGCGTCCGGTCAGCAGGTCACTCat                                                                                                              >  1:546009/1‑150 (MQ=255)
                                        cAAACTCATCGCGTACACGGTCATTGAACAGGATTGACTTCGCCCAACCTTCACCGGCACCCCGCACATCACCGTAGGAGAAACCACCGCACGCGACCAGGGCGTGGAAATCTTCCAGGCCCGTGCGTCCGGTCAGCAGGTCACTCatat                                                                                                            >  2:1196535/1‑150 (MQ=255)
                                         aaaCTCATCGCGTACACGGTCATTGAACAGGATTGACTTCGCCCAACCTTCACCGGCACCCAGCACATCACCGTAGGAGAAACCACCGCACGCGACCAGGGCGTGGAAATCTTCCAGGCCCGTGCGTCCGGTCAGCAGGTCACTCATATg                                                                                                           <  2:252755/150‑1 (MQ=255)
                                         aaaCTCATCGCGTACACGGTCATTGAACAGGATTGACTTCGCCCAACCTTCACCGGCACCCAGCACATCACCGTAGGAGAAACCACCGCACGCGACCAGGGCGTGGAAATCTTCCAGGCCCGTGCGTCCGGTCAGCAGGTCACTCATATg                                                                                                           <  2:503865/150‑1 (MQ=255)
                                          aaCTCATCGCGTACACGGTCATTGAACAGGATTGACTTCGCCCAACCTTCACCGGCACCCAGCACATCACCCTAGGAGAAACCACCGCACGCGACCAGGGCGTGGAAATCTTCCAGGCCCGTGCGTCCGGTCAGCAGGTCACTCATATGc                                                                                                          >  2:578685/1‑150 (MQ=255)
                                               aTCGCGTACACGGTCATTGAACAGGATTGACTTCGCCCAACCTTCACCGGCACCCAGCACATCACCGTAGGAGAAACCACCGCACGCGACCAGGGCGTGGAAATCTTCCAGGCCCGTGCGTCCGGTCAGCAGGTCACTCATATGCACGTc                                                                                                     >  1:650939/1‑150 (MQ=255)
                                                tCGCGTACACGGTCATTGAACAGGATTGACTTCGCCCAACCTTCACCGGCACCCAGCACATCACCGTAGGAGAAACCACCGCACGCGACCAGGGCGTGGAAATCTTCCAGGCCCGTGCGTCCGGTCAGCAGGTCACTCATATGCACGTCg                                                                                                    <  1:512824/150‑1 (MQ=255)
                                                 cgcgTACACGGTCATTGAACAGGATTGACTTCGCCCAACCTTCACCGGCACCCAGCACATCACCGTAGGAGAAACCACCGCACGCGACCAGGGCGTGGAAATCTTCCAGGCCCGTGCGTCCGGTCAGCAGGTCACTCATATGCACGTCGa                                                                                                   >  2:732800/1‑150 (MQ=255)
                                                 cgcgTACACGGTCATTGAACAGGATTGACTTCGCCCAACCTTCACCGGCACCCAGCACATCACCGTAGGAGAAACCACCGCACGCGACCAGGGCGTGGAAATCTTCCAGGCCCGTGCGTCCGGTCAGCAGGTCACTCATATGCACGTCGa                                                                                                   >  1:1186419/1‑150 (MQ=255)
                                                   cgTACACGGTCATTGAACAGGATTGACTTCGCCCAACCTTCACCGGCACCCAGCACATCACCGTAGGAGAAACCACCGCACGCGACCAGGGCGTGGAAATCTTCCAGGCCCGTGCGTCCGGTCAGCAGGTCACTCATATGCACGTCGAtg                                                                                                 <  1:1104975/150‑2 (MQ=255)
                                                   cgTACACGGTCATTGAACAGGATTGACTTCGCCCAACCTTCACCGGCACCCAGCACATCACCGTAGGAGAAACCACCGCACGCGACCAGGGCGTGGAAATCTTCCAGGCCCGTGCGTCCGGTCAGCAGGTCACTCATATGCACGTCGATa                                                                                                 >  1:1611759/1‑150 (MQ=255)
                                                      acacGGTCATTGAACAGGATTGACTTCGCCCAACCTTCACCGGCACCCAGCACATCACCGTAGGAGAAACCACCGCACGCGACCAGGGCGTGGAAATCTTCCAGGCCCGTGCGTCCGGTCAGCAGGTCACTCATATGCACGTCGATAGCa                                                                                              <  1:1160709/150‑1 (MQ=255)
                                                      acacGGTCATTGAACAGGATTGACTTCGCCCAACCTTCACCGGCACCCAGCACATCACCGTAGGAGAAACCACCGCACGCGACCAGGGCGTGGAAATCTTCCAGGCCCGTGCGTCCGGTCAGCAGGTCACTCATATGCACGTCGATAGCa                                                                                              <  1:821429/150‑1 (MQ=255)
                                                       cacGGTCATTGAACAGGATTGACTTCGCCCAACCTTCACCGGCACCCAGCACATCACCGTAGGAGAAACCACCGCACGCGACCAGGGCGTGGAAATCTTCCAGGCCCGTGCGTCCGGTCAGCAGGTCACTCATATGCACGTCGATAGCAt                                                                                             <  2:1186419/150‑1 (MQ=255)
                                                         cGGTCATTGAACAGGATTGACTTCGCCCAACCTTCACCGGCACCCAGCACATCACCGTAGGAGAAACCACCGCACGCGACCAGGGCGTGGAAATCTTCCAGGCCCGTGCGTCCGGTCAGCAgg                                                                                                                      <  1:231079/123‑1 (MQ=255)
                                                         cGGTCATTGAACAGGATTGACTTCGCCCAACCTTCACCGGCACCCAGCACATCACCGTAGGAGAAACCACCGCACGCGACCAGGGCGTGGAAATCTTCCAGGCCCGTGCGTCCGGTCAGCAgg                                                                                                                      >  2:231079/1‑123 (MQ=255)
                                                         cGGTCATTGAACAGGATTGACTTCGCCCAACCTTCACCGGCACCCAGCACATCACCGTAGGAGAAACCACCGCACGCGACCAGGGCGTGGAAATCTTCCAGGCCCGTGCGTCCGGTCAGCAGGTCACTCATATGCACGTCGATAGCATCa                                                                                           >  1:581136/1‑150 (MQ=255)
                                                             cATTGAACAGGATTGACTTCGCCCAACCTTCACCGGCACCCAGCACATCACCGTAGGAGAAACCACCGCACGCGACCAGGGCGTGGAAATCTTCCAGGCCCGTGCGTCCGGTCAGCAGGTCACTCATATGCACGTCGATAGCATCAAAGc                                                                                       >  1:417976/1‑150 (MQ=255)
                                                                 gAACAGGATTGACTTCGCCCAACCTTCACCGGCACCCAGCACATCACCGTAGGAGAAACCACCGCACGCGACCAGGGCGTGGAAATCTTCCAGGCCCGTGCGTCCGGTCAGCAGGTCACTCATATGCACGTCGATAGCATCAAAGCCTGc                                                                                   >  2:370923/1‑150 (MQ=255)
                                                                 gAACAGGATTGACTTCGCCCAACCTTCACCGGCACCCAGCACATCACCGTAGGAGAAACCACCGCACGCGACCAGGGCGTGGAAATCTTCCAGGCCCGTGCGTCCGGTCAGCAGGTCACTCATATGCACGTCGATAGCATCAAAGCCTGc                                                                                   >  1:5259/1‑150 (MQ=255)
                                                                 gAACAGGATTGACTTCGCCCAACCTTCACCGGCACCCAGCACATCACCGTAGGAGAAACCACCGCACGCGACCAGGGCGTGGAAATCTTCCAGGCCCGTGCGTCCGGTCAGCAGGTCACTCATATGCACGTCGATAGCATCAAAGCCTGc                                                                                   >  2:1277434/1‑150 (MQ=255)
                                                                  aaCAGGATTGACTTCGCCCAACCTTCACCGGCACCCAGCACATCACCGTAGGAGAAACCACCGCACGCGACCAGGGCGTGGAAATCTTCCAGGCCCGTGCGTCCGGTCAGCAGGTCACTCATATGCACGTCGATAGCATCAAAGCCTGCa                                                                                  >  2:306210/1‑150 (MQ=255)
                                                                  aaCAGGATTGACTTCGCCCAACCTTCACCGGCACCCAGCACATCACCGTAGGAGAAACCACCGCACGCGACCAGGGCGTGGAAATCTTCCAGGCCCGTGCGTCCGGTCAGCAGGTCACTCATATGCACGTCGATAGCATCAAAGCCTGCa                                                                                  <  1:864896/150‑1 (MQ=255)
                                                                     aGGATTGACTTCGCCCAACCTTCACCGGCACCCAGCACATCACCGTAGGAGAAACCACCGCACGCGACCAGGGCGTGGAAATCTTCCAGGCCCGTGCGTCCGGTCAGCAGGTCACTCATATGCACGTCGATAGCATCAAAGCCTGCACgg                                                                               <  2:1529270/150‑1 (MQ=255)
                                                                     aGGATTGACTTCGCCCAACCTTCACCGGCACCCAGCACATCACCGTAGGAGAAACCACCGCACGCGACCAGGGCGTGGAAATCTTCCAGGCCCGTGCGTCCGGTCAGCAGGTCACTCATATGCACGTCGATAGCATCAAAGCCTGCACgg                                                                               <  2:1523613/150‑1 (MQ=255)
                                                                      ggATTGACTTCGCCCAACCTTCACCGGCACCCAGCACATCACCGTAGGAGAAACCACCGCACGCGACCAGGGCATGGAAATCTTCCAGGCCCGTGCGTCCGGTCAGCAGGTCACTCATATGCACGTCGATAGCATCAAACCCTGCACggg                                                                              >  1:1555203/1‑149 (MQ=255)
                                                                        aTTGACTTCGCCCAACCTTCACCGGCACCCAGCACATCACCGTAGGAGAAACCACCGCACGCGACCAGGGCGTGGAAATCTTCCAGGCCCGTGCGTCCGGTCAGCAGGTCACTCATATGCACGTCGATAGCATCAAAGCCTGCACGGTgg                                                                            <  1:829109/150‑1 (MQ=255)
                                                                               tCGCCCAACCTTCACCGGCACCCAGCACATCACCGTAGGAGAAACCACCGCACGCGACCAGGGCGTGGAAATCTTCCAGGCCCGTGCGTCCGGTCAGCAGGTCACTCATATGCACGTCGATAGCATCAAAGCCTGCACGGTGGAAAGCTg                                                                     >  2:991805/1‑150 (MQ=255)
                                                                                 gCCCAACCTTCACCGGCACCCAGCACATCACCGTAGGAGAAACCACCGCACGCGACCAGGGCGTGGAAATCTTCCAGGCCCGTGCGTCCGGTCAGCAGGTCACTCATATGCACGTCGATAGCATCAAAGCCTGCACGGTGGAAAGCTgcc                                                                   >  1:396332/1‑150 (MQ=255)
                                                                                   ccAACCTTCACCGGCACCCAGCACATCACCGTAGGAGAAACCACCGCACGCGACCAGGGCGTGGAAATCTTCCAGGCCCGTGCGTCCGGTCAGCAGGTCACTCATATGCACGTCGATAGCATCAAAGCCTGCACGGTGGAAAGCTgccgc                                                                 >  1:554129/1‑150 (MQ=255)
                                                                                        cTTCACCGGCACCCAGCACATCACCGTAGGAGAAACCACCGCACGCGACCAGGGCGTGGAAATCTTCCAGGCCCGTGCGTCCGGTCAGCAGGTCACTCATATGCACGTCGATAGCATCAAAGCCTGCACGGTGGAAAGCTGCCGCCAttt                                                            <  1:732800/150‑1 (MQ=255)
                                                                                         ttCACCGGCACCCAGCACATCACCGTAGGAGAAACCACCGCACGCGACCAGGGCGTGGAAATCTTCCAGGCCCGTGCGTCCGGTCAGCAGGTCACTCATATGCACGTCGATAGCATCAAAGCCTGCACGGTGGAAAGCTGCCGCCATTTc                                                           <  2:241256/150‑1 (MQ=255)
                                                                                            aCCGGCACCCAGCACATCACCGTAGGAGAAACCACCGCACGCGACCAGGGCGTGGAAATCTTCCAGGCCCGTGCGTCCGGTCAGCAGGTCACTCATATGCACGTCGATAGCATCAAAGCCTGCACGGTGGAAAGCTGCCGCCATTTCAAc                                                        >  2:207476/1‑150 (MQ=255)
                                                                                              cGGCACCCAGCACATCACCGTAGGAGAAACCACCGCACGCGACCAGGGCGTGGAAATCTTCCAGGCCCGTGCGTCCGGTCAGCAGGTCACTCATATGCACGTCGATAGCATCAACGCCTGCACGGTGGAAAGCTGCCGCCATTTCACCAt                                                      >  2:1142421/1‑150 (MQ=255)
                                                                                                   cccAGCACATCACCGTAGGAGAAACCACCGCACGCGACCAGGGCGTGGAAATCTTCCAGGCCCGTGCGTCCGGTCAGCAGGTCACTCATATGCACGTCGATAGCATCAAAGCCTGCACGGTGGAAAGCTGCCGCCATTTCAACATGCGAg                                                 >  1:1134178/1‑150 (MQ=255)
                                                                                                   cccAGCACATCACCGTAGGAGAAACCACCGCACGCGACCAGGGCGTGGAAATCTTCCAGGCCCGTGCGTCCGGTCAGCAGGTCACTCATATGCACGTCGATAGCATCAAAGCCTGCACGGTGGAAAGCTGCCGCCATTTCAACATGCGAg                                                 >  2:1102249/1‑150 (MQ=255)
                                                                                                   cccAGCACATCACCGTAGGAGAAACCACCGCACGCGACCAGGGCGTGGAAATCTTCCAGGCCCGTGCGTCCGGTCAGCAGGTCACTCATATGCACGTCGATAGCATCAAAGCCTGCACGGTGGAAAGCTGCCGCCATTTCAACATGCGAg                                                 <  1:1034367/150‑1 (MQ=255)
                                                                                                   cccAGCACATCACCCTAGGAGAAACCACCGCACGCGACCAGGGCGTGGAAATCTTCCAGGCCCGTGCGTCCGGTCAGCAGGTCACTCATATGCACGTCGATAGCATCAAAGCCTGCACGGTGGAAAGCTGCCGCCATTTCAACATGCGAg                                                 <  1:578685/150‑1 (MQ=255)
                                                                                                    ccAGCACATCACCGTAGGAGAAACCACCGCACGCGACCAGGGCGTGGAAATCTTCCAGGCCCGTGCGTCCGGTCAGCAGGTCACTCATATGCACGTCGATAGCATCAAAGCCTGCACGGTGGAAAGCTGCCGCCATTTCAACATGCGAGt                                                <  1:803493/150‑1 (MQ=255)
                                                                                                         acaTCACCGTAGGAGAAACCACCGCACGCGACCAGGGCGTGGAAATCTTCCAGGCCCGTGCGTCCGGTCAGCAGGTCACTCATATGCACGTCGATAGCATCAAAGCCTGCACGGTGGAAAGCTGCCGCCATTTCAACATGCGAGTTCACg                                           <  1:891765/150‑1 (MQ=255)
                                                                                                               ccGTAGGAGAACCCACCGCACGCGACCAGGGCGTGGAAATCTTCCAGGCCCGGGCGTCCGGTCAGCAGGTCACTCATATGCCCGTCGATAGCATCAAAGCCTGCACGGTGGAAAGCTGCCGCCATTTCAACATGCGAGTTCACGCCCTGc                                     >  2:998034/1‑150 (MQ=255)
                                                                                                                  tAGGAGAAACCACCGCACGCGACCAGGGCGTGGAAATCTTCCAGGCCCGTGCGTCCGGTCAGCAGGTCACTCATATGCACGTCGATAGCATCAAAGCCTGCACGGTGGAAAGCTGCCGCCATTTCAACATGCGAGTTCACGCCCTGCTCa                                  <  2:1526294/150‑1 (MQ=255)
                                                                                                                  tAGGAGAAACCACCGCACGCGACCAGGGCGTGGAAATCTTCCAGGCCCGTGCGTCCGGTCAGCAGGTCACTCATATGCACGTCGATAGCATCAAAGCCTGCACGGTGGAAAGCTGCCGCCATTTCAACATGCGAGTTCACGCCCTGCTCa                                  <  1:727680/150‑1 (MQ=255)
                                                                                                                   aGGAGAAACCACCGCACGCGACCAGGGCGTGGAAATCTTCCAGGCCCGTGCGTCCGGTCAGCAGGTCACTCATATGCACGTCGATAGCATCAAAGCCTGCACGGTGGAAAGCTGCCGCCATTTCAACATGCGAGTTCACGCCCTGCTCAc                                 <  1:778208/150‑1 (MQ=255)
                                                                                                                    ggAGAAACCACCGCACGCGACCAGGGCGTGGAAATCTTCCAGGCCCGTGCGTCCGGTCAGCAGGTCACTCATATGCACGTCGATAGCATCAAAGCCTGCACGGTGGAAAGCTGCCGCCATTTCAACATGCGAGTTCACGCCCTGCTCACg                                >  1:807915/1‑150 (MQ=255)
                                                                                                                    ggAGAAACCACCGCACGCGACCAGGGCGTGGAAATCTTCCAGGCCCGTGCGTCCGGTCAGCAGGTCACTCATATGCACGTCGATAGCATCAAAGCCTGCACGGTGGAAAGCTGCCGCCATTTCAACATGCGAGTTCACGCCCTGCTCACg                                >  2:743909/1‑150 (MQ=255)
                                                                                                                     gagaAACCACCGCACGCGACCAGGGCGTGGAAATCTTCCAGGCCCGTGCGTCCGGTCAGCAGGTCACTCATATGCACGTCGATAGCATCAAAGCCTGCATGGTGGAAAGCTGCCGCCATTTCAACATGCGAGTTCACGCCCTGCTCACGc                               <  1:1388929/150‑1 (MQ=255)
                                                                                                                      agaAACCACCGCACGCGACCAGGGCGTGGAAATCTTCCAGGCCCGTGCGTCCGGTCAGCAGGTCACTCATATGCACGTCGATAGCATCAAAGCCTGCACGGTGGAAAGCTGCCGCCATTTCAACATGCGAGTTCACGCCCTGCTCACGCa                              <  1:548600/150‑1 (MQ=255)
                                                                                                                      agaAACCACCGCACGCGACCAGGGCGTGGAAATCTTCCAGGCCCGTGCGTCCGGTCAGCAGGTCACTCATATGCACGTCGATAGCATCAAAGCCTGCACGGTGGAAAGCTGCCGCCATTTCAACATGCGAGTTCACGCCCTGCTCACGCa                              >  2:504756/1‑150 (MQ=255)
                                                                                                                         aaCCCCCGCACGCGACCAGGGCATTGAAATCTTCCAGGCCCGTGCGTCGGGTCAGCAGGTCACTCATATGCACGTCGATAGCATCAAAGCCTGCACGGGGGAAAGCTGCCGCCATTTCAACATGCGAGTTCCCGCCCTGCTCACGCAGtt                           <  2:1555203/150‑2 (MQ=255)
                                                                                                                            caccGCACGCGACCAGGGCGTGGAAATCTTCCAGGCCCGTGCGTCCGGTCAGCAGGTCACTCATATGCACGTCGATAGCATCAAAGCCTGCACGGTGGAAAGCTGCCGCCATTTCAACATGCGAGTTCACGCCCTGCTCACGCAGTACgg                        >  1:1266682/1‑150 (MQ=255)
                                                                                                                             accGCACGCGACCAGGGCGTGGAAATCTTCCAGGCCCGTGCGTCCGGTCAGCAGGTCACTCATATGCACGTCGATAGCATCAAAGCCTGCACGGTGGAAAGCTGCCGCCATTTCAACATGCGAGTTCACGCCCTGCTCACGCAGTACGGc                       <  1:743909/150‑1 (MQ=255)
                                                                                                                                    gcgACCAGGGCGTGGAAATCTTCCAGGCCCGTGCGTCCGGTCAGCAGGTCACTCATATGCACGTCGATAGCATCAAAGCCTGCACGGTGGAAAGCTGCCGCCATTTCAACATGCGAGTTCACGCCCTGCTCACGCAGTACGGCAACTTTc                <  1:60610/150‑1 (MQ=255)
                                                                                                                                     cgACCAGGGCGTGGAAATCTTCCAGGCCCGTGCGTCCGGTCAGCAGGTCACTCATATGCACGTCGATAGCATCAAAGCCTGCACGGTGGAAAGCTGCCGCCATTTCAACATGCGAGTTCACg                                           <  1:592085/122‑1 (MQ=255)
                                                                                                                                     cgACCAGGGCGTGGAAATCTTCCAGGCCCGTGCGTCCGGTCAGCAGGTCACTCATATGCACGTCGATAGCATCAAAGCCTGCACGGTGGAAAGCTGCCGCCATTTCAACATGCGAGTTCACg                                           >  2:592085/1‑122 (MQ=255)
                                                                                                                                     cgACCAGGGCGTGGAAATCTTCCAGGCCCGTGCGTCCGGTCAGCAGGTCACTCATATGCACGTCGATACCATCAAAGCCTGCACGGTGGAAAGCTGCCGCCATTTCAACATGCGAGTTCACGCCCTGCTCACGCAGTACGGCAACTTTCg               >  2:682401/1‑150 (MQ=255)
                                                                                                                                      cgcacgggCCTGGAAGATTTCCACGCCCGTGCGTCCGGTCAGCAGGTCACTCATATGCACGTCGATAGCATCAAAGCCTGCACGGTGGAAAGCTGCCGCCATTTCAACATGCGAGTTCACGCCCTGCTCACGCAGTACGGCAACTTTCgg              >  2:1385809/6‑150 (MQ=255)
                                                                                                                                        cacgggCCTGGAAGATTTCCACGCCCGTGCGTCCGGTCAGCAGGTCACTCATATGCACGTCGATAGCATCAAAGCCTGCACGGTGGAAAGCTGCCGCCATTTCAACATGCGAGTTCACGCCCTGCTCACGCAGTACGGCAACTTTCGGAc            >  1:1575861/4‑150 (MQ=255)
                                                                                                                                         cAGGGCGTGGAAATCTTCCAGGCCCGTGCGTCCGGTCAGCAGGTCACTCATATGCACGTCGATAGCATCAAAGCCTGCACGGTGGAAAGCTGCCGCCATTTCAACATGCGAGTTCACGCCCTGCTCACGCAGTACGGCAACTTTCGGACg           <  1:133937/150‑1 (MQ=255)
                                                                                                                                         acgggCCTGGAAGATTTCCACGCCCGTGCGTCCGGTCAGCAGGTCACTCATATGCACGTCGATAGCATCAAAGCCTGCACGGTGGAAAGCTGCCGCCATTTCAACATGCGAGTTCACGCCCTGCTCACGCAGTACGGCAACTTTCGGACg           >  1:434057/3‑150 (MQ=255)
                                                                                                                                          aGGGCGTGGAAATCTTCCAGGCCCGTGCGTCCGGTCAGCAGGTCACTCATATGCACGTCGATAGCATCAAAGCCTGCACGGTGGAAAGCTGCCGCCATTTCAACATGCGAGTTCACGCCCTGCTCACGCAGTACGGCAACTTTCGGACGt          >  2:119307/1‑150 (MQ=255)
                                                                                                                                               gTGGAAATCTTCCAGGCCCGTGCGTCCGGTCAGCAGGTCACTCATATGCACGTCGATAGCATCAAAGCCTGCACGGTGGAAAGCTGCCGCCATTTCAACATGCGAGTTCACGCCCTGCTCACGCAGTACGGCAACTTTCGGACGTGCGcc     <  1:64148/150‑1 (MQ=255)
                                                                                                                                                  gAAATCTTCCAGGCCCGTGCGTCCGGTCAGCAGGTCACTCATATGCACGTCGATAGCATCAAAGCCTGCACGGTGGAAAGCTGCCGCCATTTCAACATGCGAGTTCACGCCCTGCTCACGCAGTACGGCAACTTTCGGACGTGCGCCAGt  >  2:1190349/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                  
ACCCCCAGCGCCAGCGTTTGCGGACGGTGGAAGAAGGTTGCAAACTCATCGCGTACACGGTCATTGAACAGGATTGACTTCGCCCAACCTTCACCGGCACCCAGCACATCACCGTAGGAGAAACCACCGCACGCGACCAGGGCGTGGAAATCTTCCAGGCCCGTGCGTCCGGTCAGCAGGTCACTCATATGCACGTCGATAGCATCAAAGCCTGCACGGTGGAAAGCTGCCGCCATTTCAACATGCGAGTTCACGCCCTGCTCACGCAGTACGGCAACTTTCGGACGTGCGCCAGT  >  NC_000913/2692143‑2692438

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