Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,621,555 C→A 6.1% S120* (TCA→TAA)  purM → phosphoribosylformylglycinamide cyclo‑ligase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,621,5550CA 6.1% 250.3 / 3.2 75S120* (TCA→TAA) purMphosphoribosylformylglycinamide cyclo‑ligase
Reads supporting (aligned to +/- strand):  ref base C (36/34);  new base A (2/3);  total (38/37)
Fisher's exact test for biased strand distribution p-value = 6.74e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCGTCTGGCAATGGACTTAAAACGTCACGACACCATTGGTATTGATCTGGTCGCCATGTGCGTTAATGACCTGGTGGTGCAAGGTGCAGAGCCGCTGTTTTTCCTCGACTATTACGCAACCGGAAAACTGGATGTTGATACCGCTTCAGCGGTGATCAGCGGCATTGCGGAAGGTTGTCTGCAATCAGGCTGTTCACTGGTGGGTGGCGAAACGGCAGAAATGCCGGGGATGTATCACGGTGAGGATTACGATGTCGCGGGTTTCTGCGTTGGCGTGGTAGAAAAATCAG  >  NC_000913/2621409‑2621698
                                                                                                                                                  |                                                                                                                                               
gCGTCTGGCAATGGACTTAAAACGTCACGACACCATTGGTATTGATCTGGTCGCCATGTGCGTTAATGACCTGGTGGTGCAAGGTGCAGAGCCGCTGTTTTTCCTCGACTATTACGCAACCGGAAAACTGGATGTTGATACCGCTTCAGc                                                                                                                                              >  2:1071185/1‑150 (MQ=255)
   tCTGGCAATGGACTTAAAACGTCACGACACCATTGGTATTGATCTGGTCGCCATGTGCGTTAATGACCTGGTGGTGCAAGGTGCAGAGCCGCTGTTTTTCCTCGACTATTACGCAACCGGAAAACTGGATGTTGATACCGCTTCAGCGGt                                                                                                                                           <  2:1409526/150‑1 (MQ=255)
      ggCAATGGACTTAAAACGTCACGACACCATTGGTATTGATCTGGTCGCCATGTGCGTTAATGACCTGGTGGTGCAAGGTGCAGAGCCGCTGTTTTTCCTCGACTATTACGCAACCGGAAAACTGGATGTTGATACCGCTGAAGCGGtatc                                                                                                                                        <  2:1002637/150‑4 (MQ=255)
      ggCAATGGACTTAAAACGTCACGACACCATTGGTATTGATCTGGTCGCCATGTGCGTTAATGACCTGGTGGTGCAAGGTGCAGAGCCGCTGTTTTTCCTCGACTATTACGCAACCGGAAAACTGGATGTTGATACCGCTGAAGCGGtatc                                                                                                                                        <  2:1342207/150‑4 (MQ=255)
       gCAATGGACTTAAAACGTCACGACACCATTGGTATTGATCTGGTCGCCATGTGCGTTAATGACCTGGTGGTGCAAGGTGCAGAGCCGCTGTTTTTCCTCGACTATTACGCAACCGGAAAACTGGATGTTGATACCGCTGAAGCGGtatca                                                                                                                                       <  2:721503/150‑5 (MQ=255)
           tGGACTTAAAACGTCACGACACCATTGGTATTGATCTGGTCGCCATGTGCGTTAATGACCTGGTGGTGCAAGGTGCAGAGCCGCTGTTTTTCCTCGACTATTACGCAACCGGAAAACTGGATGTTGATACCGCTTCAGCGGTGATCAGCg                                                                                                                                   >  1:999588/1‑150 (MQ=255)
            ggACTTAAAACGTCACGACACCATTGGTATTGATCTGGTCGCCATGTGCGTTAATGACCTGGTGGTGCAAGGTGCAGAGCCGCTGTTTTTCCTCGACTATTACGCAACCGGAAAACTGGATGTTGATACCGCTTCAGCGGTGATCAGCgg                                                                                                                                  <  1:921513/150‑1 (MQ=255)
             gACTTAAAACGTCACGACACCATTGGTATTGATCTGGTCGCCATGTGCGTTAATGACCTGGTGGTGCAAGGTGCAGAGCCGCTGTTTTTCCTCGACTATTACGCAACCGGAAAACTGGATGTTGATACCGCTTCAGCGGTGATCAGCGGc                                                                                                                                 >  2:719583/1‑150 (MQ=255)
               cTTAAAACGTCACGACACCATTGGTATTGATCTGGTCGCCATGTGCGTTAATGACCTGGTGGTGCAAGGTGCAGAGCCGCTGTTTTTCCTCGACTATTACGCAACCGGAAAACTGGATGTTGATACCGCTTCAGCGGTGATCAGCGGCAt                                                                                                                               >  2:397294/1‑150 (MQ=255)
               cTTAAAACGTCACGACACCATTGGTATTGATCTAGTAGAAATCTGAATTAATTAAATAATTATTCAAGGTGCAGAGCCGCTGTTATTAATCGACTATTACGCCACAGGAAAACTGGATGTTGATACCGCTTCAGCGGTGATCAGCGGCAt                                                                                                                               >  2:1269642/1‑150 (MQ=255)
                 tAAAACGTCACGACACCATTGGTATTGATCTGGTCGCCATGTGCGTTAATGACCTGGTGGTGCAAGGTGCAGAGCCGCTGTTTTTCCTCGACTATTACGCAACCGGAAAACTGGATGTTGATACCGCTTCAGCGGTGATCAGCGGCATTg                                                                                                                             <  2:6378/150‑1 (MQ=255)
                  aaaaCGTCACGACACCATTGGTATTGATCTGGTCGCCATGTGCGTTAATGACCTGGTGGTGCAAGGTGCAGAGCCGCTGTTTTTCCTCGACTATTACGCAACCGGAAAACTGGATGTTGATACCGCTTCAGCGGTGATCAGCGGCATTGc                                                                                                                            >  2:1261441/1‑150 (MQ=255)
                    aaCGTCACGCCACCATTGGTATTGATCTGGTCGCCATGTGCGTTAATGACCTGGTGGTGCAAGGTGCAGAGCCGCTGTTTTTCCTCGACTATTACGCAACCGGAAAACTGGATGTTGATACCGCTTCAGCGGTGATCAGCGGCATTGCgg                                                                                                                          >  1:1175420/1‑150 (MQ=255)
                    aaCGTCACGACACCATTGGTATTGATCTGGTCGCCATGTGCGTTAATGACCTGGTGGTGCAAGGTGCAGAGCCGCTGTTTTTCCTCGACTATTACGCAACCGGAAAACTGGATGTTGATACCGCTTCAGCGGTGATCAGCGGCATTGCgg                                                                                                                          <  1:894362/150‑1 (MQ=255)
                      cGTCACGACACCATTGGTATTGATCTGGTCGCCATGTGCGTTAATGACCTGGTGGTGCAAGGTGCAGAGCCGCTGTTTTTCCTCGACTATTACGCAACCGGAAAACTGGATGTTGATACCGCTTCAGCGGTGATCAGCGGCATTGCGGaa                                                                                                                        <  1:1287330/150‑1 (MQ=255)
                        tCACGACACCATTGGTATTGATCTGGTCGCCATGTGCGTTAATGACCTGGTGGTGCAAGGTGCAGAGCCGCTGTTTTTCCTCGACTATTACGCAACCGGAAAACTGGATGTTGATACCGCTTCAGCGGTGATCAGCGGCATTGCGGAAgg                                                                                                                      <  2:441884/150‑1 (MQ=255)
                               acCATTGGTATTGATCTGGTCGCCATGTGCGTTAATGACCTGGTGGTGCAAGGTGCAGAGCCGCTGTTTTTCCTCGACTATTACGCAACCGGAAAACTGGATGTTGATACCGCTTCAGCGGTGATCAGCGGCATTGCGGAAGGTTGTCTg                                                                                                               <  1:1391592/150‑1 (MQ=255)
                                ccATTGGTATTGATCTGGTCGCCATGTGCGTTAATGACCTGGTGGTGCAAGGTGCAGAGCCGCTGTTTTTCCTCGACTATTACGCAACCGGAAAACTGGATGTTGATACCGCTTCAGCGGTGATCAGCGGCATTGCGGAAGGTTGTCTGc                                                                                                              <  1:1603197/150‑1 (MQ=255)
                                       tATTGATCTGGTCGCCATGTGCGTTAATGACCTGGTGGTGCAAGGTGCAGAGCCGCTGTTTTTCCTCGACTATTACGCAACCGGAAAACTGGATGTTGATACCGCTTCAGCGGTGATCAGCGGCATTGCGGAAGGTTGTCTGCAATCAgg                                                                                                       >  1:321150/1‑150 (MQ=255)
                                         ttGATCTGGTCGCCATGTGCGTTAATGACCTGGTGGTGCAAGGTGCAGAGCCGCTGTTTTTCCTCGACTATTACGCAACCGGAAAACTGGATGTTGATACCGCTTCAGCGGTGATCAGCGGCATTGCGGAAGGTTGTCTGCAATCAGGCt                                                                                                     >  2:610824/1‑150 (MQ=255)
                                           gATCTGGTCGCCATGTGCGTTAATGACCTGGTGGTGCAAGGTGCAGAGCCGCTGTTTTTCCTCGACTATTACGCAACCGGAAAACTGGATGTTGATACCGCTTCAGCGGTGATCAGCGGCATTGCGGAAGGTTGTCTGCAATCAGGCTGt                                                                                                   >  1:1614626/1‑150 (MQ=255)
                                           gATCTGGTCGCCATGTGCGTTAATGACCTGGTGGTGCAAGGTGCAGAGCCGCTGTTTTTCCTCGACTATTACGCAACCGGAAAACTGGATGTTGATACCGCTTCAGCGGGGATCAGCGGCATTGCGGAAGGTTGTCTGCAATCAGGCTGt                                                                                                   >  1:555288/1‑150 (MQ=255)
                                            aTCTGGTCGCCATGTGCGTTAATGACCTGGTGGTGCAAGGTGCAGAGCCGCTGTTTTTCCTCGACTATTACGCAACCGGAAAACTGGATGTTGATACCGCTTCAGCGGTGATCAGCGGCATTGCGGAAGGTTGTCTGCAATCAGGCTGtt                                                                                                  >  1:133414/1‑150 (MQ=255)
                                                   cGCCATGTGCGTTAATGACCTGGTGGTGCAAGGTGCAGAGCCGCTGTTTTTCCTCGACTATTACGCAACCGGAAAACTGGATGTTGATACCGCTTCAGCGGTGATCAGCGGCATTGCGGAAGGTTGTCTGCCATCAGGCTGTTCACtggt                                                                                           >  1:114300/1‑150 (MQ=255)
                                                    gCCATGTGCGTTAATGACCTGGTGGTGCAAGGTGCAGAGCCGCTGTTTTTCCTCGACTATTACGCAACCGGAAAACTGGATGTTGATACCGCTTCAGCGGTGATCAGCGGCATTGCGGAAGGTTGTCTGCAATCAGGCTGTTCACtggtg                                                                                          >  2:165625/1‑150 (MQ=255)
                                                       aTGTGCGTTAATGACCTGGTGGTGCAAGGTGCAGAGCCGCTGTTTTTCCTCGACTATTACGCAACCGGAAAACTGGATGTTGATACCGCTTCAGCGGTGATCAGCGGCATTGCGGAAGGTTGTCTGCAATCAGGCTGTTCACTggtgggt                                                                                       <  1:1014342/150‑1 (MQ=255)
                                                          tgCGTTAATGACCTGGTGGTGCAAGGTGCAGAGCCGCTGTTTTTCCTCGACTATTACGCAACCGGAAAACTGGATGTTGATACCGCTTCAGCGGTGATCAGCGGCATTGCGGAAGGTTGTCTGCAATCAGGCTGTTCACTGGTGGGTGGc                                                                                    <  1:397294/150‑1 (MQ=255)
                                                            cGTTAATGACCTGGTGGTGCAAGGTGCAGAGCCGCTGTTTTTCCTCGACTATTACGCAACCGGAAAACTGGATGTTGATACCGCTTCAGCGGTGATCAGCGGCATTGCGGAAGGTTGTCTGCAATCAGGCTGTTCACTGGTGGGTGGCGa                                                                                  <  1:1078545/150‑1 (MQ=255)
                                                              ttAATGACCTGGTGGTGCAAGGTGCAGAGCCGCTGTTTTTCCTCGAATATTACGCAACCGGAAAACTGGATGTTGATACCGATTAAGCGGTGATCAGCGGCATTGCGGAACGTTGTCTGCAATCAGGCTGTTCCCTGGTGGGTGGCGaaa                                                                                >  1:1605812/1‑150 (MQ=255)
                                                               tAATGACCTGGTGGTGCAAGGTGCAGAGCCGCTGTTTTTCCTCGACTATTACGCAACCGGAAAACTGGATGTTGATACCGCTTCAGCGGTGATCAGCGGCATTGCGGAAGGTTGTCTGCAATCAGGCTGTTCACTGGTGGGTGGCGAAAc                                                                               >  1:1283055/1‑150 (MQ=255)
                                                                 aTGACCTGGTGGTGCAAGGTGCAGAGCCGCTGTTTTTCCTCGACTATTGCGCAACCGGAAAACTGGATGTTGATACCGCTTCAGCGGTGATCAGCGGCATTGc                                                                                                                            <  2:772580/103‑1 (MQ=255)
                                                                 aTGACCTGGTGGTGCAAGGTGCAGAGCCGCTGTTTTTCCTCGACTATTGCGCAACCGGAAAACTGGATGTTGATACCGCTTCAGCGGTGATCAGCGGCATTGc                                                                                                                            >  1:772580/1‑103 (MQ=255)
                                                                  tGACCTGGTGGTGCAAGGTGCAGAGCCGCTGTTTTTCCTCGACTATTACGCAACCGGAAAACTGGATGTTGATACCGCTTCAGCGGTGATCAGCGGCATTGCGGAAGGTTGTCTGCAATCAGGCTGTTCACTGGTGGGTGGCGAAACGGc                                                                            <  2:1010976/150‑1 (MQ=255)
                                                                  tGACCTGGTGGTGCAAGGTGCAGAGCCGCTGTTTTTCCTCGACTATTACGCAACCGGAAAACTGGATGTTGATACCGCTTCAGCGGTGATCAGCGGCATTGCGGAAGGTTGTCTGCAATCAGGCTGTTCACTGGTGGGTGGCGAAACGGc                                                                            >  1:345255/1‑150 (MQ=255)
                                                                   gACCTGGTGGTGCAAGGTGCAGAGCCGCTGTTTTTCCTCGACTATTACGCAACCGGAAAACTGGATGTTGATACCGCTTCAGCGGTGATCAGCGGCATTGCGGAAGGTTGTCTGCAATCAGGCTGTTCACTGGTGGGTGGCGAAACGGCa                                                                           <  1:1232316/150‑1 (MQ=255)
                                                                    aCCTGGTGGTGCAAGGTGCAGAGCCGCTGTTTTTCCTCGACTATTACGCAACCGGAAAACTGGATGTTGATACCGCTTCAGCGGTGATCAGCGGCATTGCGGAAGGTTGTCTGCAATCAGGCTGTTCACTGGTGGGTGGCGAAACGGCAg                                                                          <  2:379804/150‑1 (MQ=255)
                                                                     ccTGGTGGTGCAAGGTGCAGAGCCGCTGTTTTTCCTCGACTATTACGCAACCGGAAAACTGGATGTTGATACCGCTTCAGCGGTGATCAGCGGCATTGCGGAAGGTTGTCTGCAATCAGGCTGTTCACTGGTGGGTGGCGAAACGGCAGa                                                                         <  2:256382/150‑1 (MQ=255)
                                                                      cTGGTGGTGCAAGGTGCAGAGCCGCTGTTTTTCCTCGACTATTACGCAACCGGAAAACTGGATGTTGATACCGCTTCAGCGGTGATCAGCGGCATTGCGGAAGGTTGTCTGCAATCAGGCTGTTCACTGGTGGGTGGCGAAACGGCAGaa                                                                        >  2:555428/1‑150 (MQ=255)
                                                                        ggtggtGCAAGGTGCAGAGCCGCTGTTTTTCCTCGACTATTACGCAACCGGAAAACTGGATGTTGATACCGCTTCAGCGGTGATCAGCGGCATTGCGGAAGGTTGTCTGCAATCAGGCTGTTCACTGGTGGGTGGCGAAACGGCAGAAAt                                                                      >  1:724198/1‑150 (MQ=255)
                                                                              gCAAGGTGCAGAGCCGCTGTTTTTCCTCGACTATTACGCAACCGGAAAACTGGATGTTGATACCGCTTCAGCGGTGATCAGCGGCATTGCGGAAGGTTGTCTGCAATCAGGCTGTTCACTGGTGGGTGGCGAAACGGCAGAAATGCCggg                                                                >  1:1457084/1‑150 (MQ=255)
                                                                                  ggTGCAGAGCCGCTGTTTTTCCTCGACTATTACGCAACCGGAAAACTGGATGTTGATACCGCTTCAGCGGTGATCAGCGGCATTGCGGAAGGTTGTCTGCAATCAGGCTGTTCACTGGTGGGTGGCGAAACGGCAGAAATGCCGGGGATg                                                            <  1:727982/150‑1 (MQ=255)
                                                                                   gTGCAGAGCCGCTGTTTTTCCTCGACTATTACGCAACCGGAAAACTGGATGTTGATACCGCTTCAGCGGTGATCAGCGGCATTGCGGAAGGTTGTCTGCAATCAGGCTGTTCACTGGTGGGTGGCGAAACGGCAGAAATGCCGGGGATGt                                                           >  1:1520880/1‑150 (MQ=255)
                                                                                     gCAGAGCCGCTGTTTTTCCTCGACTATTACGCAACCGGAAAACTGGATGTTGATACCGCTTCAGCGGTGATCAGCGGCATTGCGGAAGGTTGTCTGCAATCAGGCTGTTCACTGGTGGGTGGCGAAACGGCAGAAATGCCGGGGATGTAt                                                         >  1:1192953/1‑150 (MQ=255)
                                                                                      cAGAGCCGCTGTTTTTCCTCGACTATTACGCAACCGGAAAACTGGATGTTGATACCGCTTCAGCGGTGATCAGCGGCATTGCGGAAGGTTGTCTGCAATCAGGCTGTTCACTGGTGGGTGGCGAAACGGCAGAAATGCCGGGGATGTATc                                                        >  1:852834/1‑150 (MQ=255)
                                                                                       agagCCGCTGTTTTTCCTCGACTATTACGCAACCGGAAAACTGGATGTTGATACCGCTTCAGCGGTGATCAGCGGCATTGCGGAAGGTTGTCTGCAATCAGGCTGTTCACTGGTGGGTGGCGAAACGGCAGAAATGCCGGGGATGTATCa                                                       <  1:719583/150‑1 (MQ=255)
                                                                                        gagCCGCTGTTTTTCCTCGACTATTACGCAACCGGAAAACTGGATGTTGATACCGCTTCAGCGGTGATCAGCGGCATTGCGGAAGGTTGTCTGCAATCAGGCTGTTCACTGGTGGGTGGCGAAACGGCAGAAATGCCGGGGATGTATCAc                                                      <  1:610824/150‑1 (MQ=255)
                                                                                          gCCGCTGTTTTTCCTCGACTATTACGCAACCGGAAAACTGGATGTTGATACCGCTTCAGCGGTGATCAGCGGCATTGCGGAAGGTTGTCTGCAATCAGGCTGTTCACTGGTGGGTGGCGAAACGGCAGAAATGCCGGGGATGTATCACgg                                                    >  1:408992/1‑150 (MQ=255)
                                                                                            cGCTGTTTTTCCTCGACTATTACGCAACCGGAAAACTGGATGTTGATACCGCTTCAGCGGTGATCAGCGGCATTGCGGAAGGTTGTCTGCAATCAGGCTGTTCACTGGTGGGTGGCGAAACGGCAGAAATGCCGGGGATGTATCACGGTg                                                  >  1:1274735/1‑150 (MQ=255)
                                                                                              cTGTTTTTCCTCGACTATTACGCAACCGGAAAACTGGATGTTGATACCGCTTCAGCGGTGATCAGCGGCATTGCGGAAGGTTGTCTGCAATCAGGCTGTTCACTGGTGGGTGGCGAAACGGCAGAAATGCCGGGGATGTTTCACGGTGAg                                                <  2:1336564/150‑1 (MQ=255)
                                                                                               tGTTTTTCCTCGACTATTACGCAACCGGAAAACTGGATGTTGATACCGCTTCAGCGGTGATCAGCGGCATTGCGGAAGGTTGTCTGCAATCAGGCTGTTCACTGGTGGGTGGCGAAACGGCAGAAATGCCGGGGATGTATCACGGTGAgg                                               <  2:1614626/150‑1 (MQ=255)
                                                                                                gTTTTTCCTCGACTATTACGCAACCGGAAAACTGGATGTTGATACCGCTTCAGCGGTGATCAGCGGCATTGCGGAAGGTTGTCTGCAATCAGGCTGTTCACTGGTGGGTGGCGAAACGGCAGAAATGCCGGGGATGTATCACGGTGAGGa                                              <  2:1415556/150‑1 (MQ=255)
                                                                                                  ttttCCTCGACTATTACGCAACCGGAAAACTGGATGTTGATACCGCTTCAGCGGTGATCAGCGGCATTGCGGAAGGTTGTCTGCAATCAGGCTGTTCACTGGTGGGTGGCGAAACGGCAGAAATGCCGGGGATGTATCACGGTGAGGAtt                                            >  1:1391843/1‑150 (MQ=255)
                                                                                                    ttCCTCGACTATTACGCAACCGGAAAACTGGATGTTGATACCGCTTCAGCGGTGATCAGCGGCATTGCGGAAGGTTGTCTGCAATCAGGCTGTTCACTGGTGGGTGGCGAAACGGCAGAAATGCCGGGGATGTATCACGGTGAGGATTAc                                          <  2:999588/150‑1 (MQ=255)
                                                                                                     tCCTCGACTATTACGCAACCGGAAAACTGGATGTTGATACCGCTTCAGCGGTGATCAGCGGCATTGCGGAAGGTTGTCTGCAATCAGGCTGTTCACTGGTGGGTGGCGAAACGGCAGAAATGCCGGGGATGTATCACGGTGAGGATTACg                                         >  1:154321/1‑150 (MQ=255)
                                                                                                        tCGCCTATTACGCAACCGGAAAACTGGATGTTGATACCGCTTCAGCGGTGATCAGCGGCATTGCGGAAGGTTGTCTGCAATCAGGCTGTTCACTGGTGGGTGGCGAAACGGCAGAAATGCCGGGGATGTATCACGGTGAGGATTACGATg                                      >  2:1126493/1‑150 (MQ=255)
                                                                                                         cGACTATTACGCAACCGGAAAACTGGATGTTGATACCGCTTCAGCGGTGATCAGCGGCATTGCGGAAGGTTGTCTGCAATCAGGCTGTTCACTGGTGGGTGGCGAAACGGCAGAAATGCCGGGGATGTATCACGGTGAGGATTACGATGt                                     >  2:1017110/1‑150 (MQ=255)
                                                                                                         cGACTATTACGCAACCGGAAAACTGGATGTTGATACCGCTTCAGCGGTGATCAGCGGCATTGCGGAAGGTTGTCTGCAATCAGGCTGTTCACTGGTGGGTGGCGAAACGGCAGAAATGCCGGGGATGTATCACGGTGAGGATTACGATGt                                     <  1:973845/150‑1 (MQ=255)
                                                                                                         cGACTATTACGCAACCGGAAAACTGGATGTTGATACCGCTTCAGCGGTGATCAGCGGCATTGCGGAAGGTTGTCTGCAATCAGGCTGTTCACTGGTGGGTGGCGAAACGGCAGAAATGCCGGGGATGTATCACGGTGAGGATTACGATGt                                     >  2:550558/1‑150 (MQ=255)
                                                                                                              aTTACGCAACCGGAAAACTGGATGTTGATACCGCTTCAGCGGTGATCAGCGGCATTGCGGAAGGTTGTCTGCAATCAGGCTGTTCACTGGTGGGTGGCGAAACGGCAGAAATGCCGGGGATGTATCACGGTGAGGATTACGATGTCGCgg                                >  1:1158264/1‑150 (MQ=255)
                                                                                                               ttACGCAACCGGAAAACTGGATGTTGATACCGCTTCAGCGGTGATCAGCGGCATTGCGGAAGGTTGTCTGCAATCAGGCTGTTCACTGGTGGGTGGCGAAACGGCAGAAATGCCGGGGATGTATCACGGTGAGGATTACGATGTCGCggg                               <  2:808000/150‑1 (MQ=255)
                                                                                                                  cGCAACCGGAAAACTGGATGTTGATACCGCTTCAGCGGTGATCAGCGGCATTGCGGAAGGTTGTCTGCAATCAGGCTGTTCACTGGTGGGTGGCGAAACGGCAGAAATGCCGGGGATGTATCACGGTGAGGATTACGATGTCGCGGGttt                            >  1:618379/1‑150 (MQ=255)
                                                                                                                   gCAACCGGAAAACTGGATGTTGATACCGCTTCAGCGGTGATCAGCGGCATTGCGGAAGGTTGTCTGCAATCAGGCTGTTCACTGGTGGGTGGCGAAACGGCAGAAATGCCGGGGATGTATCACGGTGAGGATTACGATGTCGCGGGTTTc                           <  2:1283055/150‑1 (MQ=255)
                                                                                                                         ggAAAACTGGATGTTGATACCCCTTCAGCGGTGATCAGCGGCATTGCGCAAGGTTGTCTGCAATCAGGCAGTTCACAGGTGGGTGGCGAAACGGCAGAAATGCAGGGGATGTATCACGGTGAGCATAACGATGTCGCGGGTTTCTGCGtt                     >  2:919061/1‑150 (MQ=255)
                                                                                                                          gAAAACTGGATGTTGATACCGCTTCAGCGGTGATCAGCGGCATTGCGGAAGGTTGTCTGCAATCAGGCTGTTCACTGGTGGGTGGCGAAACGGCAGAAATGCCGGGGATGTATCACGGTGAGGATTACGATGTCGCGGGTTTCTGCGTTg                    <  2:618379/150‑1 (MQ=255)
                                                                                                                           aaaaCTGGATGTTTATACCTCTTCAGCGGTGATCAGCGGCATTGCGGAAGGTTGTCTGCAATCAGGCTGTTCACTGGTGGGTGGGGAAACGGCAGAAATGACGGGGATGTATCACGGTGAGGATTACGATGTCGCGGGTTTCTGCGTTgt                   <  1:1269642/150‑2 (MQ=255)
                                                                                                                             aaCTGGATGTTGATACCGCTTCAGCGGTGATCAGCGGCATTGCGGAAGGTTGTCTGCAATCAGGCTGTTCACTGGTGGGTGGCGAAACGGCAGAAATGCCGGGGATGTATCACGGTGAGGATTACGATGTCGCGGGTTTCTGCGTTGGCg                 <  2:852834/150‑1 (MQ=255)
                                                                                                                                  gATGTTGATACCGCTTCAGCGGTGATCAGCGGCATTGCGGAAGGTTGTCTGCAATCAGGCTGTTCACTGGTGGGTGGCGAAACGGCAGAAATGCCGGGGATGTATCACGGTGAGGATTACGATGTCGCGGGTTTCTGCGTTGGCGTGGTa            <  1:581939/150‑1 (MQ=255)
                                                                                                                                      ttGATACCGCTTCAGCGGTGATCAGCGGCATTGCGGAAGGTTGTCTGCAATCAGGCTGTTCACTGGTGGGTGGCGAAACGGCAGAAATGCCGGGGATGTATCACGGTGAGGATTACGATGTCGCGGGTTTCTGCGTTGGCGTGGTAGaaa        <  1:1071185/150‑1 (MQ=255)
                                                                                                                                       gatcaCCGCTGAAGCGGTGATCAGCGGCATTGCGGAAGGTTGTCTGCAATCAGGCTGTTCACTGGTGGGTGGCGAAACGGCAGAAATGCCGGGGATGTATCACGGTGAGGATTACGATGTCGCGGGTTTCTGCGTTGGCGTGGTAGaaaa       >  1:873555/5‑150 (MQ=255)
                                                                                                                                        taTACCGCTTCAGCGGTGATCAGCGGCATTGCGGAAGGTTGTCTGCAATCAGGCTGTTCACTGGTGGGTGGCGAAACGGCAGAAATGCCGGGGATGTATCACGGTGAGGATTACGATGTCGCGGGTTTCTGCGTTGGCGTGGTAGaaaaa      <  1:1363111/149‑1 (MQ=255)
                                                                                                                                        gATACCGCTTCAGCGGTGATCAGCGGCATTGCGGAAGGTTGTCTGCAATCAGGCTGTTCACTGGTGGGTGGCGAAACGGCAGAAATGCCGGGGATGTATCACGGTGAGGATTACGATGTCGCGGGTTTCTGCGTTGGCGTGGTAGaaaaa      <  1:165625/150‑1 (MQ=255)
                                                                                                                                          tACCGCTTCAGCGGTGATCAGCGGCATTGCGGAAGGTTGTCTGCAATCAGGCTGTTGACTGCTGGGTGGCGAAACGGCAGAAATGCCGGGGATGTATCACGGTGAGGATTACGATGTCGCGGGTTTCTGCGTTGGCGTGGTAGAAAAAta    >  1:1355884/1‑149 (MQ=255)
                                                                                                                                          tACCGCTTCAGCGGTGATCAGCGGCATTGCGGAAGGTTGTCTGCAATCAGGCTGTTCACTGGTGGGTGGCGAAACGGCAGAAATGCCGGGGATGTATCACGGTGAGGATTACGATGTCGCGGGTTTCTGCGTTGGCGTGGTAGAAAAATc    >  2:462365/1‑150 (MQ=255)
                                                                                                                                          tACCGCTTCAGCGGTGATCAGCGGCATTGCGGAAGGTTGTCTGCAATCAGGCTGTTCACTGGTGGGTGGCGAAACGGCAGAAATGCCGGGGATGTATCACGGTGAGGATTACGATGTCGCGGGTTTCTGCGTTGGCGTGGTAGAAAAATc    >  2:799993/1‑150 (MQ=255)
                                                                                                                                            ccGCTTCAGCGGTGATCAGCGGCATTGCGGAAGGTTGTCTGCAATCAGGCTGTTCACTGGTGGGTGGCGAAACGGCAGAAATGCCGGGGATGTATCACGGTGAGGATTACGATGTCGCGGGTTTCTGCGTTGGCGTGGTAGAAAAATCAg  <  2:354621/150‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                               
GCGTCTGGCAATGGACTTAAAACGTCACGACACCATTGGTATTGATCTGGTCGCCATGTGCGTTAATGACCTGGTGGTGCAAGGTGCAGAGCCGCTGTTTTTCCTCGACTATTACGCAACCGGAAAACTGGATGTTGATACCGCTTCAGCGGTGATCAGCGGCATTGCGGAAGGTTGTCTGCAATCAGGCTGTTCACTGGTGGGTGGCGAAACGGCAGAAATGCCGGGGATGTATCACGGTGAGGATTACGATGTCGCGGGTTTCTGCGTTGGCGTGGTAGAAAAATCAG  >  NC_000913/2621409‑2621698

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: