Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 33,060 G→A 5.6% A748A (GCG→GCA carB → carbamoyl‑phosphate synthase large subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00091333,0600GA 5.6% 225.8 / 2.2 72A748A (GCG→GCAcarBcarbamoyl‑phosphate synthase large subunit
Reads supporting (aligned to +/- strand):  ref base G (29/38);  new base A (2/2);  total (32/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 = 6.15e-01

AAATGGCGGTAGAGAAGGCGAAAGAGATTGGCTACCCGCTGGTGGTACGTCCGTCTTACGTTCTCGGCGGTCGGGCGATGGAAATCGTCTATGACGAAGCTGACCTGCGTCGCTACTTCCAGACGGCGGTCAGCGTGTCTAACGATGCGCCAGTGTTGCTGGACCACTTCCTCGATGACGCGGTAGAAGTTGACGTGGATGCCATCTGCGACGGCGAAATGGTGCTGATTGGCGGCATCATGGAGCATATTGAGCAGGCGGGCGTGCACTCCGGTGACTCCGCATGTTCT  >  NC_000913/32912‑33201
                                                                                                                                                    |                                                                                                                                             
aaaTGGCGGTAGAGAAGGCGAAAGAGATTGGCTACCCGCTGGTGGTACGTCCGTCTTACGTTCTCGGCGGTCGGGCGATGGAAATCGTCTATGACGAAGCTGACCTGCGTCGCTACTTCGAGCCGGCGGTCAGCGTGTCTAACGCTgcg                                                                                                                                               >  8:454412/1‑149 (MQ=255)
 aaTGGCGGTAGAGAAGGCGAAAGAGATTGGCTACCCGCTGGTGGTACGTCCGTCTTACGTTCTCGGCGGTCGGGCGATGGAAATCGTCTATGACGAAGCTGACCTGCGTCGCTACTTCCAGACGGCGGTCAGCGTGTCTAACGATgcgc                                                                                                                                              <  7:122149/149‑1 (MQ=255)
       ggTAGAGAAGGCGAAAGAGATTGGCTACCCGCTGGTGGTACGTCCGTCTTACGTTCTCGGCGGTCGGGCGATGGAAATCGTCTATGACGAAGCTGACCTGCGTCGCTACTTCCAGACGGCGGTCAGCGTGTCTAACGATGCGCCAgtgt                                                                                                                                        <  3:124534/149‑1 (MQ=255)
         tAGAGAAGGCGAAAGAGATTGGCTACCCGCTCGTGGTAAGTCCGTCTTACGTTCTCGGCCGTCGGGCGATGGAAATCTTCTAACACGAAGCTGAACTGCGTCGCTACTTCCAGACGGAGGTCAGCGAGTCTAACGATGCGCCACTGATg                                                                                                                                      >  5:6954/1‑149 (MQ=255)
            agaAGGCGAAAGAGATTGGCTACCCGCTGGTGGTACGTCCGTCTTACGTTCTCGGCGGTCGGGCGATGGAAATCGTCTATGACGAAGCTGACCTGCGTCGCTACTTCCAGACGGCGGTCAGCGTGTCTAACGATgcgc                                                                                                                                              <  8:194632/138‑1 (MQ=255)
            agaAGGCGAAAGAGATTGGCTACCCGCTGGTGGTACGTCCGTCTTACGTTCTCGGCGGTCGGGCGATGGAAATCGTCTATGACGAAGCTGACCTGCGTCGCTACTTCCAGACGGCGGTCAGCGTGTCTAACGATgcgc                                                                                                                                              >  7:194632/1‑138 (MQ=255)
             gaAGGCGAAAGAGATTGGCTACCCGCTGGTGGTACGTCCGTCTTACGTTCTCGGCGGTCGGGCGATGGAAATCGTCTATGACGAAGCTGACCTGCGTCGCTACTTCCAGACGGCGGTCAGCGTGTCTAACGATGCGCCAGTGTTGCTg                                                                                                                                   >  8:463527/1‑148 (MQ=255)
             gaAGGCGAAAGAGATTGGCTACCCGCTGGTGGTACGTCCGTCTTACGTTCTCGGCGGTCGGGCGATGGAAATCGTCTATGACGAAGCTGACCTGCGTCGCTACTTCCAGACGGCGGTCAGCGTGTCTAACGATGCGCCAGTGTTGCTg                                                                                                                                   <  3:144450/148‑1 (MQ=255)
                 gCGAAAGAGATTGGCTACCCGCTGGTGGTACGTCCGTCTTACGTTCTCGGCGGTCGGGCGATGGAAATCGTCTATGACGAAGCTGACCTGCGTCGCTACTTCCAGACGGCGGTCAGCGTGTCTAACGATGCGCCAGTGTTGCTGGACCa                                                                                                                              <  4:395491/149‑1 (MQ=255)
                 gCGAAAGAGATTGGCTACCCGCTGGTGGTACGTCCGTCTTACGTTCTCGGCGGTCGGGCGATGGAAATCGTCTATGACGAAGCTGACCTGCGTCGCTACTTCCAGACGGCGGTCAGCGTGTCTAACGATGCGCCAGTGTTGCTGGACCa                                                                                                                              >  6:64112/1‑149 (MQ=255)
                   gAAAGAGATTGGCTACCCGCTGGTGGTACGTCCGTCTTACGTTCTCGGCGGTCGGGCGATGGAAATCGTCTATGACGAAGCTGACCTGCGTCGCTACTTCCAGACGGCGGTCAGCGTGTCTAACGATGCGCCAGTGTTGCTGGACCACt                                                                                                                            >  8:468101/1‑149 (MQ=255)
                      agagATTGGCTACCCGCTGGTGGTACGTCCGTCTTACGTTCTCGGCGGTCGGGCGATGGAAATCGTCTATGACGAAGCTGACCTGCGTCGCTACTTCCAGACGGCGGTCAGCGTGTCTAACGATGCGCCAGTGTTGCTGGACCACTTcc                                                                                                                         <  7:287086/149‑1 (MQ=255)
                       gagaTTGGCTACCCGCTGGTGGTACGTCCGTCTTACGTTCTCGGCGGTCGGGCGATGGAAATCGTCTATGACGAAGCTGACCTGCGTCGCTACTTCCAGACGGCGGTCAGCGTGTCTAACGATGCGCCAGTGTTGCTGGACCACTTCCt                                                                                                                        <  5:64112/149‑1 (MQ=255)
                        agaTTGGCTACCCGCTGGTGGTACGTCCGTCTTACGTTCTCGGCGGTCGGGCGATGGAAATCGTCTATGACGAAGCTGACCTGCGTCGCTACTTCCAGACGGCGGTCAGCGTGTCTAACGATGCGCCAGTGTTGCTGGACCACTTCCTc                                                                                                                       <  7:229022/149‑1 (MQ=255)
                                   ccGCTGGTGGTACGTCCGTCTTACGTTCTCGGCGGTCGGGCGATGGAAATCGTCTATGACGAAGCTGACCTGCGTCGCTACTTCCAGACGGCGGTCAGCGTGTCTAACGATGCGCCAGTGTTGCTGGACCACTTCCTc                                                                                                                       >  8:102676/1‑138 (MQ=255)
                                   ccGCTGGTGGTACGTCCGTCTTACGTTCTCGGCGGTCGGGCGATGGAAATCGTCTATGACGAAGCTGACCTGCGTCGCTACTTCCAGACGGCGGTCAGCGTGTCTAACGATGCGCCAGTGTTGCTGGACCACTTCCTc                                                                                                                       <  7:102676/138‑1 (MQ=255)
                                    cGCTGGTGGTACGTCCGTCTTACGTTCTCGGCGGTCGGGCGATGGAAATCGTCTATGACGAAGCTGACCTGCGTCGCTACTTCCAGACGGCGGTCAGCGTGTCTAACGATGCGCCAGTGTTGCTGGACCACTTCCTCGATGACGCGGTa                                                                                                           <  2:62075/149‑1 (MQ=255)
                                         gtggtACGTCCGTCTTACGTTCTCGGCGGTCGGGCGATGGAAATCGTCTATGACGAAGCTGACCTGCGTCGCTACTTCCAGACGGCGGTCAGCGTGTCTAACGATGCGCCAGTGTTGCTGGACCACTTCCTCGATGACGCGGTAGAAGt                                                                                                      <  7:195752/149‑1 (MQ=255)
                                           ggtACGTCCGTATTACGTTCTCGGCGGTCGGGCGATGGAACTCGTCTATGACGAAGATGACCTGAGTCGCTACTTCCAGACGGCTGTCAGAGTGTATAACTCTGCTCCAGTGTTGCTGGACCAGTGCCTCGATGACGCGGTAGAAGtt                                                                                                     >  3:35611/1‑148 (MQ=255)
                                                gtccgtcTTACGTTCTCGGCGGTCGGGCGATGGAAATCGTCTATGACGAAGCTGACCTGCGTCGCTACTTCCAGACGGCGGTCAGCGTGTCTAACGATGCGCCAGTGTTGCTGGACCACTTCCTCGATGACGCgg                                                                                                             <  7:236536/135‑1 (MQ=255)
                                                gtccgtcTTACGTTCTCGGCGGTCGGGCGATGGAAATCGTCTATGACGAAGCTGACCTGCGTCGCTACTTCCAGACGGCGGTCAGCGTGTCTAACGATGCGCCAGTGTTGCTGGACCACTTCCTCGATGACGCgg                                                                                                             >  8:236536/1‑135 (MQ=255)
                                                gtccgtcTTACGTTCTAGGCGGTCGGGCGATGGAAATAGTATATTACGAAGCTGACCTGCGTCGCTACTTCCAGACGGCGGTCATCGTGTCTAACGATGCGCCAGTGTTGCTGTACCAATTCCTCTATGACGCGGTAGAAGTTGACGTg                                                                                               >  4:370839/1‑149 (MQ=255)
                                                 tccgtcTTACGTTCTCGGCGGTCGGGCGATGGAAATCGTCTATGACGAAGCTGACCTGCGTCGCTACTTCCAGACGGCGGTCAGCGTGTCTAACGATGCGCCAGTGTTGCTGGACCACTTCCTCGATGACGCGGTAGAAGTTGa                                                                                                   <  2:194753/144‑1 (MQ=255)
                                                 tccgtcTTACGTTCTCGGCGGTCGGGCGATGGAAATCGTCTATGACGAAGCTGACCTGCGTCGCTACTTCCAGACGGCGGTCAGCGTGTCTAACGATGCGCCAGTGTTGCTGGACCACTTCCTCGATGACGCGGTAGAAGTTGa                                                                                                   >  1:194753/1‑144 (MQ=255)
                                                 tccgtcTTACGTTCTCGGCGGTCGGGCGATGGAAATCGTCTATGACGAAGCTGACCTGCGTCGCTACTTCCAGACGGCGGTCAGCGTGTCTAACGATGCGCCAGTGTTGCTGGACCACTTCCTCGATGACGCGGTAGAAGTTGACGTg                                                                                               >  5:149278/1‑148 (MQ=255)
                                                  ccgtcTTACGTTCTCGGCGGTCGGGCGATGGAAATCGTCTATGACGAAGCTGACCTGCGTCGCTACTTCCAGACGGCGGTCAGCGTGTCTAACGATGCGCCAGTGTTGCTGGACCACTTCCTCGATGACGCGGTAGAAGTTGACGTGGa                                                                                             >  2:304414/1‑149 (MQ=255)
                                                    gtcTTACGTTCTCGGCGGTCGGGCGATGGAAATCGTCTATGACGAAGCTGACCTGCGTCGCTACTTCCAGACGGCGGTCAGCGTGTCTAACGATGCGCCAGTGTTGCTGGACCACTTCCTCGATGACGCGGTAGAAGTTGACGTGGATg                                                                                           >  8:331702/1‑149 (MQ=255)
                                                         aCGTTCTCGGCGGTCGGGCGATGGAAATCGTCTATGACGAAGCTGACCTGCGTCGCTACTTCCAGACGGCGGTCAGCGTGTCTAACGATGCGCCAGTGTTGCTGGACCACTTCCTCGATGACGCGGTAGAAGTTGACGTGGATGCCAt                                                                                       <  3:150928/148‑1 (MQ=255)
                                                           gTTCTCGGCGGTCGGGCGATGGAAATCGTCTATGACGAAGCTGACCTGCGTCGCTACTTCCAGACGGCGGTCAGCGTGTCTAACGATGCGCCAGTGTTGCTGGACCACTTCCTCGATGACGCGGTAGAAGTTGACGTGGATGCCATCTg                                                                                    <  7:468101/149‑1 (MQ=255)
                                                               tcGGCGGTCGGGCGATGGAAATCGTCTATGACGAAGCTGACCTGCGTCGCTACTTCCAGACGGCGGTCAGCGTGTCTAACGATGCGCCAGTGTTGCTGGACCACTTCCTCGATGACGCGGTAGAAGTTGACGTGGATGCCATCTGCGa                                                                                 >  7:184262/1‑148 (MQ=255)
                                                               tcGGCGGTCGGGCGATGGAAATCGTCTATGACGAAGCTGACCTGCGTCGCTACTTCCAGACGGCGGTCAGCGTGTCTAACGATGCGCCAGTGTTGCTGGACCACTTCCTCGATGACGCGGTAGAAGTTGACGTGGATGCCATCTGCGAc                                                                                <  7:331702/149‑1 (MQ=255)
                                                                 ggcggTCGGGCGATGGAAATCGTCTATGACGAAGCTGACCTGCGTCGCTACTTCCAGACGGCGGTCAGCGTGTCTAACGATGCGCCAGTGTTGCTGGACCACTTCCTCGATGACGCGGTAGAAGTTGACGTGGATGCCATCTGCGa                                                                                 <  8:184262/146‑1 (MQ=255)
                                                                 ggcggTCGGGCGATGGAAATCGTCTATGACGAAGCTGACCTGCGTCGCTACTTCCAGACGGCGGTCAGCGTGTCTAACGATGCGCCAGTGTTGCTGGACCACTTCCTCGATGACGCGGTAGAAGTTGACGTGGATGCCATCTGCGACg                                                                               <  3:25203/148‑1 (MQ=255)
                                                                    ggTCGGGCGATGGAAATCGTCTATGACGAAGCTGACCTGCGTCGCTACTTCCAGACGGCGGTCAGCGTGTCTAACGATGCGCCAgtg                                                                                                                                         <  6:136369/87‑1 (MQ=255)
                                                                    ggTCGGGCGATGGAAATCGTCTATGACGAAGCTGACCTGCGTCGCTACTTCCAGACGGCGGTCAGCGTGTCTAACGATGCGCCAgtg                                                                                                                                         >  5:136369/1‑87 (MQ=255)
                                                                         ggCGATGGAAATCGTCTATGACGAAGCTGACCTGCGTCGCTACTTCCAGACGGCGGTCAGCGTGTCTAACGATGCGCCAGTGTTGCTGGACCACTTCCTCGATGACGCGGTAGAAGTTGACGTGGATGCCATCTGCGACGGCGAAATgg                                                                      <  7:454412/149‑1 (MQ=255)
                                                                                    tCGTCTATGACGAAGCTGACCTGCGTCGCTACTTCCAGACGGCGGTCAGCGTGTCTAACGATGCGCCAGTGTTGCTGGACCACTTCCTCGATGACGCGGTAGAAGTTGACGTGGATGCCATCTGCGACGGCGAAATGGTGCTGATTggc                                                           >  6:197293/1‑149 (MQ=255)
                                                                                    tCGTCTATGACGAAGCTGACCTGCGTCGCTACTTCCAGACGGCGGTCAGCGTGTCTAACGATGCGCCAGTGTTGCTGGACCACTTCCTCGATGACGCGGTAGAAGTTGACGTGGATGCCATCTGCGACGGCGAAATGGTGCTGATTggc                                                           <  7:86247/149‑1 (MQ=255)
                                                                                         tATGACGAAGCTGACCTGCGTCGCTACTTCCAGACGGCGGTCAGCGTGTCTAACGATGCGCCAGTGTTGCTGGACCACTTCCTCGATGACGCGGTAGAAGTTGACGTGGATGCCATCTGCGACGGCGAAATGGTGCTGATTGGCGGcat                                                      >  4:291412/1‑149 (MQ=255)
                                                                                         tATGACGAAGCTGACCTGCGTCGCTACTTCCAGACGGCGGTCAGCGTGTCTAACGATGCGCCAGTGTTGCTGGACCACTTCCTCGATGACGCGGTAGAAGTTGACGTGGATGCCATCTGCGACGGCGAAATGGTGCTGATTGGCGGcat                                                      >  3:56183/1‑149 (MQ=255)
                                                                                           tGACGAAGCTGACCTGCGTCGCTACTTCCAGACGGCGGGCAGCGTGTCTAACGATGCGCCAGTGTTGCTGGACCACTTCCTCGATGACGCGGTAGAAGTTGACGTGGATGCCATCTGCGACGGCGAAATGGTGCTGATTGGCGGcagc                                                     <  3:235078/148‑3 (MQ=255)
                                                                                             aCGTAGCTGACCTGCGTCGCTACTTCAAGTCGGCGGTCAGCGTGTCTAACGCTGCGCCAGTGTTTCTGGTCCACTTCCTCGATGACGCGGTAGAAGTTGACGTGGATGCCATCTGCTACGGCGAAATTGTGCTGATTGGCGGCATCATg                                                  <  6:6954/149‑1 (MQ=255)
                                                                                               gAAGCTGACCTGCGTCGCTACTTCCAGACGGCGGTCAGCGTGTCTAACGATGCGCCAGTGTTGCTGGACCACTTCCTCGATGACGCGGTAGAAGTTGACGTGGATGCCATCTGCGACGGCGAAATGGTGCTGATTGGCGGCATCATGGa                                                <  8:3556/149‑1 (MQ=255)
                                                                                                      aCCTGCGTCGCTACTTCCAGACGGCGGTCAGCGTGTCTAACGATGCGCCAGTGTTGCTGGACCACTTCCTCGATGACGCGGTAGAAGTTGACGTGGATGCCATCTGCGACGGCGAAATGGTGCTGATTGGCGGCATCATGGAGCATAtt                                         <  4:56183/149‑1 (MQ=255)
                                                                                                      aCCTGCGTCGCTACTTCCAGACGGCGGTCAGCGTGTCTAACGATGCGCCAGTGTTGCTGGACCACTTCCTCGATGACGCGGTAGAAGTTGACGTGGATGCCATCTGCGACGGCGAAATGGTGCTGATTGGCGGCATCATGGAGCATAtt                                         <  7:463527/149‑1 (MQ=255)
                                                                                                        cTGCGTCGCTACTTCCAGACGGCGGTCAGCGTGTCTAACGATGCGCCAGTGTTGCTGGACCACTTCCTCGATGACGCGGTAGAAGTTGACGTGGATGCCATCTGCGACGGCGAAATGGTg                                                                    <  3:296222/120‑1 (MQ=255)
                                                                                                        cTGCGTCGCTACTTCCAGACGGCGGTCAGCGTGTCTAACGATGCGCCAGTGTTGCTGGACCACTTCCTCGATGACGCGGTAGAAGTTGACGTGGATGCCATCTGCGACGGCGAAATGGTg                                                                    >  4:296222/1‑120 (MQ=255)
                                                                                                        cTGCGTCGCTACTTCCAGACGGCGGTCAGCGTGTCTAACGATGCGCCAGTGTTGCTGGACCACTTCCTCGATGACGCGGTAGAAGTTGACGTGGATGCCATCTGCGACGGCGAAATGGTGCTGATTGGCGGCATCATGGa                                                >  3:20109/1‑140 (MQ=255)
                                                                                                        cTGCGTCGCTACTTCCAGACGGCGGTCAGCGTGTCTAACGATGCGCCAGTGTTGCTGGACCACTTCCTCGATGACGCGGTAGAAGTTGACGTGGATGCCATCTGCGACGGCGAAATGGTGCTGATTGGCGGCATCATGGa                                                <  4:20109/140‑1 (MQ=255)
                                                                                                        cTGCGTCGCTACTTCCAGACGGCGGTCAGCGTGTCTAACGATGCGCCAGTGTTGCTGGACCACTTCCTCGATGACGCGGTAGAAGTTGACGTGGATGCCATCTGCGACGGCGAAATGGTGCTGATTGGCGGCATCATGGAGCATATTGa                                       <  2:19805/149‑1 (MQ=255)
                                                                                                                    ttCCAGACGGCGGTCAGCGTGTCTAACGATGCGCCAGTGTTGCTGGACCACTTCCTCGATGACGCGGTAGAAGTTg                                                                                                    >  5:361784/1‑76 (MQ=255)
                                                                                                                    ttCCAGACGGCGGTCAGCGTGTCTAACGATGCGCCAGTGTTGCTGGACCACTTCCTCGATGACGCGGTAGAAGTTg                                                                                                    <  6:361784/76‑1 (MQ=255)
                                                                                                                    ttCAAGACGGCGGTCAGCGTGTCTAACGAGGAGCCAGTTTTGCTGGACCACTTAATCGATGAAGCGGTAGAAGTTGACGTGGATGCAATCTGCGAAGGCTAAATGGTGCTGATTGGCGGCATAATGGAGCATATTTAGCAGGCTGGCGt                           >  8:469899/1‑149 (MQ=255)
                                                                                                                      ccAGACGGCGGTCAGCGTGTCTAACGATGCGCCAGTGTTGCTGGACCACTTCCTCGATGACGCGGTAGAAGTTGACGTGGATGCCATCTGCGACGGCGAAATGGTGCTGATTGGCGGCATCATGGAGCATATTGAGCAGGCGGGCGt                           <  5:197293/147‑1 (MQ=255)
                                                                                                                      ccAGACGGCGGTCAGCGTGTCTAACGATGCGCCAGTGTTGCTGGACCACTTCCTCGATGACGCGGTAGAAGTTGACGTGGATGCCATCTGCGACGGCGAAATGGTGCTGATTGGCGGCATCATGGAGCATATTGAGCAGGCGGGCGt                           <  3:291412/147‑1 (MQ=255)
                                                                                                                        aGACGGCGGTCAGCGTGTCTAACGATGCGCCAGTGTTGCTGGACCACTTCCTCGATGACGCGGTAGAAGTTGACGTGGATGCCATCTGCGACGGCGAAATGGTGCTGATTGGCGGCATCATGGAGCATATTGAGCAGGCGGGCGTGCa                        >  4:222480/1‑148 (MQ=255)
                                                                                                                        aGACGGCGGTCAGCGTGTCTAACGATGCGCCAGTGTTGCTGGACCACTTCCTCGATGACGCGGTAGAAGTTGACGTGGATGCCATCTGCGACGGCGAAATGGTGCTGATTGGCGGCATCATGGAGCATATTGAGCAGGCGGGCGTGCAc                       >  3:364682/1‑149 (MQ=255)
                                                                                                                          aCGGCGGTCAGCGTGTCTAACGATGCGCCAGTGTTGCTGGACCACTTCCTCGATGACGCGGTAGAAGTTGACGTGGATGCCATCTGCGACGGCGAAATGGTGCTGATTGGCGGCATCATGGAGCATATTGAGCAGGCGGGCGTGCACTc                     >  8:218586/1‑149 (MQ=255)
                                                                                                                             gcggTCAGCGTGTCTAACGATGCGCCAGTGTTGCTGGACCACTTCCTCGATGACGCGGTAGAAGTTGACGTGGATGCCATCTGCGACGGCGAAATGGTGCTGATTGGCGGCATCATGGAGCATATTGAGCAGGCGGGCGTGCACTCCgg                  <  3:222480/149‑1 (MQ=255)
                                                                                                                              cggTCAGCGTGTCTAACGATGCGCCAGTGTTGCTGGACCACTTCCTCGATGACGCGGTAGAAGTTGACGTGGATGCCATCTGCGACGGCGAAATGGTGCTGATTGGCGGCATCATGGAGCATATTGAGCAGGCGGGCGTGCACTCCGGt                 <  3:370839/149‑1 (MQ=255)
                                                                                                                                    gCGTGTCTAACGATGCGCCAGTGTTGCTGGACCACTTCCTCGATGACGCGGTAGAAGTTGACGTGGATGCCATCTGCGACGGCGAAATGGTGCTGATTGGCGGCATCATGGAGCATATTGAGCAGGCGGGCGTGCACTCCGGTGACTcc           <  1:302696/149‑1 (MQ=255)
                                                                                                                                    gCGTGTCTAACGATGCGCCAGTGTTGCTGGACCACTTCCTCGATGACGCGGTAGAAGTTGACGTGGATGCCATCTGCGACGGCGAAATGGTGCTGATTGGCGGCATCATGGAGCATATTGAGCAGGCGGGCGTGCACTCCGGTGACTcc           >  2:302696/1‑149 (MQ=255)
                                                                                                                                    gCGTGTCTAACGATGCACCAGTGTTGCTGGACCACTTCCTCGATGACGCGGTAGAAGTTGACGTGGATGCCATCTGCGACGGCGAAATGGTGCTGATTGGCGGCATCATGGAGCATATTGAGCAggc                                 <  4:206242/127‑1 (MQ=255)
                                                                                                                                    gCGTGTCTAACGATGCACCAGTGTTGCTGGACCACTTCCTCGATGACGCGGTAGAAGTTGACGTGGATGCCATCTGCGACGGCGAAATGGTGCTGATTGGCGGCATCATGGAGCATATTGAGCAggc                                 <  4:212892/127‑1 (MQ=255)
                                                                                                                                    gCGTGTCTAACGATGCACCAGTGTTGCTGGACCACTTCCTCGATGACGCGGTAGAAGTTGACGTGGATGCCATCTGCGACGGCGAAATGGTGCTGATTGGCGGCATCATGGAGCATATTGAGCAggc                                 >  3:206242/1‑127 (MQ=255)
                                                                                                                                    gCGTGTCTAACGATGCACCAGTGTTGCTGGACCACTTCCTCGATGACGCGGTAGAAGTTGACGTGGATGCCATCTGCGACGGCGAAATGGTGCTGATTGGCGGCATCATGGAGCATATTGAGCAggc                                 >  3:212892/1‑127 (MQ=255)
                                                                                                                                        gtCTAACGATGCGCCAGTGTTGCTGGACCACTTCCTCGATGACGCGGTAGAAGTTGACGTGGATGCCATCTGCGACGGCGAAATGGTGCTGATTGGCGGCATCATGGAGCATAtt                                         >  2:468932/1‑115 (MQ=255)
                                                                                                                                        gtCTAACGATGCGCCAGTGTTGCTGGACCACTTCCTCGATGACGCGGTAGAAGTTGACGTGGATGCCATCTGCGACGGCGAAATGGTGCTGATTGGCGGCATCATGGAGCATAtt                                         <  1:468932/115‑1 (MQ=255)
                                                                                                                                         tCTAACGATGCGCCAGTGTTGCTGGACCACTTCCTCGATGACGCGGTAGAAGTTGACGTGGATGCCATCTGCGACGGCGAAATGGTGCTGATTGGCGGCATCATGGAGCATATTGAGCAGGCGGGCGTGCACTCCGGTGACTCCGCAt       <  7:218586/148‑1 (MQ=255)
                                                                                                                                          cTAACGATGCGCCAGTGTTGCTGGACCACTTCCTCGATGACGCGGTAGAAGTTGACGTGGATGCCATCTGCGACGGCGAAATGGTGCTGATTGGCGGCATCATGGAGCATATTGAGCAGGCGGGCGTGCACTCCGGTGACTCCGCATGt     >  3:188016/1‑149 (MQ=255)
                                                                                                                                           tAACGATGCGCCAGTGTGGCTGGACCCCTTCCTCGATGAGGCGGTAGAAGTTGACGTGGATGCCCTCTGCGACGGCGAAATGGGGCTGCTTGGCGGCATCATGGAGCATATAGAGCAGGCAGGCGTGCACTCCGGTGACTCCGCATGtt    >  6:128340/1‑149 (MQ=255)
                                                                                                                                             aCGATGCGCCAGTGTTGCTGGACCACTTCCTCGATGACGCGGTAGAAGTTGACGTGGATGCCATCTTCGACGGCGAAATGGTGCTGATTGGCGGCATCATGGAGCATATTGAGCAGGCGTGCGTGCACTCCGGTGACTCCGCATGTtc   >  7:25181/1‑148 (MQ=255)
                                                                                                                                             aCGATGCGCCAGTGTTGCTGGACCACTTCCTCGATGACGCGGTAGAAGTTGACGTGGATGCCATCTGCGACGGCGAAATGGTGCTGATTGGCGGCATCATGGAGCATATTGAGCAGGCGGGCGTGCACTCCGGTGACTCCGCATGtt    <  4:149737/147‑1 (MQ=255)
                                                                                                                                              gtaTGCGCCATTGTTTCTTTACCAGTTCCTCGATTACGCGGTAGAAGTTGAAGTTGATGCCTTCTGCGACTGCGAAATGGTGCTGATTGGAGGCTTCATGGAGCATATTGTGCAGGCGGGAGTGCACTCCGGTGACTCCGCATGTtct  <  6:149278/146‑1 (MQ=255)
                                                                                                                                              cGATGCGCCAGTGTTGCTGGACCACTTCCTCGATGACGCGGTAGAAGTTGACGTGGATGCCATCTGCGACGGCGAAATGGTGCTGATTGGCGGCATCATGGAGCATATTGAGCAGGCGGGCGt                           >  4:175331/1‑123 (MQ=255)
                                                                                                                                              cGATGCGCCAGTGTTGCTGGACCACTTCCTCGATGACGCGGTAGAAGTTGACGTGGATGCCATCTGCGACGGCGAAATGGTGCTGATTGGCGGCATCATGGAGCATATTGAGCAGGCGGGCGt                           <  3:175331/123‑1 (MQ=255)
                                                                                                                                                    gcCAGTGTTGCTGGACCACTTCCTCGATGACGCGGTAGAAGTTGACGTGGATGCCATCTGCGACGGCGAAATGGTGCTGATTggc                                                           >  1:180967/1‑85 (MQ=255)
                                                                                                                                                    gcCAGTGTTGCTGGACCACTTCCTCGATGACGCGGTAGAAGTTGACGTGGATGCCATCTGCGACGGCGAAATGGTGCTGATTggc                                                           <  2:180967/85‑1 (MQ=255)
                                                                                                                                                    gcCAGTGTTGCTGGACCACTTCCTCGATGACGCGGTAGAAGTTGACGTGGATGCCATCTGCGACGGCGAAATGGTGCTGATTGGCGGCATCATgg                                                 <  7:448320/95‑1 (MQ=255)
                                                                                                                                                    gcCAGTGTTGCTGGACCACTTCCTCGATGACGCGGTAGAAGTTGACGTGGATGCCATCTGCGACGGCGAAATGGTGCTGATTGGCGGCATCATgg                                                 >  8:448320/1‑95 (MQ=255)
                                                                                                                                                    |                                                                                                                                             
AAATGGCGGTAGAGAAGGCGAAAGAGATTGGCTACCCGCTGGTGGTACGTCCGTCTTACGTTCTCGGCGGTCGGGCGATGGAAATCGTCTATGACGAAGCTGACCTGCGTCGCTACTTCCAGACGGCGGTCAGCGTGTCTAACGATGCGCCAGTGTTGCTGGACCACTTCCTCGATGACGCGGTAGAAGTTGACGTGGATGCCATCTGCGACGGCGAAATGGTGCTGATTGGCGGCATCATGGAGCATATTGAGCAGGCGGGCGTGCACTCCGGTGACTCCGCATGTTCT  >  NC_000913/32912‑33201

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: