Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 32,912 A→C 8.2% E699A (GAA→GCA) ‡ carB → carbamoyl phosphate synthetase subunit beta

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00091332,9120AC 8.2% 293.4 / 16.3 98E699A (GAA→GCA) ‡carBcarbamoyl phosphate synthetase subunit beta
Reads supporting (aligned to +/- strand):  ref base A (46/44);  new base C (3/5);  total (49/49)
Fisher's exact test for biased strand distribution p-value = 7.15e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.63e-01

CGCGCGCGCTGGAAGCTGCTGGCGTACCGGTTATCGGCACCAGCCCGGATGCTATCGACCGTGCAGAAGACCGTGAACGCTTCCAGCATGCGGTTGAGCGTCTGAAACTGAAACAACCGGCGAACGCCACCGTTACCGCTATTGAAATGGCGGTAGAGAAGGCGAAAGAGATTGGCTACCCGCTGGTGGTACGTCCGTCTTACGTTCTCGGCGGTCGGGCGATGGAAATCGTCTATGACGAAGCTGACCTGCGTCGCTACTTCCAGACGGCGGTCAGCGTGTCTAACGATG  >  NC_000913/32768‑33058
                                                                                                                                                |                                                                                                                                                  
cgcgcgcgcTGGAAGCTGCTGGCGTACCGGTTATCGGCACCAGCCCGGATGCTATCGACCGTGCAGAAGACCGTGAACGCTTCCAGCATGCGGTTGAGCGTCTGAAACTGAAACAACCGGCGAACGCCACCGTTACCGCTATTGAAATgg                                                                                                                                               >  1:326332/1‑150 (MQ=255)
     gcgcTGGAAGCTGCTGGCGTACCGGTTATCGGCACCAGCCCGGATGCTATCGACCGTGCAGAAGACCGTGAACGCTTCCAGCATGCGGTTGAGCGTCTGAAACTGAAACAACCGGCGAACGCCACCGTTACCGCTATTGAAATGGCGGTa                                                                                                                                          <  1:591282/150‑1 (MQ=255)
      cgcTGGAAGCTGCTGGCGTACCGGTTATCGGCACCAGCCCGGATGCTATCGACCGTGCAGAAGACCGTGAACGCTTCCAGCATGCGGTTGAGCGTCTGAAACTGAAACAACCGGCGAACGCCACCGTTACCGCTATTGAAATGGCGGTag                                                                                                                                         >  1:1011717/1‑150 (MQ=255)
       gcTGGAAGCTGCTGGCGTACCGGTTATCGGCACCAGCCCGGATGCTATCGACCGTGCAGAAGACCGTGAACGCTTCCAGCATGCGGTTGAGCGTCTGAAACTGAAACAACCGGCGAACGCCACCGTTACCGCCATTTCAATAGCGGTaac                                                                                                                                        <  2:1068672/150‑3 (MQ=255)
         tGGAAGCTGCTGGCGTACCGGTTATCGGCACCAGCCCGGATGCTATCGACCGTGCAGAAGACCGTGAACGCTTCCAGCATGCGGTTGAGCGTCTGAAACTGAAACAACCGGCGAACGCCACCGTTACCGCCATTTCAATAGCGGTaacgg                                                                                                                                      <  1:1265546/150‑5 (MQ=255)
         gggAAGCTGCTGGCGTACCGGTTATCGGCACCAGCCCGGATGCTATCGACCGTGCAGAAGACCGTGAACGCTTCCAGCATGCGGTTGAGCGTCTGAAACTGAAACAACCGGCGAACGCCACCGTTACCGCCATTTCAATAGCGGTaacgg                                                                                                                                      <  1:1259521/149‑5 (MQ=255)
          ggAAGCTGCTGGCGTACCGGTTATCGGCACCAGCCCGGATGCTATCGACCGTGCAGAAGACCGTGAACGCTTCCAGCATGCGGTTGAGCGTCTGAAACTGAAACAACCGGCGAACGCCACCGTTACCGCTATTGAAATGGCGGTAGAGaa                                                                                                                                     >  2:65158/1‑150 (MQ=255)
          ggAAGCTGCTGGCGTACCGGTTATCGGCACCAGCCCGGATGCTATCGACCGTGCAGAAGACCGTGAACGCTTCCAGCATGCGGTTGAGCGTCTGAAACTGAAACAACCGGCGAACGCCACCGTTACCGCCATTTCAATAGCGGTa                                                                                                                                          <  1:1039731/145‑1 (MQ=255)
          ggAAGCTGCTGGCGTACCGGTTATCGGCACCAGCCCGGATGCTATCGACCGTGCAGAAGACCGTGAACGCTTCCAGCATGCGGTTGAGCGTCTGAAACTGAAACAACCGGCGAACGCCACCGTTACCGCCATTTCAATAGCGGTa                                                                                                                                          >  2:1039731/1‑145 (MQ=255)
           gAAGCTGCTGGCGTACCGGTTATCGGCACCAGCCCGGATGCTATCGACCGTGCAGAAGACCGTGAACGCTTCCAGCATGCGGTTGAGCGTCTGAAACTGAAACAACCGGCGAACGCCACCGTTACCGCCATTTCAATAGCGGTaacggtg                                                                                                                                    <  1:348457/150‑7 (MQ=255)
            aaGCTGCTGGCGTACCGGTTATCGGCACCAGCCCGGATGCTATCGACCGTGCAGAAGACCGTGAACGCTTCCAGCATGAGGTTGAGCGTCTGAAACTGAAACAACCGGCGAACGCCACCGTTACCGCTATTGAAATGGCGGTAGAGAAgg                                                                                                                                   >  2:473587/1‑150 (MQ=255)
             aGCTGCTGGCGTACCGGTTATCGGCACCAGCCGGGATGCTATCGACCGTGCAGAAGACCGTGAACGCTTCCAGCATGCGGTTGAGCGTCTGAAACTGAAACAACCGGCGAACGCCACCGTTACCGCTATTGAAATGGCGGTAGAGAAGGc                                                                                                                                  <  2:1381477/150‑1 (MQ=255)
             aGCTGCTGGCGTACCGGTTATCGGCACCAGCCCGGATGCTATCGACCGTGCAGAAGACCGTGAACGCTTCCAGCATGCGGTTGAGCGTCTGAAACTGAAACAACCGGCGAACGCCACCGTTACCGCTATTGAAATGGCGGTAGAGAAGGc                                                                                                                                  <  2:327192/150‑1 (MQ=255)
             aGCTGCTGGCGTACCGGTTATCGGCACCAGCCCGGATGCTATCGACCGTGCAGAAGACCGTGAACGCTTCCAGCATGCGGTTGAGCGTCTGAAACTGAAACAACCGGCGAACGCCACCGTTACCGCTATTGAAATGGCGGTAGAGAAGGc                                                                                                                                  <  2:1008291/150‑1 (MQ=255)
             aGCTGCTGGCGTACCGGTTATCGGCACCAGCCCGGATGCTATCGACCGTGCAGAAGACCGTGAACGCTTCCAGCATGCGGTTGAGCGTCTGAAACTGAAACAACCGGCGAACGCCACCGTTACCGCTATTGAAATGGCGGTAGAGAAGGc                                                                                                                                  <  2:1008294/150‑1 (MQ=255)
                tgctgGCGTACCGGTTATCGGCACCAGCCCGGATGCTATCGACCGTGCAGAAGACCGTGAACGCTTCCAGCATGCGGTTGAGCGTCTGAAACTGAAACAACCGGCGAACGCCACCGTTACCGCTATTGAAATGGCGGTAGAGAAGGCGaa                                                                                                                               <  2:270359/150‑1 (MQ=255)
                  ctgGCGTACCGGTTATCGGCACCAGCCCGGATGCTATCGACCGTGCAGAAGACCGTGAACGCTTCCAGCATGCGGTTGAGCGTCTGAAACTGAAACAACCGGCGAACGCCACCGTTACCGCTATTGAAATGGCGGTAGAGAAGGCGAAag                                                                                                                             >  1:6949/1‑150 (MQ=255)
                  ctgGCGTACCGGTTATCGGCACCAGCCCGGATGCTATCGACCGTGCAGAAGACCGTGAACGCTTCCAGCATGCGGTTGAGCGTCTGAAACTGAAACAACCGGCGAACGCCACCGTTACCGCTATTGAAATGGCGGTAGAGAAGGCGAAag                                                                                                                             <  1:674852/150‑1 (MQ=255)
                  ctgGCGTACCGGTTATCGGCACCAGCCCGGATGCTATCGACCGTGCAGAAGACCGTGAACGCTTCCAGCATGCGGTTGAGCGTCTGAAACTGAAACAACCGGCGAACGCCACCGTTACCGCTATTGAAATGGCGGTAGAGAAGGCGAAag                                                                                                                             <  1:270379/150‑1 (MQ=255)
                    gacGTACCGGTTATCGGCACCAGCCCGGATGCTATCGACCGTGCAGAAGACCGTGAACGCTTCCAGCATGCGGTTGAGCGTCTGAAACTGAAACAACCGGCGAACGCCACCGTTACCGCTATTGAAATGGCGGTAGAGAAGGCGAAagag                                                                                                                           >  1:954315/3‑150 (MQ=255)
                        tACCGGTTATCGGCACCAGCCCGGATGCTATCGACCGTGCAGAAGACCGTGAACGCTTCCAGCATGCGGTTGAGCGTCTGAAACTGAAACAACCGGCGAACGCCACCGTTACCGCTATTGAAATGGCGGTAGAGAAGGCGAAAGAGATTg                                                                                                                       >  1:593428/1‑150 (MQ=255)
                          ccGGTTATCGGCACCAGCCCGGATGCTATCGACCGTGCAGAAGACCGTGAACGCTTCCAGCATGCGGTTGAGCGTCTGAAACTGAAACAACCGGCGAACGCCACCGTTACCGCTATTGAAATGGCGGTAGAGAAGGCGAAAGAGATTGGc                                                                                                                     >  1:287092/1‑150 (MQ=255)
                           cGGTTATCGGCACCAGCCCGGATGCTATCGACCGTGCAGAAGACCGTGAACGCTTCCAGCATGCGGTTGAGCGTCTGAAACTGAAACAACCGGCGAACGCCACCGTTACCGCTATTGAAATGGCGGTAGAGAAGGCGAAAGAGATTGGCt                                                                                                                    >  1:571052/1‑150 (MQ=255)
                            ggTTATCGGCACCAGCCCGGATGCTATCGACCGTGCAGAAGACCGTGAACGCTTCCAGCATGCGGTTGAGCGTCTGAAACTGAAACAACCGGCGAACGCCACCGTTACCGCTATTGAAATGGCGGTAGAGAAGGCGAAAGAGATTGGCTa                                                                                                                   >  1:1541610/1‑150 (MQ=255)
                            ggTTATCGGCACCAGCCCGGATGCCATCGACCGTGCAGAAGACCGTGAACGCTTCCAGCATGCGGTTGAGCGTCTGAAACTGAAACAACCGGCGAACGCCACCGTTACCGCTATTGAAATGGCGGTAGAGAAGGCGAAAGAGATTGGCTa                                                                                                                   >  2:1242557/1‑150 (MQ=255)
                             gTTATCGGCACCAGCCCGGATGCTATCGACCGTGCAGAAGACCGTGAACGCTTCCAGCATGCGGTTGAGCGTCTGAAACTGAAACAACCGGCGAACGCCACCGTTACCGCTATTGAAATGGCGGTAGAGAAGGCGAAAGAGATTGGCTAc                                                                                                                  >  1:459780/1‑150 (MQ=255)
                               tATCGGCACCAGCCCGGATGCTATCGACCGTGCAGAAGACCGTGAACGCTTCCAGCATGCGGTTGAGCGTCTGAAACTGAAACAACCGGCGAACGCCACCGTTACCGCTATTGAAATGGCGGTAGAGAAGGCGAAAGAGATTGGCTAccc                                                                                                                >  1:1516244/1‑150 (MQ=255)
                                  cGGCACCAGCCCGGATGCTATCGACCGTGCAGAAGACCGTGAACGCTTCCAGCATGCGGTTGAGCGTCTGAAACTGAAACAACCGGCGAACGCCACCGTTACCGCTATTGAAATGGCGGTAGAGAAGGCGAAAGAGATTGGCTACCCGCt                                                                                                             <  1:683114/150‑1 (MQ=255)
                                   ggCACCAGCCCGGATGCTATCGACCGTGCAGAAGACCGTGAACGCTTCCAGCATGCGGTTGAGCGTCTGAAACTGAAACAACCGGCGAACGCCACCGTTACCGCTATTGAAATGGCGGTAGAGAAGGCGAAAGAGATTGGCTACCCGCtg                                                                                                            >  1:1199509/1‑150 (MQ=255)
                                   ggCACCAGCCCGGATGCTATCGACCGTGCAGAAGACCGTGAACGCTTCCAGCATGCGGTTGAGCGTCTGAAACTGAAACAACCGGCGAACGCCACCGTTACCGCTATTGAAATGGCGGTAGAGAAGGCGAAAGAGATTGGCTACCCGCtg                                                                                                            >  1:328561/1‑150 (MQ=255)
                                     cACCAGCCCGGATGCTATCGACCGTGCAGAAGACCGTGAACGCTTCCAGCATGCGGTTGAGCGTCTGAAACTGAAACAACCGGCGAACGCCACCGTTACCGCTATTGAAATGGCGGTAGAGAAGGCGAAAGAGATTGGCTACCCGCtggt                                                                                                          >  1:823355/1‑150 (MQ=255)
                                     cACCAGCCCGGATGCTATCGACCGTGCAGAAGACCGTGAACGCTTCCAGCATGCGGTTGAGCGTCTGAAACTGAAACAACCGGCGAACGCCACCGTTACCGCTATTGAAATGGCGGTAGAGAAGGCGAAAGAGATTGGCTACCCGCtggt                                                                                                          <  2:286131/150‑1 (MQ=255)
                                      aCCAGCCAGGATGCTATCGACCGTGCAGAAGACCGTGAACGCTTCCAGCATGCGGTTGAGCGTCTGAAACTGAAACAACCGGCGAACGCCACCGTTACCGCTATTGAAATGGCGGTAGAGAAGGCGAAAGAGATTGGCTACCCGCtggtg                                                                                                         >  2:671326/1‑150 (MQ=255)
                                       ccAGCCCGGATGCTATCGACCGTGCAGAAGACCGTGAACGCTTCCAGCATGCGGTTGAGCGTCTGAAACTGAAACAACCGGCGAACGCCACCGTTACCGCTATTGAAATGGCGGTAGAGAAGGCGAAAGAGATTGGCTACCCGCtggtgg                                                                                                        >  1:385052/1‑150 (MQ=255)
                                         aGCCCGGATGCTATCGACCGTGCAGAAGACCGTGAACGCTTCCAGCATGCGGTTGAGCGTCTGAAACTGAAACAACCGGCGAACGCCACCGTTACCGCTATTGAAATGGCGGTAGAGAAGGCGAAAGAGATTGGCTACCCGCTGGTGGTa                                                                                                      <  2:1541610/150‑1 (MQ=255)
                                              ggATGCTATCGACCGTGCAGAAGACCGTGAACGCTTCCAGCATGCGGTTGAGCGTCTGAAACTGAAACAACCGGCGAACGCCACCGTTACCGCTATTGAAATGGCGGTAGAGAAGGCGAAAGAGATTGGCTACCCGCTGGTGGTAcgtcc                                                                                                 <  2:385052/150‑1 (MQ=255)
                                                 tGCTATCGACCGTGCAGAAGACCGTGAACGCTTCCAGCATGCGGTTGAGCGTCTGAAACTGAAACAACCGGCGAACGCCACCGTTACCGCTATTGAAATGGCGGTAGAGAAGGCGAAAGAGATTGGCTACCCGCTGGTGGTAcgtccgtc                                                                                              <  2:398085/150‑1 (MQ=255)
                                                 tGCTATCGACCGTGCAGAAGACCGTGAACGCTTCCAGCATGCGGTTGAGCGTCTGAAACTGAAACAACCGGCGAACGCCACCGTTACCGCTATTGAAATGGCGGTAGAGAAGGCGAAAGAGATTGGCTACCCGCTGGTGGTAcgtccgtc                                                                                              <  2:185973/150‑1 (MQ=255)
                                                 tGCTATCGACCGTGCAGAAGACCGTGAACGCTTCCAGCATGCGGTTGAGCGTCTGAAACTGAAACAACCGGCGAACGCCACCGTTACCGCTATTGAAATGGCGGTAGAGAAGGCGAAAGAGATTGGCTACCCGCTGGTGGTAcgtccgtc                                                                                              <  2:616321/150‑1 (MQ=255)
                                                 tGCTATCGACCGTGCAGAAGACCGTGAACGCTTCCAGCATGCGGTTGAGCGTCTGAAACTGAAACAACCGGCGAACGCCACCGTTACCGCTATTGAAATGGCGGTAGAGAAGGCGAAAGAGATTGGCTACCCGCTGGTGGTAcgcccgtc                                                                                              <  1:76438/150‑1 (MQ=255)
                                                  gCTATCGACCGTGCAGAAGACCGTGAACGCTTCCAGCATGCGGTTGAGCGTCTGAAACTGAAACAACCGGCGAACGCCACCGTTACCGCTATTGAAATGGCGGTAGAGAAGGCGAAAGAGATTGGCTACCCGCTGGTGGTACGTCCGTCt                                                                                             <  1:473587/150‑1 (MQ=255)
                                                   cTATCGACCCTGCAGAAGACCGTGAACGCTTCCAGCATGCGCTTGAGCGTCTCAAACTGAAACAACCGGCGAACGCCACCGTTACCGCTATTGAATTCGCGGTAGAGAAGGCGAAAGAGATTGGCTACCCGCTGGTGGTAcgtccgtcat                                                                                            >  2:32534/1‑148 (MQ=255)
                                                    tATCGACCGTGCAGAAGACCGTGAACGCTTCCAGCATGCGGTTGAGCGTCTGAAACTGAAACAACCGGCGAACGCCACCGTTACCGCTATTGAAATGGCGGTAGAGAAGGCGAAAGAGATTGGCTACCCGCTGGTGGTACGTCCGTCTTa                                                                                           <  2:1059849/150‑1 (MQ=255)
                                                    tATCGACCGTGCAGAAGACCGTGAACGCTTCCAGCATGCGGTTGAGCGTCTGAAACTGAAACAACCGGCGAACGCCACCGTTACCGCTATTGAAATGGCGGTAGAGAAGGCGAAAGAGATTGGCTACCCGCTGGTGGTACGTCCGTCTTa                                                                                           <  2:1515610/150‑1 (MQ=255)
                                                     aTCGACCGTGCAGAAGACCGTGAACGCTTCCAGCATGCGGTTGAGCGTCTGAAACTGAAACAACCGGCGAACGCCACCGTTACCGCTATTGAAATGGCGGTAGAGAAGGCGAAAGAGATTGGCTACCCGCTGGTGGTACGTCCGTCTTAc                                                                                          >  2:248891/1‑150 (MQ=255)
                                                     aTCGACCGTGCAGAAGACCGTGAACGCTTCCAGCATGCGGTTGAGCGTCTGAAACTGAAACAACCGGCGAACGCCACCGTTACCGCTATTGAAATGGCGGTAGAGAAGGCGAAAGAGATTGGCTACCCGCTGGTGGTACGTCCGTCTTAc                                                                                          <  2:1557779/150‑1 (MQ=255)
                                                          ccGTGCAGAAGACCGTGAACGCTTCCAGCATGCGGTTGAGCGTCTGAAACTGAAACAACCGGCGAACGCCACCGTTACCGCTATTGAAATGGCGGTAGAGAAGGCGAAAGAGATTGGCTACCCGCTGGTGGTACGTCCGTCTTACGTtct                                                                                     <  2:741183/150‑1 (MQ=255)
                                                          ccGTGCAGAAGACCGTGAACGCTTCCAGCATGCGGTTGAGCGTCTGAAACTGAAACAACCGGCGAACGCCACCGTTACCGCTATTGAAATGGCGGTAGAGAAGGCGAAAGAGATTGGCTACACGCTGGTGGTACGTCCGTCTTACGTtct                                                                                     >  2:1611673/1‑150 (MQ=255)
                                                              gCAGAAGACCGTGAACGCTTCCAGCATGCGGTTGAGCGTCTGAAACTGAAACAACCGGCGAACGCCACCGTTACCGCTATTGAAATGGCGGTAGAGAAGGCGAAAGAGATTGGCTACCCGCTGGTGGTACGTCCGTCTTACGTTCTcggc                                                                                 <  1:1251844/150‑1 (MQ=255)
                                                                agaagaCCGTGAACGCTTCCAGCATGCGGTTGAGCGTCTGAAACTGAAACAACCGGCGAACGCCACCGTTACCGCTATTGAAATGGCGGTAGAGAAGGCGAAAGAGATTGGCTACCCGCTGGTGGTACGTCCGTCTTACGTTCTcggcgg                                                                               <  1:122087/150‑1 (MQ=255)
                                                                 gaagaCCGTGAACGCTTCCAGCATGCGGTTGAGCGTCTGAAACTGAAACAACCGGCGAACGCCACCGTTACCGCTATTGAAATGGCGGTAGAGAAGGCGAAAGAGATTGGCTACCCGCTGGTGGTACGTCCGTCTTACGTTCTCGGCGGt                                                                              >  2:1081530/1‑150 (MQ=255)
                                                                 gaagaCCGTGAACGCTTCCAGCATGCGGTTGAGCGTCTGAAACTGAAACAACCGGCGAACGCCACCGTTACCGCTATTGAAATGGCGGTAGAGAAGGCGAAAGAGATTGGCTACCCGCTGGTGGTACGTCCGTCTTACGTTCTCGGCGGt                                                                              >  2:123032/1‑150 (MQ=255)
                                                                   agaCCGTGAACGCTTCCAGCATGCGGTTGAGCGTCTGAAACTGAAACAACCGGCGAACGCCACCGTTACCGCTATTGAAATGGCGGTAGAGAAGGCGAAAGAGATTGGCTACCCGCTGGTGGTACGTCCGTCTTACGTTCTCGGCGGTCg                                                                            >  2:50428/1‑150 (MQ=255)
                                                                      ccGTGAACGCTTCCAGCATGCGGTTGAGCGTCTGAAACTGAAACAACCGGCGAACGCCACCGTTACCGCTATTGAAATGGCGGTAGAGAAGGCGAAAGAGATTGGCTACCCGCTGGTGGTACGTCCGTCTTACGTTCTCGGCGGTCGGGc                                                                         >  2:557403/1‑150 (MQ=255)
                                                                      ccGTGAACGCTTCCAGCATGCGGCTGAGCGTCTGAAACTGAAACAACCGGCGAACGCCACCGTTACCGCTATTGAAATGGCGGTAGAGAAGGCGAAAGAGATTGGCTACCCGCTGGTGGTACGTCCGTCTTACGTTCTCGGCGGTCGGGc                                                                         >  1:176990/1‑150 (MQ=255)
                                                                       cGTGAACGCTTCCAGCATGCGGTTGAGCGTCTGAAACTGAAACAACCGGCGAACGCCACCGTTACCGCTATTGAAATGGCGGTAGAGAAGGCGAAAGAGATTGGCTACCCGCTGGTGGTACGTCCGTCTTACGTTCTCGGCGGTCGGGCg                                                                        >  2:967979/1‑150 (MQ=255)
                                                                           aaCGCTTCCAGCATGCGGTTGAGCGTCTGAAACTGAAACAACCGGCGAACGCCACCGTTACCGCTATTGAAATGGCGGTAGAGAAGGCGAAAGAGATTGGCTACCCGCTGGTGGTACGTCCGTCTTACGTTCTCGGCGGTCGGGCGATgg                                                                    <  2:684831/150‑1 (MQ=255)
                                                                             cGCTTCCAGCATGCGGTTGAGCGTCTGAAACTGAAACAACCGGCGAACGCCACCGTTACCGCTATTGAAATGGCGGTAGAGAAGGCGAAAGAGATTGGCTACCCGCTGGTGGTACGTCCGTCTTACGTTCTCGGCGGTCGGGCGATGGaa                                                                  <  1:254876/150‑1 (MQ=255)
                                                                             cGCTTCCAGCATGCGGTTGAGCGTCTGAAACTGAAACAACCGGCGAACGCCACCGTTACCGCTATTGAAATGGCGGTAGAGAAGGCGAAAGAGATTGGCTACCCGCTGGTGGTACGTCCGTCTTACGTTCTCGGCGGTCGGGCGATGGaa                                                                  >  2:1625285/1‑150 (MQ=255)
                                                                              gCTTCCAGCATGCGGTTGAGCGTCTGAAACTGAAACAACCGGCGAACGCCACCGTTACCGCTATTGAAATGGCGGTAGAGAAGGCGAAAGAGATTGGCTACCCGCTGGTGGTACGTCCGTCTTACGTTCTCGGCGGTCGGGCGATGGaaa                                                                 <  1:1298840/150‑1 (MQ=255)
                                                                                  ccAGCATGCGGTTGAGCGTCTGAAACTGAAACAACCGGCGAACGCCACCGTTACCGCTATTGAAATGGCGGTAGAGAAGGCGAAAGAGATTGGCTACCCGCTGGTGGTACGTCCGTCTTACGTTCTCGGCGGTCGGGCGATGGAAATCGt                                                             <  2:516848/150‑1 (MQ=255)
                                                                                    aGCATGCGGTTGAGCGTCTGAAACTGAAACAACCGGCGAACGCCACCGTTACCGCTATTGAAATGGCGGTAGAGAAGGCGAAAGAGATTGGCTACCCGCTGGTGGTACGTCCGTCTTACGTTCTCGGCGGTCGGGCGATGGAAATCGTCt                                                           <  2:593428/150‑1 (MQ=255)
                                                                                    aGCATGCGGTTGAGCGTCTGAAACTGAAACAACCGGCGAACGCCACCGTTACCGCTATTGAAATGGCGGTAGAGAAGGCGAAAGAGATTGGCTACCCGCTGGTGGTACGTCCGTCTTACGTTCTCGGCGGTCGGGCGATGGAAATCGTCt                                                           <  2:159222/150‑1 (MQ=255)
                                                                                     gCATGCGGTTGAGCGTCTGAAACTGAAACAACCGGCGAACGCCACCGTTACCGCTATTGAAATGGCGGTAGAGAAGGCGAAAGAGATTGGCTACCCGCTGGTGGTACGTCCGTCTTACGTTCTCGGCGGTCGGGCGATGGAAATCGTCTa                                                          <  1:1625285/150‑1 (MQ=255)
                                                                                     gCATGCGGTTGAGCGTCTGAAACTGAAACAACCGGCGAACGCCACCGTTACCGCTATTGAAATGGCGGTAGAGAAGGCGAAAGAGATTGGCTACCCGCTGGTGGTACGTCCGTCTTACGTTCTCGGCGGTCGGGCGATGGAAATCGTCTa                                                          <  2:6949/150‑1 (MQ=255)
                                                                                       aTGCGGTTGAGCGTCTGAAACTGAAACAACCGGCGAACGCCACCGTTACCGCTATTGAAATGGCGGTAGAGAAGGCGAAAGAGATTGGCTACCCGCTGGTGGTACGTCCGTCTTACGTTCTCGGCGGTCGGGCGATGGAAATCGTCTATg                                                        >  1:699988/1‑150 (MQ=255)
                                                                                        tGCGGTTGAGCGTCTGAAACTGAAACAACCGGCGAACGCCACCGTTACCGCTATTGAAATGGCGGTAGAGAAGGCGAAAGAGATTGGCTACCCGCTGGTGGTACGTCCGTCTTACGTTCTCGGCGGTCGGGCGATGGAAATCGTCTATGa                                                       >  2:849122/1‑150 (MQ=255)
                                                                                        tGCGGTTGAGCGTCTGAAACTGAAACAACCGGCGAACGCCACCGTTACCGCTATTGAAATGGCGGTAGAGAAGGCGAAAGAGATTGGCTACCCGCTGGTGGTACGTCCGTCTTACGTTCTCGGCGGTCGGGCGATGGAAATCGTCTATGa                                                       >  1:234130/1‑150 (MQ=255)
                                                                                        tGCGGTTGAGCGTCTGAAACTGAAACAACCGGCGAACGCCACCGTTACCGCTATTGAAATGGCGGTAGAGAAGGCGAAAGAGATTGGCTACCCGCTGGTGGTACGTCCGTCTTACGTTCTCGGCGGTCGGGCGATGGAAATCGTCTATGa                                                       <  1:1611673/150‑1 (MQ=255)
                                                                                              tGAGCGTCTGAAACTGAAACAACCGGCGAACGCCACCGTTACCGCTATTGAAATGGCGGTAGAGAAGGCGAAAGAGATTGGCTACCCGCTGGGGGGACGTCCGTCTTACGTTCTCGGCGGTCGGGCGATGGAAATCGGCTATGACGAAGc                                                 <  1:1152572/150‑1 (MQ=255)
                                                                                                    tCTGAAACTGAAACAACCGGCGAACGCCACCGTTACCGCTATTGAAATGGCGGTAGAGAAGGCGAAAGAGATTGGCTACCCGCTGGTGGTACGTCCGTCTTACGTTCTCGGCGGTCGGGCTATGGAAATCGTCTATGACGAAGCTGACCt                                           >  2:62514/1‑150 (MQ=255)
                                                                                                     ctgaaactgaaaCAACCGGCGAACGCCACCGTTACCGCTATTGAAATGGCGGTAGAGAAGGCGAAAGAGATTGGCTACCCGCTGGTGGTACGTCCGTCTTACGTTCTCGGCGGTCGGGCGATGGAAATCGTCTATGACGAAGCTGACCTg                                          >  2:1612420/1‑150 (MQ=255)
                                                                                                       gaaactgaaacAACCGGCGAACGCCACCGTTACCGCTATTGAAATGGCGGTAGAGAAGGCGAAAGAGATTGGCTACCCGCTGGTGGTACGTCCGTCTTACGTTCTCGGCGGTCGGGCGATGGAAATCGTCTATGACGAAGCTGACCTGCg                                        <  1:914957/150‑1 (MQ=255)
                                                                                                          actgaaacAACCGGCGAACGCCACCGTTACCGCTATTGAAATGGCGGTAGAGAAGGCGAAAGAGATTGGCTACCCGCTGGTGGTACGTCCGTCTTACGTTCTCGGCGGTCGGGCGATGGAAATCGTCTATGACGAAGCTGACCTGCGTCg                                     >  1:397319/1‑150 (MQ=255)
                                                                                                             gaaacAACCGGCGAACGCCACCGTTACCGCTATTGAAATGGCGGTAGAGAAGGCGAAAGAGATTGGCTACCCGCTGGTGGTACGTCCGTCTTACGTTCTCGGCGGTCGGGCGATGGAAATCGTCTATGACGAAGCTGACCTGCGTCGCTa                                  >  1:1537634/1‑150 (MQ=255)
                                                                                                              aaacAACCGGCGAACGCCACCGTTACCGCTATTGAAATGGCGGTAGAGAAGGCGAAAGAGATTGGCTACCCGCTGGTGGTACGTCCGTCTTACGTTCTCGGCGGTCGGGCGATGGAAATCGTCTATGACGAAGCTGACCTGCGTCGCTAc                                 <  2:667715/150‑1 (MQ=255)
                                                                                                               aacaacCGGCGAACGCCACCGTTACCGCTATTGAAATGGCGGTAGAGAAGGCGAAAGAGATTGGCTACCCGCTGGTGGTACGTCCGTCTTACGTTCTCGGCGGTCGGGCGATGGAAATCGTCTATGACGAAGCTGACCTGCGTCGCTACt                                >  2:1619045/1‑150 (MQ=255)
                                                                                                                  aacCGGCGAACGCCACCGTTACCGCTATTGAAATGGCGGTAGAGAAGGCGAAAGAGATTGGCTACCCGCTGGTGGTACGTCCGTCTTACGTTCTCGGCGGTCGGGCGATGGAAATCGTCTATGACGAAGCTGACCTGCGTCGCTACTTcc                             <  2:848575/150‑1 (MQ=255)
                                                                                                                  aacCGGCGAACGCCACCGTTACCGCTATTGAAATGGCGGTAGAGAAGGCGAAAGAGATTGGCTACCCGCTGGTGGTACGTCCGTCTTACGTTCTCGGCGGTCGGGCGATGGAAATCGTCTATGACGAAGCTGACCTGCGTCGCTACTTcc                             <  1:119580/150‑1 (MQ=255)
                                                                                                                    ccGGCGAACGCCACCGTTACCGCTATTGAAATGGCGGTAGAGAAGGCGAAAGAGATTGGCTACCCGCTGGTGGTACGTCCGTCTTACGTTCTCGGCGGTCGGGCGATGGAAATCGTCTATGACGAAGCTGACCTGCGTCGCTACTTCCAg                           >  2:1259563/1‑150 (MQ=255)
                                                                                                                     cGGCGAACGCCACCGTTACCGCTATTGAAATGGCGGTAGAGAAGGCGAAAGAGATTGGCTACCCGCTGGTGGTACGTCCGTCTTACGTTCTCGGCGGTCGGGCGATGGAAATCGTCTATGACGAAGCTGACCTGCGTCGCTACTTCCAGa                          >  1:490926/1‑150 (MQ=255)
                                                                                                                      ggCGAACGCCACCGTTACCGCTATTGAAATGGCGGTAGAGAAGGCGAAAGAGATTGGCTACCCGCTGGTGGTACGTCCGTCTTACGTTCTCGGCGGTCGGGCGATGGAAATCGTCTATGACGAAGCTGACCTGCGTCGCTACTTCCAGAc                         <  2:326332/150‑1 (MQ=255)
                                                                                                                        cGAACGCCACCGTTACCGCTATTGAAATGGCGGTAGAGAAGGCGAAAGAGATTGGCTACCCGCTGGTGGTACGTCCGTCTTACGTTCTCGGCGGTCGGGCGATGGAAATCGTCTATGACGAAGCTGACCTGCGTCGCTACTTCCAGAcgg                       >  2:1614128/1‑150 (MQ=255)
                                                                                                                         gAACGCCACCGTTACCGCTATTGAAATGGCGGTAGAGAAGGCGAAAGAGATTGGCTACCCGCTGGTGGTACGTCCGTCTTACGTTCTCGGCGGTCGGGCGATGGAAATCGTCTATGACGAAGCTGACCTGCGTCGCTACTTCCAGAcggc                      <  1:671326/150‑1 (MQ=255)
                                                                                                                              ccACCGTTACCGCTATTGAAATGGCGGTAGAGAAGGCGAAAGAGATTGGCTACCCGCTGGTGGTACGTCCGTCTTACGTTCTCGGCGGTCGGGCGATGGAAATCGTCTATGACGAAGCTGAcc                                            <  1:181420/123‑1 (MQ=255)
                                                                                                                              ccACCGTTACCGCTATTGAAATGGCGGTAGAGAAGGCGAAAGAGATTGGCTACCCGCTGGTGGTACGTCCGTCTTACGTTCTCGGCGGTCGGGCGATGGAAATCGTCTATGACGAAGCTGAcc                                            >  2:181420/1‑123 (MQ=255)
                                                                                                                              ccACCGTTACCGCTATTGAAATGGCGGTAGAGAAGGCGAAAGAGATTGGCTACCCGCTGGTGGTACGTCCGTCTTACGTTCTCGGCGGTCGGGCGATGGAAATCGTCTATGACGAAGCTGACCTGCGTCGCTACTTCCAGACGGCGGTCa                 >  2:390792/1‑150 (MQ=255)
                                                                                                                                 tctctACCGCCATTTCAATAGCGGTAGAGAAGGCGAAAGAGATTGGCTACCCGCTGGTGGTACGTCCGTCTTACGTTCTCGGCGGTCGGGCGATGGAAATCGTCTATGACGAAGCTGACCTGCGTCGCTACTTCCAGACGGCGGTCAGCg              >  1:948287/5‑150 (MQ=255)
                                                                                                                                 ccGTTACCGCTATTGAAATGGCGGTAGAGAAGGCGAAAGAGATTGGCTACCCGCTGGTGGTACGTCCGTCTTACGTTCTCGGCGGTCGGGCGATGGAAATCGTCTATGACGAAGCTGACCTGCGTCGCTACTTCCAGACGGCGGTCAGCg              <  1:1242557/150‑1 (MQ=255)
                                                                                                                                   tctACCGCCATTTCAATAGCGGTAGAGAAGGCGAAAGAGATTGGCTACCCGCTGGTGGTACGTCCGTCTTACGTTCTCGGCGGTCGGGCGATGGAAATCGTCTATGACGAAGCTGACCTGCGTCGCTACTTCCAGACGGCGGTCAGCgtg            >  2:1468967/3‑150 (MQ=255)
                                                                                                                                   gTTACCGCTATTGAAATGGCGGTAGAGAAGGCGAAAGAGATTGGCTACCCGCTGGTGGTACGTCCGTCTTACGTTCTCGGCGGTCGGGCGATGGAAATCGTCTATGACGAAGCTGACCTGCGTCGCTACTTCCAGACGGCGGTCAGCgtg            <  2:823355/150‑1 (MQ=255)
                                                                                                                                     tACCGCTATTGAAATGGCGGTAGAGAAGGCGAAAGAGATTGGCTACCCGCTGGTGGTACGTCCGTCTTACGTTCTCGGCGGTCGGGCGATGGAAATCGTCTATGACGAAGCTGACCTGCGTCGCTACTTCCAGACGGCGGTCAGCGTGTc          <  2:1516244/150‑1 (MQ=255)
                                                                                                                                     tACCGCTATTGAAATGGCGGTAGAGAAGGCGAAAGAGATTGGCTACCCGCTGGTGGTACGTCCGTCTTACGTTCTCGGCGGTCGGGCGATGGAAATCGTCTATGACGAAGCTGACCTGCGTCGCTACTTCCAGACGGCGGTCAGCGTGTc          <  2:459780/150‑1 (MQ=255)
                                                                                                                                      aCCGCTATTGAAATGGCGGTAGAGAAGGCGAAAGAGATTGGCTACCCGCTGGTGGTACGTCCGTCTTACGTTCTCGGCGGTCGGGCGATGGAAATCGTCTATGACGAAGCTGACCTGCGTCGCTACTTCCAGACGGCGGTCAGCGTGTCt         >  2:884922/1‑150 (MQ=255)
                                                                                                                                       ccGCTATTGAAATGGCGGTAGAGAAGGCGAAAGAGATTGGCTACCCGCTGGTGGTACGTCCGTCTTACGTTCTCGGCGGTCGGGCGATGGAAATCGTCTATGACGAAGCTGACCTGCGTCGCTACTTCCAGACGGCGGTCAGCGTGTCTa        >  2:421358/1‑150 (MQ=255)
                                                                                                                                        cGCTATTGAAATGGCGGTAGAGAAGGCGAAAGAGATTGGCTACCCGCTGGTGGTACGTCCGTCTTACGTTCTCGGCGGTCGGGCGATGGAAATCGTCTATGACGAAGCTGACCTGCGTCGCTACTTCCAGACGGCGGTCAGCGTGTCTaa       >  1:467487/1‑150 (MQ=255)
                                                                                                                                          cTATTGAAATGGCGGTAGAGAAGGCGAAAGAGATTGGCTACCCGCTGGTGGTACGTCCGTCTTACGTTCTCGGCGGTCGGGCGATGGAAATCGTCTATGACGAAGCTGACCTGCGTCGCTACTTCCAGACGGCGGTCAGCGTGTCTAACg     >  2:1295646/1‑150 (MQ=255)
                                                                                                                                             ttGAAATGGCGGTAGAGAAGGCGAAAGAGATTGGCTACCCGCTGGTGGTACGTCCGTCTTACGTTCTCGGCGGTCGGGCGATGGAAATCGTCTATGACGAAGCTGACCTGCGTCGCTACTTCCAGACGGCGGTCAGCGTGTCTAACGATg  >  2:458170/1‑150 (MQ=255)
                                                                                                                                                |                                                                                                                                                  
CGCGCGCGCTGGAAGCTGCTGGCGTACCGGTTATCGGCACCAGCCCGGATGCTATCGACCGTGCAGAAGACCGTGAACGCTTCCAGCATGCGGTTGAGCGTCTGAAACTGAAACAACCGGCGAACGCCACCGTTACCGCTATTGAAATGGCGGTAGAGAAGGCGAAAGAGATTGGCTACCCGCTGGTGGTACGTCCGTCTTACGTTCTCGGCGGTCGGGCGATGGAAATCGTCTATGACGAAGCTGACCTGCGTCGCTACTTCCAGACGGCGGTCAGCGTGTCTAACGATG  >  NC_000913/32768‑33058

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: