Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,692,286 G→C 7.1% H1086Q (CAC→CAG)  purL ← phosphoribosylformylglycinamide synthetase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,692,2860GC 7.1% 275.6 / 9.8 86H1086Q (CAC→CAG) purLphosphoribosylformylglycinamide synthetase
Reads supporting (aligned to +/- strand):  ref base G (41/37);  new base C (3/3);  total (45/41)
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.62e-01

CATACCCCCAGCGCCAGCGTTTGCGGACGGTGGAAGAAGGTTGCAAACTCATCGCGTACACGGTCATTGAACAGGATTGACTTCGCCCAACCTTCACCGGCACCCAGCACATCACCGTAGGAGAAACCACCGCACGCGACCAGGGCGTGGAAATCTTCCAGGCCCGTGCGTCCGGTCAGCAGGTCACTCATATGCACGTCGATAGCATCAAAGCCTGCACGGTGGAAAGCTGCCGCCATTTCAACATGCGAGTTCACGCCCTGCTCACGCAGTACGGCAACTTTCGGACGTGCGCC  >  NC_000913/2692140‑2692435
                                                                                                                                                  |                                                                                                                                                     
cATACCCCCAGCGCCAGCGTTTGCGGACGGTGGAAGAAGGTTGCAAACTCATCGCGTACACGGTCATTGAACAGGATTGACTTCGCCCAACCTTCACCGGCACCCAGCACATCACCGTAGGAGAAACCACCGCACGCGACCAGGGCGTgg                                                                                                                                                    >  1:354026/1‑150 (MQ=255)
cATACCCCCAGCGCCAGCGTTTGCGGACGGTGGAAGAAGGTTGCAAACTCATCGCGTACACGGTCATTGAACAGGATTGACTTCGCCCAACCTTCACCGGCACCCAGCACATCACCGTAGGAGAAACCACCGCACGCGACCAGGGCGTgg                                                                                                                                                    >  1:1526294/1‑150 (MQ=255)
   aCCCCCAGCGCCAGCGTTTGCGGACGGTGGAAGAAGGTTGCAAACTCATCGCGTACACGGTCATTGAACAGGATTGACTTCGCCCAACCTTCACCGGCACCCAGCACATCACCGTAGGAGAAACCACCGCACGCGACCAGGGCGTGGaaa                                                                                                                                                 >  2:1160709/1‑150 (MQ=255)
      cccAGCGCCAGCGTTTGCGGACGGTGGAAGAAGGTTGCAAACTCATCGCGTACACGGTCATTGAACAGGATTGACTTCGCCCAACCTTCACCGGCACCCAGCACATCACCGTAGGAGAAACCACCGCACGCGACCAGGGCGTGGAAATCt                                                                                                                                              >  2:842485/1‑150 (MQ=255)
      cccAGCGCCAGCGTTTGCGGACGGTGGAAGAAGGTTGCAAACTCATCGCGTACACGGTCATTGAACAGGATTGACTTCGCCCAACCTTCACCGGCACCCAGCACATCACCGTAGGAGAAACCACCGCACGCGACCAGGGCGTGGAAATCt                                                                                                                                              <  2:1382/150‑1 (MQ=255)
      cccAGCGCCAGCGTTTGCGGACGGTGGAAGAAGGTTGCAAACTCATCGCGTACACGGTCATTGAACAGGATTGACTTCGCCCAACCTTCACCGGCACCCAGCACATCACCGTAGGAGAAACCACCGCACGCGACCAGGGCGTGGAAATCt                                                                                                                                              >  1:1529270/1‑150 (MQ=255)
      cccAGCGCCAGCGTTTGCGGACGGTGGAAGAAGGTTGCAAACTCATCGCGTACACGGTCATTGAACAGGATTGACTTCGCCCAACCTTCACCGGCACCCAGCACATCACCGTAGGAGAAACCACCGCACGCGACCAGGGCGTGGAAATCt                                                                                                                                              >  1:1523613/1‑150 (MQ=255)
         agcgccagcgTTTGCGGACGGTGGAAGAAGGTTGCAAACTCATCGCGTACACGGTCATTGAACAGGATTGACTTCGCCCAACCTTCACCGGCACCCAGCACATCACCGTAGGAGAAACCACCGCACGCGACCAGGGCGTGGAAATCTTcc                                                                                                                                           <  1:152957/150‑1 (MQ=255)
         agcgccagcgTTTGCGGACGGTGGAAGAAGGTTGCAAACTCATCGCGTACACGGTCATTGAACAGGATTGACTTCGCCCAACCTTCACCGGCACCCAGCACATCACCGTAGGAGAAACCACCGCACGCGACCAGGGCGTGGAAATCTTcc                                                                                                                                           <  2:809441/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:503865/150‑1 (MQ=255)
                                            aaaCTCATCGCGTACACGGTCATTGAACAGGATTGACTTCGCCCAACCTTCACCGGCACCCAGCACATCACCGTAGGAGAAACCACCGCACGCGACCAGGGCGTGGAAATCTTCCAGGCCCGTGCGTCCGGTCAGCAGGTCACTCATATg                                                                                                        <  2:252755/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:821429/150‑1 (MQ=255)
                                                         acacGGTCATTGAACAGGATTGACTTCGCCCAACCTTCACCGGCACCCAGCACATCACCGTAGGAGAAACCACCGCACGCGACCAGGGCGTGGAAATCTTCCAGGCCCGTGCGTCCGGTCAGCAGGTCACTCATATGCACGTCGATAGCa                                                                                           <  1:1160709/150‑1 (MQ=255)
                                                          cacGGTCATTGAACAGGATTGACTTCGCCCAACCTTCACCGGCACCCAGCACATCACCGTAGGAGAAACCACCGCACGCGACCAGGGCGTGGAAATCTTCCAGGCCCGTGCGTCCGGTCAGCAGGTCACTCATATGCACGTCGATAGCAt                                                                                          <  2:1186419/150‑1 (MQ=255)
                                                            cGGTCATTGAACAGGATTGACTTCGCCCAACCTTCACCGGCACCCAGCACATCACCGTAGGAGAAACCACCGCACGCGACCAGGGCGTGGAAATCTTCCAGGCCCGTGCGTCCGGTCAGCAgg                                                                                                                   >  2:231079/1‑123 (MQ=255)
                                                            cGGTCATTGAACAGGATTGACTTCGCCCAACCTTCACCGGCACCCAGCACATCACCGTAGGAGAAACCACCGCACGCGACCAGGGCGTGGAAATCTTCCAGGCCCGTGCGTCCGGTCAGCAgg                                                                                                                   <  1:231079/123‑1 (MQ=255)
                                                            cGGTCATTGAACAGGATTGACTTCGCCCAACCTTCACCGGCACCCAGCACATCACCGTAGGAGAAACCACCGCACGCGACCAGGGCGTGGAAATCTTCCAGGCCCGTGCGTCCGGTCAGCAGGTCACTCATATGCACGTCGATAGCATCa                                                                                        >  1:581136/1‑150 (MQ=255)
                                                                cATTGAACAGGATTGACTTCGCCCAACCTTCACCGGCACCCAGCACATCACCGTAGGAGAAACCACCGCACGCGACCAGGGCGTGGAAATCTTCCAGGCCCGTGCGTCCGGTCAGCAGGTCACTCATATGCACGTCGATAGCATCAAAGc                                                                                    >  1:417976/1‑150 (MQ=255)
                                                                    gAACAGGATTGACTTCGCCCAACCTTCACCGGCACCCAGCACATCACCGTAGGAGAAACCACCGCACGCGACCAGGGCGTGGAAATCTTCCAGGCCCGTGCGTCCGGTCAGCAGGTCACTCATATGCACGTCGATAGCATCAAAGCCTGc                                                                                >  2:370923/1‑150 (MQ=255)
                                                                    gAACAGGATTGACTTCGCCCAACCTTCACCGGCACCCAGCACATCACCGTAGGAGAAACCACCGCACGCGACCAGGGCGTGGAAATCTTCCAGGCCCGTGCGTCCGGTCAGCAGGTCACTCATATGCACGTCGATAGCATCAAAGCCTGc                                                                                >  2:1277434/1‑150 (MQ=255)
                                                                    gAACAGGATTGACTTCGCCCAACCTTCACCGGCACCCAGCACATCACCGTAGGAGAAACCACCGCACGCGACCAGGGCGTGGAAATCTTCCAGGCCCGTGCGTCCGGTCAGCAGGTCACTCATATGCACGTCGATAGCATCAAAGCCTGc                                                                                >  1:5259/1‑150 (MQ=255)
                                                                     aaCAGGATTGACTTCGCCCAACCTTCACCGGCACCCAGCACATCACCGTAGGAGAAACCACCGCACGCGACCAGGGCGTGGAAATCTTCCAGGCCCGTGCGTCCGGTCAGCAGGTCACTCATATGCACGTCGATAGCATCAAAGCCTGCa                                                                               <  1:864896/150‑1 (MQ=255)
                                                                     aaCAGGATTGACTTCGCCCAACCTTCACCGGCACCCAGCACATCACCGTAGGAGAAACCACCGCACGCGACCAGGGCGTGGAAATCTTCCAGGCCCGTGCGTCCGGTCAGCAGGTCACTCATATGCACGTCGATAGCATCAAAGCCTGCa                                                                               >  2:306210/1‑150 (MQ=255)
                                                                        aGGATTGACTTCGCCCAACCTTCACCGGCACCCAGCACATCACCGTAGGAGAAACCACCGCACGCGACCAGGGCGTGGAAATCTTCCAGGCCCGTGCGTCCGGTCAGCAGGTCACTCATATGCACGTCGATAGCATCAAAGCCTGCACgg                                                                            <  2:1523613/150‑1 (MQ=255)
                                                                        aGGATTGACTTCGCCCAACCTTCACCGGCACCCAGCACATCACCGTAGGAGAAACCACCGCACGCGACCAGGGCGTGGAAATCTTCCAGGCCCGTGCGTCCGGTCAGCAGGTCACTCATATGCACGTCGATAGCATCAAAGCCTGCACgg                                                                            <  2:1529270/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                             >  2:743909/1‑150 (MQ=255)
                                                                                                                       ggAGAAACCACCGCACGCGACCAGGGCGTGGAAATCTTCCAGGCCCGTGCGTCCGGTCAGCAGGTCACTCATATGCACGTCGATAGCATCAAAGCCTGCACGGTGGAAAGCTGCCGCCATTTCAACATGCGAGTTCACGCCCTGCTCACg                             >  1:807915/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                                        >  2:592085/1‑122 (MQ=255)
                                                                                                                                        cgACCAGGGCGTGGAAATCTTCCAGGCCCGTGCGTCCGGTCAGCAGGTCACTCATATGCACGTCGATAGCATCAAAGCCTGCACGGTGGAAAGCTGCCGCCATTTCAACATGCGAGTTCACg                                        <  1:592085/122‑1 (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)
                                                                                                                                                  |                                                                                                                                                     
CATACCCCCAGCGCCAGCGTTTGCGGACGGTGGAAGAAGGTTGCAAACTCATCGCGTACACGGTCATTGAACAGGATTGACTTCGCCCAACCTTCACCGGCACCCAGCACATCACCGTAGGAGAAACCACCGCACGCGACCAGGGCGTGGAAATCTTCCAGGCCCGTGCGTCCGGTCAGCAGGTCACTCATATGCACGTCGATAGCATCAAAGCCTGCACGGTGGAAAGCTGCCGCCATTTCAACATGCGAGTTCACGCCCTGCTCACGCAGTACGGCAACTTTCGGACGTGCGCC  >  NC_000913/2692140‑2692435

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: ‑