Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 30,053 G→C 5.7% G135R (GGC→CGC)  carA → carbamoyl phosphate synthetase subunit alpha

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00091330,0530GC 5.7% 295.5 / 7.0 87G135R (GGC→CGC) carAcarbamoyl phosphate synthetase subunit alpha
Reads supporting (aligned to +/- strand):  ref base G (44/38);  new base C (3/2);  total (47/40)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.33e-01

GCTGATTGCCAGCAACTTCCGTAATACCGAAGACCTCTCTTCTTACCTGAAACGCCATAACATCGTGGCGATTGCCGATATCGATACCCGTAAGCTGACGCGTTTACTGCGCGAGAAAGGCGCACAGAATGGCTGCATTATCGCGGGCGATAACCCGGATGCGGCGCTGGCGTTAGAAAAAGCCCGCGCGTTCCCAGGTCTGAATGGCATGGATCTGGCAAAAGAAGTGACCACCGCAGAAGCCTATAGCTGGACACAAGGGAGCTGGACGTTGACCGGTGGCCTGCCAGAAGCG  >  NC_000913/29908‑30202
                                                                                                                                                 |                                                                                                                                                     
gCTGATTGCCAGCAACTTCCGTAATACCGAAGACCTCTCTTCTTACCTGAAACGCCATAACATCGTGGCGATTGCCGATATCGATACCCGTAAGCTGACGCGTTTACTGCGCGAGAAAGGCGCACAGAATGGCTGCATTATCGCGGGCGa                                                                                                                                                   >  1:191914/1‑150 (MQ=255)
 cTGATTGCCAGCAACTTCCGTAATACCGAAGACCTCTCTTCTTACCTGAAACGCCATAACATCGTGGCGATTGCCGATATCGATACCCGTAAGCTGACGCGTTTACTGCGCGAGAAAGGCGCACAGAATGGCTGCATTATCGCGGGCGAt                                                                                                                                                  <  1:722880/150‑1 (MQ=255)
  tGATTGCCAGCAACTTCCGTAATACCGAAGACCTCTCTTCTTACCTGAAACGCCATAACATCGTGGCGATTGCCGATATCGATACCCGTAAGCTGACGCGTTTACTGCGCGAGAAAGGCGCACAGAATGGCTGCATTATCGCGGGCGATa                                                                                                                                                 >  2:1597133/1‑150 (MQ=255)
  tGATTGCCAGCAACTTCCGTAATACCGAAGACCTCTCTTCTTACCTGAAACGCCATAACATCGTGGCGATTGCCGATATCGATACCCGTAAGATGACGCGTTTACTGCGCGAGAAAGGCGCACAGAATGGCTGCATTATCGCGGGCGATa                                                                                                                                                 >  2:800741/1‑150 (MQ=255)
    aTTGCCAGCAACTTCCGTAATACCGAAGACCTCTCTTCTTACCTGAAACGCCATAACATCGTGGCGATTGCCGATATCGATACCCGTAAGCTGACGCGTTTACTGCGCGAGAAAGGCGCACAGAATGGCTGCATTATCGCGGGCGATAAc                                                                                                                                               >  2:1346485/1‑150 (MQ=255)
      tGCCAGCAACTTCCGTAATACCGAAGACCTCTCTTCTTACCTGAAACGCCATAACATCGTGGCGATTGCCGATATCGATACCCGTAAGCTGACGCGTTTACTGCGCGAGAAAGGCGCACAGAATGGCTGCATTATCGCCCGCGATaatgc                                                                                                                                             <  1:479258/150‑4 (MQ=255)
       gCCAGCAACTTCCGTAATACCGAAGACCTCTCTTCTTACCTGAAACGCCATAACATCGTGGCGATTGCCGATATCGATACCCGTAAGCTGACGCGTTTACTGCGCGAGAAAGGCGCACAGAATGGCTGCATTATCGCGGGCGATAACCCg                                                                                                                                            >  1:910023/1‑150 (MQ=255)
       gCCAGCAACTTCCGTAATACCGAAGACCTCTCTTCTTACCTGAAACGCCATAACATCGTGGCGATTGCCGATATCGATACCCGTAAGCTGACGCGTTTACTGCGCGAGAAAGGCGCACAGAATGGCTGCATTATCGCGGGCGATAACCCg                                                                                                                                            >  1:903988/1‑150 (MQ=255)
       gCCAGCAACTTCCGTAATACCGAAGACCTCTCTTCTTACCTGAAACGCCATAACATCGTGGCGATTGCCGATATCGATACCCGTAAGCTGACGCGTTTACTGCGCGAGAAAGGCGCACAGAATGGCTGCATTATCGCGGGCGACAACCCg                                                                                                                                            >  2:1443345/1‑150 (MQ=255)
        ccAGCAACTTCCGTAATACCGAAGACCTCTCTTCTTACCTGAAACGCCATAACATCGTGGCGATTGCCGATATCGATACCCGTAAGCTGACGCGTTTACTGCGCGAGAAAGGCGCACAGAATGGCTGCATTATCGCGGGCGATAACCCgg                                                                                                                                           >  2:273643/1‑150 (MQ=255)
         cAGCAACTTCCGTAATACCGAAGACCTCTCTTCTTACCTGAAACGCCATAACATCGTGGCGATTGCCGATATCGATACCCGTAAGCTGACGCGTTTACTGCGCGAGAAAGGCGCACAGAATGGCTGCATTATCGCCCGCGATaatgcagc                                                                                                                                          <  1:808652/150‑7 (MQ=255)
         cAGCAACTTCCGTAAAACCGAAGACCTCTCTTCTAACCTGAAACGCCATAACATCGTGGCGATTGCCGAAATCGAAACCCGTAAGCTGACGCCTTTACTGCGCGAGAAAGGCGCACAGAATGGCTGCATTAACGCGGGCGATAACCCGGa                                                                                                                                          >  2:1499070/1‑150 (MQ=255)
           gCAACTTCCGTAATACCGAAGACCTCTCTTCTTACCTGAAACGCCATAACATCGTGGCGATTGCCGATATCGATACCCGTAAGCTGACGCGTTTACTGCGCGAGAAAGGCGCACAGAATGGCTGCATTATCGCGGGCGATAACCCGGATg                                                                                                                                        >  2:399943/1‑150 (MQ=255)
             aaCTTCCGTAATACCGAAGACCTCTCTTCTTACCTGAAACGCCATAACATCGTGGCGATTGCCGATATCGATACCCGTAAGCTGACGCGTTTACTGCGCGAGAAAGGCGCACAGAATGGCTGCATTATCGCGGGCGATAACCCGGATgcg                                                                                                                                      <  2:495707/150‑1 (MQ=255)
             aaCTTCCGTAATACCGAAGACCTCTCTTCTTACCTGAAACGCCATAACATCGTGGCGATTGCCGATATCGATACCCGTAAGCTGACGCGTTTACTGCGCGAGAAAGGCGCACAGAATGGCTGCATTATCGCGGGCGATAACCCGGATgcg                                                                                                                                      <  1:22427/150‑1 (MQ=255)
               cTTCCGTAATACCGAAGACCTCTCTTCTTACCTGAAACGCCATAACATCGTGGCGATTGCCGATATCGATACCCGTAAGCTGACGCGTTTACTGCGCGAGAAAGGCGCACAGAATGGCTGCATTATCGCGGGCGATAACCCGGATgcggc                                                                                                                                    >  2:775547/1‑150 (MQ=255)
               cTTCCGTAATACCGAAGACCTCTCTTCTTACCTGAAACGCCATAACATCGTGGCGATTGCCGATATCGATACCCGTAAGCTGACGCGTTTACTGCGCGAGAAAGGCGCACAGAATGGCTGCATTATCGCGGGCGATAACCCGGATgcggc                                                                                                                                    >  1:836695/1‑150 (MQ=255)
                   cGTAATACCGAAGACCTCTCTTCTTACCTGAAACGCCATAACATCGTGGCGATTGCCGATATCGATACCCGTAAGCTGACGCGTTTACTGCGCGAGAACGGCGCACAGAATGGCTGCATTATCGCGGGCGATAACCCGGATGCGGCGCTg                                                                                                                                <  1:696606/150‑1 (MQ=255)
                   cGTAATACCGAAGACCTCTCTTCTTACCTGAAACGCCATAACATCGTGGCGATTGCCGATATCGATACCCGTAAGCTGACGCGTTTACTGCGCGAGAAAGGCGCACAGAATGGCTGCATTATCGCGGGCGATAACCCGGATGCGGCGCTg                                                                                                                                <  1:399943/150‑1 (MQ=255)
                     tcaTACCGAAGACCTCTCTTCTTACCTGAAACGCCATAACATCGTGGCGATTGCCGATAGCGATACCCGTAAGCTGACGCGTTTACTGCGCGAGAAAGGCGCACAGAATGGCTGCATTATCGCGGGCGATAACCCGGATGCGGCGCGGGc                                                                                                                              <  1:715753/148‑1 (MQ=255)
                      aaTACCGAAGACCTCTCTTCTTACCTGAAACGCCATAACATCGTGGCGATTGCCGATATCGATACCCGTAAGCTGACGCGTTTACTGCGCGAGAAAGGCGCACAGAATGGCTGCATTATCGCGGGCGATAACCCGGATGCGGCGCTGGCg                                                                                                                             <  2:585786/150‑1 (MQ=255)
                        tACCGAAGACCTCTCTTCTTACCTGAAACGCCATAACATCGTGGCGATTGCCGATATCGATACCCGTAAGCTGACGCGTTTACTGCGCGAGAAAGGCGCACAGAATGGCTGCATTATCGCGGGCGATAACCCGGATGCGGCGCTGGCGtt                                                                                                                           <  2:560631/150‑1 (MQ=255)
                          ccGAAGACCTCTCTTCTTACCTGAAACGCCATAACATCGTGGCGATTGCCGATATCGATACCCGTAAGCTGACGCGTTTACTGCGCGAGAAAGGCGCACAGAATGGCTGCATTATCGCGGGCGATAACCCGGATGCGGCGCTGGCGTTAg                                                                                                                         >  2:592483/1‑150 (MQ=255)
                            gAAGACCTCTCTTCTTACCTGAAACGCCATAACATCGTGGCGATTGCCGATATCGATACCCGTAAGCTGACGCGTTTACTGCGCGAGAAAGGCGCACAGAATGGCTGCATTATCGCGGGCGATAACCCGGATGCGGCGCTGGCGTTAGaa                                                                                                                       >  1:221618/1‑150 (MQ=255)
                              aGACCTCTCTTCTTACCTGAAACGCCATAACATCGTGGCGATTGCCGATATCGATACCCGTAAGCTGACGCGTTTACTGCGCGAGAAAGGCGCACAGAATGGCTGCATTATCGCGGGCGATAACCCGGATGCGGCGCTGGCGTTAGaaaa                                                                                                                     <  1:1009977/150‑1 (MQ=255)
                                 ccTCTCTTCTTACCTGAAACGCCATAACATCGTGGCGATTGCCGATATCGATACCCGTAAGCTGACGCGTTTACTGCGCGAGAAAGGCGCACAGAATGGCTGCATTATCGCGGGCGATAACCCGGATGCGGCGCTGGCGTTAGAAAAAGc                                                                                                                  >  1:712135/1‑150 (MQ=255)
                                 ccTCTCTTCTTACCTGAAACGCCATAACATCGTGGCGATTGCCGATATCGATACCCGTAAGCTGACGCGTTTACTGCGCGAGAAAGGCGCACAGAATGGCTGCATTATCGCGGGCGATAACCCGGATGCGGCGCTGGCGTTAGAAAAAGc                                                                                                                  >  1:1086057/1‑150 (MQ=255)
                                  ctctctTCTTACCTGAAACGCCATAACATCGTGGCGATTGCCGATATCGATACCCGTAAGCTGACGCGTTTACTGCGCGAGAAAGGCGCACAGAATGGCTGCATTATCGCGGGCGATAACCAGGATGCGGCGCTGGCGTTAGAAAAAGcc                                                                                                                 >  1:494914/1‑150 (MQ=255)
                                    ctctTCTTACCTGAAACGCCATAACATCGTGGCGATTGCCGATATCGATACCCGTAAGCTGACGCGTTTACTGCGCGAGAAAGGCGCACAGAATGGCTGCATTATCGCGGGCGATAACCCGGATGCGGCGCTGGCGTTAGAAAAAGCCcg                                                                                                               >  2:649865/1‑150 (MQ=255)
                                    ctctTCTTACCTGAAACGCCATAACATCGTGGCGATTGCCGATATCGATACCCGTAAGCTGACGCGTTTACTGCGCGAGAAAGGCGCACAGAATGGCTGCATTATCGCGGGCGATAACCCGGATGCGGCGCTGGCGTTAGAAAAAGCCcg                                                                                                               <  1:1659047/150‑1 (MQ=255)
                                      cttcttACCTGAAACGCCATAACATCGTGGCGATTGCCGATATCGATACCCGTAAGCTGACGCGTTTACTGCGCGAGAAAGGCGCACAGAATGGCTGCATTATCGCGGGCGATAACCCGGATGCGGCGCTGGCGTTAGAAAAAGCCcgcg                                                                                                             <  2:421015/150‑1 (MQ=255)
                                        tcttACCTGAAACGCCATAACATCGTGGCGATTGCCGATATCGATACCCGTAAGCTGACGCGTTTACTGCGCGAGAAAGGCGCACAGAATGGCTGCATTATCGCGGGCGATAACCCGGATGCGGCGCTGGCGTTAGAAAAAGCCcgcgcg                                                                                                           <  1:649865/150‑1 (MQ=255)
                                             ccTGAAACGCCATAACATCGTGGCGATTGCCGATATCGATACCCGTAAGCTGACGCGTTTACTGCGCGAGAAAGGCGCACAGAATGGCTGCATTATCGCGGGCGATAACCCGGATGCGGCGCTGGCGTTAGAAAAAGCCCGCGCGTTccc                                                                                                      >  1:1078104/1‑150 (MQ=255)
                                              cTGAAACGCCATAACATCGTGGCGATTGCCGATATCGATACCCGTAAGCTGACGCGTTTACTGCGCGAGAAAGGCGCACAGAATGGCTGCATTATCGCGGGCGATAACCCGGATGCGGCGCTGGCGTTAGAAAAAGCCCGCGCGTTCCCa                                                                                                     <  1:142845/150‑1 (MQ=255)
                                                gAAACGCCATAACATCGTGGCGATTGCCGATATCGATACCCGTAAGCTGACGCGTTTACTGCGCGAGAAAGGCGCACAGAATGGCTGCATTATCGCGGGCGATAACCCGGATGCGGCGCTGGCGTTAGAAAAAGCCCGCGCGTTCCCAgg                                                                                                   >  2:1461538/1‑150 (MQ=255)
                                                  aaCGCCATAACATCGTGGCGATTGCCGATATCGATACCCGTAAGCTGACGCGTTTACTGCGCGAGAAAGGCGCACAGAATGGCTGCATTATCGCGGGCGATAACCCGGATGCGGCGCTGGCGTTAGAAAAAGCCCGCGCGTTCCCAGGTc                                                                                                 <  1:502304/150‑1 (MQ=255)
                                                   aCGCCATAACATCGTGGCGATTGCCGATATCGATACCCGTAAGCTGACGCGTTTACTGCGCGAGAAAGGCGCACAGAATGGCTGCATTATCGCGGGCGATAACCCGGATGCGGCGCTGGCGTTAGAAAAAGCCCGCGCGTTCCCAGGTCt                                                                                                <  2:749061/150‑1 (MQ=255)
                                                   aCGCCATAACATCGTGGCGATTGCCGATATCGATACCCGTAAGCTGACGCGTTTACTGCGCGAGAAAGGCGCACAGAATGGCTGCATTATCGCGGGCGATAACCCGGATGCGGCGCTGGCGTTAGAAAAAGCCCGCGCGTTCCCAGGTCt                                                                                                <  2:569160/150‑1 (MQ=255)
                                                    cGCCATAACATCGTGGCGATTGCCGATATCGATACCCGTAAGCTGACGCGTTTACTGCGCGAGAAAGGCGCACAGAATGGCTGCATTATCGCGGGCGATAACCCGGATGCGGCGCTGGCGTTAGAAAAAGCCCGCGCGTTCCCAGGTCTg                                                                                               <  1:539059/150‑1 (MQ=255)
                                                             aTCGTGGCGATTGCCGATATCGATACCCGTAAGCTGACGCGTTTACTGCGCGAGAAAGGCGCACAGAATGGCTGCATTATCGCGGGCGATAACCCGGATGCGGCGCTGGCGTTAGAAAAAGCCCGCGCGTTCCCAGGTCTGAATGGCATg                                                                                      <  2:221618/150‑1 (MQ=255)
                                                               cGTGGCGATTGCCGATATCGATACCCGTAAGCTGACGCGTTTACTGCGCGAGAAAGGCGCACAGAATGGCTGCATTATCGCGGGCGATAACCCGGATGCGGCGCTGGCGTTAGAAAAAGCCCGCGCGTTCCCAGGTCTGAATGGCATGGa                                                                                    >  2:1074517/1‑150 (MQ=255)
                                                                    cGATTGCCGATATCGATACCCGTAAGCTGACGCGTTTACTGCGCGAGAAAGGCGCACAGAATGGCTGCATTATCGCGGGCGATAACCCGGATGCGGCGCTGGCGTTAGAAAAAGCCCGCGCGTTCCCAGGTCTGAATGGCATGGATCTgg                                                                               <  2:324246/150‑1 (MQ=255)
                                                                     gATTGCCGATATCGATACGCGTAAGCTGACGCGTTTACTGCGCGAGAAAGGCGCACAGAATGGCTTCATTATCGCGGGCGATAACCCGGTTGCGGCGCTGGCGTTAGAAAAAGCCCGCGCGTTCCCAGGTCTGAATGGCATGGATCTGGc                                                                              <  1:695751/150‑1 (MQ=255)
                                                                     gATTGCCGATATCGATACCCGTAAGCTGACGCGTTTACTGCGCGAGAAAGGCGCACAGAATGGCTGCATTATCGCGGGCGATAACCCGGATGCGGCGCTGGCGTTAGAAAAAGCCCGCGCGTTCCCAGGTCTGAATGGCATGGATCTGGc                                                                              <  1:150827/150‑1 (MQ=255)
                                                                     gATTGCCGATATCGATACCCGTAAGCTGACGCGTTTACTGCGCGAGAAAGGCGCACAGAATGGCTGCATTATCGCGGGCGATAACCCGGATGCGGCGCTGGCGTTAGAAAAAGCCCGCGCGTTCCCAGGTCTGAATGGCATGGATCTGGc                                                                              <  1:273643/150‑1 (MQ=255)
                                                                     gATTGCCGATATCGATACCCGTAAGCTGACGCGTTTACTGCGCGAGAAAGGCGCACAGAATGGCTGCATTATCGCGGGCGATAACCCGGATGCGGCGCTGGCGTTAGAAAAAGCCCGCGCGTTCCCAGGTCTGAATGGCATGGATCTGGc                                                                              <  2:1350876/150‑1 (MQ=255)
                                                                       ttGCCGATATCGATACCCGTAAGCTGACGCGTTTACTGCGCGAGAAAGGCGCACAGAATGGCTGCATTATCGCGGGCGATAACCCGGATGCGGCGCTGGCGTTAGAAAAAGCCCGCGCGTTCCCAGGTCTGAATGGCATGGATCTGGCaa                                                                            >  1:184903/1‑150 (MQ=255)
                                                                         gCCGATATCGATACCCGTAAGCTGACGCGTTTACTGCGCGAGAAAGGCGCACAGAATGGCTGCATTATCGCGGGCGATAACCCGGATGCGGCGCTGGCGTTAGAAAAAGCCCGCGCGTTCCCAGGTCTGAATGGCATGGATCTGGCaaaa                                                                          >  1:223091/1‑150 (MQ=255)
                                                                          ccGATATCGATACCCGTAAGCTGACGCGTTTACTGCGCGAGAAAGGCGCACAGAATGGCTGCATTATCGCGGGCGATAACCCGGATGCGGCGCTGGCGTTAGAAAAAGCCCGCGCGTTCCCAGGTCTGAATGGCATGGATCTGGCAAaag                                                                         <  2:1513673/150‑1 (MQ=255)
                                                                           cgttatCGATACCCGTAAGCTGACGCGTTTACTGCGCGAGAAAGTCGCACAGAATGGCTGCATTATCGCGGGCGATAACCCGGATGCGGCGCTTGCGTTAGAAAAAGCCCGCGCGTTCCCAGTTCTGAATGGCATGGATCTGGCAAaaga                                                                        <  2:575314/147‑1 (MQ=255)
                                                                            gATATCGATACCCGTAAGCTGACGCGTTTACTGCGCGAGAAAGGCGCACAGAATGGCTGCATTATCGCGGGCGATAACCCGGATGCGGCGCTGGCGTTAGAAAAAGCCCGCGCGTTCCCAGGTCTGAATGGCATGGATCTGGCAAaagaa                                                                       >  1:300215/1‑150 (MQ=255)
                                                                               atCGATACCCGTAAGCTGACGCGTTTACTGCGCGAGAAAGGCGCACAGAATGGCTGCATTATCGCGGGCGATAACCCGGATGCGGCGCTGGCGTTAGAAAAAGCCCGCGCGTTCCCAGGTCTGAATGGCATGGATCTGGCAAAAGAAGTg                                                                    <  1:1514305/150‑1 (MQ=255)
                                                                                 cGATACCCGTAAGCTGACGCGTTTACTGCGCGAGAAAGGCGCACAGAATGGCTGCATTATCGCGGGCGATAACCCGGATGCGGCGCTGGCGTTAGAAAAAGCCCGCGCGTTCCCAGGTCTGAATGGCATGGATCTGGCAAAAGAAGTGac                                                                  >  2:72496/1‑150 (MQ=255)
                                                                                  gATACCCGTAAGCTGACGCGTTTACTGCGCGAGAAAGGCGCACAGAATGGCTGCATTATCGCGGGCGATAACCCGGATGCGGCGCTGGCGTTAGAAAAAGCCCGCGCGTTCCCAGGTCTGAATGGCATGGATCTGGCAAAAGAAGTGacc                                                                 <  2:440287/150‑1 (MQ=255)
                                                                                      cccGTAAGCTGACGCGTTTACTGCGCGAGAAAGGCGCACAGAATGGCTGCATTATCGCGGGCGATAACCCGGATGCGGCGCTGGCGTTAGAAAAAGCCCGCGCGTTCCCAGGTCTGAATGGCATGGATCTGGCAAAAGAAGTGACCACCg                                                             >  2:25220/1‑150 (MQ=255)
                                                                                      cccGTAAGCTGACGCGTTTACTGCGCGAGAAAGGCGCACAGAATGGCTGCATTATCGCGGGCGATAACCCGGATGCGGCGCTGGCGTTAGAAAAAGCCCGCGCGTTCCCAGGTCTGAATGGCATGGATCTGGCAAAAGAAGTGACCACCg                                                             >  1:762936/1‑150 (MQ=255)
                                                                                       ccGTAAGCTGACGCGTTTACTGCGCGAGAAAGGCGCACAGAATGGCTGCATTATCGCGGGCGATAACCCGGATGCGGCGCTGGCGTTAGAAAAAGCCCGCGCGTTCCCAGGTCTGaa                                                                                             >  2:1359801/1‑117 (MQ=255)
                                                                                       ccGTAAGCTGACGCGTTTACTGCGCGAGAAAGGCGCACAGAATGGCTGCATTATCGCGGGCGATAACCCGGATGCGGCGCTGGCGTTAGAAAAAGCCCGCGCGTTCCCAGGTCTGaa                                                                                             <  1:1359801/117‑1 (MQ=255)
                                                                                            aGCTGACGCGTTTACTGCGCGAGAAAGGCGCACAGAATGGCTGCATTATCGCGGGCGATAACCCGGATGCGGCGCTGGCGTTAGAAAAAGCCCGCGCGTTCCCAGGTCTGAATGGCATGGATCTGGCAAAAGAAGTGACCACCGCAGAAg                                                       >  2:1280526/1‑150 (MQ=255)
                                                                                             gCTGACGCGTTTACTGCGCGAGAAAGGCGCACAGAATGGCTGCATTATCGCGGGCGATAACCCGGATGCGGCGCTGGCGTTAGAAAAAGCCCGCGCGTTCCCAGGTCTGAATGGCATGGATCTGGCAAAAGAAGTGACCACCGCAGAAGc                                                      >  1:542273/1‑150 (MQ=255)
                                                                                              cTGACGCGTTTACTGCGCGAGAAAGGCGCACAGAATGGCTGCATTATCGCGGGCGATAACCCGGATGCGGCGCTGGCGTTAGAAAAAGCCCGCGCGTTCCCAGGTCTGAATGGCATGGATCTGGCAAAAGAAGTGACCACCGCAGAAGcc                                                     <  1:949818/150‑1 (MQ=255)
                                                                                              cTGACGCGTTTACTGCGCGAGAAAGGCGCACAGAATGGCTGCATTATCGCGGGCGATAACCCGGATGCGGCGCTGGCGTTAGAAAAAGCCCGCGCGTTCCCAGGTCTGAATGGCATGGATCTGGCAAAAGAAGTGACCACCGCAGAAGcc                                                     <  1:592483/150‑1 (MQ=255)
                                                                                               tGACGCGTTTACTGCTCGAGAAAGGCGCACATAATGGCTGCATTATCGCGGGCGATAACCCGGATGCGGCGCTGCCGTTAGAAAAAGCCCGCGCGATCCCAGGTCTGAATCCCATCGATCTGCCAAAAGAACTCACCACCGCAGAAccca                                                    >  2:1522595/1‑149 (MQ=255)
                                                                                                     gTTTACTGCGCGAGAAAGGCGCACAGAATGGCTGCATTATCGCGGGCGATAACCCGGATGCGGCGCTGGCGTTAGAAAAAGCCCGCGCGTTCCCAGGTCTGAATGGCATGGATCTGGCAAAAGAAGTGACCACCGCAGAAGCCTATAGCt                                              <  1:775547/150‑1 (MQ=255)
                                                                                                         aCTGCGCGAGAAAGGCGCACAGAATGGCTGCATTATCGCGGGCGATAACCCGGATGCGGCGCTGGCGTTAGAAAAAGCCCGCGCGTTCCCAGGTCTGAATGGCATGGATCTGGCAAAAGAAGTGACCACCGCAGAAGCCTATAGCTGGac                                          <  2:1535975/150‑1 (MQ=255)
                                                                                                          cTGCGCGAGAAAGGCGCACAGAATGGCTGCATTATCGCGGGCGATAACCCGGATGCGGCGCTGGCGTTAGAAAAAGCCCGCGCGTTCCCAGGTCTGAATGGCATGGATCTGGCAAAAGAAGTGACCACCGCAGAAGCCTATAGCTGGaca                                         >  1:1283984/1‑150 (MQ=255)
                                                                                                           tGCGCGAGAAAGGCGCACAGAATGGCTGCATTATCGCGGGCGATAACCCGGATGCGGCGCTGGCGTTAGAAAAAGCCCGCGCGTTCCCAGGTCTGAATGGCATGGATCTGGCAAAAGAAGTGACCACCGCAGAAGCCTATAGCTGGacac                                        >  2:1516913/1‑150 (MQ=255)
                                                                                                            gcgcGAGAAAGGCGCACAGAATGGCTGCATTATCGCGGGCGATAACCCGGATGCGGCGCTGGCGTTAGAAAAAGCCCGCGCGTTCCCAGGTCTGAATGGCATGGATCTGGCAAAAGAAGTGACCACCGCAGAAGCCTATAGCTGGAcaca                                       >  1:665245/1‑150 (MQ=255)
                                                                                                              gcgAGAAAGGCGCACAGAATGGCTGCATTATCGCGGGCGATAACCCGGATGCGGCGCTGGCGTTAGAAAAAGCCCGCGCGTTCCCAGGTCTGAATGGCATGGATCTGGCAAAAGAAGTGACCACCGCAGAAGCCTATAGCTGGACACAAg                                     <  1:1597133/150‑1 (MQ=255)
                                                                                                              gcgAGAAAGGCGCACAGAATGGCTGCATTATCGCGGGCGATAACCAGGATGCGGCGCTGGCGTTAGAAAAAGCCCGCGCGTTCGCAGGTCTGAATGGCATGGATCTGGCAAAATAAGTGACCACCTCAGAAGCCTATAGCTGGACACAAg                                     <  2:494914/150‑1 (MQ=255)
                                                                                                                gagaAAGGCGCACAGAATGGCTGCATTATCGCGGGCGATAACCCGGATGCGGCGCTGGCGTTAGAAAAAGCCCGCGCGTTCCCAGGTCTGAATGGCATGGATCTGGCAAAAGAAGTGACCACCGCAGAAGCCTATAGCTGGACACAAggg                                   >  1:571964/1‑150 (MQ=255)
                                                                                                                 agaAAGGCGCACAGAATGGCTGCATTATCGCGGGCGATAACCCGGATGCGGCGCTGGCGTTAGAAAAAGCCCGCGCGTTCCCAGGTCTGAATGGCATGGATCTGGCAAAAGAAGTGACCACCGCAGAAGCCTATAGCTGGACACAAGGGa                                  >  1:592998/1‑150 (MQ=255)
                                                                                                                     aGGCGCACAGAATGGCTGCATTATCGCGGGCGATAACCCGGATGCGGCGCTGGCGTTAGAAAAAGCCCGCGCGTTCCCAGGTCTGAATGGCATGGATCTGGCAAAAGAAGTGACCACCGCAGAAGCCTATAGCTGGACACAAGGGAGCTg                              <  2:1078104/150‑1 (MQ=255)
                                                                                                                       gcgcACAGAATGGCTGCATTATCGCGGGCGATAACCCGGATGCGGCGCTGGCGTTAGAAAAAGCCCGCGCGTTCCCAGGTCTGAATGGCATGGATCTGGCAAAAGAAGTGACCACCGCAGAAGCCTATAGCTGGACACAAGGGAGCTGGa                            >  1:7636/1‑150 (MQ=255)
                                                                                                                              gAATGGCTGCATTATCGCGGGCGATAACCCGGATGCGGCGCTGGCGTTAGAAAAAGCCCGCGCGTTCCCAGGTCTGAATGGCATGGATCTGGCAAAAGAAGTGACCACCGCAGAAGCCTATAGCTGGACACAAGGGAGCTGGACGTTGAc                     >  1:265450/1‑150 (MQ=255)
                                                                                                                              gAATGGCTGCATTATCGCGGGCGATAACCCGGATGCGGCGCTGGCGTTAGAAAAAGCCCGCGCGTTCCCAGGTCTGAATGGCATGGATCTGGCAAAAGAAGTGACCACCGCAGAAGCCTATAGCTGGACACAAGGGAGCTGGACGTTGAc                     >  1:1218519/1‑150 (MQ=255)
                                                                                                                              gAATGGCTGCATTATCGCGGGCGATAACCCGGATGCGGCGCTGGCGTTAGAAAAAGCCCGCGCGTTCCCAGGTCTGAATGGCATGGATCTGGCAAAAGAAGTGACCACCGCAGAAGCCTATAGCTGGACACAAGGGAGCTGGACGTTGAc                     >  1:1301504/1‑150 (MQ=255)
                                                                                                                               aaTGGCTGCATTATCGCGGGCGATAACCCGGATGCGGCGCTGGCGTTAGAAAAAGCCCGCGCGTTCCCAGGTCTGAATGGCATGGATCTGGCAAAAGAAGTGACCACCGCAGAAGCCTATAGCTGGACACAAGGGAGCTGGACGTTGAcc                    <  1:1522595/150‑1 (MQ=255)
                                                                                                                                  ggCTGCATTATCGCGGGCGATAACCCGGATGCGGCGCTGGCGTTAGAAAAAGCCCGCGCGTTCCCAGGTCTGAATGGCATGGATCTGGCAAAAGAAGTGACCACCGCAGAAGCCTATAGCTGGACACAAGGGAGCTGGACGTTGACCGGt                 <  1:1346485/150‑1 (MQ=255)
                                                                                                                                    ccgggttATCGCCCGCGATAACCCGGATGCGGCGCTGGCGTTAGAAAAAGCCCGCGCGTTCCCAGGTCTGAATGGCATGGATCTGGCAAAAGAAGTGACCACCGCAGAAGCCTATAGCTGGACACAAGGGAGCTGGACGTTGACCGGTgg               >  1:1212212/6‑150 (MQ=255)
                                                                                                                                    ccgggttATCGCCCGCGATAACCCGGATGCGGCGCTGGCGTTAGAAAAAGCCCGCGCGTTCCCAGGTCTGAATGGCATGGATCTGGCAAAAGAAGTGACCACCGCAGAAGCCTATAGCTGGACACAAGGGAGCTGGACGTTGACCGGTgg               >  2:94786/6‑150 (MQ=255)
                                                                                                                                     tGCATTATCGCGGGCGATAACCCGGATGCGGCGCTGGCGTTAGAAAAAGCCCGCGCGTTCCCAGGTCTGAATGGCATGGATCTGGCAAAAGAAGTGACCACCGCAGAAGCCTATAGCTGGACACAAGGGAGCTGGACGTTGACCGGTGGc              <  2:542273/150‑1 (MQ=255)
                                                                                                                                      gggttATCGCCCGCGATAACCCGGATGCGGCGCTGGCGTTAGAAAAAGCCCGCGCGTTCCCAGGTCTGAATGGCATGGATCTGGCAAAAGAAGTGACCACCGCAGAAGCCTATAGCTGGACACAAGGGAGCTGGACGTTGACCGGTGGcc             >  1:1163082/4‑150 (MQ=255)
                                                                                                                                      gCATTATCGCGGGCGATAACCCGGATGCGGCGCTGGCGTTAGAAAAAGCCCGCGCGTTCCCAGGTCTGAATGGCATGGATCTGGCAAAAGAAGTGACCACCGCAGAACCCTATAGCTGGACACAAGGGAGCTGGACGTTGACCGGTGGcc             >  2:113470/1‑150 (MQ=255)
                                                                                                                                         ttATCGCGGGCGATAACCCGGATGCGGCGCTGGCGTTAGAAAAAGCCCGCGCGTTCCCAGGTCTGAATGGCATGGATCTGGCAAAAGAAGTGACCACCGCAGAAGCCTATAGCTGGACACAAGGGAGCTGGACGTTGACCGGTGGCCTGc          >  2:1179349/1‑150 (MQ=255)
                                                                                                                                             cgcgGGCGATAACCCGGATGCGGCGCTGGCGTTAGAAAAAGCCCGCGCGTTCCCAGGTCTGAATGGCATGGATCTGGCAAAAGAAGTGACCACCGCAGAAGCCTATAGCTGGACACAAGGGAGCTGGACGTTGACCGGTGGCCTGCCAGa      >  2:689606/1‑150 (MQ=255)
                                                                                                                                              gcgGGCGATAACCCGGATGCGGCGCTGGCGTTAGAAAAAGCCCGCGCGTTCCCAGGTCTGAATGGCATGGATCTGGCAAAAGAAGTGACCACCGCAGAAGCCTATAGCTGGACACAAGGGAGCTGGACGTTGACCGGTGGCCTGCCAGaa     >  1:326645/1‑150 (MQ=255)
                                                                                                                                                 ggCGATAACCCGGATGCGGCGCTGGCGTTAGAAAAAGCCCGCGCGTTCCCAGGTCTGAATGGCATGGATCTGGCAAAAGAAGTGACCACCGCAGAAGCCTATAGCTGGACACAAGGGAGCTGGACGTTGACCGGTGGCCTGCCAGAAGCg  >  1:1169116/1‑150 (MQ=255)
                                                                                                                                                 |                                                                                                                                                     
GCTGATTGCCAGCAACTTCCGTAATACCGAAGACCTCTCTTCTTACCTGAAACGCCATAACATCGTGGCGATTGCCGATATCGATACCCGTAAGCTGACGCGTTTACTGCGCGAGAAAGGCGCACAGAATGGCTGCATTATCGCGGGCGATAACCCGGATGCGGCGCTGGCGTTAGAAAAAGCCCGCGCGTTCCCAGGTCTGAATGGCATGGATCTGGCAAAAGAAGTGACCACCGCAGAAGCCTATAGCTGGACACAAGGGAGCTGGACGTTGACCGGTGGCCTGCCAGAAGCG  >  NC_000913/29908‑30202

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: