Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 30,052 G→C 5.9% A134A (GCG→GCC carA → carbamoyl phosphate synthetase subunit alpha

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00091330,0520GC 5.9% 291.1 / 7.1 85A134A (GCG→GCCcarAcarbamoyl phosphate synthetase subunit alpha
Reads supporting (aligned to +/- strand):  ref base G (42/38);  new base C (3/2);  total (45/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.99e-01

CTGCCGCTGATTGCCAGCAACTTCCGTAATACCGAAGACCTCTCTTCTTACCTGAAACGCCATAACATCGTGGCGATTGCCGATATCGATACCCGTAAGCTGACGCGTTTACTGCGCGAGAAAGGCGCACAGAATGGCTGCATTATCGCGGGCGATAACCCGGATGCGGCGCTGGCGTTAGAAAAAGCCCGCGCGTTCCCAGGTCTGAATGGCATGGATCTGGCAAAAGAAGTGACCACCGCAGAAGCCTATAGCTGGACACAAGGGAGCTGGACGTTGACCGGTGGCCTGCCAGAA  >  NC_000913/29903‑30199
                                                                                                                                                     |                                                                                                                                                   
cTGCCGCTGATTGCCAGCAACTTCCGTAATACCGAAGACCTCTCTTCTTACCTGAAACGCCATAACATCGTGGCGATTGCCGATATCGATACCCGTAAGCTGACGCGTTTACTGCGCGAGAAAGGCGCACAGAATGGCTGCATTATcgcg                                                                                                                                                     <  1:937177/150‑1 (MQ=255)
     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:903988/1‑150 (MQ=255)
            gCCAGCAACTTCCGTAATACCGAAGACCTCTCTTCTTACCTGAAACGCCATAACATCGTGGCGATTGCCGATATCGATACCCGTAAGCTGACGCGTTTACTGCGCGAGAAAGGCGCACAGAATGGCTGCATTATCGCGGGCGATAACCCg                                                                                                                                         >  1:910023/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:1086057/1‑150 (MQ=255)
                                      ccTCTCTTCTTACCTGAAACGCCATAACATCGTGGCGATTGCCGATATCGATACCCGTAAGCTGACGCGTTTACTGCGCGAGAAAGGCGCACAGAATGGCTGCATTATCGCGGGCGATAACCCGGATGCGGCGCTGGCGTTAGAAAAAGc                                                                                                               >  1:712135/1‑150 (MQ=255)
                                       ctctctTCTTACCTGAAACGCCATAACATCGTGGCGATTGCCGATATCGATACCCGTAAGCTGACGCGTTTACTGCGCGAGAAAGGCGCACAGAATGGCTGCATTATCGCGGGCGATAACCAGGATGCGGCGCTGGCGTTAGAAAAAGcc                                                                                                              >  1:494914/1‑150 (MQ=255)
                                         ctctTCTTACCTGAAACGCCATAACATCGTGGCGATTGCCGATATCGATACCCGTAAGCTGACGCGTTTACTGCGCGAGAAAGGCGCACAGAATGGCTGCATTATCGCGGGCGATAACCCGGATGCGGCGCTGGCGTTAGAAAAAGCCcg                                                                                                            <  1:1659047/150‑1 (MQ=255)
                                         ctctTCTTACCTGAAACGCCATAACATCGTGGCGATTGCCGATATCGATACCCGTAAGCTGACGCGTTTACTGCGCGAGAAAGGCGCACAGAATGGCTGCATTATCGCGGGCGATAACCCGGATGCGGCGCTGGCGTTAGAAAAAGCCcg                                                                                                            >  2:649865/1‑150 (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                                                          >  1:762936/1‑150 (MQ=255)
                                                                                           cccGTAAGCTGACGCGTTTACTGCGCGAGAAAGGCGCACAGAATGGCTGCATTATCGCGGGCGATAACCCGGATGCGGCGCTGGCGTTAGAAAAAGCCCGCGCGTTCCCAGGTCTGAATGGCATGGATCTGGCAAAAGAAGTGACCACCg                                                          >  2:25220/1‑150 (MQ=255)
                                                                                            ccGTAAGCTGACGCGTTTACTGCGCGAGAAAGGCGCACAGAATGGCTGCATTATCGCGGGCGATAACCCGGATGCGGCGCTGGCGTTAGAAAAAGCCCGCGCGTTCCCAGGTCTGaa                                                                                          <  1:1359801/117‑1 (MQ=255)
                                                                                            ccGTAAGCTGACGCGTTTACTGCGCGAGAAAGGCGCACAGAATGGCTGCATTATCGCGGGCGATAACCCGGATGCGGCGCTGGCGTTAGAAAAAGCCCGCGCGTTCCCAGGTCTGaa                                                                                          >  2:1359801/1‑117 (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:1301504/1‑150 (MQ=255)
                                                                                                                                   gAATGGCTGCATTATCGCGGGCGATAACCCGGATGCGGCGCTGGCGTTAGAAAAAGCCCGCGCGTTCCCAGGTCTGAATGGCATGGATCTGGCAAAAGAAGTGACCACCGCAGAAGCCTATAGCTGGACACAAGGGAGCTGGACGTTGAc                  >  1:1218519/1‑150 (MQ=255)
                                                                                                                                    aaTGGCTGCATTATCGCGGGCGATAACCCGGATGCGGCGCTGGCGTTAGAAAAAGCCCGCGCGTTCCCAGGTCTGAATGGCATGGATCTGGCAAAAGAAGTGACCACCGCAGAAGCCTATAGCTGGACACAAGGGAGCTGGACGTTGAcc                 <  1:1522595/150‑1 (MQ=255)
                                                                                                                                       ggCTGCATTATCGCGGGCGATAACCCGGATGCGGCGCTGGCGTTAGAAAAAGCCCGCGCGTTCCCAGGTCTGAATGGCATGGATCTGGCAAAAGAAGTGACCACCGCAGAAGCCTATAGCTGGACACAAGGGAGCTGGACGTTGACCGGt              <  1:1346485/150‑1 (MQ=255)
                                                                                                                                         ccgggttATCGCCCGCGATAACCCGGATGCGGCGCTGGCGTTAGAAAAAGCCCGCGCGTTCCCAGGTCTGAATGGCATGGATCTGGCAAAAGAAGTGACCACCGCAGAAGCCTATAGCTGGACACAAGGGAGCTGGACGTTGACCGGTgg            >  2:94786/6‑150 (MQ=255)
                                                                                                                                         ccgggttATCGCCCGCGATAACCCGGATGCGGCGCTGGCGTTAGAAAAAGCCCGCGCGTTCCCAGGTCTGAATGGCATGGATCTGGCAAAAGAAGTGACCACCGCAGAAGCCTATAGCTGGACACAAGGGAGCTGGACGTTGACCGGTgg            >  1:1212212/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)
                                                                                                                                                     |                                                                                                                                                   
CTGCCGCTGATTGCCAGCAACTTCCGTAATACCGAAGACCTCTCTTCTTACCTGAAACGCCATAACATCGTGGCGATTGCCGATATCGATACCCGTAAGCTGACGCGTTTACTGCGCGAGAAAGGCGCACAGAATGGCTGCATTATCGCGGGCGATAACCCGGATGCGGCGCTGGCGTTAGAAAAAGCCCGCGCGTTCCCAGGTCTGAATGGCATGGATCTGGCAAAAGAAGTGACCACCGCAGAAGCCTATAGCTGGACACAAGGGAGCTGGACGTTGACCGGTGGCCTGCCAGAA  >  NC_000913/29903‑30199

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: