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

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00091330,0530GC 5.1% 284.9 / 5.8 78G135R (GGC→CGC) carAcarbamoyl phosphate synthetase small subunit, glutamine amidotransferase
Reads supporting (aligned to +/- strand):  ref base G (37/37);  new base C (2/2);  total (39/39)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.38e-01

CCAGCAACTTCCGTAATACCGAAGACCTCTCTTCTTACCTGAAACGCCATAACATCGTGGCGATTGCCGATATCGATACCCGTAAGCTGACGCGTTTACTGCGCGAGAAAGGCGCACAGAATGGCTGCATTATCGCGGGCGATAACCCGGATGCGGCGCTGGCGTTAGAAAAAGCCCGCGCGTTCCCAGGTCTGAATGGCATGGATCTGGCAAAAGAAGTGACCACCGCAGAAGCCTATAGCTGGACACAAGGGAGCTGGACGTTGACCGGTGGCCTGCCAGA  >  NC_000913/29916‑30198
                                                                                                                                         |                                                                                                                                                 
ccAGCAACTTCCGTAATACCGAAGACCTCTCTTCTTACCTGAAACGCCATAACATCGTGGCGATTGCCGATATCGATACCCGTAAGCTGACGCGTTTACTGCGCGAGAAAGGCGCACAGAATGGCTGCATTATCGCGGGCGATAAccc                                                                                                                                         <  3:955738/148‑1 (MQ=255)
    cAACTTCCGTAATACCGAAGACCTCTCTTCTTACCTGAAACGCCATAACATCGTGGCGATTGCCGATATCGATACCCGTAAGCTGACGCGTTTACTGCGCGAGAAAGGCGCACAGAATGGCTGCATTATCGCGGGCGATAAccc                                                                                                                                         >  2:985732/1‑144 (MQ=255)
    cAACTTCCGTAATACCGAAGACCTCTCTTCTTACCTGAAACGCCATAACATCGTGGCGATTGCCGATATCGATACCCGTAAGCTGACGCGTTTACTGCGCGAGAAAGGCGCACAGAATGGCTGCATTATCGCGGGCGATAAccc                                                                                                                                         <  1:985732/144‑1 (MQ=255)
      aCTTCCGTAATACCGAAGACCTCTCTTCTTACCTGAAACGCCATAACATCGTGGCGATTGCCGATATCGATACCCGTAAGCTGACGCGTTTACTGCGCGAGAAAGGCGCACAGAATGGCTGCATTATCGCGGGCGATAACCCGGATgc                                                                                                                                   >  4:583547/1‑148 (MQ=255)
       cTTCCGTAATACCGAAGACCTCTCTTCTTACCTGAAACGCCATAACATCGTGGCGATTGCCGATATCGATACCCGTAAGCTGACGCGTTTACTGCGCGAGAAAGGCGCACAGAATGGCTGCATTATCGCGGGCGATAACCCGGATgcg                                                                                                                                  <  3:583547/148‑1 (MQ=255)
          ccGTAATACCGAAGACCTCTCTTCTTACCTGAAACGCCATAACATCGTGGCGATTGCCGATATCGATACCCGTAAGCTGACGCGTTTACTGCGCGAGAAAGGCGCACAGAATGGCTGCATTATCGCGGGCGATAACCCGGATgcggcg                                                                                                                               >  3:646104/1‑148 (MQ=255)
                        aCCTCTCTTCTTACCTGAAACGCCATAACATCGTGGCGATTGCCGATATCGATACCCGTAAGCTGACGCGTTTACTGCGCGAGAAAGGCGCACAGAATGGCTGCATTATCGCGGGCGATAACCCGGATGCGGCGCTGGCGTTAGaaaa                                                                                                                 <  4:646104/148‑1 (MQ=255)
                           tctctTCTTACCTGAAACGCCATAACATCGTGGCGATTGCCGATATCGATACCCGTAAGCTGACGCGTTTACTGCGCGAGAAAGGCGCACAGAATGGCTGCATTATCGCGGGCGATa                                                                                                                                             <  3:515371/117‑1 (MQ=255)
                           tctctTCTTACCTGAAACGCCATAACATCGTGGCGATTGCCGATATCGATACCCGTAAGCTGACGCGTTTACTGCGCGAGAAAGGCGCACAGAATGGCTGCATTATCGCGGGCGATa                                                                                                                                             >  4:515371/1‑117 (MQ=255)
                             tcttctTACCTGAAACGCCATAACATCGTGGCGATTGCCGATATCGATACCCGTAAGCTGACGCGTTTACTGCGCGAGAAAGGCGCACAGAATGGCTGCATTATCGCGGGCGATAACCCGGATGCGGCGCTGGCGTTAGaaaa                                                                                                                 <  4:1131823/143‑1 (MQ=255)
                             tcttctTACCTGAAACGCCATAACATCGTGGCGATTGCCGATATCGATACCCGTAAGCTGACGCGTTTACTGCGCGAGAAAGGCGCACAGAATGGCTGCATTATCGCGGGCGATAACCCGGATGCGGCGCTGGCGTTAGaaaa                                                                                                                 >  3:1131823/1‑143 (MQ=255)
                             tcttctTACCTGAAACGCCATAACATCGTGGCGATTGCCGATATCGATACCCGTAAGCTGACGCGTTACTGCGCGAGAAAGGCGCACAGAATGGCTGCATTATCGCGGGCGATAACCCGGATGCGGCGCTGGCGTTAGAAAAAGcc                                                                                                             <  1:860493/146‑1 (MQ=255)
                             tcttctTACCTGAAACGCCATAACATCGTGGCGATTGCCGATATCGATACCCGTAAGCTGACGCGTTACTGCGCGAGAAAGGCGCACAGAATGGCTGCATTATCGCGGGCGATAACCCGGATGCGGCGCTGGCGTTAGAAAAAGcc                                                                                                             >  2:860493/1‑146 (MQ=255)
                                  ttACCTGAAACGCCATAACATCGTGGCGATTGCCGATATCGATACCCGTAAGCTGACGCGTTTACTGCGCGAGAAAGGCGCACAGAATGGCTGCATTATCGCGGGCGATAACCCGGATGCGGCGCTGGCGTTAGaaa                                                                                                                  <  4:576761/137‑1 (MQ=255)
                                  ttACCTGAAACGCCATAACATCGTGGCGATTGCCGATATCGATACCCGTAAGCTGACGCGTTTACTGCGCGAGAAAGGCGCACAGAATGGCTGCATTATCGCGGGCGATAACCCGGATGCGGCGCTGGCGTTAGaaa                                                                                                                  >  3:576761/1‑137 (MQ=255)
                                       tGAAACGCCATAACATCGTGGCGATTGCCGATATCGATACCCGTAAGCTGACGCGTTTACTGCGCGAGAAAGGCGCACAGAATGGCTGCATTATCGCGGGCGATAACCCGGATgcg                                                                                                                                  >  3:586878/1‑116 (MQ=255)
                                       tGAAACGCCATAACATCGTGGCGATTGCCGATATCGATACCCGTAAGCTGACGCGTTTACTGCGCGAGAAAGGCGCACAGAATGGCTGCATTATCGCGGGCGATAACCCGGATgcg                                                                                                                                  <  4:586878/116‑1 (MQ=255)
                                        gAAACGCCATAACATCGTGGCGATTGCCGATATCGATACCCGTAAGCTGACGCGTTTACTGCGCGAGAAAGGCGCACAGAATGGCTGCATTATCGCGGGCGATAACCCGGATGCGGCGCTGGCGTTAGAAAAAGCCCGCGCGTTCCCa                                                                                                 >  4:634095/1‑148 (MQ=255)
                                        gAAACGCCATAACATCGTGGCGATTGCCGATATCGATACCCGTAAGCTGACGCGTTTACTGCGCGAGAAAGGCGCACAGAATGGCTGCATTATCGCGGGCGATAACCCGGATGCGGCGCTGGCGTTAGAAAAAGCCCGCGCGTTCCCa                                                                                                 >  4:692792/1‑148 (MQ=255)
                                                aTAACATCGTGGCGATTGCCGATATCGATACCCGTAAGCTGACGCGTTTACTGCGCGAGAAAGGCGCACAGAATGGCTGCATTATCGCGGGCGATAAcc                                                                                                                                          <  4:398239/99‑1 (MQ=255)
                                                aTAACATCGTGGCGATTGCCGATATCGATACCCGTAAGCTGACGCGTTTACTGCGCGAGAAAGGCGCACAGAATGGCTGCATTATCGCGGGCGATAAcc                                                                                                                                          >  3:398239/1‑99 (MQ=255)
                                                             gATTGCCGATATCGATACCCGTAAGCTGACGCGTTTACTGCGCGAGAAAGGCGCACAGAATGGCTGCATTATCGCGGGCGATAACCCGGATGCGGCGCTGGCGTTAGAAAAAGCCCGCGCGTTCCCAGGTCTGAATGGCATGGATCTg                                                                            <  3:692792/148‑1 (MQ=255)
                                                              aTTGCCGATATCGATACCCGTAAGCTGACGCGTTTACTGCGCGAGAAAGGCGCACAGAATGGCTGCATTATCGCGGGCGATAACCCGGATGCGGCGCTgg                                                                                                                           >  3:829307/1‑100 (MQ=255)
                                                              aTTGCCGATATCGATACCCGTAAGCTGACGCGTTTACTGCGCGAGAAAGGCGCACAGAATGGCTGCATTATCGCGGGCGATAACCCGGATGCGGCGCTgg                                                                                                                           <  4:829307/100‑1 (MQ=255)
                                                                         cGATACCCGTAAGCTGACGCGTTTACTGCGCGAGAAAGGCGCACAGAATGGCTGCATTATCGCGGGCGATAACCCGGATGCGGCGCTGGCGTTAGAAAAAGCCCGCGCGTTCCCAGGTCTGAATGGCATGGATCTGGCAAAAGAAGTg                                                                >  2:742160/1‑148 (MQ=255)
                                                                                cGTAAGCTGACGCGTTTACTGCGCGAGAAAGGCGCACAGAATGGCTGCATTATCGCGGGCGATAACCCGGa                                                                                                                                      <  3:466728/71‑1 (MQ=255)
                                                                                cGTAAGCTGACGCGTTTACTGCGCGAGAAAGGCGCACAGAATGGCTGCATTATCGCGGGCGATAACCCGGa                                                                                                                                      >  4:466728/1‑71 (MQ=255)
                                                                                cGTAAGCTGACGCGTTTACTGCGCGAGAAAGGCGCACAGAATGGCTGCATTATCGCGGGCGATAACCCGGATGCGGCGCTGGCGTTAGAAAAAGCCCGCGCGTTCCCAGGTCTGAATGGCATGGATCTGGCaaaa                                                                      <  4:755275/135‑1 (MQ=255)
                                                                                cGTAAGCTGACGCGTTTACTGCGCGAGAAAGGCGCACAGAATGGCTGCATTATCGCGGGCGATAACCCGGATGCGGCGCTGGCGTTAGAAAAAGCCCGCGCGTTCCCAGGTCTGAATGGCATGGATCTGGCaaaa                                                                      <  2:1085550/135‑1 (MQ=255)
                                                                                cGTAAGCTGACGCGTTTACTGCGCGAGAAAGGCGCACAGAATGGCTGCATTATCGCGGGCGATAACCCGGATGCGGCGCTGGCGTTAGAAAAAGCCCGCGCGTTCCCAGGTCTGAATGGCATGGATCTGGCaaaa                                                                      >  3:755275/1‑135 (MQ=255)
                                                                                cGTAAGCTGACGCGTTTACTGCGCGAGAAAGGCGCACAGAATGGCTGCATTATCGCGGGCGATAACCCGGATGCGGCGCTGGCGTTAGAAAAAGCCCGCGCGTTCCCAGGTCTGAATGGCATGGATCTGGCaaaa                                                                      >  1:1085550/1‑135 (MQ=255)
                                                                                 gTAAGCTGACGCGTTTACTGCGCGAGAAAGGCGCACAGAATGGCTGCATTATCGCGGGCGATAACCCGGATGCGGCGCTGGCGTTAGAAAAAGCCCGCGCGTTCCCAGGTCTGAATGGCATGGATCTGGCAAAAGAAGTGACCACCGc                                                        >  1:89845/1‑148 (MQ=255)
                                                                                  tAAGCTGACGCGTTTACTGCGCGAGAAAGGCGCACAGAATGGCTGCATTATCGCGGGCGATAACCCGGATGCGGCGCTGGCGTTAGAAAAAGCCCGCGCGTTCCCAGGTCTGAATgg                                                                                      >  3:339408/1‑117 (MQ=255)
                                                                                  tAAGCTGACGCGTTTACTGCGCGAGAAAGGCGCACAGAATGGCTGCATTATCGCGGGCGATAACCCGGATGCGGCGCTGGCGTTAGAAAAAGCCCGCGCGTTCCCAGGTCTGAATgg                                                                                      <  4:339408/117‑1 (MQ=255)
                                                                                   aaGCTGACGCGTTTACTGCGCGAGAAAGGCGCACAGAATGGCTGCATTATCGCGGGCGATAACCCGGATGCGGCGCTGGCGTTAGAAAAAGCCCGCGCGTTCCCAGGTCTGAATGGCATGGATCTGGCAAaagaag                                                                  >  3:503984/1‑136 (MQ=255)
                                                                                   aaGCTGACGCGTTTACTGCGCGAGAAAGGCGCACAGAATGGCTGCATTATCGCGGGCGATAACCCGGATGCGGCGCTGGCGTTAGAAAAAGCCCGCGCGTTCCCAGGTCTGAATGGCATGGATCTGGCAAaagaag                                                                  <  4:503984/136‑1 (MQ=255)
                                                                                          cgcgTTTACTGCGCGAGAAAGGCGCACAGAATGGCTGCATTATCGCGGGCGATAACCCGGATGCGGCGCTGGCGt                                                                                                                        >  1:924226/1‑75 (MQ=255)
                                                                                          cgcgTTTACTGCGCGAGAAAGGCGCACAGAATGGCTGCATTATCGCGGGCGATAACCCGGATGCGGCGCTGGCGt                                                                                                                        <  2:924226/75‑1 (MQ=255)
                                                                                            cgTTTACTGCGCGAGAAAGGCGCACAGAATGGCTGCATTATCGCGGGCGATAACCCGGATGCGGCGCTGGCGTTAGAAAAAGCCCGCGCGTTc                                                                                                    >  1:174071/1‑93 (MQ=255)
                                                                                            cgTTTACTGCGCGAGAAAGGCGCACAGAATGGCTGCATTATCGCGGGCGATAACCCGGATGCGGCGCTGGCGTTAGAAAAAGCCCGCGCGTTc                                                                                                    <  2:174071/93‑1 (MQ=255)
                                                                                              tttACTGCGCGAGAAAGGCGCACAGAATGGCTGCATTATCGCGGGCGATAACCCGGATGCGGCGCTGGCGTTAGAAAAAGCCCGCGCGTTCCCAGGTCTGAATGGCATGGATCTGGCAAAAGAAGTGACCACCGCAGAAGCCTATAGc                                           >  1:1109923/1‑148 (MQ=255)
                                                                                               ttACTGCGCGAGAAAGGCGCACAGAATGGCTGCATTATCGCGGGCGATAACCCGGATGCGGCGCTGGCg                                                                                                                         <  1:848109/69‑1 (MQ=255)
                                                                                               ttACTGCGCGAGAAAGGCGCACAGAATGGCTGCATTATCGCGGGCGATAACCCGGATGCGGCGCTGGCg                                                                                                                         >  2:848109/1‑69 (MQ=255)
                                                                                                  cTGCGCGAGAAAGGCGCACAGAATGGCTGCATTATCGCGGGCGATAACCCGGATGCGGCGCTGGCGTTAGAAAAAGCCCGCGCGTTCCCAGGTCTGAATGGCATGGATCTGGCAAAAGAAGTGACCACCGCAGAAGCCTATAGCt                                          >  3:240759/1‑145 (MQ=255)
                                                                                                  cTGCGCGAGAAAGGCGCACAGAATGGCTGCATTATCGCGGGCGATAACCCGGATGCGGCGCTGGCGTTAGAAAAAGCCCGCGCGTTCCCAGGTCTGAATGGCATGGATCTGGCAAAAGAAGTGACCACCGCAGAAGCCTATAGCt                                          <  4:240759/145‑1 (MQ=255)
                                                                                                      gcgAGAAAGGCGCACAGAATGGCTGCATTATCGCGGGCGATAACCCGGATGCGGCGCTGGCGTTAGAAAAAGCCCGCGCGTTCCCAGGTCTGAATGGCATGGATCTGGCAAAAGAAGTGACCACCGCAGAAGCCTATAGCTGGAcaca                                   <  3:634095/148‑1 (MQ=255)
                                                                                                         agaAAGGCGCACAGAATGGCTGCATTATCGCGGGCGATAACCCGGATGCGGCGCTGGCGTTAGAAAAAGCCCGCGCGTTCCCAGGTCTGAATGGCATGGATCTGGCAAAAGAAGTGACCACCGCAGAAGCCTATAGCTGGACACAAgg                                >  2:1026964/1‑148 (MQ=255)
                                                                                                         agaAAGGCGCACAGAATGGCTGCATTATCGCGGGCGATAACCCGGATGCGGCGCTGGCGTTAGAAAAAGCCCGCGCGTTCCCAGGTCTGAATGGCATGGATCTGGCAAAAGAAGTGACCACCGCAGAAGCCTATAGCTGGACACAAgg                                >  4:465937/1‑148 (MQ=255)
                                                                                                          gaAAGGCGCACAGAATGGCTGCATTATCGCGGGCGATAACCCGGATGCGGCGCTGGCGTTAGaaaa                                                                                                                 <  1:627331/66‑1 (MQ=255)
                                                                                                          gaAAGGCGCACAGAATGGCTGCATTATCGCGGGCGATAACCCGGATGCGGCGCTGGCGTTAGaaaa                                                                                                                 >  2:627331/1‑66 (MQ=255)
                                                                                                            aaGGCGCACAGAATGGCTGCATTATCGCGGGCGATAACCCGGATGCGGCGCTGGCGTTAGAAAAAGCCCGCGCGTTCCCAGGTCTGAATGGCATGGATCTGGCAAAAGAAGTGACCACCGCAGAAGCCTATAGCTGGACACAAGGGAg                             >  3:1086621/1‑148 (MQ=255)
                                                                                                             aGGCGCACAGAATGGCTGCATTATCGCGGGCGATAACCCGGATGCGGCGCTGGCGTTAGAAAAAGCCCGCGCGTTCCCAGGTCTGAATGGCATGGATCTGGCAAAAGAAGTGACCACCGCAGAAGCCTATAGCTGGACACAAGGGAGc                            <  4:1086621/148‑1 (MQ=255)
                                                                                                              ggCGCACAGAATGGCTGCATTATCGCGGGCGATAACCCGGATGCGGCGCTGGCGTTAGAAAAAGCCCGCGCGTTCCCAGGTCTGAATGGCATGGATCTGGCAAAAGAAGTGACCACCGCAGAAGCCTATAGCTGGACACAAGGGAGCt                           >  2:445770/1‑148 (MQ=255)
                                                                                                                cgcACAGAATGGCTGCATTATCGCGGGCGATAACCCGGATGCGGCGCTGGCGTTAGAAAAAGCCCGCGCGTTCCCa                                                                                                 >  4:552410/1‑76 (MQ=255)
                                                                                                                cgcACAGAATGGCTGCATTATCGCGGGCGATAACCCGGATGCGGCGCTGGCGTTAGAAAAAGCCCGCGCGTTCCCa                                                                                                 <  3:552410/76‑1 (MQ=255)
                                                                                                                 gcACAGAATGGCTGCATTATCGCGGGCGATAACCCGGATGCGGCGCTGGCGTTAGAAAAAGCCCGCGCGTTCCCAGGTCTGAATGGCATGGATCTGGCAAAAGAAGTGACCACCGCAGAAGCCTATAGCTGGACACAAGGGAGCTGGa                        <  1:742160/148‑1 (MQ=255)
                                                                                                                  cacaGAATGGCTGCATTATCGCGGGCGATAACCCGGATGCGGCGCTGGCGTTAGAAAAAGCCCGCGCGTTCCCAGGTCTGAATGGCATGGATCTGGCAAaagaag                                                                  <  3:1101410/105‑1 (MQ=255)
                                                                                                                  cacaGAATGGCTGCATTATCGCGGGCGATAACCCGGATGCGGCGCTGGCGTTAGAAAAAGCCCGCGCGTTCCCAGGTCTGAATGGCATGGATCTGGCAAaagaag                                                                  >  4:1101410/1‑105 (MQ=255)
                                                                                                                    caGAATGGCTGCATTATCGCGGGCGATAACCCGGATGCGGCGCTGGCGTTAGAAAAAGCCCGCGCGTTCCCAGGTCTGAATGGCATGGATCTGGCAAAAGAAGTGACCACCGc                                                        <  1:1077920/113‑1 (MQ=255)
                                                                                                                    caGAATGGCTGCATTATCGCGGGCGATAACCCGGATGCGGCGCTGGCGTTAGAAAAAGCCCGCGCGTTCCCAGGTCTGAATGGCATGGATCTGGCAAAAGAAGTGACCACCGc                                                        >  2:1077920/1‑113 (MQ=255)
                                                                                                                           gCTGCATTATCGCGGGCGATAACCCGGATGCGGCGCTGGCGTTAGAAAAAGCCCGCGCGTTCCCAGGTCTGAATGGCATGGATCTGGCAAAAGAAGTGACCACCGCAGAAGCCTATAGCTGGACACAAGGGAGCTGGACGTTGACCgg              <  3:465937/148‑1 (MQ=255)
                                                                                                                            cTGCATTATCGCGGGCGATAACCCGGATGCGGCGCTGGCGTTAGAAAAAGCCCGCGCGTTCCCAGGTCTGAATGGCATGGATCTGGCAAAAGAAGTGACCACCGCAGAAGCCTATAGCt                                          <  4:21/119‑1 (MQ=255)
                                                                                                                            cTGCATTATCGCGGGCGATAACCCGGATGCGGCGCTGGCGTTAGAAAAAGCCCGCGCGTTCCCAGGTCTGAATGGCATGGATCTGGCAAAAGAAGTGACCACCGCAGAAGCCTATAGCt                                          >  3:21/1‑119 (MQ=255)
                                                                                                                             tGCATTATCGCGGGCGATAACCCGGATGCGGCGCTGGCGTTAGAAAAAGCCCGCGCGTTCCCAGGTCTGAATGGCATGGATCTGGCAAaagaa                                                                   <  3:167885/93‑1 (MQ=255)
                                                                                                                             tGCATTATCGCGGGCGATAACCCGGATGCGGCGCTGGCGTTAGAAAAAGCCCGCGCGTTCCCAGGTCTGAATGGCATGGATCTGGCAAaagaa                                                                   >  4:167885/1‑93 (MQ=255)
                                                                                                                             tGCATTATCGCGGGCGATAACCCGGATGCGGCGCTGGCGTTAGAAAAAGCCCGCGCGTTCCCAGGTCTGAATGGCATGGATCTGGCAAAAGAAGTGACCACCGCAGAAg                                                   <  4:59748/109‑1 (MQ=255)
                                                                                                                             tGCATTATCGCGGGCGATAACCCGGATGCGGCGCTGGCGTTAGAAAAAGCCCGCGCGTTCCCAGGTCTGAATGGCATGGATCTGGCAAAAGAAGTGACCACCGCAGAAg                                                   >  3:59748/1‑109 (MQ=255)
                                                                                                                              gCATTATCGCGGGCGATAACCCGGATGCGGCGCTGGCGTTAGAAAAAGCCCGCGCGTTCCCAGGTCTGAATGGCATGGATCTGGCAAAAGAAGTGACCACCGCAGAAGCCTATAGCTGGACACAAGGGAGCTGGACGTTGACCGGTgg           <  2:89845/148‑1 (MQ=255)
                                                                                                                                 ttATCGCCCGCGATAACCCGGATGCGGCGCTGGCGTTAGAAAAAGCCCGCGCGTTCCCAGGTCTGAATGGCATGGATCTGGCAAAAGAAGTGACCACCGCAGa                                                     >  2:297359/1‑103 (MQ=255)
                                                                                                                                 ttATCGCCCGCGATAACCCGGATGCGGCGCTGGCGTTAGAAAAAGCCCGCGCGTTCCCAGGTCTGAATGGCATGGATCTGGCAAAAGAAGTGACCACCGCAGa                                                     <  1:297359/103‑1 (MQ=255)
                                                                                                                                 ttATCGCCCGCGATAACCCGGATGCGGCGCTGGCGTTAGAAAAAGCCCGCGCGTTCCCAGGTCTGAATGGCATGGATCTGGCAAAAGAAGTGACCACCGCAGAAGCCTATAGCTGGAcaca                                   <  1:1036544/121‑1 (MQ=255)
                                                                                                                                 ttATCGCCCGCGATAACCCGGATGCGGCGCTGGCGTTAGAAAAAGCCCGCGCGTTCCCAGGTCTGAATGGCATGGATCTGGCAAAAGAAGTGACCACCGCAGAAGCCTATAGCTGGAcaca                                   >  2:1036544/1‑121 (MQ=255)
                                                                                                                                  tATCGCGGGCGATAACCCGGATGCGGCGCTGGCGTTAGAAAAAGCCCGCGCGTTCCCAGGTCTGAATGGCATGGATCTGGCAAaaga                                                                    <  1:841927/87‑1 (MQ=255)
                                                                                                                                  tATCGCGGGCGATAACCCGGATGCGGCGCTGGCGTTAGAAAAAGCCCGCGCGTTCCCAGGTCTGAATGGCATGGATCTGGCAAaaga                                                                    >  2:841927/1‑87 (MQ=255)
                                                                                                                                   aTCGCGGGCGATAACCCGGATGCGGCGCTGGCGTTAGAAAAAGCCCGCGCGTTCCTAGGTCTGAATGGCATGGATCTGGCAAAAGAAGTGACCACCGCAGAAGCCTATAGCTGGACACAAGGGAGCTGGACGTTGACCGGTGGCCTGc      <  2:1109923/148‑1 (MQ=255)
                                                                                                                                       cgGGCGATAACCCGGATGCGGCGCTGGCGTTAGAAAAAGCCCGCGCGTTCCCAGGTCTGAATGGCATGGATCTGGCAAaagaa                                                                   >  2:1059461/1‑83 (MQ=255)
                                                                                                                                       cgGGCGATAACCCGGATGCGGCGCTGGCGTTAGAAAAAGCCCGCGCGTTCCCAGGTCTGAATGGCATGGATCTGGCAAaagaa                                                                   <  1:1059461/83‑1 (MQ=255)
                                                                                                                                       cgGGCGATAACCCGGATGCGGCGCTGGCGTTAGAAAAAGCCCGCGCGTTCCCAGGTCTGAATGGCATGGATCTGGCAAAAGAAGTGACCACCGCAGAAGCCTATAGCTGGACACAAGGGAGCTGGACGTTGACCGGTGGCCTGCCAGa  <  1:1026964/148‑1 (MQ=255)
                                                                                                                                         ggCGATAACCCGGATGCGGCGCTGGCGTTAGAAAAAGCCCGCGCGTTCCCAGGTCTGAATGGCATGGATCTGGCAAAAGAAGTGACCACCGCAGAAGCCtata                                             >  2:1075007/1‑103 (MQ=255)
                                                                                                                                         ggCGATAACCCGGATGCGGCGCTGGCGTTAGAAAAAGCCCGCGCGTTCCCAGGTCTGAATGGCATGGATCTGGCAAAAGAAGTGACCACCGCAGAAGCCtata                                             <  1:1075007/103‑1 (MQ=255)
                                                                                                                                         ggCGATAACCCGGATGCGGCGCTGGCGTTAGAAAAAGCCCGCGCGTTCCCAGGTCTGAATGGCATGGATCTGGCAAAAGAAGTGACCACCGCAGAAGCCTATAGCTGGACACAAg                                 >  4:1080172/1‑115 (MQ=255)
                                                                                                                                         ggCGATAACCCGGATGCGGCGCTGGCGTTAGAAAAAGCCCGCGCGTTCCCAGGTCTGAATGGCATGGATCTGGCAAAAGAAGTGACCACCGCAGAAGCCTATAGCTGGACACAAg                                 <  3:1080172/115‑1 (MQ=255)
                                                                                                                                         |                                                                                                                                                 
CCAGCAACTTCCGTAATACCGAAGACCTCTCTTCTTACCTGAAACGCCATAACATCGTGGCGATTGCCGATATCGATACCCGTAAGCTGACGCGTTTACTGCGCGAGAAAGGCGCACAGAATGGCTGCATTATCGCGGGCGATAACCCGGATGCGGCGCTGGCGTTAGAAAAAGCCCGCGCGTTCCCAGGTCTGAATGGCATGGATCTGGCAAAAGAAGTGACCACCGCAGAAGCCTATAGCTGGACACAAGGGAGCTGGACGTTGACCGGTGGCCTGCCAGA  >  NC_000913/29916‑30198

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: